Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913884,1840AG36.2% 65.7 / 12.1 59E171G (GAG→GGG) mdfAmultidrug efflux pump MdfA/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base A (20/17);  new base G (17/4);  total (38/21)
Fisher's exact test for biased strand distribution p-value = 5.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAG  >  NC_000913/883939‑884425
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
ggtgtggtTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGAAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATAAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTTGGCGCGGCGGGGGTCCATGTGCTGCCCTggggggggg                                                                                                                                                                                                                                               >  1:542926/1‑250 (MQ=255)
 gtgtggtTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAt                                                                                                                                                                                                                                              <  1:440137/251‑1 (MQ=255)
 gtgtggtTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAt                                                                                                                                                                                                                                              <  1:439838/251‑1 (MQ=255)
     ggtTTATCGTCAACTGTCTGGCTATATGTCTTGCGTAAAACATGGAGCCATTCACCCTTTTGCGCTGCTTGCAGGGCATATGCCTCTGTTTCACTTTCGCTGTGGGATACTCCGCAATTCAGGCATCCTTCGAAGGGGCGTTTTGTTTCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTCAGTTGTCCTCTGGTGGGCGCGGCGTTTATCTATGTGCTGACCGGGGAGGGGAtgtt                                                                                                                                                                                                                                           <  2:497521/250‑1 (MQ=255)
         tATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGAtgtttgtt                                                                                                                                                                                                                                       <  1:345232/250‑1 (MQ=255)
               cACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTTTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTtttgtttg                                                                                                                                                                                                                                <  1:282720/251‑1 (MQ=255)
                  cTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATCGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCCTGTGCTGCCCTGGGAGGGGATGGTTGTTTTGTTTGCCg                                                                                                                                                                                                                             >  1:489359/1‑251 (MQ=255)
                        ggCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGTGGGTCCATGTGCTGCCCTGGGAGGGGGTGTTTGGTTTGTTTGCCGCATTgg                                                                                                                                                                                                                       >  1:506385/1‑251 (MQ=255)
                         gCAATATTGCTGGCGCACAACATTGAACCATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGATAAAAGGTCACCCCGCTGGTTGCGCAAGGGGCGCTTGTTGTCCCGCTACTTGGGTCGCTGGGTGGCGACGCGGGGATCCATGGGCTGCCCTGGNGGGGGGGGGTTGTTTTGGTTATCGCATTggg                                                                                                                                                                                                                      >  2:248893/1‑250 (MQ=255)
                         gCAATATTGATGGCTCAAAACTTTGAATAATTCTCCCTGTTTCGCTGTTTGCAGGGCATAAGCATCTGTTTCATTGTCGCTTTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTTGTACGCTGGTGGGCGCGGCGTGGATCCATGTGCTGTCCCGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGc                                                                                                                                                                                                                      <  2:517248/251‑1 (MQ=255)
                              atTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGGGCGGCGGGGATCCCTGTGCTGCCCTGGGAGGGGGTGTTTGTTTTGTTTGCCGCATTGGCAggg                                                                                                                                                                                                                  >  1:89163/1‑248 (MQ=255)
                              atTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGTGGGTACATGTGGTGCCCTGGGGGGGGGTGGTTGTTTTGGTTTCCGCAtttgccgcgg                                                                                                                                                                                                                 >  2:433401/1‑243 (MQ=255)
                                 gCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGAtct                                                                                                                                                                                                              <  1:311767/251‑1 (MQ=255)
                                        cAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCttttt                                                                                                                                                                                                       <  1:190788/251‑1 (MQ=255)
                                           aaCATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGGGCGGCGTGGGTCCATGGGGTGCCCTGGGCGGGGGTGTTTGTTTTTNTTTCCCCATTGGAGGCGGGCTCCtttttacg                                                                                                                                                                                                    >  2:256385/1‑248 (MQ=255)
                                             cATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGGATCAAGATCACCGCGCTGATGGCGAAAGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGGTTTTTTTGCCGCATTGGCAGCGATCTTCTTTTTCGGTc                                                                                                                                                                                                 >  1:15016/1‑251 (MQ=255)
                                             cATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGAGCTGTGGGCTACGCCGCACTTCAGGAATCCTTCGAAGAGCCTGTTTTTCTAAAGATCCCCGCGCTCATGGCTGACGAGGCGATGATTGCTCCGCTTCTTGGTCCGCTGGATGGCGCGGCGGCGGGTCATGTCCTGCCCGGGGAGGAGATGTTTGGTTTGGTTGCCGCACTGGCAGCGGACTTCTTTTTccgac                                                                                                                                                                                                  >  1:257757/1‑247 (MQ=255)
                                                     aaTTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTTCCGCATTGGCAGCGATCTCCTTTTTAGGTCTGTAACg                                                                                                                                                                                          >  1:324845/1‑250 (MQ=255)
                                                        tCACCCTGTTGCGCTTCTTGCAGGGCATAATCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTGTTCGGTCTGCAACGAGc                                                                                                                                                                                       <  2:68596/250‑1 (MQ=255)
                                                         cACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGCTGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGTGGATCCCTGTGCTGCCCTGGGAGGGGGTGTTTGTTTTGTTTGCCGCATTGGCCGCGCTCTCCTTTTTCGGTCTGCAACGGGccc                                                                                                                                                                                     >  2:174741/1‑250 (MQ=255)
                                                         cACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCAACGCGCTGATGGCGGACGTGGCGCTGGTTGGTCCGCTACTTGGTCCGCCGGTGGGCGCGGCGTGGGTCCATGTGTTGCCCTGGGAGGGGGTGGTTGTTTTGTTTGGCCCCTTGGGAGCGATTTCCttttttgggttgcaacaggcc                                                                                                                                                                                       >  2:376736/1‑234 (MQ=255)
                                                               gTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATAACCGCGCTTATGGCGCACGTGGCGCGGATTGCTCTGATACCTGGGTCGCTGGGGGGCGCGGCGTGGGTACCAGTGGTGGCCTGGGAGGGGGTATTTGTTTTGTTTGCCGCATTGGCCGCGAGCTCCTTTTTCGGCCTGCAAcaaggcctgccgg                                                                                                                                                                               >  2:439376/1‑239 (MQ=255)
                                                                  gcgcTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGaa                                                                                                                                                                             >  1:120479/1‑250 (MQ=255)
                                                                      ttcttGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTCTCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCGTTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGc                                                                                                                                                                        >  1:329289/1‑251 (MQ=255)
                                                                       tcttGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGGTGTTTGGTTTGTTTGCCGCATTGGCAGCGAACTCCTTTTTCGGTCTGCAACGAGCCAAGCCTGAAACCGcc                                                                                                                                                                       >  1:136436/1‑251 (MQ=255)
                                                                            cAGGGCATCAGACTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCCAGAACACCGCGCTGGTGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTACGCTGGTGGGCGCGGGGCGGGTCCATGGGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTGGCGCATTGGCCGCGCTCTCCTTTTTCGGTTTGCCACGGGCCCTGCCTGCAACCGccccccg                                                                                                                                                                  >  2:305173/1‑246 (MQ=255)
                                                                            cAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCAcgcg                                                                                                                                                                  >  1:78049/1‑251 (MQ=255)
                                                                             gggtcATAAGCCTCTGTTTAATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTGCGAAGAGGCGGTTTGTAGCAAGATCCCCGCGCGGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCGGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCACCGAGCCATGCCTGAAACCGCCACGCGt                                                                                                                                                                 <  1:437627/247‑1 (MQ=255)
                                                                             aGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATCAAGATCACCGCCCTGATGGCGAACGGGGCGCTGGTTGCTCCGCTACTTGGTCCGCTGGGGGGCGGGGCGGGGACCCAAGTTCTGGCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTTGGTCTTCAACTAGCCCTGCCTTAAACCGCCACCCGt                                                                                                                                                                 >  2:224457/1‑251 (MQ=255)
                                                                              gggCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCCAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGGTCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCAATGGCCGGGGTCTCCCTTTTCTGGCTGCAACGAGCCCTGCCTGAAACCGCCACGCGtg                                                                                                                                                                >  2:458608/1‑250 (MQ=255)
                                                                                 cATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATAAAACTCCACGCGCTGAAGGCGAAAGTGGCGCCGAGTGCCCCGCTACTTTGGCCGCCGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGCTGTTTGTTTTGTGTGCCGCATTTGCAGCGTTCTTCTTTTTCGTTCTGCAACGACCCATggcctataccgcatcgctaaacg                                                                                                                                                             >  2:22947/1‑229 (MQ=255)
                                                                                            gTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGCTTGTATCAAGATCACCGCGCTGCTCGCGAACGCGGCGCTGTTTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGCGGATGCCTGCGCTGCCCCGGGGGGGGATGTTTGGTTTGTTTGCCGCAATGGTAGCGATCTCCTTTTTTGGACTTCAACGCGCCCTGCCTGAAACCGCCACGCGGATAGGGGAGAAAccg                                                                                                                                                  >  1:127898/1‑249 (MQ=255)
                                                                                              ttCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGAAGGCGAACGTGGCGCTGATTGTTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGGGGATCCATGTGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTGCCGCATTGGCAGCGGTCTCCTTTTTTGGTCTGGAACGAGCCATGCCCGAAACCGCCACGCGGAAAGGCgggacaatggc                                                                                                                                                >  1:502186/1‑241 (MQ=255)
                                                                                                 aTTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGTGGATCCCTGTGCTGCCCTGGGGGGGGGTGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCCTGCCTGAANCCGCCCCGCGTATTGGCGCGAAACTGTccct                                                                                                                                             >  2:235225/1‑248 (MQ=255)
                                                                                                       gcTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGCTGGCGAACGTGGCGCTGATTGGTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGGCCTGGGAGGGGGAGGTTGTTTTGTTTTCCGCCTTGGGCGGGGTCTCCTTTTTTGGTCTGCAAAGAGGCATGCCTGACACCGCCACGCGGATAAGCGAGGAAATGGCACTGaaaa                                                                                                                                        >  1:52342/1‑249 (MQ=255)
                                                                                                                     cgcAATTCAGGAATCCTTCGAAGAGGCGGTTTTTATCAAGATCACCGCGCTGATGGCGAACGGGGCGCTGATTGCTCCGCTACTTTGTCCGCTGGTGGGCGCGGCGGGGATCCATGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCCCATTGGCAGCGATATCCTTTTTCGGTCTGCAACGAGCCCTGCCTGAAACCGCCCCGCGTATAGGCGGGAAACTGTCACTGAAAGAAATCGGTCGGGAcc                                                                                                                         >  2:45434/1‑250 (MQ=255)
                                                                                                                      gcAATTCAGGAATCCTTCGAAGAGGCGGATTGTATCAAGATCACCGCGCTGATTGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGGGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGGGGGGATGTGTGGTTTGTTTGCCGCATTGGCAGCGGTCTCCTTTTTCGGTCTGCAACTAGCCATGCCTGAAACCGCCACGCGTATATGCGAGAAACTGTCACTGTAAGAAATCTGTCGTGACt                                                                                                                         >  1:391779/1‑250 (MQ=255)
                                                                                                                        aaTTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAAATCACCGCGCTGATTGCGAACGTGGCGCTGGTTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCCTGTGCTGCCCTGGGGGGGGATGTTTTTTTTGTTTGCCGCCTTGGTCGCGGTCTCCTTTTTTGGTCTGGCACGAGCCCTGCCTGAAACCGCCACGCGTATAGGCGGGAAACTGTCCCTGTAAGAACCCGGTCGTGGCtat                                                                                                                       >  2:220550/1‑250 (MQ=255)
                                                                                                                            cAGGAATCCTTCGAAGAGGCGGCTTGTAACAAAATCACCGCGCTGATGGGGAACGTGGCGCTGGTTGCTCCGCTACTTTGTACGCTGGTGGGCGCGGCGTGGATCCCTGTGCTGCCCTGGGAGGGGATGGTTGTTTTGGTTGCCGCCTTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGACACCGCCACGCGTATAGGCGCGAAACTGTCACTGAAAGAACCCGGTCGTGACTATAAGCt                                                                                                                  >  2:354644/1‑251 (MQ=255)
                                                                                                                              ggAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCCAACGGGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATTCATGTGCTGCCCTGGGGGGGGAGGGTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGGCTGCAACGAGCCACGCCTGTAAACGGCACGCGTATTGGGGGGGAACTGTCACTGGAAGAACTCGGGGGTGACGATAAGcggg                                                                                                                >  2:97144/1‑248 (MQ=255)
                                                                                                                                          ccgttctgTTTGTACCAAGATCCCCGCGCTGNGGGCGAACGTGGCCCTGATTGCTCCGTTCCTTGGTCCGCTGGTGGGCGCGGCGTGGACCCAGGTGCTGCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACg                                                                                                    <  2:239430/244‑1 (MQ=255)
                                                                                                                                             ggCGTTTTGTATAAAGATCCCCGCTCTGATGGCGACCGGGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGGGGATCCATGTGTGGCCCTGGGAGGGGATGTTTTTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCACCGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGcc                                                                                                 <  1:363697/251‑1 (MQ=255)
                                                                                                                                                 gTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACg                                                                                                    >  1:239430/1‑244 (MQ=255)
                                                                                                                                                    tgtacaaaGAACACCGCGCTGAGGGGGACCGTGGCGCTGATTGCCCCGCTCCGTGGTCCGCTGGTGGGCGGGGCGTGGATCCATGTGCGGCCCGGGGAGGGGGTGTTTGTTTTGTTTGCCGCGTGGGCGGCGTTCTCCTTTTTCGGTCGGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTtg                                                                                           <  2:138112/244‑1 (MQ=255)
                                                                                                                                                        tCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCg                                                                                      >  1:264078/1‑251 (MQ=255)
                                                                                                                                                              aTCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGACACCGCCACGCGTATAGGCGAGCAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCggaggg                                                                                 >  1:39710/1‑247 (MQ=255)
                                                                                                                                                              aTCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCgggggg                                                                                 >  2:195881/1‑250 (MQ=255)
                                                                                                                                                                    gcgcTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGCCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCggggggggcggg                                                                           >  1:517581/1‑244 (MQ=255)
                                                                                                                                                                     cgcTGATGGAGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGCTGGGCGCGGCGTGGATCCATGTGATGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGTCTATAAGCTGGTGCTGTAGAACGGGCGCTTTGTTGCggggggggggggg                                                                          >  1:350542/1‑243 (MQ=255)
                                                                                                                                                                               gAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCgggggggcgggggcggggggtgg                                                                >  1:139741/1‑233 (MQ=255)
                                                                                                                                                                                 aggTGGCGCTGATTGCTCCGCTTCTGGGCCCGCTGGGGGGCGGGGCGTGGATCCATGTGCTCCCCTGGGGGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCACCGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAGCGGCCGCTTTGTGGCGGTGGCGCTGGCGCTGGGATTCGtt                                                              <  1:435563/249‑1 (MQ=255)
                                                                                                                                                                                   gTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCgggggggccggggccgggggtggtt                                                              >  1:326762/1‑229 (MQ=255)
                                                                                                                                                                                     ggggcTGATTGCTCAGCTACTTGGTCGGTTGGTGGGCGCGGGTGGGATCCTTGTGCTGCCCGGGGGGGGGATGTTTGTTTTTTTTGCCGCATTGGGAGCTATCGCCTTTTTCGGTCGGCAACGAGCCATGCATGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTc                                                          <  1:495331/248‑1 (MQ=255)
                                                                                                                                                                                               gCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCAGGTGCTGCCCTGGGAGGGGATGTTTGTTTTTTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGc                                                 <  1:45434/250‑1 (MQ=255)
                                                                                                                                                                                                              ccctcTGGTGGGCGCGGCGTGGTTCCATGTGCTGCCCTGGGAGGGAATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTGTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGccc                                  <  1:489090/247‑1 (MQ=255)
                                                                                                                                                                                                                 gCTGGTGGGCTCGGCTTGGATCCTTGTTCTGCCTGGGTAGTGGGTTTTTTTTTTTTTTGCCGCATTGGCAGCTTCCTCCTTTTCGGGTCGGCAGCGGGCCAGGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACACGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGTCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGt                               <  2:1870/251‑1 (MQ=255)
                                                                                                                                                                                                                        ggggcggcgTGGATCCATGTGCTGCCCTGGGAGGGGATTTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCACCGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGAt                        <  1:195881/248‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgctttGGGGGGGGTTGTTTGTTTTGTTGGCCGCTGTGGCCGCGATCCCTTTTTCCGGTCTACAGCGGGCCAGGCCGGAAACCGCCTCGCGTAAAGGCGAGAAACTGTCACTGAAAGAACTCGGGCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGGTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTTGCCGATTATCATCATTACCGGCGAGCAg  <  2:504079/246‑1 (MQ=255)
                                                                                                                                                                                                                                                       gggAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCa                                                                                                                                                                      >  1:236458/1‑81 (MQ=255)
                                                                                                                                                                                                                                                       gggAGGGGATGTTTGTTTTGTTGTCCGCATTGGCAGCGACCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCa                                                                                                                                                                      <  2:236458/81‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                  
GGTGTGGTTTATCGTCACCTGTCTGGCAATATTGCTGGCGCAAAACATTGAACAATTCACCCTGTTGCGCTTCTTGCAGGGCATAAGCCTCTGTTTCATTGGCGCTGTGGGATACGCCGCAATTCAGGAATCCTTCGAAGAGGCGGTTTGTATCAAGATCACCGCGCTGATGGCGAACGTGGCGCTGATTGCTCCGCTACTTGGTCCGCTGGTGGGCGCGGCGTGGATCCATGTGCTGCCCTGGGAGGGGATGTTTGTTTTGTTTGCCGCATTGGCAGCGATCTCCTTTTTCGGTCTGCAACGAGCCATGCCTGAAACCGCCACGCGTATAGGCGAGAAACTGTCACTGAAAGAACTCGGTCGTGACTATAAGCTGGTGCTGAAGAACGGCCGCTTTGTGGCGGGGGCGCTGGCGCTGGGATTCGTTAGTCTGCCGTTGCTGGCGTGGATCGCCCAGTCGCCGATTATCATCATTACCGGCGAGCAG  >  NC_000913/883939‑884425

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: