Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913867,5860AT28.3% 122.3 / 13.3 60M23L (ATG→TTG) gsiAglutathione ABC transporter ATP binding subunit GsiA
Reads supporting (aligned to +/- strand):  ref base A (8/35);  new base T (17/0);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 2.79e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATAT  >  NC_000913/867338‑867836
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
ggAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATg                                                                                                                                                                                                                                                             <  1:54001/251‑1 (MQ=255)
           aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCCACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTCATGTGCTGGCGGTTGGAACTCTAGATATTTCCTTTTTGCAGGCCcagcg                                                                                                                                                                                                                                                 >  2:247475/1‑250 (MQ=255)
                                     gcgctgtgCGGTAGCGGTGGCTTAATCGCTATCGCCCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTAGCGCGCCGGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAAtctct                                                                                                                                                                                                                       <  2:131171/244‑1 (MQ=255)
                                       gcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATctctct                                                                                                                                                                                                                     <  1:153754/251‑1 (MQ=255)
                                       gcTTGGCGGTAACGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGACTGGGGCTACGCAGGCGGTACGCCAACCCACGGTATCTACCGTGAAAAAGGGGACACCGCTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTTAAAAACTGAATATTTGCTTTTTTCAGGACCAGCAGAAAACAGCTCCgggcggcacaccctcc                                                                                                                                                                                                                     >  2:452105/1‑237 (MQ=255)
                                          tGGCGGTAGCTGTTGCTTAATCGCTATCGACCATGAAGTGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCttt                                                                                                                                                                                                                   <  2:14515/250‑1 (MQ=255)
                                                   cgcgggCTTAATCGCTATCGACCAGGAAGTGAATGTCGCGCTCCCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGCTACGCAAACCACCGGTATCTACCGGGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAAGCTCTCTTTTAGTCTGCaa                                                                                                                                                                                                        <  2:108473/247‑1 (MQ=255)
                                                    ggTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGGATATTGCCTTTTTGCAGGACCAGCAGAAAATAGCTGGGGTCCGCCATCTCTCTTTTTGTCTGCAAc                                                                                                                                                                                                       >  1:459264/1‑251 (MQ=255)
                                                    ggTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAAc                                                                                                                                                                                                       <  2:459264/251‑1 (MQ=255)
                                                     gTGGCTTAATCGCTATGGACCATGAAGGGAATGTCGCGATCCCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTATGCCGGACCAGCAGAAAATAGCTGCGGTCCGCACTCTCTCTTTTTGTCTGCAAcg                                                                                                                                                                                                      >  1:167937/1‑251 (MQ=255)
                                                       ggCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAAcgc                                                                                                                                                                                                     <  1:338747/250‑1 (MQ=255)
                                                           tAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGGAAATCTGAATAATGCCTTTTTGCAGGACCAGCAGAAAATAGCTGCGGTCCGCCATCTCCCTTTTTGTCTGCAAAGCggggg                                                                                                                                                                                                >  2:365100/1‑248 (MQ=255)
                                                             aTCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCAGGACCAGCAGAAAATAGCTGCGGGCCCCAATTTCTCTTTTTGTCTGCAACGCgggg                                                                                                                                                                                                 >  2:136202/1‑246 (MQ=255)
                                                               cGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACg                                                                                                                                                                                            <  1:135604/251‑1 (MQ=255)
                                                                    tCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCAGGACCAGCAGAAAATAGCTGCGGGCCGCAATCTCTCTTTTAGGCTGCAACGCGGGGAGACGCTggg                                                                                                                                                                                       >  1:234665/1‑250 (MQ=255)
                                                                       aCCAGGAAGGGAATGTCGCGCTACCGGTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAAt                                                                                                                                                                                    <  2:149118/251‑1 (MQ=255)
                                                                        ccATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAAtt                                                                                                                                                                                   <  1:41460/251‑1 (MQ=255)
                                                                          aTGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGTAAATCTGAATATTGCCTTTTTGCAGGACCAGCAGAAAATAGCTCCGGGCCGCAATCTCTCTTTTTGTCTGCAACGCGGGgagaagatgggcattgt                                                                                                                                                                                 >  2:256108/1‑238 (MQ=255)
                                                                            gAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAAttgtt                                                                                                                                                                                <  1:120241/250‑1 (MQ=255)
                                                                             aaGGGAATGTCGCGCTACCGTTTACCCCCGAAGGAATGAACCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACCGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGCAGGACCAGCACAAAATAGATGCGGTCCGCAAACTTTCTTTTCATCTGCAAAGCgggtggcagccggccattgttgg                                                                                                                                                                              >  2:212592/1‑230 (MQ=255)
                                                                                 gAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTTAAAAACTGGAAATTTTCTTTATGCAGGGACAGCAGAAAATAAGTGCGGGCCGCAAACTCTCTTTTAATTGTCAAAGCGGGGAGACGCTGGCCATTGTTGGCGa                                                                                                                                                                           >  2:53926/1‑250 (MQ=255)
                                                                                      tCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGGTGAAAAACTGAATATTGCCTTTTTGCAGGACCAGCAGAAAATAGACGCGGTCCCCAATCTCTCTTTTAGACTGCAAAGCGGGGGGACGCTGGCAATTGTTGGCGAATCCgg                                                                                                                                                                     >  2:155838/1‑251 (MQ=255)
                                                                                                  ttAACTCCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGt                                                                                                                                                          <  1:247475/250‑1 (MQ=255)
                                                                                                    aaCACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAg                                                                                                                                                       <  1:299637/251‑1 (MQ=255)
                                                                                                        ccGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCAGGACCAGCCGAAAATAGCTGCGGTCCGCCATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTTGCGAATCCCGCTCCGGTAAGTCAGGGAc                                                                                                                                                  >  1:325222/1‑251 (MQ=255)
                                                                                                          gAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGGAAATCTTTATATTGCCTTTTTGGCGGGCCCGCAAAAAATAGCTGCGGGCCCGAATCTCTTTTTTAGTCTGCAACGCGGGGAGACGCATGCCATTGTTTGCGAATCCGGCCCCCGTAAGTCAggggccg                                                                                                                                                 >  2:181338/1‑245 (MQ=255)
                                                                                                               aaTGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGGCCGCAATCTCTCTTTTAGTCTTCCACGCGGTGAGACGCTGGCAATTGTTGGGGGATCCGGCTCCGGTAAGGCAGGGACTGCGTTg                                                                                                                                           >  1:155693/1‑251 (MQ=255)
                                                                                                                aTGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTGCCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCCGGACCAGGAGAAAAAAACGGCGGGCCGCAAACCCTTTTTTATTCCTCAACGCCGGGCGACGCCGGCAATTGTTTGCGGATCCGGCCCCGGTAAGTCCGTggaggcgtttg                                                                                                                                           >  2:431123/1‑241 (MQ=255)
                                                                                                                  gTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGCAGGGCCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTTCAAAGCGGGGAGACGCTGGCAATTTTTTGCGAATACGGCGCCGGGCAGGCAGTGACCGCGTTgggg                                                                                                                                        >  2:363126/1‑249 (MQ=255)
                                                                                                                    aTCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTGAGTCTGCAACGCGGTGAGACGCTGGCAAGTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCAtt                                                                                                                                      <  1:408019/251‑1 (MQ=255)
                                                                                                                                   aCGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTgg                                                                                                                        <  2:93589/250‑1 (MQ=255)
                                                                                                                                   aCGCAGGCGATACGCCAACCACCGGGATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATAGTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATGGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCTTTGGCATTGATGCGCCTGTTg                                                                                                                         <  2:420239/249‑1 (MQ=255)
                                                                                                                                        ggCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAgg                                                                                                                  <  1:457408/251‑1 (MQ=255)
                                                                                                                                          cGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGGATATTTCCTTTTTGCCGGACCCGCAGAAAATAGCTGCGGGCCGCCATCTATCTTTTAGTCTGCCAAGCGGGGGGGCGCTGGCAATTGTTGGCGAATCCGGCTCCGGGAAGGCAGCGGCTGCGTTGGCCTTTATGGGCCTGTTGGGACAgggg                                                                                                                >  1:205419/1‑249 (MQ=255)
                                                                                                                                          cGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcg                                                                                                                <  1:305507/251‑1 (MQ=255)
                                                                                                                                            aTACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcggg                                                                                                              <  1:141097/251‑1 (MQ=255)
                                                                                                                                              aCGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggc                                                                                                             <  1:256108/250‑1 (MQ=255)
                                                                                                                                              aCGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAAACTGAATATTGCCTTTATGCAGTACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTGTAGTCTGCAACGCGGTGAGACGCGGGCAATTGGTGGCGAATCCGGCTCCGGTAAGTCAGGGACTGCGTTGGCATGGATGCGCCTGTTGGAACAggcgggg                                                                                                             <  2:454138/250‑2 (MQ=255)
                                                                                                                                                   aaCCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGTAGGGCCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGGCTGCAACGCGGTGAGACGCCGGTAATTGTTTGCGAATTCGGCTCCGGTAAGTCAGGGGACGCGGTGGGCTTTCTGCGCCTGTTGGGACCGgcgggcgggttt                                                                                                       >  2:147747/1‑247 (MQ=255)
                                                                                                                                                    cccaccGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                      <  2:17831/250‑1 (MQ=255)
                                                                                                                                                           gTATCTACCGTGAAAAAGGGGACACCGTTGCCACACCGTGGTAAACTTGGTGCCGGTAATGTGCTGGCGGTTGAAAATCTAAATATTTCCTTTTTGCATGAACCACAGAAAATAACAGCGGGACGCAAACTCTCTTTTAGTTTGCAACGCGGTGAGACGCTAGCAATTGTTGGCGAATCCGGCTTCGGTTAGTCAGGGACCGCGTTGGCATTGATGCGCCTTTTGGGAAAggcgggggggttttggacgtg                                                                                               >  1:427328/1‑236 (MQ=255)
                                                                                                                                                             aTCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGCATATTGCCTTTATGCAGGACCCGCAGAAAATAGCTGGGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGGGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGGAAGTCAGTGGCTGCGTTGGCCTTGATGCGCCTGTTGGAAAAggcgggggggtttggacagtggg                                                                                             >  1:284583/1‑234 (MQ=255)
                                                                                                                                                                       aaaaaGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCt                                                                                   <  1:109947/251‑1 (MQ=255)
                                                                                                                                                                           aGGGGACAGCGTTGCCACACTGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATAGTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTTGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTg                                                                               <  2:330840/251‑1 (MQ=255)
                                                                                                                                                                                acacCGTTGCCACACAGTGATGAACTTGAGGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATGGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcg                                                                          <  1:442132/251‑1 (MQ=255)
                                                                                                                                                                                       ggCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATTCGCCTGTTGGAACAGGCGGGCGGTTTATTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAg                                                                    <  2:167937/249‑1 (MQ=255)
                                                                                                                                                                                        gCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTc                                                                  <  1:136202/251‑1 (MQ=255)
                                                                                                                                                                                               aGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGt                                                           <  1:175568/251‑1 (MQ=255)
                                                                                                                                                                                                      aaCTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGTCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCCCAATCTCTCTTTTAGTCTGCAAAGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGGGGGTTTGTAAAGTGCGGTAAAATGCTGTTGGCGCGGGGCGGTCGCGGAGGGAttta                                                     >  1:153231/1‑248 (MQ=255)
                                                                                                                                                                                                      aaCTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGaa                                                    <  2:325222/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaTGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGGGACGCTGGCAATTGTTGGCGGATCCGGCCCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTTGAACAGGCGGGGGGGTTAGTAAAGTGCGGTAAAAAACTGTTGCAACGGGGCCGGCGCGAAGTGGTTTGACTTAGCGAGCAGa                                       >  1:91221/1‑250 (MQ=255)
                                                                                                                                                                                                                              cgggtGAAAATCGGAATATGGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCTGCAATCTCTCTTTTAGGCTGCAACGCGGTGAGACGCTGGCAATGGTTGGCGAATCCGGCTCCGGTAAGTCAGTGTCTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGTCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACGTAGCGAGCAGAACGCTGCACaa                            <  2:419131/247‑1 (MQ=255)
                                                                                                                                                                                                                                gTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAAt                          <  2:153231/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTGAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATgc                        <  1:295225/250‑1 (MQ=255)
                                                                                                                                                                                                                                               atTGCCTTTATGCAGGTCCAGCAGAAAATAGCTGCGGTCCTCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGGTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTGTAGTATAGGGCTATAAAATGCTGTTGCAGCGGCGCAGTCGCTATTTGATGGAACTTAGCGAGCATAACGCTGCACAAATGCGCCATGTTCGCgg           <  1:421127/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ttGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTAGTACAGTGGGGTAAAATGCTGTTGCAGCGGCGCTGTCGCGAAGTGATTGAACTTAGCGGGCAGCAAGCTGCACAAATGTGCCCTGTTCGCggg          >  2:44754/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGaa                                                             <  1:275723/198‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcgcgga                                                             >  2:275723/1‑196 (MQ=255)
                                                                                                                                                                                                                                                       tATGCAGGACCAGCAGAAAAGACCGGCGGTCAGCAATCGCTCGTTTAGTCTGCAACTCGTTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGATTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTAGTACAGGGCGATAAAATGCGGTTGCAGCAGCGCAGTCGCGAAGTGATTGAACGTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCgtat    <  2:190616/250‑4 (MQ=255)
                                                                                                                                                                                                                                                       tATGCAGAACAAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCACTTGTTGGCGAATCAGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGTGGGTTTGTAAAGTGGGATACAATGCTGTTTCCGCGGCGCGGGCGCGGTGTGAATGAACTTAACGAGCAGGACGCTGCACAAATGCACCATGGTCGCgggtcgggtt   >  2:387170/1‑243 (MQ=255)
                                                                                                                                                                                                                                                        aTGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTAAAGTGCGATAAAATGCTGTTGCAGCGGCGCAGACGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACCAATGGGCCATGTTCGCGGTGCGGgtat  >  1:415097/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
GGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCAGGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTGCAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTGACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATAT  >  NC_000913/867338‑867836

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