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.

TTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCC  >  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)
                    atatTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCCTGTGCTGGCCTGGGGGGGGGGGTTTGGt‑ttgtttg                                                                                                                                                                                                        >  1:1419376/1‑238 (MQ=255)
                    atatTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGc                                                                                                                                                                                                       <  1:236594/239‑1 (MQ=255)
                       ttGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCCAGATCACCGCGCTGATGGCGAACGTGGCGCTGGTTGCTCCGCTACTTGGTCCGCTGGTGGGCGAGGCGCGGAGCCATGTTCTGCCCTGGGGGGGGagggt                                                                                                                                                                                                                    >  2:1224427/1‑220 (MQ=255)
                         gCTGGCGCAAAACATTGAACAATTCGCCCTTTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAt                                                                                                                                                                                                  <  1:335721/239‑1 (MQ=255)
                            ggCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑Tgg                                                                                                                                                                                              <  1:992225/239‑1 (MQ=255)
                                   aaCATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGAt                                                                                                                                                                                       <  1:2060824/239‑1 (MQ=255)
                                    aCATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGAtc                                                                                                                                                                                      <  1:1276809/239‑1 (MQ=255)
                                        tGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCCTGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCa                                                                                                                                                                                            >  1:583988/1‑229 (MQ=255)
                                           aCAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCtttt                                                                                                                                                                               <  2:1682162/239‑1 (MQ=255)
                                              aTTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCg                                                                                                                                                                            <  1:2248125/239‑1 (MQ=255)
                                               ttCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAAAGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGGGCGGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGGCGCATGTGGC‑GCGAGCTCCtttt                                                                                                                                                                               >  2:1960769/1‑235 (MQ=255)
                                                tCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGt                                                                                                                                                                          <  1:351218/239‑1 (MQ=255)
                                                       gTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCCAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGGGCTGCCCTGGGGGGGGGTGTTTGTT‑TTGTTTGGCGCAT‑TGGGAGCGGTCTCCTTTTTCGGTCTG‑CAAc                                                                                                                                                                  >  2:2090010/1‑239 (MQ=255)
                                                       gTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGACCGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGGGCGGCGCGGAGCCATGTGCTGCCCTGGGAGGGGGTGATTGTT‑TTGTTTGACGCAT‑TGGCAGCGAGTTCCTTTTTCGGTGTG‑Cacc                                                                                                                                                                  >  2:5089/1‑237 (MQ=255)
                                                            gcTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGcc                                                                                                                                                             <  1:1619252/239‑1 (MQ=255)
                                                                cttGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTTTTCAGATCACCGCGCTGATGGCGAAAGTGGCGCCGATTGCTCCGCTACTTGGTCCGCCGGTGGGGGCGGCGTGGATCCATGTGCTGCCCTTGGGGGGGATGTTTGTT‑TTTTTTGCCGCAT‑TGGCTGCGATCTCCTTTTTCGGTCTG‑CAACCAgacctc                                                                                                                                                          >  2:1377976/1‑233 (MQ=255)
                                                                  tGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCTAGTTCACCGCGCTGATGGCGAAAGTGGCGCTGATTTCTTCCCTAATTGGGCCCATGGTGGGCGCGGGCTGGATCCAAGTGGTGGCCTCGGGGGGGTAGTTTGTT‑TTGTTTGCCGCAT‑AGGCAACGCGTTCCCTTTTTGGTCTT‑Caaagtaccaag                                                                                                                                                          >  2:1731521/1‑227 (MQ=255)
                                                                          aTAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGa                                                                                                                                                                                                                             >  2:750023/1‑164 (MQ=255)
                                                                          aTAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGa                                                                                                                                                                                                                             <  1:750023/164‑1 (MQ=255)
                                                                                   tGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATTAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTTCCGCAT‑TGGCAGCGAGCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGcccaccgtat                                                                                                                                      >  2:1725395/1‑231 (MQ=255)
                                                                                   tGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTTTCAAGATCCCCGCGCCGATGGCGAACGGGGCGCCGATTGCTCCGCTACATGGGCCGCTTGTGGGGGCGGCGTGGGACCATGTGTTGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGCCGCTT‑TGGCAGGGATCTCCTTTTTTGGCTTG‑CAACTAACCCTGCCTGAAACCGccccccgctc                                                                                                                                      >  2:354112/1‑231 (MQ=255)
                                                                                   tGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGTTTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGGGCTGCCCTGGGGGGGGATGGTTGTT‑TTGTTTGCCGCCT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACTAGCCATTCCTTGAACCGcccaccgtat                                                                                                                                      >  2:1495873/1‑231 (MQ=255)
                                                                                        cATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATCAAAATCACCGCGCTGATGGCGAACGCGGCGCCGATTGCTCCGCTTCTTGGGCCGCTTGTGGGCGCGGCCTGGATCCATTTGCTGCCCTTGGGGGGGATGTTTTGT‑TTTTTTTCCCCAT‑TGGCCGGCACCTCCCTTTTCGGGCTGCCACCGGCCCAGTCCTGAAACC‑CCCCGCGtattaggc                                                                                                                                 >  2:2081821/1‑234 (MQ=255)
                                                                                                ttgtgGTATCCGCTGCAATCAAGGATTCCTCGGAAGAGGCGTTTTGTATAAAGATCACGGCGCGGAGGGCGAGCGTGGCGCTGTTTCCTCCGCTTTTTGGTCCGCTGGGGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TGTTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCGG‑CAACGGGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACt                                                                                                                          <  2:1356989/237‑1 (MQ=255)
                                                                                                           gccgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGAACACCGCCCTGGAGGCGAAACCGGCGCTGATTGCCCCGCCACTTTGTCCCCTGGGGGGCGCGGCGGGGAGCCATTGCCTGCCCTGGGGGGGGGTGCTTTTT‑TTTGTTGCCGCCT‑TGGCCGCCGACTCCTTTTTTGGGCTG‑CCACGAGCCAACCCGGAAAACCCCCCGCGTACATGCGGGCAACAGtacacggaaga                                                                                                              >  2:1919067/1‑229 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATCAAGACCCCCGCGCCGGTGGCGAACGGGGCGCCTGTTGCTCCGCTACTTGGTCCGCCGGGGGGCGCGGGGGGGATCCATGTGCTGCCCTGGGGGGGGGGGTTTGTT‑TTGTTTTCCGCAT‑TGGTAGCGGTCTCCTTTTTTGGGCTG‑CACCCAGCCATTCCTGAAACCGCCCCCCGTATAGTCGGGGAACAGttctctaaaaaa                                                                                                             >  2:874830/1‑227 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGaaaaa                                                                                                              >  1:1609011/1‑235 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGGGGCGCTGATTGCCCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGGTTGCCGCAT‑TGGGAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGaa                                                                                                                             >  1:379305/1‑222 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCCCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGGTGTTTGTT‑TTGTTTGCCGCCT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCCGAAACCGCCCCCCGTATAGGCGAGAAAATGTCACTGAAAGaaa                                                                                                            >  2:2237855/1‑238 (MQ=255)
                                                                                                             cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAAATCACCGCGCCGATGGCGAACGTGGCGCTGATTGCTCCGCGACTTGGGCCGCCGGGGGGCGCCGCGCGGGCTCCTGTGCTGCCCTGGGGGGGGGTGGTTGTT‑TTGTTTGCCCCAA‑TTGGAGCCATCTTCTTTTTCGGTCTG‑TCACGAGCCATGCCTGTAACCGCCCCGCGCATAGGCGCGAAACAGTcgacgtaaga                                                                                                              >  2:121284/1‑228 (MQ=255)
                                                                                                             ccaatttCAGGAACCCTTGGAAGGGGCGGTTTGTATCAAGATCACCGCGCGGTTGGGGAACGGGGCGCTGATTGCCCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCa                                                                                                                      <  2:857740/224‑1 (MQ=255)
                                                                                                               cAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAAAGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGGGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGGTGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTTGGGCTG‑CAACGGGCCATGCCTGGAACGGCCACCCGCATAGGCGGGAAACAGTTCCCGAAAAAAccc                                                                                                          >  2:922987/1‑237 (MQ=255)
                                                                                                                       gAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGGTCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTTTCACTGAAAGAAc                                                                                                            >  1:857740/1‑229 (MQ=255)
                                                                                                                        aacccTGCGGAGAGGCGGTTTGTATCAATATCCCGGCGCTGAGGGCGAACGTGGCGCTGATGGTGCGGCTACTTGGTCCGCTGGTGGACGCGGCGTGGATCCAGGTGCTGCCCTGGGGGGGGATTTTTGTT‑TTGTTTGCCGCAA‑TGGCAGCGTTCTCTTTTTTCGGTTTG‑CAACGACCAAGGCCGGAAACCGCCACGCGTATAGGCGAGAATCTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:1558270/233‑1 (MQ=255)
                                                                                                                        aacccTGCGAAGAGGCGGTTTGTTTCAAGACCCCCGCGCTGAGGGCGACCGGGGCGTGGATTGCTCGGCTACTGGGTCCGCTGGTGGGCGGGGCGTGGATCCATGTGCTGCCCGGGGAGGGGATGTTTTTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:1984240/233‑1 (MQ=255)
                                                                                                                        aaTCTTGCGAGGAGGCGGTTGGTATCAAGATCCCCGGGCGGTGGGGGAACGGGGCGCTGATTGCCCCGCTACTTGGCCGGCTGGTGGGCGGGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATTTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCg                                                                                                     <  2:1003909/235‑1 (MQ=255)
                                                                                                                        aaTCCTTCAAAGAGGCGGTTTTGATCAAGACACCCGCGCTGATGGCGGCGGGGGCACTGATTGCGCCGCTACTGGGTCCGGTGGTGGGCGCGGCGTGGATCCATTGGCTCCCCTGGGAGGGAAT‑TTTGTTGTTGTGTGCGGCTT‑GGGCAGCGATCTCCTTTTTCGGTCGG‑CAACGAGCCGTGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:1614255/236‑1 (MQ=255)
                                                                                                                            cTTCGAAGAGGCGGTTTGTATCAAGATCCCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACt                                                                                                                          <  2:1482424/210‑1 (MQ=255)
                                                                                                                            cTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACtata                                                                                             >  1:1984240/1‑239 (MQ=255)
                                                                                                                            cTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACt                                                                                                                          >  1:1482424/1‑210 (MQ=255)
                                                                                                                               cGAAGAGGCTGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGc                                                                                          <  1:1960769/239‑1 (MQ=255)
                                                                                                                               cGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGTGCTGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGGCGTGACTATAAGc                                                                                          >  1:1670223/1‑239 (MQ=255)
                                                                                                                                 aaGAGGCGGTTTTTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAAATGTCACTGAAAGAAATCGGTCGGGACTATaa                                                                                            >  1:476430/1‑235 (MQ=255)
                                                                                                                                      gCGGTTTTTATCAAGATCACCGCGCTGATGGCGGACGTTGCGCTGATTGCTCCGCTACTTGGTACGCTGGGGGGCGCGGCGTGGATCCATGTTCTGCCCTGGGGGGGGGTGTTTGTG‑TTGTTTGCCGCGT‑TGGCAGCGATATCCTTTTTTGGTCTG‑CAACGCGCCATGCCTGAAACCGCCACGCGTATATGCGGGAATCTGTGCCTGAAAGAATCCGGGGGTGCCTATAAGCTggggtc                                                                                   >  2:949042/1‑235 (MQ=255)
                                                                                                                                      gCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGGGACTATAAGCTGGt                                                                                      >  1:126200/1‑236 (MQ=255)
                                                                                                                                        ggTTTGTATCAAGATCACCGCGCTGATGGCGACCGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCCTGTGCTGCCCTGGGGGGGGGTGGTTGTT‑TTGTTTGCCGCCT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCCTGCCCGAAAACGCCACGCGTATAGGCGAGAAAATGTCCCTGAAAGAACTCGGGTGTGAATACAAACTgggggtgg                                                                                 >  2:493487/1‑233 (MQ=255)
                                                                                                                                               aTCAAGATCCCCGCGCTGAGGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACgg                                                                          <  1:2237855/239‑1 (MQ=255)
                                                                                                                                                       ccccGCGCTGAGGGCACCCGGGGCGCTGTTGGCTCCGCTCCTGGGTCCGCTGGTGGGCGCGGCGTGGATCCTTGTGTTGCCCGGGGGGGGGTGTTTTGTT‑TTGTTTGCCGCTT‑GGGCGGCGATCTCCTTTTTCGGTCGG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCttt                                                                   <  2:1419376/236‑1 (MQ=255)
                                                                                                                                                        aCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑Caa                                                                                                                                                                   >  2:1400108/1‑141 (MQ=255)
                                                                                                                                                        aCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑Caa                                                                                                                                                                   <  1:1400108/141‑1 (MQ=255)
                                                                                                                                                        aCCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCATGGATCCATTAGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCGGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGGATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTtgt                                                                 >  2:1310138/1‑239 (MQ=255)
                                                                                                                                                                    gggacggTGGCGCTGTTTGCTCCGCTCCTGGGTCGGCTGGTGGGCGCGGGGTGGATCCTTGTGCTGCCCTGGGGGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGt                                                                                                    <  2:534926/186‑1 (MQ=255)
                                                                                                                                                                           tGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGTTCCATGTGCTGCCCTGGGAGGGGTTGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCATGCGTCTAGGCGAGaa                                                                                                                             >  1:727936/1‑160 (MQ=255)
                                                                                                                                                                              cgcTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGc                                                                                          >  1:534926/1‑192 (MQ=255)
                                                                                                                                                                                   tttGCCCCGCTACTGGGTCCGCTGGTGGGCGCGGCGTGGATCCAGGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGgcgctggcgctg                                              <  1:502838/230‑1 (MQ=255)
                                                                                                                                                                                             tACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGGGCTGGGGCGGGGGTTGGtgt                                    >  2:502838/1‑229 (MQ=255)
                                                                                                                                                                                                  ggTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTg                                                                                        <  1:46205/174‑1 (MQ=255)
                                                                                                                                                                                                  ggTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTg                                                                                        >  2:46205/1‑174 (MQ=255)
                                                                                                                                                                                                    tCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATTTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTg                      <  1:1698786/238‑1 (MQ=255)
                                                                                                                                                                                                       ggtggtggGCGGGGGGGGGATCCAGGTGCGGCCGTGGGAGGGGATGTTTGTT‑TTGTGTGCCGCTT‑TGGCAGCGATCTCCTTTTTCGGTCGG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGt                  <  2:227940/237‑1 (MQ=255)
                                                                                                                                                                                                        cTGGTGGGGGGGGTGGGGTTCCAGGTGCTGCCCGGGGAGGGGGTGTTTGTT‑TTGTTGGCCGCAT‑TGGCAGGGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGcc                                                                                                                                              <  2:1441711/114‑1 (MQ=255)
                                                                                                                                                                                                         tggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACtat                                                                                              <  1:763980/161‑1 (MQ=255)
                                                                                                                                                                                                         tggtggGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACtat                                                                                              >  2:763980/1‑161 (MQ=255)
                                                                                                                                                                                                                 gcggtgTGGATCCGTGTGCTGCCCTGGGAGTGTATGTTTGTT‑TTGTTTGCCGCAT‑TGGCGGCGTCCTCCTTTTCCGTTCTG‑AAGCGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCa        <  1:1087952/239‑1 (MQ=255)
                                                                                                                                                                                                                  cggcgTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAg                                                                                                                                    >  1:1441711/1‑114 (MQ=255)
                                                                                                                                                                                                                     cggggATCATTGTGCTGCCTGGGGAGGGTAGTTTTGTT‑TTTTTGGCCGCAT‑GGGCAGCGACCTCCTTTTTCGGTCTG‑CAACGACCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGggcgc                                                      <  2:1139708/186‑1 (MQ=255)
                                                                                                                                                                                                                        ggATCCTTGTGCTGCCCTGGGGGGGGATGTTTGTT‑TTGTTGGCCGATT‑GGGCAGGGATCTCTTTTTCCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGcc  <  1:2210810/238‑1 (MQ=255)
                                                                                                                                                                                                                                    tgcCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGcc                                                                                                                                              <  2:1075732/86‑1 (MQ=255)
                                                                                                                                                                                                                                    tgcCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGcc                                                                                                                                              >  1:1075732/1‑86 (MQ=255)
                                                                                                                                                                                                                                      cccTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCg                                                                                                                                                                            >  1:1468559/1‑55 (MQ=255)
                                                                                                                                                                                                                                      cccTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCg                                                                                                                                                                            <  2:1468559/55‑1 (MQ=255)
                                                                                                                                                                                                                                      cccTGGGA‑GGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTg                         <  1:1523464/200‑1 (MQ=255)
                                                                                                                                                                                                                                        cTGGGAGGGGATGTTTTTT‑TTGTTGGCCCCAT‑TGGCAGCGATCTCCTTTTTCGGGCTG‑CAACGAGCCATGCCTAAAACCGcc                                                                                                                                              <  2:2086851/82‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAGGGGATTTTTGTT‑TTGTTTCCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGt                           <  1:556670/194‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                           
TTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTT‑TTGTTTGCCGCAT‑TGGCAGCGATCTCCTTTTTCGGTCTG‑CAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCC  >  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: ‑