Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,1840AG33.3% 107.3 / 10.6 75E171G (GAG→GGG) mdfAmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (20/30);  new base G (19/6);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 3.75e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCC  >  NC_000913/883947‑884400
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                           
ttATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCGCCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGa                                                                                                                                                                                                                             <  1:1599857/238‑1 (MQ=255)
ttATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCGCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAg                                                                                                                                                                                                                            <  1:993318/239‑1 (MQ=255)
        aCCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGAggtt                                                                                                                                                                                                                    >  1:1681482/1‑239 (MQ=255)
        aCCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAtgtt                                                                                                                                                                                                                    >  1:1356989/1‑239 (MQ=255)
                    atatTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCCTGTGCTGGCCTGGGGGGGGGGGTTTGGtttgtttg                                                                                                                                                                                                        >  1:1419376/1‑238 (MQ=255)
                    atatTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGc                                                                                                                                                                                                       <  1:236594/239‑1 (MQ=255)
                       ttGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCCAGATCACCGCGCTGATGGCGAACGTGGCGCTGGTTGCTCCGCTACTTGGTCCGCTGGTGGGCGAGGCGCGGAGCCATGTTCTGCCCTGGGGGGGGagggt                                                                                                                                                                                                                    >  2:1224427/1‑220 (MQ=255)
                         gCTGGCGCAAAACATTGAACAATTCGCCCTTTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCAt                                                                                                                                                                                                  <  1:335721/239‑1 (MQ=255)
                            ggCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTgg                                                                                                                                                                                              <  1:992225/239‑1 (MQ=255)
                                   aaCATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGAt                                                                                                                                                                                       <  1:2060824/239‑1 (MQ=255)
                                    aCATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGAtc                                                                                                                                                                                      <  1:1276809/239‑1 (MQ=255)
                                        tGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCCTGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCa                                                                                                                                                                                            >  1:583988/1‑229 (MQ=255)
                                           aCAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCtttt                                                                                                                                                                               <  2:1682162/239‑1 (MQ=255)
                                              aTTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCg                                                                                                                                                                            <  1:2248125/239‑1 (MQ=255)
                                               ttCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAAAGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGGGCGGCCCTGGGGGGGGATGTTTGTTTTGTTTGGCGCATGTGGCGCGAGCTCCtttt                                                                                                                                                                               >  2:1960769/1‑235 (MQ=255)
                                                tCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGt                                                                                                                                                                          <  1:351218/239‑1 (MQ=255)
                                                       gTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCCAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGGGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTGGCGCATTGGGAGCGGTCTCCTTTTTCGGTCTGCAAc                                                                                                                                                                  >  2:2090010/1‑239 (MQ=255)
                                                       gTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGACCGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGGGCGGCGCGGAGCCATGTGCTGCCCTGGGAGGGGGTGATTGTTTTGTTTGACGCATTGGCAGCGAGTTCCTTTTTCGGTGTGCacc                                                                                                                                                                  >  2:5089/1‑237 (MQ=255)
                                                            gcTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGcc                                                                                                                                                             <  1:1619252/239‑1 (MQ=255)
                                                                cttGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTTTTCAGATCACCGCGCTGATGGCGAAAGTGGCGCCGATTGCTCCGCTACTTGGTCCGCCGGTGGGGGCGGCGTGGATCCATGTGCTGCCCTTGGGGGGGATGTTTGTTTTTTTTGCCGCATTGGCTGCGATCTCCTTTTTCGGTCTGCAACCAgacctc                                                                                                                                                          >  2:1377976/1‑233 (MQ=255)
                                                                  tGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCTAGTTCACCGCGCTGATGGCGAAAGTGGCGCTGATTTCTTCCCTAATTGGGCCCATGGTGGGCGCGGGCTGGATCCAAGTGGTGGCCTCGGGGGGGTAGTTTGTTTTGTTTGCCGCATAGGCAACGCGTTCCCTTTTTGGTCTTCaaagtaccaag                                                                                                                                                          >  2:1731521/1‑227 (MQ=255)
                                                                          aTAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGa                                                                                                                                                                                                                             >  2:750023/1‑164 (MQ=255)
                                                                          aTAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGa                                                                                                                                                                                                                             <  1:750023/164‑1 (MQ=255)
                                                                                   tGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATTAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTTCCGCATTGGCAGCGAGCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGcccaccgtat                                                                                                                                      >  2:1725395/1‑231 (MQ=255)
                                                                                   tGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTTTCAAGATCCCCGCGCCGATGGCGAACGGGGCGCCGATTGCTCCGCTACATGGGCCGCTTGTGGGGGCGGCGTGGGACCATGTGTTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCTTTGGCAGGGATCTCCTTTTTTGGCTTGCAACTAACCCTGCCTGAAACCGccccccgctc                                                                                                                                      >  2:354112/1‑231 (MQ=255)
                                                                                   tGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGTTTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGGGCTGCCCTGGGGGGGGATGGTTGTTTTGTTTGCCGCCTTGGCAGCGATCTCCTTTTTCGGTCTGCAACTAGCCATTCCTTGAACCGcccaccgtat                                                                                                                                      >  2:1495873/1‑231 (MQ=255)
                                                                                        cATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATCAAAATCACCGCGCTGATGGCGAACGCGGCGCCGATTGCTCCGCTTCTTGGGCCGCTTGTGGGCGCGGCCTGGATCCATTTGCTGCCCTTGGGGGGGATGTTTTGTTTTTTTTCCCCATTGGCCGGCACCTCCCTTTTCGGGCTGCCACCGGCCCAGTCCTGAAACCCCCCGCGtattaggc                                                                                                                                 >  2:2081821/1‑234 (MQ=255)
                                                                                                ttgtgGTATCCGCTGCAATCAAGGATTCCTCGGAAGAGGCGTTTTGTATAAAGATCACGGCGCGGAGGGCGAGCGTGGCGCTGTTTCCTCCGCTTTTTGGTCCGCTGGGGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTGTTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCGGCAACGGGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACt                                                                                                                          <  2:1356989/237‑1 (MQ=255)
                                                                                                           gccgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGAACACCGCCCTGGAGGCGAAACCGGCGCTGATTGCCCCGCCACTTTGTCCCCTGGGGGGCGCGGCGGGGAGCCATTGCCTGCCCTGGGGGGGGGTGCTTTTTTTTGTTGCCGCCTTGGCCGCCGACTCCTTTTTTGGGCTGCCACGAGCCAACCCGGAAAACCCCCCGCGTACATGCGGGCAACAGtacacggaaga                                                                                                              >  2:1919067/1‑229 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATCAAGACCCCCGCGCCGGTGGCGAACGGGGCGCCTGTTGCTCCGCTACTTGGTCCGCCGGGGGGCGCGGGGGGGATCCATGTGCTGCCCTGGGGGGGGGGGTTTGTTTTGTTTTCCGCATTGGTAGCGGTCTCCTTTTTTGGGCTGCACCCAGCCATTCCTGAAACCGCCCCCCGTATAGTCGGGGAACAGttctctaaaaaa                                                                                                             >  2:874830/1‑227 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGaaaaa                                                                                                              >  1:1609011/1‑235 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGGGGCGCTGATTGCCCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGGTTGCCGCATTGGGAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGaa                                                                                                                             >  1:379305/1‑222 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCCCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTGCCGCCTTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCCGAAACCGCCCCCCGTATAGGCGAGAAAATGTCACTGAAAGaaa                                                                                                            >  2:2237855/1‑238 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAAATCACCGCGCCGATGGCGAACGTGGCGCTGATTGCTCCGCGACTTGGGCCGCCGGGGGGCGCCGCGCGGGCTCCTGTGCTGCCCTGGGGGGGGGTGGTTGTTTTGTTTGCCCCAATTGGAGCCATCTTCTTTTTCGGTCTGTCACGAGCCATGCCTGTAACCGCCCCGCGCATAGGCGCGAAACAGTcgacgtaaga                                                                                                              >  2:121284/1‑228 (MQ=255)
                                                                                                             ccaatttCAGGAACCCTTGGAAGGGGCGGTTTGTATCAAGATCACCGCGCGGTTGGGGAACGGGGCGCTGATTGCCCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCa                                                                                                                      <  2:857740/224‑1 (MQ=255)
                                                                                                               cAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAAAGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGGGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGGTGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTTGGGCTGCAACGGGCCATGCCTGGAACGGCCACCCGCATAGGCGGGAAACAGTTCCCGAAAAAAccc                                                                                                          >  2:922987/1‑237 (MQ=255)
                                                                                                                       gAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGGTCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTTTCACTGAAAGAAc                                                                                                            >  1:857740/1‑229 (MQ=255)
                                                                                                                        aacccTGCGGAGAGGCGGTTTGTATCAATATCCCGGCGCTGAGGGCGAACGTGGCGCTGATGGTGCGGCTACTTGGTCCGCTGGTGGACGCGGCGTGGATCCAGGTGCTGCCCTGGGGGGGGATTTTTGTTTTGTTTGCCGCAATGGCAGCGTTCTCTTTTTTCGGTTTGCAACGACCAAGGCCGGAAACCGCCACGCGTATAGGCGAGAATCTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:1558270/233‑1 (MQ=255)
                                                                                                                        aacccTGCGAAGAGGCGGTTTGTTTCAAGACCCCCGCGCTGAGGGCGACCGGGGCGTGGATTGCTCGGCTACTGGGTCCGCTGGTGGGCGGGGCGTGGATCCATGTGCTGCCCGGGGAGGGGATGTTTTTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:1984240/233‑1 (MQ=255)
                                                                                                                        aaTCTTGCGAGGAGGCGGTTGGTATCAAGATCCCCGGGCGGTGGGGGAACGGGGCGCTGATTGCCCCGCTACTTGGCCGGCTGGTGGGCGGGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATTTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCg                                                                                                     <  2:1003909/235‑1 (MQ=255)
                                                                                                                        aaTCCTTCAAAGAGGCGGTTTTGATCAAGACACCCGCGCTGATGGCGGCGGGGGCACTGATTGCGCCGCTACTGGGTCCGGTGGTGGGCGCGGCGTGGATCCATTGGCTCCCCTGGGAGGGAATTTTGTTGTTGTGTGCGGCTTGGGCAGCGATCTCCTTTTTCGGTCGGCAACGAGCCGTGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:1614255/236‑1 (MQ=255)
                                                                                                                            cTTCGAAGAGGCGGTTTGTATCAAGATCCCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACt                                                                                                                          <  2:1482424/210‑1 (MQ=255)
                                                                                                                            cTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACtata                                                                                             >  1:1984240/1‑239 (MQ=255)
                                                                                                                            cTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACt                                                                                                                          >  1:1482424/1‑210 (MQ=255)
                                                                                                                               cGAAGAGGCTGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGc                                                                                          <  1:1960769/239‑1 (MQ=255)
                                                                                                                               cGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGGCGTGACTATAAGc                                                                                          >  1:1670223/1‑239 (MQ=255)
                                                                                                                                 aaGAGGCGGTTTTTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAAATGTCACTGAAAGAAATCGGTCGGGACTATaa                                                                                            >  1:476430/1‑235 (MQ=255)
                                                                                                                                      gCGGTTTTTATCAAGATCACCGCGCTGATGGCGGACGTTGCGCTGATTGCTCCGCTACTTGGTACGCTGGGGGGCGCGGCGTGGATCCATGTTCTGCCCTGGGGGGGGGTGTTTGTGTTGTTTGCCGCGTTGGCAGCGATATCCTTTTTTGGTCTGCAACGCGCCATGCCTGAAACCGCCACGCGTATATGCGGGAATCTGTGCCTGAAAGAATCCGGGGGTGCCTATAAGCTggggtc                                                                                   >  2:949042/1‑235 (MQ=255)
                                                                                                                                      gCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGGGACTATAAGCTGGt                                                                                      >  1:126200/1‑236 (MQ=255)
                                                                                                                                        ggTTTGTATCAAGATCACCGCGCTGATGGCGACCGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCCTGTGCTGCCCTGGGGGGGGGTGGTTGTTTTGTTTGCCGCCTTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCCTGCCCGAAAACGCCACGCGTATAGGCGAGAAAATGTCCCTGAAAGAACTCGGGTGTGAATACAAACTgggggtgg                                                                                 >  2:493487/1‑233 (MQ=255)
                                                                                                                                               aTCAAGATCCCCGCGCTGAGGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACgg                                                                          <  1:2237855/239‑1 (MQ=255)
                                                                                                                                                       ccccGCGCTGAGGGCACCCGGGGCGCTGTTGGCTCCGCTCCTGGGTCCGCTGGTGGGCGCGGCGTGGATCCTTGTGTTGCCCGGGGGGGGGTGTTTTGTTTTGTTTGCCGCTTGGGCGGCGATCTCCTTTTTCGGTCGGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCttt                                                                   <  2:1419376/236‑1 (MQ=255)
                                                                                                                                                        aCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCaa                                                                                                                                                                   >  2:1400108/1‑141 (MQ=255)
                                                                                                                                                        aCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCaa                                                                                                                                                                   <  1:1400108/141‑1 (MQ=255)
                                                                                                                                                        aCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCATGGATCCATTAGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCGGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGGATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTtgt                                                                 >  2:1310138/1‑239 (MQ=255)
                                                                                                                                                                    gggacggTGGCGCTGTTTGCTCCGCTCCTGGGTCGGCTGGTGGGCGCGGGGTGGATCCTTGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:534926/186‑1 (MQ=255)
                                                                                                                                                                           tGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGTTCCATGTGCTGCCCTGGGAGGGGTTGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCATGCGTCTAGGCGAGaa                                                                                                                             >  1:727936/1‑160 (MQ=255)
                                                                                                                                                                              cgcTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGc                                                                                          >  1:534926/1‑192 (MQ=255)
                                                                                                                                                                                   tttGCCCCGCTACTGGGTCCGCTGGTGGGCGCGGCGTGGATCCAGGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGgcgctggcgctg                                              <  1:502838/230‑1 (MQ=255)
                                                                                                                                                                                             tACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGGGCGGGGGTTGGtgt                                    >  2:502838/1‑229 (MQ=255)
                                                                                                                                                                                                  ggTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTg                                                                                        <  1:46205/174‑1 (MQ=255)
                                                                                                                                                                                                  ggTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTg                                                                                        >  2:46205/1‑174 (MQ=255)
                                                                                                                                                                                                    tCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATTTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTg                      <  1:1698786/238‑1 (MQ=255)
                                                                                                                                                                                                       ggtggtggGCGGGGGGGGGATCCAGGTGCGGCCGTGGGAGGGGATGTTTGTTTTGTGTGCCGCTTTGGCAGCGATCTCCTTTTTCGGTCGGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGt                  <  2:227940/237‑1 (MQ=255)
                                                                                                                                                                                                        cTGGTGGGGGGGGTGGGGTTCCAGGTGCTGCCCGGGGAGGGGGTGTTTGTTTTGTTGGCCGCATTGGCAGGGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGcc                                                                                                                                              <  2:1441711/114‑1 (MQ=255)
                                                                                                                                                                                                         tggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACtat                                                                                              <  1:763980/161‑1 (MQ=255)
                                                                                                                                                                                                         tggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACtat                                                                                              >  2:763980/1‑161 (MQ=255)
                                                                                                                                                                                                                 gcggtgTGGATCCGTGTGCTGCCCTGGGAGTGTATGTTTGTTTTGTTTGCCGCATTGGCGGCGTCCTCCTTTTCCGTTCTGAAGCGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCa        <  1:1087952/239‑1 (MQ=255)
                                                                                                                                                                                                                  cggcgTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAg                                                                                                                                    >  1:1441711/1‑114 (MQ=255)
                                                                                                                                                                                                                     cggggATCATTGTGCTGCCTGGGGAGGGTAGTTTTGTTTTTTTGGCCGCATGGGCAGCGACCTCCTTTTTCGGTCTGCAACGACCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgc                                                      <  2:1139708/186‑1 (MQ=255)
                                                                                                                                                                                                                        ggATCCTTGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTGGCCGATTGGGCAGGGATCTCTTTTTCCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGcc  <  1:2210810/238‑1 (MQ=255)
                                                                                                                                                                                                                                    tgcCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGcc                                                                                                                                              <  2:1075732/86‑1 (MQ=255)
                                                                                                                                                                                                                                    tgcCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGcc                                                                                                                                              >  1:1075732/1‑86 (MQ=255)
                                                                                                                                                                                                                                      cccTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCg                                                                                                                                                                            >  1:1468559/1‑55 (MQ=255)
                                                                                                                                                                                                                                      cccTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCg                                                                                                                                                                            <  2:1468559/55‑1 (MQ=255)
                                                                                                                                                                                                                                      cccTGGGAGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTg                         <  1:1523464/200‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGGGAGGGGATGTTTTTTTTGTTGGCCCCATTGGCAGCGATCTCCTTTTTCGGGCTGCAACGAGCCATGCCTAAAACCGcc                                                                                                                                              <  2:2086851/82‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAGGGGATTTTTGTTTTGTTTCCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGt                           <  1:556670/194‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                           
TTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCC  >  NC_000913/883947‑884400

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