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

AACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAT  >  NC_000913/4518854‑4519125
                                                                                                                                          |                                                                                                                                     
aaCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                       <  2:136165/139‑1 (MQ=255)
 aCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                      >  2:1016722/1‑139 (MQ=255)
  cAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCgg                                                                                                                                     >  1:817314/1‑139 (MQ=255)
       gggCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCttt                                                                                                                                >  2:767955/1‑139 (MQ=255)
        ggCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCttga                                                                                                                               <  1:825229/139‑3 (MQ=255)
             cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  2:33903/1‑139 (MQ=255)
             cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:906386/139‑1 (MQ=255)
              gTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                         <  2:627785/139‑1 (MQ=255)
                             cGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACgg                                                                                                          >  1:479989/1‑139 (MQ=255)
                               tCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAggtggt                                                                                                                        >  1:747003/1‑123 (MQ=255)
                               tCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                                                                                                                        <  2:747003/123‑1 (MQ=255)
                                 aGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtgg                                                                                                      <  1:33903/139‑1 (MQ=255)
                                       ggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                <  1:236989/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               <  2:684881/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               <  1:767955/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               >  1:709443/1‑139 (MQ=255)
                                                  gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  1:497483/1‑102 (MQ=255)
                                                  gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:497483/102‑1 (MQ=255)
                                                      aTGACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc                                                                                >  2:546014/1‑139 (MQ=1)
                                                             gctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                         >  1:924793/5‑139 (MQ=11)
                                                             gctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                         <  1:945116/135‑1 (MQ=11)
                                                               tccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                       <  1:394742/137‑1 (MQ=11)
                                                               tACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAggg                                                                        >  1:925281/1‑139 (MQ=255)
                                                                aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaat                                                                                                               <  2:6455/99‑4 (MQ=255)
                                                                aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaat                                                                                                               >  1:6455/1‑96 (MQ=255)
                                                                aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                       <  2:709443/139‑1 (MQ=255)
                                                                  ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                    <  2:311946/139‑1 (MQ=14)
                                                                   gaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGc                                                                   <  1:546014/136‑1 (MQ=11)
                                                                       aCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                               >  1:587382/1‑139 (MQ=14)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  1:239380/81‑1 (MQ=255)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  2:239380/1‑81 (MQ=255)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGTCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                >  1:524361/1‑139 (MQ=255)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                >  1:1409/1‑139 (MQ=255)
                                                                           ttgttgCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               <  2:101395/103‑1 (MQ=14)
                                                                           ttgttgCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               >  1:101395/1‑103 (MQ=14)
                                                                           ttgttgCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTGGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCCGGAGGGAAAGCTGCACGta                                                           >  1:524889/1‑138 (MQ=14)
                                                                           ttGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGt                                                            >  2:318056/1‑139 (MQ=255)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                     <  2:924793/139‑1 (MQ=14)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                     <  2:625030/139‑1 (MQ=14)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAg                                                  >  1:273691/1‑139 (MQ=255)
                                                                                                  cGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAgc                                     >  2:378856/1‑139 (MQ=255)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTc                                                                                                                          >  2:395510/1‑53 (MQ=14)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTc                                                                                                                          <  1:395510/53‑1 (MQ=14)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGt                                                                                                                     >  1:513338/1‑58 (MQ=17)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGt                                                                                                                     <  2:513338/58‑1 (MQ=17)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                   >  2:342056/1‑92 (MQ=17)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCa                                                                                   <  1:342056/92‑1 (MQ=18)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                 >  1:864559/1‑94 (MQ=17)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                 <  2:864559/94‑1 (MQ=17)
                                                                                                  cGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAgcg                                    >  2:109081/1‑139 (MQ=11)
                                                                                                       tatGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                    <  2:209692/68‑1 (MQ=255)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                >  1:214375/1‑74 (MQ=255)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGcc                                                                                                <  2:214375/74‑1 (MQ=255)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTg                               >  2:951296/1‑139 (MQ=255)
                                                                                                          ggCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAg                             >  2:985606/1‑139 (MQ=255)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTc                                                                                                                          >  2:883640/3‑40 (MQ=37)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTc                                                                                                                          <  1:883640/38‑1 (MQ=37)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               >  2:43441/3‑67 (MQ=21)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               <  1:43441/65‑1 (MQ=21)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtat                       >  2:774838/3‑135 (MQ=11)
                                                                                                               cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          <  2:996487/41‑1 (MQ=255)
                                                                                                               cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                          >  1:996487/1‑41 (MQ=255)
                                                                                                               cacTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACa                        <  2:925281/139‑1 (MQ=255)
                                                                                                                   ggAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATatacgc                                        <  1:1084177/118‑6 (MQ=17)
                                                                                                                   ggAGCGTCTTCTGAGCCTGCTGTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATatacgc                                        >  2:1084177/1‑116 (MQ=14)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                            >  2:452317/1‑129 (MQ=255)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                            <  1:452317/129‑1 (MQ=255)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                  <  1:242464/139‑1 (MQ=255)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                  >  1:251854/1‑139 (MQ=255)
                                                                                                                         tcttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaat              <  2:985194/139‑1 (MQ=255)
                                                                                                                            tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                            <  1:979617/122‑1 (MQ=255)
                                                                                                                            tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                            >  2:979617/1‑122 (MQ=255)
                                                                                                                                aGCCTGCTGTGGGCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 <  2:143565/48‑1 (MQ=25)
                                                                                                                                aGCCTGCTGTGGGCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                 >  1:143565/1‑48 (MQ=25)
                                                                                                                                   ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaa                                                                      >  1:35387/1‑73 (MQ=255)
                                                                                                                                   ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaa                                                                      <  2:35387/73‑1 (MQ=255)
                                                                                                                                    tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTa                                                          <  1:710626/84‑1 (MQ=255)
                                                                                                                                    tgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTa                                                          >  2:710626/1‑84 (MQ=255)
                                                                                                                                    tgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAg                                                      >  2:925073/1‑87 (MQ=17)
                                                                                                                                    tgctgTGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAg                                                      <  1:925073/87‑1 (MQ=17)
                                                                                                                                     gctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAt  >  2:805093/1‑139 (MQ=255)
                                                                                                                                      ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               <  2:182662/45‑1 (MQ=255)
                                                                                                                                      ctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                               >  1:182662/1‑45 (MQ=255)
                                                                                                                                        gTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                  <  2:124657/40‑1 (MQ=38)
                                                                                                                                        gTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                  >  1:124657/1‑40 (MQ=255)
                                                                                                                                          cGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                              >  1:272489/1‑106 (MQ=255)
                                                                                                                                          cGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                              <  2:272489/106‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                     
AACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACAT  >  NC_000913/4518854‑4519125

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