Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,882,564 G→T 5.9% Q159H (CAG→CAT hyi → Hydroxypyruvate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,882,5640GT 5.9% 236.7 / 2.4 68Q159H (CAG→CAThyiHydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base G (33/31);  new base T (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

TGCGCCACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAAC  >  ACB143/4882430‑4882690
                                                                                                                                      |                                                                                                                              
tGCGCCACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCcag                                                                                                                                >  5:155180/1‑135 (MQ=255)
     cACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACCTCCCCGGCTTCTACCTGAACACCACCCAGCAGGcc                                                                                                                          >  5:146886/1‑136 (MQ=255)
     cACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCATCAGGcc                                                                                                                          >  4:388259/1‑136 (MQ=255)
     cACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCATCAGGcc                                                                                                                          >  4:23220/1‑136 (MQ=255)
     cACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGcc                                                                                                                          >  6:302012/1‑136 (MQ=255)
     cACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGcc                                                                                                                          >  6:253813/1‑136 (MQ=255)
     cACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACACGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGcc                                                                                                                          >  7:342652/1‑136 (MQ=255)
      aCCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGccc                                                                                                                         >  2:393000/1‑136 (MQ=255)
           gCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGa                                                                                                                     >  7:44444/1‑135 (MQ=255)
               aaGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGATCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATc                                                                                                                 >  8:371635/1‑135 (MQ=255)
                    cTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGcc                                                                                                           >  7:33580/1‑136 (MQ=255)
                    cTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGcc                                                                                                           <  8:246968/136‑1 (MQ=255)
                    cTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGcc                                                                                                           >  3:65273/1‑136 (MQ=255)
                       cgtcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCCCCCAGCAGGCCCTGGAAATCCAGGCCg                                                                                                          >  6:348827/1‑134 (MQ=255)
                         tcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGt                                                                                                      <  1:365239/136‑1 (MQ=255)
                         tcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGt                                                                                                      <  2:198211/136‑1 (MQ=255)
                         tcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGt                                                                                                      <  5:354103/136‑1 (MQ=255)
                         tcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGt                                                                                                      >  1:81301/1‑136 (MQ=255)
                         tcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGt                                                                                                      <  3:340124/136‑1 (MQ=255)
                         tcgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGt                                                                                                      >  3:90807/1‑136 (MQ=255)
                          cgAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTg                                                                                                     >  2:180042/1‑136 (MQ=255)
                                 cTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGc                                                                                               >  7:358412/1‑135 (MQ=255)
                                    cgcTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGa                                                                                             <  6:84575/134‑1 (MQ=255)
                                         cgccgcTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTg                                                                                      <  6:131815/136‑1 (MQ=255)
                                            cgcTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTc                                                                                   <  2:242658/136‑1 (MQ=255)
                                            cgcTGACAAGCTCAAGGACGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTc                                                                                   >  5:123960/1‑136 (MQ=255)
                                                 aCAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGc                                                                               >  3:366112/1‑135 (MQ=255)
                                                 aCAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGc                                                                               <  1:17590/135‑1 (MQ=255)
                                                  cAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCaa                                                                             <  2:311647/136‑1 (MQ=255)
                                                   aaGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAAt                                                                            <  2:363165/136‑1 (MQ=255)
                                                         aaGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGa                                                                        <  5:282408/134‑1 (MQ=255)
                                                         aaGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTACTGCAATACTACa                                                                      >  4:109424/1‑136 (MQ=255)
                                                            ggCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCt                                                                   <  8:358412/136‑1 (MQ=255)
                                                                 cGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCATCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACcac                                                              <  3:388259/136‑1 (MQ=255)
                                                                 cGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCATCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACcac                                                              <  3:23220/136‑1 (MQ=255)
                                                                 cGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACcac                                                              <  4:49747/136‑1 (MQ=255)
                                                                 cGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACcac                                                              <  8:282267/136‑1 (MQ=255)
                                                                 cGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACcac                                                              <  5:344728/136‑1 (MQ=255)
                                                                 cGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACcac                                                              <  1:17392/136‑1 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  6:198633/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  7:235355/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  6:177949/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  6:157508/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  3:403016/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  7:74777/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  8:124854/1‑136 (MQ=255)
                                                                   ggATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACAt                                                            >  1:72196/1‑136 (MQ=255)
                                                                    gATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATg                                                           <  8:364965/136‑1 (MQ=255)
                                                                    gATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATg                                                           <  2:30853/136‑1 (MQ=255)
                                                                              gTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATgg                                                 <  6:199931/136‑1 (MQ=255)
                                                                                    gAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGt                                            <  3:153104/135‑1 (MQ=255)
                                                                                    gAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGt                                            <  5:302012/135‑1 (MQ=255)
                                                                                    gAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGt                                            <  6:7829/135‑1 (MQ=255)
                                                                                    gAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGt                                            <  7:163995/135‑1 (MQ=255)
                                                                                      aaTGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGAc                                         >  6:94846/1‑136 (MQ=255)
                                                                                           tCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGc                                    >  7:62356/1‑136 (MQ=255)
                                                                                                  ccGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCAc                              >  3:183418/1‑135 (MQ=255)
                                                                                                       aCATCCCCGGCTTCTACCTGATCACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGa                        <  1:105961/136‑1 (MQ=255)
                                                                                                       aCATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGa                        <  3:9613/136‑1 (MQ=255)
                                                                                                       aCATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGa                        <  2:39122/136‑1 (MQ=255)
                                                                                                        cATCCCAGGCTTCTACCTGAACACCAACCAGCAGGCCCTGGAAATCCATGCCGAAGTTGGCAGCGACAACCTGTTCCTTCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGa                        <  6:47873/135‑1 (MQ=255)
                                                                                                             ccGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCaa                  >  3:36999/1‑136 (MQ=255)
                                                                                                                   tCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTg              >  4:309395/1‑134 (MQ=255)
                                                                                                                   tCTACCTGAAAACCACCCAGAAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACATGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGCCCTGGCTCGCACCATGGAAGCCAACCTg              >  5:251867/1‑134 (MQ=255)
                                                                                                                       ccTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCt         <  5:77093/135‑1 (MQ=255)
                                                                                                                       ccTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCt         <  5:200956/135‑1 (MQ=255)
                                                                                                                              accaccCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAAc  >  1:60843/1‑135 (MQ=255)
                                                                                                                                   ccAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGa                                                                        <  1:291762/60‑1 (MQ=255)
                                                                                                                                   ccAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGa                                                                        >  2:291762/1‑60 (MQ=255)
                                                                                                                                      |                                                                                                                              
TGCGCCACCGTGCGTAAGACCTTCGTCGAGAACCTGCGCTACGCCGCTGACAAGCTCAAGGGCGCCGGGATCCGCCTGGTCATGGAAATGATCAACACCCGCGACATCCCCGGCTTCTACCTGAACACCACCCAGCAGGCCCTGGAAATCCAGGCCGAAGTGGGCAGCGACAACCTGTTCCTGCAATACGACATCTACCACATGCAGATCATGGAAGGTGACCTGGCTCGCACCATGGAAGCCAACCTGAAGCTGATCAAC  >  ACB143/4882430‑4882690

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: