Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,9920C.34.4% 284.7 / 149.7 122pseudogene (272/504 nt)insB1IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base . (23/19);  total (67/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 127 total aligned reads displayed.

AACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGG  >  NC_000913/4518854‑4519129
                                                                                                                                          |                                                                                                                                         
 aCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                          >  2:139914/1‑139 (MQ=255)
           tACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAgg                                                                                                                                <  2:884470/139‑1 (MQ=255)
              gTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                             <  2:1503786/139‑1 (MQ=255)
                 ggTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTgg                                                                                                                          >  1:1525017/1‑139 (MQ=255)
                  gTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                          >  1:157387/1‑122 (MQ=255)
                  gTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                          <  2:157387/122‑1 (MQ=255)
                  gTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGt                                                                                                                         >  2:1327718/1‑139 (MQ=255)
                  gTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGt                                                                                                                         >  2:417787/1‑139 (MQ=255)
                       aaaTCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                    <  2:333978/139‑3 (MQ=255)
                               tCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                              <  1:198315/121‑1 (MQ=255)
                               tCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                              >  2:198315/1‑121 (MQ=255)
                                  gcgcTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                           >  1:1789384/1‑105 (MQ=255)
                                  gcgcTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                    <  2:1475727/128‑3 (MQ=255)
                                  gcgcTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                    >  1:1475727/1‑126 (MQ=255)
                                      tggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                     <  1:1327718/139‑1 (MQ=255)
                                       ggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                    <  2:825523/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGGTACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGGCTGATGGCTGGCCg                                                                                                   <  1:804841/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                             >  2:250152/1‑113 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                             <  1:250152/113‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                   <  2:1385737/139‑1 (MQ=255)
                                                  gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                           >  1:1692448/1‑89 (MQ=255)
                                                  gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                           <  2:1692448/89‑1 (MQ=255)
                                                      aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                    >  1:559476/1‑106 (MQ=255)
                                                                  ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGa                                                                                                             >  2:1426072/1‑103 (MQ=255)
                                                                  ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGa                                                                                                             <  1:1426072/103‑1 (MQ=255)
                                                                  ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTg                                                                                 >  1:1127529/1‑131 (MQ=255)
                                                                  ggaggaCGGATGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGCTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTg                                                                                 <  2:1127529/131‑1 (MQ=255)
                                                                  ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                    >  1:258125/5‑111 (MQ=11)
                                                                       aCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                   >  1:921273/1‑139 (MQ=14)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTg                                                                                                                           >  2:555879/1‑84 (MQ=255)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTg                                                                                                                           <  1:555879/84‑1 (MQ=255)
                                                                         ggTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTa                                                                                              <  1:1256780/110‑1 (MQ=14)
                                                                           ttGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                   <  2:1020048/104‑1 (MQ=255)
                                                                           ttGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                   >  1:1020048/1‑104 (MQ=255)
                                                                              ttgcgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATc                                                            >  1:884924/3‑139 (MQ=14)
                                                                              tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTAt                                                             >  2:178027/1‑139 (MQ=255)
                                                                              tgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTAt                                                             >  2:1044022/1‑139 (MQ=255)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                        <  2:1469127/92‑1 (MQ=17)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                        >  1:1469127/1‑92 (MQ=14)
                                                                                   cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                       >  1:1590311/1‑108 (MQ=255)
                                                                                   cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                       <  2:1590311/108‑1 (MQ=255)
                                                                                   cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCa                                                        >  1:1012315/1‑139 (MQ=255)
                                                                                    aCGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa                                                       >  1:1078298/1‑139 (MQ=255)
                                                                                     cgtattCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTc                                                                                                                              >  1:496511/5‑66 (MQ=1)
                                                                                     cgtattCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTc                                                                                                                              <  2:496511/62‑1 (MQ=1)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                                                                                                               <  1:1521303/66‑1 (MQ=255)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                                                                                                               >  2:1521303/1‑66 (MQ=255)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                   <  2:122659/94‑1 (MQ=255)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGtt                                                               <  2:548984/130‑1 (MQ=255)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg                                                      >  1:761214/1‑139 (MQ=255)
                                                                                     cGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCg                                                                                                                             >  2:208854/1‑67 (MQ=11)
                                                                                     cGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                    >  2:1727848/1‑92 (MQ=17)
                                                                                         ttCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                        <  1:740158/84‑1 (MQ=17)
                                                                                         ttCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                        >  2:740158/1‑84 (MQ=17)
                                                                                             gTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTacac                                              <  2:1012315/139‑1 (MQ=255)
                                                                                                aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggc                                                                                                         <  2:481976/76‑1 (MQ=18)
                                                                                                aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggc                                                                                                         >  1:481976/1‑76 (MQ=18)
                                                                                                  cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgc                                         >  2:146944/1‑139 (MQ=255)
                                                                                                  cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCGGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgc                                         >  1:185977/1‑139 (MQ=255)
                                                                                                       tctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCAtt                                    >  2:330289/1‑139 (MQ=255)
                                                                                                       tctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCAtt                                    >  1:56546/1‑139 (MQ=255)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTg                                   >  1:290610/1‑139 (MQ=255)
                                                                                                        atGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTg                                                                     <  1:369900/103‑1 (MQ=18)
                                                                                                         tGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                  <  2:761214/139‑1 (MQ=255)
                                                                                                          ggcgacgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                     >  2:597660/8‑70 (MQ=21)
                                                                                                          ggcgacgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                     <  1:597660/63‑1 (MQ=21)
                                                                                                          ggcgacgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                 <  2:407767/131‑1 (MQ=11)
                                                                                                          ggCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                 >  1:407767/1‑138 (MQ=11)
                                                                                                            cgacgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCgg                              >  1:359819/6‑138 (MQ=11)
                                                                                                               cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                            <  1:417787/139‑1 (MQ=255)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                             >  2:318810/1‑40 (MQ=255)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTg                                                                                                                           >  1:759431/1‑41 (MQ=255)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                       >  1:1128529/1‑109 (MQ=18)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCaa                                                       <  2:1128529/109‑1 (MQ=18)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                  <  1:136617/130‑1 (MQ=14)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataa                         >  1:1106433/1‑133 (MQ=11)
                                                                                                                   ggAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaat                        <  1:146944/139‑1 (MQ=255)
                                                                                                                   ggAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                        <  1:1691200/58‑1 (MQ=255)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                      >  1:244993/1‑139 (MQ=255)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGTATaacctg                     >  1:1762587/1‑139 (MQ=11)
                                                                                                                         tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGCAGGGAAAGCTGCACGTTATCAGCCAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat                  >  1:531386/1‑139 (MQ=255)
                                                                                                                         tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat                  <  2:244993/139‑1 (MQ=255)
                                                                                                                          cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTc                                            >  2:742543/1‑112 (MQ=255)
                                                                                                                          cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTc                                            <  1:742543/112‑1 (MQ=255)
                                                                                                                          cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatc                 >  2:344542/1‑139 (MQ=255)
                                                                                                                                aGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGac           <  1:139914/139‑1 (MQ=255)
                                                                                                                                    tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca       >  1:12376/1‑139 (MQ=255)
                                                                                                                                    tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca       >  1:1671985/1‑139 (MQ=255)
                                                                                                                                    tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca       >  2:397652/1‑139 (MQ=255)
                                                                                                                                    tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGacaaca       <  2:1351862/139‑1 (MQ=255)
                                                                                                                                    tgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                           <  1:1015048/70‑1 (MQ=21)
                                                                                                                                    tgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                           >  2:1015048/1‑70 (MQ=21)
                                                                                                                                    tgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAgcgaa                                      <  1:1466624/107‑3 (MQ=11)
                                                                                                                                     gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                     <  1:1308760/44‑1 (MQ=255)
                                                                                                                                     gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                     >  2:1308760/1‑44 (MQ=255)
                                                                                                                                     gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                            >  2:1443558/1‑117 (MQ=255)
                                                                                                                                      ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTg                                                                     <  2:1644924/75‑1 (MQ=255)
                                                                                                                                       tgTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                     >  2:724842/1‑42 (MQ=255)
                                                                                                                                       tgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCt    >  1:1157378/1‑139 (MQ=255)
                                                                                                                                         tCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTgg  <  2:1157378/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGG  >  NC_000913/4518854‑4519129

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