Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,191,9760AG53.1% 111.9 / inf 581pseudogene (626/788 nt)phnDphosphonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (151/121);  new base G (156/152);  total (308/273)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-03

Only 100 of 588 total aligned reads displayed.

GGTGGAAGCGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACGGTTGCGTAAG  >  DH1/4191826‑4192126
                                                                                                                                                      |                                                                                                                                                      
  tGGAAGCGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCCCCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAAc                                                                                                                                                      >  1:103878/1‑151 (MQ=255)
  gggAAGCGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCAGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGAc                                                                                                                                                      <  7:1204056/150‑1 (MQ=255)
     aaGCGATTGAGCCCGGCGCGGTGGATGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACcgt                                                                                                                                                   <  4:1699886/151‑1 (MQ=255)
     aaGCGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACcgt                                                                                                                                                   >  5:123747/1‑151 (MQ=255)
      aGCGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACcgtc                                                                                                                                                  >  8:395163/1‑151 (MQ=255)
        cGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGCCcgtcgt                                                                                                                                                >  2:53676/1‑151 (MQ=255)
        cGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACcgtcgt                                                                                                                                                >  8:519885/1‑151 (MQ=255)
         gATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCgtcgtc                                                                                                                                               <  1:1960130/151‑1 (MQ=255)
            tGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGc                                                                                                                                            >  7:530997/1‑151 (MQ=255)
             gAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGCTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGc                                                                                                                                            >  4:278992/1‑150 (MQ=255)
                 ccggcccggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGACACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCAAACGTCCGCTCGGCGACCGTCGTCGGCATGGt                                                                                                                                       >  1:1607842/1‑151 (MQ=255)
                  cggcccggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTc                                                                                                                                      >  8:1402563/1‑151 (MQ=255)
                  cggcccggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTc                                                                                                                                      <  8:1688227/151‑1 (MQ=255)
                   ggcccggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCg                                                                                                                                     >  1:1869197/1‑151 (MQ=255)
                   ggcccggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCg                                                                                                                                     >  2:4192/1‑151 (MQ=255)
                      ccggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGgcg                                                                                                                                  >  5:169851/1‑151 (MQ=255)
                        ggTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCAAAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGgcgc                                                                                                                                 >  8:1177426/1‑150 (MQ=255)
                                 cATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATc                                                                                                                       <  5:82771/151‑1 (MQ=255)
                                 cATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATc                                                                                                                       >  7:1956803/1‑151 (MQ=255)
                                 cATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGAtg                                                                                                                       >  1:1459004/1‑150 (MQ=255)
                                 cATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGGGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATc                                                                                                                       <  3:962798/151‑1 (MQ=255)
                                  aTTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCg                                                                                                                      <  8:1553966/151‑1 (MQ=255)
                                     cgcgCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGc                                                                                                                    <  7:271928/150‑1 (MQ=255)
                                         ccACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCAc                                                                                                               >  5:328619/1‑151 (MQ=255)
                                           aCCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCAccc                                                                                                             >  1:553675/1‑151 (MQ=255)
                                              ggTGCCAACAAGCTCGAAGAGATGCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCtgt                                                                                                          >  6:803465/1‑151 (MQ=255)
                                                  ccAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTggg                                                                                                       >  4:635381/1‑150 (MQ=255)
                                                   caacaaGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGa                                                                                                      >  7:667315/1‑150 (MQ=255)
                                                    aacaaGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAg                                                                                                    <  1:1247199/151‑1 (MQ=255)
                                                     acaaGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGc                                                                                                   <  5:1763857/151‑1 (MQ=255)
                                                       aaGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGa                                                                                                 <  2:1799372/151‑1 (MQ=255)
                                                       aaGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATAGGCGTCACCCTGTGGGAAGCg                                                                                                  >  5:1040524/1‑150 (MQ=255)
                                                        aGCTCGAAGAGATCCTGTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGAt                                                                                                <  3:1708047/151‑1 (MQ=255)
                                                        aGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGATTCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGAt                                                                                                <  4:394048/151‑1 (MQ=255)
                                                        aGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGAt                                                                                                <  1:925181/151‑1 (MQ=255)
                                                          gtCGTAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCAGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATGGGCGTCACCCTGTGGGAAGCGATTc                                                                                              <  6:2012337/150‑1 (MQ=255)
                                                          cTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTc                                                                                              >  8:938682/1‑151 (MQ=255)
                                                           gcGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTcg                                                                                             <  6:1162760/150‑1 (MQ=255)
                                                            cGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTc                                                                                              <  4:1424568/149‑1 (MQ=255)
                                                             gAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTcgcg                                                                                           >  1:749700/1‑151 (MQ=255)
                                                             gAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTcgcg                                                                                           <  6:946380/151‑1 (MQ=255)
                                                              aaGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCgg                                                                                          >  8:1907370/1‑151 (MQ=255)
                                                               agagATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGt                                                                                         >  6:1528447/1‑151 (MQ=255)
                                                               agagATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGt                                                                                         <  5:1972170/151‑1 (MQ=255)
                                                                 agaTCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGttt                                                                                       >  6:1776606/1‑151 (MQ=255)
                                                                   aTCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTcc                                                                                     >  3:1088067/1‑151 (MQ=255)
                                                                     ccTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTcc                                                                                     >  8:281153/1‑149 (MQ=255)
                                                                     ccTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTcc                                                                                     >  2:2000149/1‑149 (MQ=255)
                                                                      ctctACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGt                                                                                  >  7:43731/1‑151 (MQ=255)
                                                                       tctACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGtt                                                                                 >  3:1258474/1‑151 (MQ=255)
                                                                        ctACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGt                                                                                  <  7:1836364/149‑1 (MQ=255)
                                                                        ctACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGCTCGGCGCAGGCGGGATCGGCGTCCCCCTGTGGGAAGCGATTCGCGGTTTCCAGTTc                                                                                >  3:617512/1‑151 (MQ=255)
                                                                         tACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTcc                                                                               >  1:27043/1‑151 (MQ=255)
                                                                          aCGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCca                                                                              <  5:1271755/151‑1 (MQ=255)
                                                                          aCGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCca                                                                              <  3:327012/151‑1 (MQ=255)
                                                                          aCGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGtt                                                                                 <  6:583168/148‑1 (MQ=255)
                                                                           cGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCcaa                                                                             <  3:1435946/151‑1 (MQ=255)
                                                                              cGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCcaacaa                                                                          <  2:1573952/151‑1 (MQ=255)
                                                                              cGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCcaacaa                                                                          <  8:1758168/151‑1 (MQ=255)
                                                                                   tgcCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTg                                                                     >  3:414648/1‑151 (MQ=255)
                                                                                     ccACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTgcg                                                                   >  4:837516/1‑151 (MQ=255)
                                                                                      cacaGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTgcgc                                                                  >  5:224536/1‑151 (MQ=255)
                                                                                            tGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGAt                                                            >  1:1538511/1‑151 (MQ=255)
                                                                                                gCCACTGCTGATCTCCTACTCCGTCTATCGCTTCGATTCGTTCGTCCGCTCGGCAAGCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAAGAAAGCTGCGCCCTGATGGt                                                         <  2:864051/150‑1 (MQ=255)
                                                                                                 ccACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGc                                                       <  5:935007/151‑1 (MQ=255)
                                                                                                  cACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCt                                                      >  2:731231/1‑151 (MQ=255)
                                                                                                  cACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCt                                                      >  8:1251966/1‑151 (MQ=255)
                                                                                                   aCTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCtt                                                     >  1:1315414/1‑151 (MQ=255)
                                                                                                   aCTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCtt                                                     >  8:1070499/1‑151 (MQ=255)
                                                                                                        tgATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTatcat                                                >  1:1797615/1‑151 (MQ=255)
                                                                                                         gATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTatcatc                                               >  3:782580/1‑151 (MQ=255)
                                                                                                         gATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTatcat                                                <  3:1911422/150‑1 (MQ=255)
                                                                                                           tctcCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATCGTCGGCGCAGGCGGGATCGCCGTCACCCTCTGGGAAGCGATTCGCCGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTAATCATCg                                              >  4:1242545/1‑150 (MQ=255)
                                                                                                               cTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTAGCATCGTCACg                                         >  2:1139551/1‑151 (MQ=255)
                                                                                                                 aCTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGt                                       >  2:1400569/1‑151 (MQ=255)
                                                                                                                 aCTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACgg                                        <  6:42503/150‑1 (MQ=255)
                                                                                                                    cccTCTATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTc                                      <  5:1353715/149‑1 (MQ=255)
                                                                                                                       tctATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGccg                                 >  8:796000/1‑150 (MQ=255)
                                                                                                                        ctATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGcc                                  <  8:667177/149‑1 (MQ=255)
                                                                                                                         tATCGCTTCGAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCctgc                               >  7:1414444/1‑151 (MQ=255)
                                                                                                                          aTCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCct                                 <  4:1559768/148‑1 (MQ=255)
                                                                                                                          aTCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCct                                 <  5:469383/148‑1 (MQ=255)
                                                                                                                            cGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTgg                            >  4:1340637/1‑151 (MQ=255)
                                                                                                                                tcgaatcgaaGGTCCGCGCGGCGTCGGTAGTCGTCATGGTCGGCGCAGGCTGGATCGTCGTCACCCTCTGTGAAGCGATTCGCGGTTTCCAGTTCCAACATACCTCCGCCCTGATGGTGCTTATCATCGTCACTTTCAGCCTGCTTGAtt                         <  6:1389712/150‑1 (MQ=255)
                                                                                                                                 cGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGAttt                        <  7:1264712/150‑1 (MQ=255)
                                                                                                                                  gAATCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTcc                      >  6:1935816/1‑151 (MQ=255)
                                                                                                                                  gAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTcg                      >  2:29931/1‑150 (MQ=255)
                                                                                                                                  gAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTcc                      >  3:857514/1‑151 (MQ=255)
                                                                                                                                   aaTCCAACGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCct                     >  1:864406/1‑151 (MQ=255)
                                                                                                                                     tCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCctct                   >  7:987220/1‑151 (MQ=255)
                                                                                                                                         aCGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCa               <  4:1988146/151‑1 (MQ=255)
                                                                                                                                          cGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCaa              <  3:1154976/151‑1 (MQ=255)
                                                                                                                                          cGTCCGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCaa              <  3:1620213/151‑1 (MQ=255)
                                                                                                                                          cGTCCGCTCGGCGACCGGCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCTCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCaa              <  7:155753/151‑1 (MQ=255)
                                                                                                                                          cGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCaa              <  5:578009/151‑1 (MQ=255)
                                                                                                                                             ccGCTCGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACgg           >  6:1271861/1‑151 (MQ=255)
                                                                                                                                               gcgcGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACGGtt         <  4:1899679/148‑1 (MQ=255)
                                                                                                                                                 gcGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACGGTTGc       <  5:1550411/150‑1 (MQ=255)
                                                                                                                                                  cGGCGACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACGGTTGCg      <  6:276791/151‑1 (MQ=255)
                                                                                                                                                  cGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACGGTTg        <  7:804137/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                      
GGTGGAAGCGATTGAGCCCGGCCCGGTGGAAGGCATTCGCGCCACCGGTGCCAACAAGCTCGAAGAGATCCTCTACGGCGTGCTGCCACAGGTGATGCCACTGCTGATCTCCTACTCCCTCTATCGCTTCGAATCCAACGTCCGCTCGGCAACCGTCGTCGGCATGGTCGGCGCAGGCGGGATCGGCGTCACCCTGTGGGAAGCGATTCGCGGTTTCCAGTTCCAACAAACCTGCGCCCTGATGGTGCTTATCATCGTCACGGTCAGCCTGCTGGATTTCCTCTCTCAACGGTTGCGTAAG  >  DH1/4191826‑4192126

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: