Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,518,9920C.34.2% 188.9 / 95.3 79pseudogene (272/504 nt)insB1IS1 transposase B,IS, phage, Tn, Transposon‑related functions, extrachromosomal, transposon related
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base . (16/11);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 6.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGC  >  NC_000913/4518854‑4519130
                                                                                                                                          |                                                                                                                                          
aaCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTc                                                                                                                                            <  1:1018591/139‑1 (MQ=255)
 aCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCg                                                                                                                                           <  2:227071/139‑1 (MQ=255)
      ggggCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCtt                                                                                                                                      >  2:274011/1‑139 (MQ=255)
       gggCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCttt                                                                                                                                     <  1:858037/139‑1 (MQ=255)
            aCGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                                                                                                                <  1:335706/139‑1 (MQ=255)
             cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                                                                                                                <  2:732032/138‑1 (MQ=255)
             cGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                                                                                                                >  1:732032/1‑138 (MQ=255)
              gTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                              <  1:670276/139‑1 (MQ=255)
              gTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                              <  2:1052413/139‑1 (MQ=255)
               tCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                                                                                                                             >  2:846835/1‑139 (MQ=255)
                    gCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtat                                                                                                                        <  1:731009/139‑1 (MQ=255)
                       aaaTCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGtataaa                                                                                                                     <  1:814490/139‑3 (MQ=255)
                              gTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                              <  1:263387/123‑1 (MQ=255)
                              gTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                              >  2:263387/1‑123 (MQ=255)
                              gTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTACAAATGACGGa                                                                                                              <  2:707309/139‑1 (MQ=255)
                                   cgcTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggct                                                                                                         <  1:586090/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                    >  1:287394/1‑139 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                    <  1:274011/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                    <  1:875128/139‑1 (MQ=255)
                                        gctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCg                                                                                                    <  1:846835/139‑1 (MQ=255)
                                             tttACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTa                                                                                               <  1:454230/139‑1 (MQ=255)
                                                      aTGACAGGCTCCGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGc                                                                                     >  1:218951/1‑139 (MQ=11)
                                                      aTGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACg                                                                                      >  2:597301/1‑139 (MQ=255)
                                                        gACAGGCTCCGGAAGACGGTCGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCt                                                                                   >  1:181639/1‑139 (MQ=11)
                                                            ggATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGt                                                                                                                                <  1:679358/91‑1 (MQ=255)
                                                            ggAAACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCAGTCGGCCTTTGAGGt                                                                                                                                >  2:679358/1‑91 (MQ=255)
                                                             gctccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAgg                                                                              >  1:341556/5‑139 (MQ=11)
                                                             gATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                              <  1:607295/92‑1 (MQ=255)
                                                             gATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCg                                                                                                                              >  2:607295/1‑92 (MQ=255)
                                                              aTACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                       >  1:594956/1‑114 (MQ=255)
                                                              aTACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGAtggctgg                                                                                                       <  2:594956/114‑1 (MQ=255)
                                                                ccGGAAGACGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGaa                                                                           >  2:997019/2‑139 (MQ=11)
                                                                aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGa                                                                            >  2:460156/1‑139 (MQ=255)
                                                                  ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                                                                                                                             >  2:938432/1‑88 (MQ=255)
                                                                  ggaggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                                                                                                                             <  1:938432/88‑1 (MQ=255)
                                                                  ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                         >  1:465987/5‑139 (MQ=14)
                                                                  ggaagaCGGTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAg                                                                         >  1:888494/1‑139 (MQ=14)
                                                                    aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGc                                                                        >  1:965920/1‑139 (MQ=255)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                     >  1:583896/1‑139 (MQ=255)
                                                                       aCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGc                                                                     >  1:117343/1‑139 (MQ=255)
                                                                         ggTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACg                                                                  >  1:445884/1‑139 (MQ=14)
                                                                         ggTTGTTGCGCACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACg                                                                  >  2:987271/1‑139 (MQ=14)
                                                                             gtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTa                                                               >  2:369210/1‑139 (MQ=255)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                          <  1:973450/139‑1 (MQ=14)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                          >  2:973450/1‑139 (MQ=14)
                                                                                 cgcACGTATTCGGTGAACGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                          <  2:181639/139‑1 (MQ=14)
                                                                                  gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                          <  2:965920/139‑1 (MQ=255)
                                                                                   cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCa                                                         >  2:594948/1‑139 (MQ=255)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                               <  2:597273/67‑1 (MQ=255)
                                                                                     cGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTc                                                                                                                               >  1:597273/1‑67 (MQ=255)
                                                                                             gTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTacac                                               <  1:460156/139‑1 (MQ=255)
                                                                                                aaCGCACTATGGCGACGCTGGAGCGTCTTCTGCGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                      >  2:628590/1‑80 (MQ=18)
                                                                                                aaCGCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGAtggctggc                                                                                                      <  1:628590/80‑1 (MQ=18)
                                                                                                   gCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGc                                                      <  2:133321/125‑1 (MQ=14)
                                                                                                   gCACTATGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGc                                                      >  1:133321/1‑125 (MQ=17)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                                                                                                                             >  2:147256/1‑50 (MQ=255)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGt                                                                                                                             <  1:147256/50‑1 (MQ=255)
                                                                                                        ctGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTg                                    >  2:992542/1‑139 (MQ=255)
                                                                                                        atGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGa                                   <  1:891650/138‑1 (MQ=11)
                                                                                                         tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                  <  2:465987/139‑1 (MQ=11)
                                                                                                         tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                  <  1:235115/139‑1 (MQ=11)
                                                                                                         tGGCGACGCTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAg                                  <  1:658400/139‑1 (MQ=11)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                    >  1:32/3‑99 (MQ=21)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCa                                                                    <  2:32/97‑1 (MQ=21)
                                                                                                               cgcTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtat                            >  1:811915/3‑135 (MQ=11)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATa                           >  2:880329/1‑139 (MQ=255)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCg                                                                                                                              <  1:1030833/39‑1 (MQ=255)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCg                                                                                                                              >  2:1030833/1‑39 (MQ=255)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                          >  1:669375/1‑107 (MQ=18)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCggtataa                          >  2:928528/1‑133 (MQ=14)
                                                                                                                 cTGGAGCGTCTTCTGAGCCTGCTGT‑GGCCATTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTAATCAGc                                                          <  2:669375/107‑1 (MQ=17)
                                                                                                                     aGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatct                       >  1:802655/1‑139 (MQ=255)
                                                                                                                         tcttctGAGCCTGCTGTCGACCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAggg                                                                             >  1:526159/1‑81 (MQ=255)
                                                                                                                         tcttctGAGCCTGCTGTCGACCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAggg                                                                             <  2:526159/81‑1 (MQ=255)
                                                                                                                          cttctGAGCCTGCTGTCGGCCTTTGAGGACGTGGTATAAATGACTGATGGCTGGCCGCTGTATGAATCACGCCTGACGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatc                  >  1:564506/1‑139 (MQ=255)
                                                                                                                            tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatctg                >  2:723288/1‑139 (MQ=255)
                                                                                                                            tctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatctg                >  2:192289/1‑139 (MQ=255)
                                                                                                                              tGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTgag              >  1:204826/1‑139 (MQ=255)
                                                                                                                                   ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCa                                                            >  1:921977/1‑88 (MQ=255)
                                                                                                                                   ctgctgTCGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCa                                                            <  2:921977/88‑1 (MQ=255)
                                                                                                                                          cGGCCTTTGAGGTCGTGGTATAAATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGc  >  2:960770/1‑139 (MQ=255)
                                                                                                                                          |                                                                                                                                          
AACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGC  >  NC_000913/4518854‑4519130

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: ‑