Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913872,2400GT77.4% 40.2 / 13.5 31A117A (GCG→GCT) gsiDglutathione ABC transporter membrane subunit GsiD
Reads supporting (aligned to +/- strand):  ref base G (5/2);  new base T (16/8);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.78e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGACCGCCGCATTATTCGTTATTTTATTGATTGTGGTGGCCATTTTTGCACGCTGGATCGCTCCCTATGACGCCGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTC  >  NC_000913/872009‑872482
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
gATGACCGCCGCATTATTCGTTATTTTATTGCTTGTGGTGTCCATTTTTGCACGCTGGATCGCTCCCTATGACGCCGAACCTTATTTTGATTATGACAATCTGAATAACGGACCATCTTTGCAGCACTGGTTTGGCGTCGATTCCCTGGGGCGTTACATTTTCATCCGTGTCCTGGTTGGTGCGCAAATCCCGCTGGCGGCGGGCGAGTTTTCCGTGTTTATCGGGGAGGCTATCTGGACGTTGCTgaggt                                                                                                                                                                                                                                 >  2:229250/1‑247 (MQ=255)
 acgcccgccgCATTATTCGTTATTTTATTGATTGTGGTGGCCATTTTTGCACGCTGGATCGCTCCCTATGACGCCGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCt                                                                                                                                                                                                                                 <  1:229250/246‑1 (MQ=255)
               attCGTTATTTTATTGATTGTGGTGGCCATTTTTGCACGCTGGATCGCTCCCTATGACGCCGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGAtatt                                                                                                                                                                                                                  >  1:266120/1‑251 (MQ=255)
                   gTTATTTTATTGATTGTGGTGGCCATTTTTGCACGCTGGATCGCTCCCTATGACGCCGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGCTATTATg                                                                                                                                                                                                               >  1:389786/1‑250 (MQ=255)
                                                tGCACGCTGGATCGCTCCCTATGACGCCGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATgcgc                                                                                                                                                                                  <  1:78598/250‑1 (MQ=255)
                                                  cACGCTGGATCGCTCCCTATGACGACGAAAATTCTTTTGATTATGACACTCTGAATAACGGACCTTCTTAGCAGCACTGGTTTGGCGTCGATTCAATGGGGCGTGACATTTTCAGCCGTGTCCTGGATGGGGAGCAAGTATAGCTGGAGGCGGGCGTGGTTGCAGCGATTATCGGTGAGGCTAACGGGACGAGGCGGGGCTTGATCGATGGGTATAATGAAGGCTGGAGGGGACGGCTGAtgaagagcatg                                                                                                                                                                               >  2:369318/1‑241 (MQ=255)
                                                     gCTGGATCGCTCCCTATGACGACGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTCTCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGGTCGGCTGGTCATGCGCCTTTGc                                                                                                                                                                            >  2:240814/1‑251 (MQ=255)
                                                                ccatATGACGCCGATAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTGTGCCGTGTTTATCGGGGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTg                                                                                                                                                                   <  1:62966/246‑1 (MQ=255)
                                                                            gAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTAACTGGTGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGGTCGGCTATCGGGACGTTGCTGGGCTTGCTCGGTGGATATTATGAAGGATGGGGGGATCGGCTGTACATACGCATTTGCGCTGGTCTGTTTGCCTTCCCggg                                                                                                                                                     >  2:221224/1‑251 (MQ=255)
                                                                              aacttattTTGATTCTGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGCTTCACTGAGGCGCGACATTTTCAGCCGTATCCTGGTTGGCGCGCAAATCTCGCTGGCGGCGGCCCTCTTTGCCGTGTTTCTCGGTGCGCCGATCGGGACGTTGCTGTGCTTGATAGCTGGATCTTTTGAAGGCTGGTGGGTTCGGCTGATCACGCGCATTGGTGATGTGATGTTTGCCTTCCCGGGt                                                                                                                                                    >  2:118528/4‑250 (MQ=255)
                                                                               aaTTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTAt                                                                                                                                                  >  1:396677/1‑251 (MQ=255)
                                                                                                            cGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGgtgt                                                                                                                     <  1:157432/251‑1 (MQ=255)
                                                                                                             ggACCTTCTTTGCAGCACTGGTTTGGCGTCGAATCCCTGGGGCGTGACCTTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGCCGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGtt                                                                                                                    >  1:240177/1‑251 (MQ=255)
                                                                                                                     gttgaagcaCTGGTTTGGCGTCTATTCATGGGGGCGTGATATTTTCAGCCGTGTCGTGGTTGGTGCGCATATCTCGCTGGCGGCTTGCGTGGTTGCCGTGTTTATCGGTTCGGCTATCGGGACGGTGCTGGGCTTGCTCGCTGTATATTATGATTGCTTGTGGGATCTGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTGCCCGGGTATTTTACTGGCTATCGGTGAGGTTGCGGTGTTGGGAAGc                                                                                                             <  2:45467/245‑1 (MQ=255)
                                                                                                                       tGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTCTCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGc                                                                                                          >  1:418548/1‑251 (MQ=255)
                                                                                                                         cagcaCTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCAt                                                                                                        >  1:28269/1‑251 (MQ=255)
                                                                                                                               gggTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGGTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCGTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTaa                                                                                                  <  2:143558/250‑1 (MQ=255)
                                                                                                                                                                 tCAGCCGTGGCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTTGCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTAGTTTACTGGCGATCGCGGTTGTTGCGGTGTAGGGAAGCGGCATTGCTAACGTGATTAGTGCGGTCGCCTTTTTTTCCATccc                                                                 <  2:418548/250‑1 (MQ=255)
                                                                                                                                                                               gttggtGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGGATTTTACTGGCGATCGCTGTTGTTGCGGGTTTGGGGAGCGGCATTGATAACGTGGTTATTGCCGTCGCCATTTTTTCCATCCCCCCGTTTTGCCGCCt                                                  >  2:115931/1‑251 (MQ=255)
                                                                                                                                                                                ttggtGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTg                                                 >  1:273602/1‑251 (MQ=255)
                                                                                                                                                                                 tggtGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTgg                                                >  1:58753/1‑251 (MQ=255)
                                                                                                                                                                                     gcgcAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTgcg                                            >  1:409817/1‑251 (MQ=255)
                                                                                                                                                                                         aaatcgggcgggcggcGGGCGTGTTTGCCGTTTTGATGGGGGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCGTGTGGTATCGGCTGATCTTGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGAGCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTGTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGc                                        <  2:348264/241‑1 (MQ=255)
                                                                                                                                                                                         aaaTCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGAATTTTACTGGCGATCGCTGTTGTTGCGGGGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGTGTTTGACCGCCTGGTGGGCgg                                         >  1:227122/1‑250 (MQ=255)
                                                                                                                                                                                         aaaTCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGACCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGGGTTGGGAAGCGGCATTGCTAACGGGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGc                                        >  1:227115/1‑251 (MQ=255)
                                                                                                                                                                                                                 ttGCCGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCa                 <  1:240814/250‑1 (MQ=255)
                                                                                                                                                                                                                  ttccGTGTTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGCTCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCAGCGTTTGCCCGACTGTTGCGCGGCAACACGCTGGTGTTGAAACAGCaa                >  2:287626/3‑250 (MQ=255)
                                                                                                                                                                                                                   gCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGACGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGAGGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGCAACAGCAAAc              >  2:142121/1‑251 (MQ=255)
                                                                                                                                                                                                                      gtgtTTATCGGTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCAtttt                                                                          <  1:273712/188‑1 (MQ=255)
                                                                                                                                                                                                                      gtgtTTATCGGAGCGGCTATCAGGACGCTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATAGGCTGATAATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCAtttt                                                                          >  2:273712/1‑188 (MQ=255)
                                                                                                                                                                                                                                gTGCGGCTATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGGATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTTCAGTCGCCATTTTTTCCATCCCCTCGTTTGCCCGCCTGGTGCGCGGCACCACGCTGGTGTTGAAACAGCAACCCTTGATTGAGTc  >  2:73875/1‑250 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
GATGACCGCCGCATTATTCGTTATTTTATTGATTGTGGTGGCCATTTTTGCACGCTGGATCGCTCCCTATGACGCCGAAAATTATTTTGATTATGACAATCTGAATAACGGACCTTCTTTGCAGCACTGGTTTGGCGTCGATTCACTGGGGCGTGACATTTTCAGCCGTGTCCTGGTTGGTGCGCAAATCTCGCTGGCGGCGGGCGTGTTTGCCGTGTTTATCGGTGCGGCGATCGGGACGTTGCTGGGCTTGCTCGCTGGATATTATGAAGGCTGGTGGGATCGGCTGATCATGCGCATTTGCGATGTGCTGTTTGCCTTCCCGGGTATTTTACTGGCGATCGCTGTTGTTGCGGTGTTGGGAAGCGGCATTGCTAACGTGATTATTGCAGTCGCCATTTTTTCCATCCCCGCGTTTGCCCGCCTGGTGCGCGGCAACACGCTGGTGTTGAAACAGCAAACCTTTATTGAGTC  >  NC_000913/872009‑872482

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