Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,625,708 A→C F11C (TTC→TGC)  tnaA ← tryptophanase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,625,7080AC98.1% 2599.7 / 17.7 677F11C (TTC→TGC) tnaAtryptophanase
Reads supporting (aligned to +/- strand):  ref base A (3/10);  new base C (318/346);  total (321/356)
Fisher's exact test for biased strand distribution p-value = 9.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 678 total aligned reads displayed.

GCGTCACCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGAACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCATTTACTAAAACAAGTAAATA  >  DH1/4625559‑4625855
                                                                                                                                                     |                                                                                                                                                   
 cGTCACCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCAc                                                                                                                                                   >  3:725072/1‑151 (MQ=255)
  gTCACCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACg                                                                                                                                                  <  3:89978/151‑1 (MQ=255)
    cACCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGt                                                                                                                                                >  5:1465673/1‑151 (MQ=255)
     aCCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGtt                                                                                                                                               >  4:1221799/1‑151 (MQ=255)
     aCCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGtt                                                                                                                                               <  3:2237690/151‑1 (MQ=255)
       cGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCa                                                                                                                                             <  5:1916594/151‑1 (MQ=255)
       cGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCa                                                                                                                                             <  5:98478/151‑1 (MQ=255)
       cGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGATTTACTGGCTCAATAACACGAATGCGGCACGGTTc                                                                                                                                              >  5:1588488/1‑150 (MQ=255)
        gCCCCGGTGCCTCTGTCGGTCATTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACTGTTCAg                                                                                                                                            <  3:1221158/151‑1 (MQ=255)
        gCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTc                                                                                                                                              >  5:1701262/1‑149 (MQ=255)
          cccGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAA‑AACACGAATGCGGCATGGTTCAGGga                                                                                                                                         <  6:947714/151‑1 (MQ=255)
             ggTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGgaga                                                                                                                                       <  6:1330147/151‑1 (MQ=255)
                gCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATg                                                                                                                                     <  1:2272058/150‑1 (MQ=255)
                    cTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTa                                                                                                                                 >  1:1842951/1‑150 (MQ=255)
                     tGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTaa                                                                                                                                >  3:1295802/1‑150 (MQ=255)
                      gTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTaaa                                                                                                                               <  7:887049/150‑1 (MQ=255)
                      gTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTaa                                                                                                                                <  7:69255/149‑1 (MQ=255)
                        cGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGtt                                                                                                                            >  8:1957126/1‑151 (MQ=255)
                         ggTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGttt                                                                                                                           <  5:431612/151‑1 (MQ=255)
                           tCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTc                                                                                                                         >  7:1794231/1‑151 (MQ=255)
                              gTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCAt                                                                                                                      <  3:1725566/151‑1 (MQ=255)
                                aaaTCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCAtt                                                                                                                     >  8:1203942/1‑150 (MQ=255)
                                aaaTCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATAACGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGTCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTATAAAGTTTTCCAtt                                                                                                                     >  4:2095199/1‑150 (MQ=255)
                                 aaTCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTAc                                                                                                                   <  7:2120017/151‑1 (MQ=255)
                                 aaTCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTAc                                                                                                                   >  3:1878043/1‑151 (MQ=255)
                                  aTCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTAc                                                                                                                   >  4:1782776/1‑150 (MQ=255)
                                  aTCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTAc                                                                                                                   >  2:731494/1‑150 (MQ=255)
                                     gATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATa                                                                                                                >  1:846237/1‑150 (MQ=255)
                                       tAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATc                                                                                                             >  4:1997799/1‑151 (MQ=255)
                                        aaaaaCATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATcc                                                                                                            <  4:2237639/151‑1 (MQ=255)
                                        aaaaaCATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATcc                                                                                                            >  6:307180/1‑151 (MQ=255)
                                        aaaaaCATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATc                                                                                                             <  8:2189768/150‑1 (MQ=255)
                                         aaaaCATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCt                                                                                                           <  4:642169/151‑1 (MQ=255)
                                         aaaaCATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCt                                                                                                           >  4:820616/1‑151 (MQ=255)
                                            aCATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTc                                                                                                         >  7:1011870/1‑150 (MQ=255)
                                             cATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCat                                                                                                       >  2:620832/1‑151 (MQ=255)
                                              aTCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatt                                                                                                      >  1:2006536/1‑151 (MQ=255)
                                               tCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCattt                                                                                                     >  1:1267446/1‑151 (MQ=255)
                                                  tCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatttat                                                                                                   <  5:1330138/150‑1 (MQ=255)
                                                  tCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCattt                                                                                                     <  4:42902/148‑1 (MQ=255)
                                                  tCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatt                                                                                                      <  2:1314361/147‑1 (MQ=255)
                                                   cGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATTATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAACGTTTTCCATTACATAATCCTTCatttatt                                                                                                  <  7:1323496/150‑1 (MQ=255)
                                                   cGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatttattt                                                                                                 >  5:2096912/1‑151 (MQ=255)
                                                   cGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatttattt                                                                                                 <  3:1490269/151‑1 (MQ=255)
                                                   cGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatttatt                                                                                                  <  6:849859/150‑1 (MQ=255)
                                                   cGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatttatt                                                                                                  <  4:346353/150‑1 (MQ=255)
                                                   cGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGAACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCatttattt                                                                                                 >  6:710931/1‑151 (MQ=255)
                                                      tATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTa                                                                                               >  8:40791/1‑150 (MQ=255)
                                                      tATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTa                                                                                               >  7:2251157/1‑150 (MQ=255)
                                                       aTCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTaa                                                                                              <  2:118854/150‑1 (MQ=255)
                                                       aTCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTa                                                                                               <  5:260602/149‑1 (MQ=255)
                                                       aTCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAAt                                                                                             >  1:2156885/1‑151 (MQ=255)
                                                         ccAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAAtt                                                                                            >  2:958422/1‑150 (MQ=255)
                                                         ccAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTa                                                                                           <  8:1582009/151‑1 (MQ=255)
                                                           agcagGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACa                                                                                         <  6:1165942/151‑1 (MQ=255)
                                                             cagGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGt                                                                                       >  8:1345345/1‑151 (MQ=255)
                                                              agGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGtt                                                                                      >  1:1037883/1‑150 (MQ=255)
                                                              agGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTg                                                                                      <  3:2016255/151‑1 (MQ=255)
                                                               ggAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTg                                                                                      >  8:687255/1‑150 (MQ=255)
                                                               ggAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGa                                                                                     <  3:1943209/151‑1 (MQ=255)
                                                                    gggTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCt                                                                                >  2:448714/1‑151 (MQ=255)
                                                                     ggTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCtt                                                                               >  6:490095/1‑150 (MQ=255)
                                                                        tCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCtgtg                                                                             <  1:1920250/150‑1 (MQ=255)
                                                                         cATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGaa                                                                           <  4:1108909/151‑1 (MQ=255)
                                                                          aTACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAat                                                                          >  8:54508/1‑151 (MQ=255)
                                                                            aCCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAata                                                                         <  7:1184725/150‑1 (MQ=255)
                                                                                 aTTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACAt                                                                   <  2:2246889/151‑1 (MQ=255)
                                                                                 aTTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACAt                                                                   <  6:2008970/151‑1 (MQ=255)
                                                                                     aataatTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTg                                                                <  4:134864/150‑1 (MQ=255)
                                                                                        aatTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCtat                                                            >  5:601475/1‑151 (MQ=255)
                                                                                           tGCCTCTTCACGATGAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATAcc                                                         >  3:1086488/1‑151 (MQ=255)
                                                                                           tGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATAcc                                                         >  6:617183/1‑151 (MQ=255)
                                                                                                  tCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGAtt                                                     <  7:463063/148‑1 (MQ=255)
                                                                                                   cACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAAtt                                                 <  7:1602514/151‑1 (MQ=255)
                                                                                                   cACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAAt                                                  <  3:242645/150‑1 (MQ=255)
                                                                                                     cGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCg                                               <  2:1715310/151‑1 (MQ=255)
                                                                                                       aTAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGcc                                             >  8:1742913/1‑151 (MQ=255)
                                                                                                        tAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGcc                                             >  6:1916062/1‑150 (MQ=255)
                                                                                                        tAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGcc                                             >  7:139399/1‑150 (MQ=255)
                                                                                                         aaGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGcc                                             >  2:1476465/1‑149 (MQ=255)
                                                                                                          aGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGa                                          <  3:359966/151‑1 (MQ=255)
                                                                                                           gcgcGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGcc                                             >  1:1624773/1‑147 (MQ=255)
                                                                                                           gcgcGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGat                                         <  6:1743649/151‑1 (MQ=255)
                                                                                                             gcgAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGatat                                       <  5:1115333/151‑1 (MQ=255)
                                                                                                               gAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCGGatatat                                     >  2:1055909/1‑151 (MQ=255)
                                                                                                                   ggTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATg                                  <  2:2039315/150‑1 (MQ=255)
                                                                                                                    gTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCa                                <  1:753480/151‑1 (MQ=255)
                                                                                                                    gTACGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCa                                <  1:959085/151‑1 (MQ=255)
                                                                                                                      aCGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCa                                >  2:1855577/1‑149 (MQ=255)
                                                                                                                      aCGTTTTACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCa                                >  7:1957268/1‑149 (MQ=255)
                                                                                                                           ttACTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCAcc                          >  5:823584/1‑150 (MQ=255)
                                                                                                                             aCTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACcat                        <  6:1464689/150‑1 (MQ=255)
                                                                                                                             aCTGGCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACcat                        <  1:31732/150‑1 (MQ=255)
                                                                                                                                ggCTCAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCAtt                    <  8:1008286/151‑1 (MQ=255)
                                                                                                                                    cAATAACACGAATGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCATTTa                  <  6:2269733/149‑1 (MQ=255)
                                                                                                                                              aaTGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCATTTACTAAAACAAGta      <  1:145827/151‑1 (MQ=255)
                                                                                                                                               aTGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCATTTACTAAAACAAGtaa     <  8:2249646/151‑1 (MQ=255)
                                                                                                                                               aTGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTCCTATACCGATTAATTCGCCAGATATATATGCAAGCAACATCATATACTAAAACAACta      >  4:1599608/1‑150 (MQ=255)
                                                                                                                                                tGCGGCACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCATTTACTAAAACAAGtaa     >  1:51293/1‑150 (MQ=255)
                                                                                                                                                tGCGGCACGGTTCAGGGAGATGTATAAAGTTTTCCATTACATAAACATTCAATTATTTTAATTACAGACATCCAATTAAATATTAAATATAATATACAGATTAATTCGCCAGATAAATATGCAAGCACAATAATTAAATAAAAAAAAtaa     >  6:1850142/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGTCACCGCCCCGGTGCCGCTGTCGGTCAGTAAATCGATAAAAACATCTTCGCTATCCAGCAGGAACGGGTTCATACCGGATTTAATAATTGCCTCTTCACGATAAGCGCGAGTGGTACGTTTTACTGGCTCAATAACACGAATGCGGAACGGTTCAGGGAGATGTTTAAAGTTTTCCATTACATAATCCTTCATTTATTTTAATTACAGTGATCCCTGTGAATATTACATCTGCTATACCGATTAATTCGCCAGATATATATGCAAGCACCATCATTTACTAAAACAAGTAAATA  >  DH1/4625559‑4625855

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