Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,628,4350GT37.7% 256.0 / 79.5 146V37F (GTT→TTT) lysPlysine‑specific permease
Reads supporting (aligned to +/- strand):  ref base G (13/78);  new base T (55/0);  total (68/78)
Fisher's exact test for biased strand distribution p-value = 4.58e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 153 total aligned reads displayed.

GCTACAATCGCGCCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCACTTACGGTC  >  DH1/1628285‑1628585
                                                                                                                                                      |                                                                                                                                                       
gCTACAATCGCGCCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTttg                                                                                                                                                         <  2:2019432/151‑1 (MQ=255)
      aTCGCGCCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCct                                                                                                                                                   <  2:290261/151‑1 (MQ=255)
        cgcgCCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCctct                                                                                                                                                 <  2:1127343/151‑1 (MQ=255)
         gcgcCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGTACAGGGCTTTTTTTTGTctctt                                                                                                                                                >  1:756361/1‑150 (MQ=255)
         gcgcCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGGCTTTTTTTTTGctctt                                                                                                                                                >  1:1841035/1‑150 (MQ=255)
           gcCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGGACAGGTCTTTttgatcggaagagg                                                                                                                                              >  1:1384256/1‑140 (MQ=255)
           gcCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGgc                                                                                                                                              <  2:477813/151‑1 (MQ=255)
           gcCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGGCTTTTTTTTGCCTCTTgc                                                                                                                                              >  1:1403051/1‑151 (MQ=255)
             cTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGGACAGGTCTTTTTTTTTCCTCTGgcgc                                                                                                                                           >  1:1718659/1‑151 (MQ=255)
                   tttAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGTACAGGTCTTTTTTTTGCCTCTGGCGCAACgct                                                                                                                                     >  1:719507/1‑149 (MQ=255)
                   tttAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGTAAAGGGCTTTTTTTTTCCTCCTGCGCCACGAt                                                                                                                                     >  1:1150352/1‑151 (MQ=255)
                   tttAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCTGTACAGGGCTTTTTTTTGtctccgatcgggag                                                                                                                                        >  1:2170999/1‑139 (MQ=255)
                   tttAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCctc                                                                                                                                                  <  2:2171007/139‑1 (MQ=255)
                   tttAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGGCTTTTTTTTGCCTCTGGCGCAACGAt                                                                                                                                     >  2:1401647/1‑151 (MQ=255)
                    ttAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGAt                                                                                                                                     <  1:1519549/150‑1 (MQ=255)
                     tAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCCTTGGGACAGGGCTTTTTTTTTCCTCTCGCGCCACGAttt                                                                                                                                   >  1:909761/1‑151 (MQ=255)
                     tAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCCGGACCGGGCTTTTTGGTGCCTCTGGCGCCACGGttt                                                                                                                                   >  2:1263027/1‑151 (MQ=255)
                     tAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGAtt                                                                                                                                    <  1:1013163/150‑1 (MQ=255)
                       aGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGAGCAGGTCTTTTTTTTGCCTCTTGCGCAACGCttttt                                                                                                                                 >  2:226604/1‑149 (MQ=255)
                         aTGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGGCTTTTTTTTTCCTCTGGCGCAACGagatcggg                                                                                                                               >  2:323075/1‑144 (MQ=255)
                          tGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGGCTTTTTTTTGCCTCTTGCGCAACGATTtttc                                                                                                                                >  2:173529/1‑149 (MQ=255)
                              tAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTTTTTCCTCTGGCGCAACGATTTCTccggc                                                                                                                            >  2:1451593/1‑149 (MQ=255)
                               aGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTcaggcag                                                                                                                          <  2:2165224/150‑1 (MQ=255)
                                 cATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTc                                                                                                                       <  1:2472162/151‑1 (MQ=255)
                                   tttttGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGAttt                                                                                                                                   <  2:609973/137‑1 (MQ=255)
                                   tttttGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGGCTTTTTGTTGCCTCTGGCGCAACGAttttgatcgggagagcg                                                                                                                     >  1:609972/1‑137 (MQ=255)
                                    ttttGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCTTCTGCGCAGGGCTTTTTTTTGGCTCTGGGGCCACGAGTTCTcagggaggggcgg                                                                                                                    >  1:79613/1‑146 (MQ=255)
                                    ttttGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCcg                                                                                                                     <  2:1841040/150‑1 (MQ=255)
                                     tttGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGTCTTTTTTTTTCCTCTGGCGCAACGATTTTTCAGGCAGGGCcggg                                                                                                                   >  2:1261917/1‑151 (MQ=255)
                                      ttGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCcggg                                                                                                                   <  2:2556096/150‑1 (MQ=255)
                                        gTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGTACAGGGCTTTTTTTTGTCTCTGGCGCCAAGATTTTTCCGGCAGGTCcgggggg                                                                                                                >  2:1481548/1‑151 (MQ=255)
                                          aTGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCCTCGGGACAGGGCTTTTTTTTGGCTTTTGCGCCACGGTTTTTCACGCAGGGCCGGGCgggg                                                                                                              >  1:70713/1‑151 (MQ=255)
                                          aTGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCcgggcggg                                                                                                               <  2:379043/150‑1 (MQ=255)
                                           tGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGTACAGGGCTTTTTTTTTCCTCTGGCGCCACGATTTCTCAGGCAGGTCCggggggggg                                                                                                             >  2:1706286/1‑150 (MQ=255)
                                           tGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGGACAGGGCTTTTTTTTGCCTCTGGCGCAACGATTTTTCCTGCAGGTCCGGGCggggg                                                                                                             >  2:589198/1‑150 (MQ=255)
                                            ggTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCa                                                                                                            <  2:750838/151‑1 (MQ=255)
                                             gTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGGGCAGGGCTTTTTTTTGTCTCCTGGGCAACGATTTTTCCGGGAGGGCCGGGCGGGGGAt                                                                                                           >  2:1864110/1‑151 (MQ=255)
                                             gTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCAt                                                                                                           <  2:756362/151‑1 (MQ=255)
                                               ttCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCAtt                                                                                                          <  1:387702/150‑1 (MQ=255)
                                                tCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGTACAGGTCTTTTTTTTGCCTCTGGCGCAACGAATTTTCCGGCAGGTCCGGGCGGGGCATTgt                                                                                                        >  2:427413/1‑150 (MQ=255)
                                                  cGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGGACAGGGCTTTTTTTTTCCTCTGGCGCAACGATTTTTCAGGCAGGTCCGGGCGGGGGGTTGctc                                                                                                      >  1:584895/1‑151 (MQ=255)
                                                  cGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGctc                                                                                                      <  1:1369071/151‑1 (MQ=255)
                                                   gAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGGACAGGGCTTTTTTTTGTCTCTGGGGCAACGATTTTTCAGGCAGGTCCGGGCGGGGGGTTGctct                                                                                                     >  1:383858/1‑151 (MQ=255)
                                                    aaaCTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGCAAAAGGGTTTTTTTTGGCCCTGGGGGAACGATTTCTTAGGCAGGTGCGGGCGGGGGATTGcttcc                                                                                                    >  1:2400667/1‑148 (MQ=255)
                                                       cTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCCGTACCGGGCTTTTTTTTGCCTCTCGCGCCACGATTTCTCCGGCCGGTCCGGGCGGGGccgtgttttcgtt                                                                                                 >  1:1956614/1‑139 (MQ=255)
                                                         taaaCCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGGACAGGGCTTTTTTTTGGCTCTGGCGGAACGATTTTTTAGGCAGGGCCGGGCGGGGGATTTTTCTCGtatt                                                                                               >  1:92862/2‑150 (MQ=255)
                                                         aaaaCCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGtata                                                                                               <  2:406192/151‑1 (MQ=255)
                                                         aaaaCCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGtat                                                                                                <  2:70713/150‑1 (MQ=255)
                                                         aaaaCCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGATCCCTCGGGTCAGGTTTTTTTTTTGTTTCTGGCGCAACCAGTTTTAAGGCCGGGGCGGGCGGGGCAGTtctccctcga                                                                                                >  1:1193083/1‑144 (MQ=255)
                                                            aCCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGGACAGGGCTTTTTTTTTGCCCTTGCGCAACGATTTTTCAGGCAGGGCCCGGGGGGGGCTTGCTCTCGTATTTGc                                                                                            >  1:2603283/1‑151 (MQ=255)
                                                             ccACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCt                                                                                           <  1:427414/151‑1 (MQ=255)
                                                              cacaGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTg                                                                                          <  2:2425226/151‑1 (MQ=255)
                                                              cacaGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGTCTTTTTTTTTTCTCTTGCGCCACGATTTCTCAGGCAGGGCCGGGCGGGGGCTTTCTCTTGTTTATTCtt                                                                                          >  1:433431/1‑150 (MQ=255)
                                                               acaGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGa                                                                                         <  1:1594577/151‑1 (MQ=255)
                                                                caGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCACTGTACAGGGCTTTTTTTTTGCTCCGGCGCCACCATTTCTCAGGCAGGTCCGGGGGGGGCCTTGCTCTCGTATATGCTGAt                                                                                        >  2:1578674/1‑151 (MQ=255)
                                                                caGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGAt                                                                                        <  2:383857/151‑1 (MQ=255)
                                                                 aGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGTACAGGTCTTTTTTTTGCCTCTTGCGCAACGAATTTTCAGGCAGGTCCGGGGGGGGGATTTCTCTCGTATATTCTGAtt                                                                                       >  2:2423797/1‑151 (MQ=255)
                                                                 aGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCCTCGGTACAGGGCTTTTTGTTGCCTCTGGCGCAACGATTTTTTAGGCAGGGCCGGGGGGGGCCTTTCTCTCGTTTATGCTGAtt                                                                                       >  1:325409/1‑151 (MQ=255)
                                                                   aaGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATTGGGACAGGGCTTTTTTTTGCCCCTTGGGCCACGATTTTTTAGGGAGGTCCCGGGGGGGGCTTTTTCTCGTTTATGCTGGTTgg                                                                                     >  1:1657950/1‑151 (MQ=255)
                                                                    aGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGGTTTCTCGGCACCGGTCTTTTTTTTTGCTCCCGCGCAACGATTTCTCAGGCGGGGCCGGGTGGGGCCTTGCTCCCGTATATGCtcagttgg                                                                                    >  2:1924714/1‑144 (MQ=255)
                                                                        ccGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGa                                                                                <  2:1641382/151‑1 (MQ=255)
                                                                         cGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTgg                                                                                     <  2:1066646/145‑1 (MQ=255)
                                                                          gggCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGGGGAGGGCTTTTTTTTGGCTCCGGCGCAAAGCTTTCTCAGGCAGGGCCGGGGGGGGGCTTGCTCTCTTTTATTCTGAtttggcttagg                                                                              >  1:618992/1‑142 (MQ=255)
                                                                             cTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTgg                                                                                     <  1:1889292/141‑1 (MQ=255)
                                                                               tACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGtt                                                                           <  2:2370083/149‑1 (MQ=255)
                                                                                  gCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCCTCGGTTCAGGGCTTTTTTTTGCCGCTGGCGGAACGAGTTTTTAGGCAGGGCCGGGGGGGGGGTTGGTTTCGTATATGGTGATTGGGCTggggggtttctt                                                                      >  2:1302520/1‑140 (MQ=255)
                                                                                    cGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTcc                                                                    <  2:92862/151‑1 (MQ=255)
                                                                                     gTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCt                                                                   <  1:173529/151‑1 (MQ=255)
                                                                                      tGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGTACAGGGCTTTTTGTTGCCCCTGGCGCAACGATTTTTCCGGGAGGTCCGGGCGGGGGATTGCTTTCGTATATGCTGATTTGCCTGATGGTTttctttctt                                                                  >  1:1631074/1‑150 (MQ=255)
                                                                                          ttAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCATCGGGACAGGGCTTTTTTTTTCCTCTTGCGCAACGATTTTTCAGGCAGGGCCGGGCGGGGGATTTCTCTCCTATATTCTGATTTGGCTGATGGTTTTCTTCCttatt                                                               >  2:169811/1‑149 (MQ=255)
                                                                                            aaaGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGAc                                                             <  2:2603299/150‑1 (MQ=255)
                                                                                             aaGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTgg                                                                                     <  2:681875/125‑1 (MQ=255)
                                                                                             aaGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGGACAGGGCTTTTTTTTTCCTCTGGCGCCACGATTTTTCTGGCCGGTCCCGGGGGGGGGTTGTTCTCGTTTTTGTTGGTTTGGCTGGTGGGTTTCTTTCTTGTGAcc                                                            >  1:1572753/1‑150 (MQ=255)
                                                                                              aGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTTCCTCGGGACAGGGCTTTTTTTTTCCTCTGGCGCCACGATTTTTCCGGCAGGGCCCGGCGGGGCCTTGCTCTCGTATATGCTGATTTGGCTGATGGGTTTCTTCCTGATTAcc                                                            >  2:2226657/1‑149 (MQ=255)
                                                                                                 cgcgTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGtct                                                       <  2:1572758/151‑1 (MQ=255)
                                                                                                   cgTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCg                                                     <  2:162088/151‑1 (MQ=255)
                                                                                                       aCCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGa                                                 <  1:1451588/151‑1 (MQ=255)
                                                                                                        ccTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCt                                                                   <  1:2399370/132‑1 (MQ=255)
                                                                                                        ccTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGGCTTTTTTTTGCCTCTTGCGCAACGATTTTTCAGGCAGGGCCGGGGGGGGCATTGCTCTCGTATATGCTGATTTGCCTGATGGTTTACTTCCt                                                                   >  2:2399381/1‑132 (MQ=255)
                                                                                                         cTGACGATGATTGCCATTGGCGGTTTCCTCGGGACAGGTCTTTTTTTTGCCTCTGGCGCCAAGGATTTTTCGGGAGGGCCCGGCGGGGGACTTCTCCTGTTTTTGGTGATTGGCCTGGTGGTTTACTTTCTTTTGAGCCGTCTCGGGGaa                                                >  1:229116/1‑150 (MQ=255)
                                                                                                         cTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTg                                                  <  1:430477/148‑1 (MQ=255)
                                                                                                           gACGATGATTGCCATTGGCGGTTTCCTCGGGACAGGTCTTTTTTTTGCCCCTGGGGCAACGATTTTTCAGGGGGGGCCGGGGGGGGGGTTTTTTTCGTTTATTTTGTTTTGCCTGTGGGTTTTTTTTCTGTTGACCaaaaacggaaaagc                                              >  2:1367919/1‑137 (MQ=255)
                                                                                                              gatgatTGCCATTGGCGGTTCCATCGGTACAGGTCTGTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCt                                          <  1:2173986/151‑1 (MQ=255)
                                                                                                                        aTTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGc                                <  2:1175652/151‑1 (MQ=255)
                                                                                                                           ggCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCg                              <  2:429700/150‑1 (MQ=255)
                                                                                                                              ggTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGttt                          <  1:1302518/151‑1 (MQ=255)
                                                                                                                              ggTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGttt                          >  1:1060868/1‑151 (MQ=255)
                                                                                                                                  ccATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCgg                      <  2:79613/151‑1 (MQ=255)
                                                                                                                                      cGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCg                  >  2:1609753/1‑151 (MQ=255)
                                                                                                                                        gTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGtt                <  1:169811/151‑1 (MQ=255)
                                                                                                                                          aCAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTg              <  1:424794/151‑1 (MQ=255)
                                                                                                                                           cAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGtt                    <  1:1305690/144‑1 (MQ=255)
                                                                                                                                           cAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGc             <  2:909762/151‑1 (MQ=255)
                                                                                                                                              gTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCAc          <  1:1924709/151‑1 (MQ=255)
                                                                                                                                                 tttttGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGc             >  1:544960/1‑145 (MQ=255)
                                                                                                                                                 tttttGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGc             <  2:544961/145‑1 (MQ=255)
                                                                                                                                                  ttttGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCACTTa       <  1:226605/150‑1 (MQ=255)
                                                                                                                                                    ttgttgCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCACTTACg     <  1:1609748/150‑1 (MQ=255)
                                                                                                                                                      gttgCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCACTTACGGt   <  2:129106/150‑1 (MQ=255)
                                                                                                                                                      gttgCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCACTTACGGTc  <  2:1657955/151‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                       
GCTACAATCGCGCCTCATTTTTAAGATGGATAGCATTTTTGTATGGTTTCCGAAACTAAAACCACAGAAGCGCCGGGCTTACGCCGTGAATTAAAGGCGCGTCACCTGACGATGATTGCCATTGGCGGTTCCATCGGTACAGGTCTTTTTGTTGCCTCTGGCGCAACGATTTCTCAGGCAGGTCCGGGCGGGGCATTGCTCTCGTATATGCTGATTGGCCTGATGGTTTACTTCCTGATGACCAGTCTCGGTGAACTGGCTGCATATATGCCGGTTTCCGGTTCGTTTGCCACTTACGGTC  >  DH1/1628285‑1628585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: