Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,708,8440CG41.2% ‑2.5 / 10.6 18A642A (GCC→GCGrsxCSoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base C (8/2);  new base G (4/3);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTG  >  NC_000913/1708695‑1708992
                                                                                                                                                     |                                                                                                                                                    
gcagCTATTGCACGTGCCAAAGCGCGCAATCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGCCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGcc                                                                                                                                                      <  1:584561/147‑1 (MQ=255)
      aTTCCCCGGGCCCAACCGCGCAAGCTGGCACACCAACAGCCTACTCCGGAAGCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCa                                                                                                                                                <  1:896174/150‑1 (MQ=11)
                                            ggatAATTCCGAAGCAGAACCACATGTTGATCCTCGCAAACCTGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc                                                                                                          <  1:1438788/147‑1 (MQ=11)
                                              cTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGAGCAAAAGCCCTAAAACACAAAAAGCAACAGGCTCATGAAGGGCAAGAAGAAAAAGAGGATAAGCGCAAAGCGGAAGGCGAAGCCg                                                                                                        >  2:1847840/1‑150 (MQ=2)
                                                                  cAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCAAAAGCCCAAAATCACGCACAGCAAAAGGATAATTACGTAACAGCAGAACAGGTGGATACACTCAAAACAGCAGTAGAAGCCGCTATTGCCCGAGCCaaatca                                                                                    >  2:278554/1‑147 (MQ=0)
                                                                        gATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCGGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGACGACGCGGCTATTGCCCGCGCCCAGGCCCACaaa                                                                              >  1:15357/1‑150 (MQ=14)
                                                                                  aaGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGCGCCAAAGCGCGCAAGCGGGAACAGc                                                                    >  1:1515791/1‑150 (MQ=255)
                                                                                            cGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCGGAGCCAGAACAACAGGTAGATCCGCGCAACGCCGCCCTCGAAGCCGCGATTGCCCGGGCCATAGCGAGCAACCGGGACCCGCAACAGGCTaa                                                          >  1:1966346/1‑150 (MQ=11)
                                                                                                        tGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCGGCCGTAGCAGCCGCTATTGCCCGCGCCCAGGCCAAAAACCCCGCACAGCAACCGCCTAATGCAGAGCCaaa                                              >  1:707662/1‑148 (MQ=14)
                                                                                                          cccGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCagaag                                            >  1:1451200/1‑150 (MQ=255)
                                                                                                                  aaaGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTTGATCCGCGCAAAGCCGCAGTAGAAGCCGCTATTGCCCGCGCACACGCCCCAAAACACGCACAGCAACAGGATAAAACGGAGCCCGAAAAActgtta                                    >  1:76928/1‑145 (MQ=9)
                                                                                                                         gCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTAGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGCGCCCAAGCGAGCAAACCCGAACAGCAGCCGGGTAATCAGGAGCCAGAAGAACAGGTTGATAcgc                             >  2:781409/1‑150 (MQ=14)
                                                                                                                               tGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                       <  2:95927/102‑1 (MQ=255)
                                                                                                                               tGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAAc                                                                       >  1:95927/1‑102 (MQ=12)
                                                                                                                                       cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                         <  2:1565146/92‑1 (MQ=16)
                                                                                                                                       cAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                         >  1:1565146/1‑92 (MQ=16)
                                                                                                                                              cTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAACCGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTAGAACCGCGCAAAGCTGCCGTCGACGCgg        >  1:165543/1‑150 (MQ=255)
                                                                                                                                                  tGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                         <  1:349682/81‑1 (MQ=255)
                                                                                                                                                  tGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGa                                                                         >  2:349682/1‑81 (MQ=255)
                                                                                                                                                    ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCGGCTATTGCCCGCGCCAAAGCCCCCAAAGGCGACCAGCAACCGGGTAATGAGGAGCAATAAGAACAGGTAGCTCAGCCCAAATCGGCTGTAGCAGCGGCAATTg  >  2:366102/1‑150 (MQ=9)
                                                                                                                                                    ccGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCCAAGCCAACAAACGGGAACAGCAACATGCTAATACGGAGCCAGAAGAACAGGTAGATCCGCGCCAAGCTGCCGTCGACGCGGCTATTg  >  2:1002644/1‑150 (MQ=14)
                                                                                                                                                     |                                                                                                                                                    
GCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACGGGAACAGCAACCGGCTAATGCGGAGCCAGAAGAACAGGTTGATCCGCGCAAAGCTGCCGTCGAAGCGGCTATTG  >  NC_000913/1708695‑1708992

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