Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_vch_tpiA4,159,0070AC24.3% 42.1 / 14.7 38E163A (GAG→GCG) oxyRoxidative and nitrosative stress transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (16/12);  new base C (6/3);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.27e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCG  >  NC_000913_3_vch_tpiA/4158862‑4159153
                                                                                                                                                 |                                                                                                                                                   
cAGACCTTTCCAAAGCTGGCAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCCTTGAAGTGCCGTTTTTTTATGTGcc                                                                                                                                                  >  2:240587/1‑149 (MQ=255)
     cTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAATGt                                                                                                                                             >  1:528167/1‑149 (MQ=255)
        tCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTTTCCTCACAACTGGACAGCAGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTACCGTTGTTTGATGAGCCAATGTTGc                                                                                                                                          >  1:575035/1‑149 (MQ=255)
            atgttgaaaaTTTATCGGCATGAAGGACAGAACTGCCAGTGACTGTCGTAACTGTACAGCGGCAAATTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGATCGAAGTTTTCATTGAAGTCCCGTTGTTTGATAAGCCAATGTTGCT‑GGc                                                                                                                                     <  2:236818/142‑1 (MQ=255)
               cTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGctat                                                                                                                                  <  2:575035/149‑1 (MQ=255)
                   aaaTGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGctatctat                                                                                                                              <  2:14340/149‑1 (MQ=255)
                                   aGCAGAGACCCACCTGTTACTGGGGCAAATGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTggg                                                                                                              <  2:400053/149‑1 (MQ=255)
                                        aGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGAGAAAGAGAGCGAAACATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAAc                                                                                                         >  1:134413/1‑149 (MQ=255)
                                            ccaccaGTTACTGGCGCAAGTGGACAGCGGCAAACTCGATTGGGTGATCCTCGGGGTGGTGAAGGAGAGCGAAGCATGCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATGACCAGTGGGCGACCtgcg                                                                                                     >  2:296329/1‑149 (MQ=255)
                                            ccaccaGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACcgcg                                                                                                     >  1:454201/1‑149 (MQ=255)
                                             caccaGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGa                                                                                                    >  2:267408/1‑149 (MQ=255)
                                              accaGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGaa                                                                                                   >  1:554036/1‑149 (MQ=255)
                                                         gcgcAACTGGACAGGGGCAAACTCGATTGCGTGATCCTCGGGTTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGAt                                                                                        <  2:80223/149‑1 (MQ=255)
                                                               cTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGa                                                                                  >  1:213406/1‑149 (MQ=255)
                                                                        ggCAAACTCGATTGCGTGCTCCGCGCGCTGGTGAAATAGACCTAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAgg                                                                         <  1:229118/149‑1 (MQ=255)
                                                                        ggCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAgg                                                                         <  2:454201/149‑1 (MQ=255)
                                                                        ggCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTAATTAAAGTACCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAgg                                                                         >  1:391761/1‑149 (MQ=255)
                                                                         gCAAACTCGATTGCGTGATCCTCGCACGGGTGAAAGAGAGCGAAGCATTGATTTAAGTGCCGTTTTTTGATGCGCCAATGTTGCT‑GGCTATCTATGAAGATGACCCGT‑GGCGAACCGCGAACGCTTACCGATGCCCGAGCTGGCAgggg                                                                       >  2:577331/1‑149 (MQ=255)
                                                                            aaCTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTAAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGaa                                                                     >  1:137773/1‑149 (MQ=255)
                                                                              cTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAATGTTGCT‑GCCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGaaaa                                                                   >  2:496824/1‑149 (MQ=255)
                                                                              cTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGaaaa                                                                   >  2:575328/1‑149 (MQ=255)
                                                                                cGATTGCGTAAGCCGCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAAc                                                                 <  2:575809/149‑1 (MQ=255)
                                                                                    tGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAActgct                                                             <  2:528167/149‑1 (MQ=255)
                                                                                       gTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATAAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGCACCGCGACTGCGTACCGATGGCCGATCTGCCAGGGGAAAAACTACTGAt                                                          >  1:290902/1‑149 (MQ=255)
                                                                                       gTGATCCTCCCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGCGCCAGTGTTGCT‑GGCTATCTATGAAGATCACACGTGGGCGAACCGCGCATGCGGACCGATGGACGATCGCGAAGGGGAAAAActgctggg                                                          >  2:154874/1‑147 (MQ=255)
                                                                                         gATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTGGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACAGATGGCCGATCTGGCAGGGAAAACAGTGCTGGTGc                                                        >  2:26682/1‑149 (MQ=255)
                                                                                                  gcTGGTGAAAGAGGTCGAAGCATTCATTGATGTGCCGTTCTTTGATTAGCCAATGTTGCT‑TGCTATCTATGAAGATCACCCGTGGGCGATCCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATg                                               <  2:7207/149‑1 (MQ=255)
                                                                                                    tGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGTTATCTATGAAGATCACCCGTGGGCGAACCGCTAATGTGTATCGATTGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGt                                             <  1:115805/149‑1 (MQ=255)
                                                                                                                      cattcattGAAGTGCTGTTTTTTGATGAGCCAATGTTGCT‑GGCTATCTATTAAGATCACCCGTGGGCGAACCGCGAATGCGTACCCATGGCTTATCTGGCATGGTAAAAACTGCTGATCCTGGAAGATGGTCACTGTTTGCGCGATCAg                           <  1:26682/149‑1 (MQ=255)
                                                                                                                        ttcattGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTGACTGTTTGCGCGATCAGGc                         >  1:25626/1‑149 (MQ=255)
                                                                                                                        ttcattGAAGTGCCGTTGTTTGATGAGCCAATATTGCT‑AAATATCTATGAAGATCACCAGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGc                         >  1:483161/1‑149 (MQ=255)
                                                                                                                               aaGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTg                                                      >  2:109358/1‑113 (MQ=255)
                                                                                                                               aaGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTg                                                      <  1:109358/113‑1 (MQ=255)
                                                                                                                                  tGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTc               <  1:267408/149‑1 (MQ=255)
                                                                                                                                      gttgttTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGACAAACCGCAAATACGTACAGATGACCGATCTGGCAGGGGAAAAACTGCTGATGCGGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGtt           >  1:79708/1‑149 (MQ=255)
                                                                                                                                       ttgttTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGttt          <  2:213406/149‑1 (MQ=255)
                                                                                                                                          tttGATGCGCCAATGTTACTAAGCT‑GCTGTGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGGTGGCCGAACTGGCAGGGGAAAAGCTGGTGATGCAGGAAGATGATCACTGTTGGCGCGAACAGGCAATGGGTTTCTGtttgta       >  2:4386/1‑146 (MQ=255)
                                                                                                                                          tttGATGAGCCAATGTTGCT‑GGCTATATATGAAGATGACCAGTGGGCGAACCGCGAATGCGTACCAATGGCGGATCTGGCAGGGGAAAAACTGCGGATGCTGGAAGATGGGCACTGTTTGGGCGATCAGGCAAGGGGTTGCCGTTTTGa       >  2:425732/1‑149 (MQ=255)
                                                                                                                                               tgaGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTTCTGCTTTAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCg  <  1:240587/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CAGACCTTTCCAAAGCTGGAAATGTATCTGCATGAAGCACAGACCCACCAGTTACTGGCGCAACTGGACAGCGGCAAACTCGATTGCGTGATCCTCGCGCTGGTGAAAGAGAGCGAAGCATTCATTGAAGTGCCGTTGTTTGATGAGCCAATGTTGCT‑GGCTATCTATGAAGATCACCCGTGGGCGAACCGCGAATGCGTACCGATGGCCGATCTGGCAGGGGAAAAACTGCTGATGCTGGAAGATGGTCACTGTTTGCGCGATCAGGCAATGGGTTTCTGTTTTGAAGCCG  >  NC_000913_3_vch_tpiA/4158862‑4159153

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