Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,1070AG43.9% 55.0 / 10.9 41I530V (ATT→GTT) gsiAglutathione ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (3/20);  new base G (18/0);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 7.52e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAG  >  NC_000913/868973‑869253
                                                                                                                                      |                                                                                                                                                             
tCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGTCATCGTGTgg                                                                                                                                                >  1:10946/1‑150 (MQ=255)
      ccGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGTGTGGCgggg                                                                                                                                          >  1:277489/1‑147 (MQ=255)
       cgacgaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTAGGCATCGTGTGGCgggg                                                                                                                                          >  1:1010975/1‑146 (MQ=255)
           gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGGGCGGGTTATTCATCTTGGGGCGGGGgtgt                                                                                                                                      >  2:384614/1‑148 (MQ=255)
           gaAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGAGGGAGGGGGGGGTAGTTTTCGTGGGGCgggggtttt                                                                                                                                     >  2:250867/1‑142 (MQ=255)
                 gTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGc                                                                                                                       <  2:799117/151‑1 (MQ=255)
                    tCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCaa                                                                                                                     <  2:832585/150‑1 (MQ=255)
                      ggCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAAtt                                                                                                                  <  2:36984/151‑1 (MQ=255)
                              aTGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTTGTCATCGTGTGGGGGGGGTGTTTTTCCGGGGAAATTTTTTaa                                                                                                            >  1:1053015/1‑150 (MQ=255)
                                 tttCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGTGGGGTTGGCATCCTGTGGGGGGGGAGTATGTTCGGGGAAATTTTTTGAAtt                                                                                                         >  1:982797/1‑150 (MQ=255)
                                         cgcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGgc                                                                                               <  2:789895/151‑1 (MQ=255)
                                          gcgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGgcg                                                                                              <  1:1195/151‑1 (MQ=255)
                                           cgGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGGGGGTTAGTCATCGTGTGGGGGGGGGGTATGTCCGGGCAAAATTTGTAAATTGTCCCAccc                                                                                                >  1:362417/1‑147 (MQ=255)
                                            gggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGgcgcg                                                                                            <  2:416767/151‑1 (MQ=255)
                                             ggCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGgcgcgc                                                                                          <  2:32248/151‑1 (MQ=255)
                                                      tCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGGGGGGGGGGATTTCCGGGGAAATTTTTGAAATTGTGCCCCGCGCGCgg                                                                                       >  1:263637/1‑144 (MQ=255)
                                                             ggtgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGAGCGCGGTCGTTGAAAAccc                                                                        <  2:919428/148‑1 (MQ=255)
                                                              ttgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCCGGTTAGTCATCGTGTGGCGGGGGATGTTTTCCGGGGCAATTTTTTAAATTGGGCCCCGCGGCGCGGCCTCCGAAAAcc                                                                         >  2:1019531/1‑149 (MQ=255)
                                                               tgCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGTGGGGGGGTTAGTAATCGTGTGGCGGGGGGGTGTTTTCGGGCAAATTTTTTAAATTGTGCCCCGCGCGCGGGGTTCGGAAAacc                                                                        >  1:789053/1‑148 (MQ=255)
                                                                gCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCgca                                                                     <  2:697458/150‑1 (MQ=255)
                                                                   cGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCgcagca                                                                  <  2:282145/150‑1 (MQ=255)
                                                                        tcCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCATCGTGTGGCGGGGGTGTATTTCGGGGGAAATTTTTAAAATTGTTCCCGGCCCGCGGGCTTTCGAAAACCCCAGCAAcc                                                               >  1:88883/1‑148 (MQ=255)
                                                                        tcCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGGGGTAGAGGGGGTTGTTCAACGTGTGGCGGGGGGGTTTCTTCGGGGAAATGGTTTAAATTGTTCCCCGGGGCGGGGGCTTGGAAAACCCCCGCATcc                                                               >  2:10718/1‑148 (MQ=255)
                                                                         ccAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTAGTCCTCGTGTGGCGGGGGTGTTTATCCGGGGAAATGTTTTAAATTGGGCACCGCGCGCCGGCTTTGgaaaaaccccgcacccct                                                             >  1:403054/1‑135 (MQ=255)
                                                                           aGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGGGGGTTAGTCATCGGGTGGCGGGGGTGTATTTCGGGGGAAATTTTTTAAATTTTGCCCCGCGCGCGGGCTTTGCAAAAcccccccaac                                                                 >  1:601515/1‑137 (MQ=255)
                                                                             cGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATAcgc                                                       <  1:802523/151‑1 (MQ=255)
                                                                             cGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATAcgc                                                       <  2:523152/151‑1 (MQ=255)
                                                                                gATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGGTTTGTCATTTTGTGGCGGGGGTTTTATTTCGGGCAAATTGTGTGAATTGGTCACGGGGGGCCGGCCTTGGAAAAACCCAGCAACCTtttaatc                                                       >  2:679096/1‑145 (MQ=255)
                                                                                           tGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTggc                                          <  2:1055051/150‑1 (MQ=255)
                                                                                            gCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTggcg                                         <  1:650644/150‑1 (MQ=255)
                                                                                             cGTATCTGTTTATCTCCCACGATATGGCGGTGGGAGAGCGGGTTTGTCATCGTGTGGCGGGGGTTTTTTTCGGGGGAAATTTTTTAATTTGGGTCCCGCGCGCGGGCTTTTGAAAACCCCCGGAACCCTTTAACCGGTAAAttcttggg                                          >  2:233840/1‑143 (MQ=255)
                                                                                                    gTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTcc                                 <  1:805718/150‑1 (MQ=255)
                                                                                                    gTTTATCTCCCACGATATGGCGGTGGGAGAGCCGGTTCGTCATCGTGTTGCGGGCGTGGATTCCCGGGCAAATTTTTGAAATTGGTCACCGCGCGCGGGCCTTCGAAAACCCCCCCAACCCTATACCCGTAAAATaactggggccggttc                                 >  2:443930/1‑136 (MQ=255)
                                                                                                         tctcCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTc                             <  2:677724/149‑1 (MQ=255)
                                                                                                            cccACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCg                            <  1:278286/147‑1 (MQ=255)
                                                                                                                    atGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTccc                <  2:344627/151‑1 (MQ=255)
                                                                                                                     tGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTccc                >  2:489071/1‑150 (MQ=255)
                                                                                                                              tAGAGCGGATTAGTCATCGTGTGGCGGCGATGGATCCCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTTTACGCGTAAATTACTGGCGGCCGTTTCGGTCGCTGAACCGTCCCGACAaccacc      >  2:629031/1‑151 (MQ=255)
                                                                                                                               agagCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGAcc      <  2:217756/150‑1 (MQ=255)
                                                                                                                                  gCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAg  <  1:486925/151‑1 (MQ=255)
                                                                                                                                    ggATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGGCGCTGAACCGTCCCGACAACGACCGCAg  >  1:331845/1‑149 (MQ=255)
                                                                                                                                      |                                                                                                                                                             
TCATTGCCGACGAAGCCGTTTCGGCGCTGGATGTTTCTATTCGCGGGCAGATTATCAACTTGTTGCTCGATCTCCAGCGTGATTTCGGCATTGCGTATCTGTTTATCTCCCACGATATGGCGGTGGTAGAGCGGATTAGTCATCGTGTGGCGGTGATGTATCTCGGGCAAATTGTTGAAATTGGTCCACGGCGCGCGGTCTTCGAAAACCCGCAGCATCCTTATACGCGTAAATTACTGGCGGCAGTTCCGGTCGCTGAACCGTCCCGACAACGACCGCAG  >  NC_000913/868973‑869253

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