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.

GGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATAT  >  NC_000913/867338‑867836
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
ggAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATg                                                                                                                                                                                                                                                             <  1:54001/251‑1 (MQ=255)
           aGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCCACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTCATGTGCTGGCGGTTGGAACTCTAGATATTTCCTTTTTGCA‑GGCCcagcg                                                                                                                                                                                                                                                 >  2:247475/1‑250 (MQ=255)
                                     gcgctgtgCGGTAGCGGTGGCTTAATCGCTATCGCCCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTAGCGCGCCGGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAAtctct                                                                                                                                                                                                                       <  2:131171/244‑1 (MQ=255)
                                       gcTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATctctct                                                                                                                                                                                                                     <  1:153754/251‑1 (MQ=255)
                                       gcTTGGCGGTAACGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGACTGGGGCTACGCAGGCGGTACGCCAACCCACGGTATCTACCGTGAAAAAGGGGACACCGCTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTTAAAAACTGAATATTTGCTTTTTTCA‑GGACCAGCAGAAAACAGCTCCgggcggcacaccctcc                                                                                                                                                                                                                     >  2:452105/1‑237 (MQ=255)
                                          tGGCGGTAGCTGTTGCTTAATCGCTATCGACCATGAAGTGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCttt                                                                                                                                                                                                                   <  2:14515/250‑1 (MQ=255)
                                                   cgcgggCTTAATCGCTATCGACCAGGAAGTGAATGTCGCGCTCCCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGCTACGCAAACCACCGGTATCTACCGGGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAAGCTCTCTTTTAGTCTG‑Caa                                                                                                                                                                                                        <  2:108473/247‑1 (MQ=255)
                                                    ggTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGGATATTGCCTTTTTGCA‑GGACCAGCAGAAAATAGCTGGGGTCCGCCATCTCTCTTTTTGTCTG‑CAAc                                                                                                                                                                                                       >  1:459264/1‑251 (MQ=255)
                                                    ggTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAAc                                                                                                                                                                                                       <  2:459264/251‑1 (MQ=255)
                                                     gTGGCTTAATCGCTATGGACCATGAAGGGAATGTCGCGATCCCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTATGCC‑GGACCAGCAGAAAATAGCTGCGGTCCGCACTCTCTCTTTTTGTCTG‑CAAcg                                                                                                                                                                                                      >  1:167937/1‑251 (MQ=255)
                                                       ggCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAAcgc                                                                                                                                                                                                     <  1:338747/250‑1 (MQ=255)
                                                           tAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGGAAATCTGAATAATGCCTTTTTGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCCATCTCCCTTTTTGTCTG‑CAAAGCggggg                                                                                                                                                                                                >  2:365100/1‑248 (MQ=255)
                                                             aTCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCA‑GGACCAGCAGAAAATAGCTGCGGGCCCCAATTTCTCTTTTTGTCTG‑CAACGCgggg                                                                                                                                                                                                 >  2:136202/1‑246 (MQ=255)
                                                               cGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACg                                                                                                                                                                                            <  1:135604/251‑1 (MQ=255)
                                                                    tCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCA‑GGACCAGCAGAAAATAGCTGCGGGCCGCAATCTCTCTTTTAGGCTG‑CAACGCGGGGAGACGCTggg                                                                                                                                                                                       >  1:234665/1‑250 (MQ=255)
                                                                       aCCAGGAAGGGAATGTCGCGCTACCGGTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAAt                                                                                                                                                                                    <  2:149118/251‑1 (MQ=255)
                                                                        ccATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAAtt                                                                                                                                                                                   <  1:41460/251‑1 (MQ=255)
                                                                          aTGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGTAAATCTGAATATTGCCTTTTTGCA‑GGACCAGCAGAAAATAGCTCCGGGCCGCAATCTCTCTTTTTGTCTG‑CAACGCGGGgagaagatgggcattgt                                                                                                                                                                                 >  2:256108/1‑238 (MQ=255)
                                                                            gAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAAttgtt                                                                                                                                                                                <  1:120241/250‑1 (MQ=255)
                                                                             aaGGGAATGTCGCGCTACCGTTTACCCCCGAAGGAATGAACCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACCGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGCA‑GGACCAGCACAAAATAGATGCGGTCCGCAAACTTTCTTTTCATCTG‑CAAAGCgggtggcagccggccattgttgg                                                                                                                                                                              >  2:212592/1‑230 (MQ=255)
                                                                                 gAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTTAAAAACTGGAAATTTTCTTTATGCA‑GGGACAGCAGAAAATAAGTGCGGGCCGCAAACTCTCTTTTAAT‑TGTCAAAGCGGGGAGACGCTGGCCATTGTTGGCGa                                                                                                                                                                           >  2:53926/1‑250 (MQ=255)
                                                                                      tCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGGTGAAAAACTGAATATTGCCTTTTTGCA‑GGACCAGCAGAAAATAGACGCGGTCCCCAATCTCTCTTTTAGACTG‑CAAAGCGGGGGGACGCTGGCAATTGTTGGCGAATCCgg                                                                                                                                                                     >  2:155838/1‑251 (MQ=255)
                                                                                                  ttAACTCCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGt                                                                                                                                                          <  1:247475/250‑1 (MQ=255)
                                                                                                    aaCACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAg                                                                                                                                                       <  1:299637/251‑1 (MQ=255)
                                                                                                        ccGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCA‑GGACCAGCCGAAAATAGCTGCGGTCCGCCATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTTGCGAATCCCGCTCCGGTAAGTCAGGG‑Ac                                                                                                                                                  >  1:325222/1‑251 (MQ=255)
                                                                                                          gAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGGAAATCTTTATATTGCCTTTTTGGC‑GGGCCCGCAAAAAATAGCTGCGGGCCCGAATCTCTTTTTTAGTCTG‑CAACGCGGGGAGACGCATGCCATTGTTTGCGAATCCGGCCCCCGTAAGTCAggggccg                                                                                                                                                 >  2:181338/1‑245 (MQ=255)
                                                                                                               aaTGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGGCCGCAATCTCTCTTTTAGTCTT‑CCACGCGGTGAGACGCTGGCAATTGTTGGGGGATCCGGCTCCGGTAAGGCAGGG‑ACTGCGTTg                                                                                                                                           >  1:155693/1‑251 (MQ=255)
                                                                                                                aTGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTGCCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTTCCTTTTTGCC‑GGACCAGGAGAAAAAAACGGCGGGCCGCAAACCCTTTTTTATTCCT‑CAACGCCGGGCGACGCCGGCAATTGTTTGCGGATCCGGCCCCGGTAAGTCCGTggaggcgtttg                                                                                                                                           >  2:431123/1‑241 (MQ=255)
                                                                                                                  gTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGCA‑GGGCCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTT‑CAAAGCGGGGAGACGCTGGCAATTTTTTGCGAATACGGCGCCGGGCAGGCAGTG‑ACCGCGTTgggg                                                                                                                                        >  2:363126/1‑249 (MQ=255)
                                                                                                                    aTCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTGAGTCTG‑CAACGCGGTGAGACGCTGGCAAGTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCAtt                                                                                                                                      <  1:408019/251‑1 (MQ=255)
                                                                                                                                   aCGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTgg                                                                                                                        <  2:93589/250‑1 (MQ=255)
                                                                                                                                   aCGCAGGCGATACGCCAACCACCGGGATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATAGTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATGGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCTTTGGCATTGATGCGCCTGTTg                                                                                                                         <  2:420239/249‑1 (MQ=255)
                                                                                                                                        ggCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAgg                                                                                                                  <  1:457408/251‑1 (MQ=255)
                                                                                                                                          cGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGGATATTTCCTTTTTGCC‑GGACCCGCAGAAAATAGCTGCGGGCCGCCATCTATCTTTTAGTCTG‑CCAAGCGGGGGGGCGCTGGCAATTGTTGGCGAATCCGGCTCCGGGAAGGCAGCG‑GCTGCGTTGGCCTTTATGGGCCTGTTGGGACAgggg                                                                                                                >  1:205419/1‑249 (MQ=255)
                                                                                                                                          cGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcg                                                                                                                <  1:305507/251‑1 (MQ=255)
                                                                                                                                            aTACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcggg                                                                                                              <  1:141097/251‑1 (MQ=255)
                                                                                                                                              aCGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAggcgggc                                                                                                             <  1:256108/250‑1 (MQ=255)
                                                                                                                                              aCGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAAACTGAATATTGCCTTTATGCA‑GTACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTGTAGTCTG‑CAACGCGGTGAGACGCGGGCAATTGGTGGCGAATCCGGCTCCGGTAAGTCAGGG‑ACTGCGTTGGCATGGATGCGCCTGTTGGAACAggcgggg                                                                                                             <  2:454138/250‑2 (MQ=255)
                                                                                                                                                   aaCCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGTA‑GGGCCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGGCTG‑CAACGCGGTGAGACGCCGGTAATTGTTTGCGAATTCGGCTCCGGTAAGTCAGGGGAC‑GCGGTGGGCTTTCTGCGCCTGTTGGGACCGgcgggcgggttt                                                                                                       >  2:147747/1‑247 (MQ=255)
                                                                                                                                                    cccaccGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAg                                                                                                      <  2:17831/250‑1 (MQ=255)
                                                                                                                                                           gTATCTACCGTGAAAAAGGGGACACCGTTGCCACACCGTGGTAAACTTGGTGCCGGTAATGTGCTGGCGGTTGAAAATCTAAATATTTCCTTTTTGCATGAACCA‑CAGAAAATAACAGCGGGACGCAAACTCTCTTTTAGTTTG‑CAACGCGGTGAGACGCTAGCAATTGTTGGCGAATCCGGCTTCGGTTAGTCAGGG‑ACCGCGTTGGCATTGATGCGCCTTTTGGGAAAggcgggggggttttggacgtg                                                                                               >  1:427328/1‑236 (MQ=255)
                                                                                                                                                             aTCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGCATATTGCCTTTATGCA‑GGACCCGCAGAAAATAGCTGGGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGGGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGGAAGTCAGTG‑GCTGCGTTGGCCTTGATGCGCCTGTTGGAAAAggcgggggggtttggacagtggg                                                                                             >  1:284583/1‑234 (MQ=255)
                                                                                                                                                                       aaaaaGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCt                                                                                   <  1:109947/251‑1 (MQ=255)
                                                                                                                                                                           aGGGGACAGCGTTGCCACACTGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATAGTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTTGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTg                                                                               <  2:330840/251‑1 (MQ=255)
                                                                                                                                                                                acacCGTTGCCACACAGTGATGAACTTGAGGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATGGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAgcg                                                                          <  1:442132/251‑1 (MQ=255)
                                                                                                                                                                                       ggCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATTCGCCTGTTGGAACAGGCGGGCGGTTTATTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAg                                                                    <  2:167937/249‑1 (MQ=255)
                                                                                                                                                                                        gCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTc                                                                  <  1:136202/251‑1 (MQ=255)
                                                                                                                                                                                               aGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGt                                                           <  1:175568/251‑1 (MQ=255)
                                                                                                                                                                                                      aaCTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGTCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCCCAATCTCTCTTTTAGTCTG‑CAAAGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGGGGGTTTGTAAAGTGCGGTAAAATGCTGTTGGCGCGGGGCGGTCGCGGAGGGAttta                                                     >  1:153231/1‑248 (MQ=255)
                                                                                                                                                                                                      aaCTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGaa                                                    <  2:325222/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaTGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTTTGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGGGACGCTGGCAATTGTTGGCGGATCCGGCCCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTTGAACAGGCGGGGGGGTTAGTAAAGTGCGGTAAAAAACTGTTGCAACGGGGCCGGCGCGAAGTGGTTTGACTTAGCGAGCAGa                                       >  1:91221/1‑250 (MQ=255)
                                                                                                                                                                                                                              cgggtGAAAATCGGAATATGGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCTGCAATCTCTCTTTTAGGCTG‑CAACGCGGTGAGACGCTGGCAATGGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑TCTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGTCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACGTAGCGAGCAGAACGCTGCACaa                            <  2:419131/247‑1 (MQ=255)
                                                                                                                                                                                                                                gTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAAt                          <  2:153231/251‑1 (MQ=255)
                                                                                                                                                                                                                                   gAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTGAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATgc                        <  1:295225/250‑1 (MQ=255)
                                                                                                                                                                                                                                               atTGCCTTTATGCA‑GGTCCAGCAGAAAATAGCTGCGGTCCTCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGGTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTGTAGTATAGGGCTATAAAATGCTGTTGCAGCGGCGCAGTCGCTATTTGATGGAACTTAGCGAGCATAACGCTGCACAAATGCGCCATGTTCGCgg           <  1:421127/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ttGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTAGTACAGTGGGGTAAAATGCTGTTGCAGCGGCGCTGTCGCGAAGTGATTGAACTTAGCGGGCAGCAAGCTGCACAAATGTGCCCTGTTCGCggg          >  2:44754/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  gCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGaa                                                             <  1:275723/198‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTcgcgga                                                             >  2:275723/1‑196 (MQ=255)
                                                                                                                                                                                                                                                       tATGCA‑GGACCAGCAGAAAAGACCGGCGGTCAGCAATCGCTCGTTTAGTCTG‑CAACTCGTTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ATTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGGTTAGTACAGGGCGATAAAATGCGGTTGCAGCAGCGCAGTCGCGAAGTGATTGAACGTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCgtat    <  2:190616/250‑4 (MQ=255)
                                                                                                                                                                                                                                                       tATGCA‑GAACAAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCACTTGTTGGCGAATCAGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGTGGGTTTGTAAAGTGGGATACAATGCTGTTTCCGCGGCGCGGGCGCGGTGTGAATGAACTTAACGAGCAGGACGCTGCACAAATGCACCATGGTCGCgggtcgggtt   >  2:387170/1‑243 (MQ=255)
                                                                                                                                                                                                                                                        aTGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTAAAGTGCGATAAAATGCTGTTGCAGCGGCGCAGACGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACCAATGGGCCATGTTCGCGGTGCGGgtat  >  1:415097/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                             
GGAAGCCTGCGAGCGGGTAGTAATGGAAAAACTCCCTGCGCTTGGCGGTAGCGGTGGCTTAATCGCTATCGACCATGAAGGGAATGTCGCGCTACCGTTTAACACCGAAGGAATGTATCGCGCCTGGGGCTACGCAGGCGATACGCCAACCACCGGTATCTACCGTGAAAAAGGGGACACCGTTGCCACACAGTGATGAACTTGATGCCGGTAATGTGCTGGCGGTTGAAAATCTGAATATTGCCTTTATGCA‑GGACCAGCAGAAAATAGCTGCGGTCCGCAATCTCTCTTTTAGTCTG‑CAACGCGGTGAGACGCTGGCAATTGTTGGCGAATCCGGCTCCGGTAAGTCAGTG‑ACTGCGTTGGCATTGATGCGCCTGTTGGAACAGGCGGGCGGTTTAGTACAGTGCGATAAAATGCTGTTGCAGCGGCGCAGTCGCGAAGTGATTGAACTTAGCGAGCAGAACGCTGCACAAATGCGCCATGTTCGCGGTGCGGATAT  >  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: ‑