Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,808,2250AC51.3% ‑5.3 / 60.0 39intergenic (‑289/‑266)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base A (10/9);  new base C (4/16);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 4.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

AATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAC  >  NC_002947/4808082‑4808372
                                                                                                                                               |                                                                                                                                                   
aaTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgct                                                                                                                                               >  1:3931525/1‑149 (MQ=255)
aaTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgct                                                                                                                                               <  1:1685870/150‑2 (MQ=255)
  tGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttc                                                                                                                                             <  2:762843/150‑4 (MQ=255)
  tGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttc                                                                                                                                             <  2:762780/150‑4 (MQ=255)
    ggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                           >  2:1596453/1‑145 (MQ=255)
     gCGGAATTCGCGCTTTGTCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                          <  1:784678/150‑1 (MQ=38)
     gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                          <  1:2877302/150‑1 (MQ=255)
     gCGGAATTCGCGCTTTGGCGGGTAGGTGTTGGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                          <  1:731581/150‑1 (MQ=38)
      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                         <  1:1899041/150‑1 (MQ=255)
      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                         <  1:1112703/150‑1 (MQ=255)
       ggAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTTGCCGGAGGATTTTCAGCGCCTGTGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg                                                                                                                                        <  2:3554654/150‑1 (MQ=25)
        gAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggt                                                                                                                                       <  1:1633413/150‑2 (MQ=255)
          aTTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaa                                                                                                                                     <  2:1210421/150‑4 (MQ=38)
           ttCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaa                                                                                                                                    >  1:3140849/1‑146 (MQ=38)
            tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                   <  1:3301775/150‑6 (MQ=38)
            tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                   <  1:3869516/150‑6 (MQ=38)
            tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                   <  1:969151/150‑6 (MQ=38)
              gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                 >  1:1104449/1‑143 (MQ=38)
              gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                 <  1:830298/150‑8 (MQ=38)
              gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                 <  1:830314/150‑8 (MQ=38)
                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa            <  2:253149/144‑1 (MQ=255)
                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa            <  2:1213308/144‑1 (MQ=255)
                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa            <  1:946055/144‑1 (MQ=255)
                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa            <  1:2732154/144‑1 (MQ=255)
                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa            <  1:178671/144‑1 (MQ=255)
                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAAGTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa            <  1:516735/144‑1 (MQ=255)
                                                                                                                                        acccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAGCccacaaca       >  2:3934487/2‑150 (MQ=255)
                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGACACGCcaacaacaca     >  1:3328017/1‑148 (MQ=255)
                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa     >  2:3664124/1‑150 (MQ=255)
                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa     >  2:1042907/1‑150 (MQ=255)
                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa     <  2:1131202/150‑1 (MQ=255)
                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa     <  1:3425629/150‑1 (MQ=255)
                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa     <  1:3157282/150‑1 (MQ=255)
                                                                                                                                          ccgcGAAGCGACGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAAGTCCAACGCCGTCCACTCAAAATTCCAACGCCGGCCCCTCAAACATCCAACGCCGCCTCCCCTCGAAAATCAAACAACaaa     >  1:1206925/1‑150 (MQ=255)
                                                                                                                                           cgcgAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg    >  1:1279395/1‑150 (MQ=255)
                                                                                                                                           cgcgAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg    >  2:3297293/1‑150 (MQ=255)
                                                                                                                                            gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAACATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa   >  2:2910273/1‑150 (MQ=255)
                                                                                                                                            gcgAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa   >  2:1386899/1‑150 (MQ=255)
                                                                                                                                             cgAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAc  >  2:2363463/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAC  >  NC_002947/4808082‑4808372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑