Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,808,1740C.53.3% ‑3.3 / 31.1 15intergenic (‑283/‑272)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base C (0/7);  new base . (2/6);  total (2/13)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCC  >  NC_002947/4808052‑4808297
                                                                                                                          |                                                                                                                           
cgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGcgccgcgc                                                                                                                           <  1:3807507/125‑3 (MQ=255)
   tttGGCGGGTAGGTGTTTGGCGGGAGTTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgcgcggg                                                                                                                        <  2:3807508/125‑3 (MQ=255)
    ttGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccctc                                                                                                                       <  1:3807510/125‑5 (MQ=255)
    ttGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccctc                                                                                                                       <  2:3807509/125‑5 (MQ=255)
     tGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGccgccccctca                                                                                                                      <  2:3807511/125‑6 (MQ=255)
         gggTAGGTGTTTGGCGGGAGGTATTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgc                                                                                                                  <  2:3807513/125‑1 (MQ=255)
          ggTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCcgcgcgg                                                                                                                     <  2:3807515/121‑4 (MQ=255)
           gTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgctt                                                                                                                <  2:3807516/125‑3 (MQ=255)
              ggTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCgccgcc                                                                                                                            <  2:3807518/110‑1 (MQ=255)
               gtgtTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                            <  1:3807519/125‑1 (MQ=255)
                tgtTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttcg                                                                                                              <  2:3807521/122‑5 (MQ=255)
                tgtTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                           <  1:3807520/125‑1 (MQ=255)
                                                    cGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTTGCCGGAGGATTTTCAGCGCCTGTGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCGTCCGGCAAACACCTCaa                                                          >  1:3807565/1‑124 (MQ=2)
                                                     gCCGGGGGATTTGCAGCGCCTGGGAGGTGTTTGCCGGAGGATTTTCAGCGCCTGTGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                             <  2:3807532/120‑1 (MQ=2)
                                                        gggggATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTc                                                     >  1:3807603/1‑125 (MQ=17)
                                                         ggggATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGCGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTcc                                                    <  2:3807534/125‑1 (MQ=17)
                                                           ggATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCaa                                                  <  2:3807535/125‑1 (MQ=17)
                                                              ttttCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGc                                               <  2:3807536/125‑1 (MQ=17)
                                                                 tCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑GCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGt                                            <  1:3807537/125‑1 (MQ=18)
                                                                   aGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGt                                            <  2:3807538/123‑1 (MQ=17)
                                                                    gcgcCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑AGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTcc                                          <  2:3807539/124‑1 (MQ=17)
                                                                                                                    ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCa       >  1:3807611/7‑125 (MQ=255)
                                                                                                                    ggtttacccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCa       >  2:3807627/7‑125 (MQ=255)
                                                                                                                          cccgcGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcgcc  >  1:3807594/1‑124 (MQ=255)
                                                                                                                          |                                                                                                                           
CGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCC  >  NC_002947/4808052‑4808297

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