Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,4170CG25.0% 37.4 / 10.6 20W30S (TGG→TCG) nuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base C (8/7);  new base G (2/3);  total (10/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTG  >  NC_000913/2391279‑2391553
                                                                                                                                          |                                                                                                                                        
tAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCACCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAtca                                                                                                                                 <  4:94846/148‑3 (MQ=255)
                 aTCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAaccacca                                                                                                               >  7:60351/1‑149 (MQ=255)
                                  gcgTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATc                                                                                              <  8:60351/149‑1 (MQ=255)
                                            cATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACgg                                                                                      >  8:23077/1‑147 (MQ=255)
                                                              aacaaACCAGAGCGAAACCAGCGCCGCGTGAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCCCCGTCAAGCCGACGATCAGCAACGGTAGCAGGGCGATCAGGgttt                                                                  >  3:10196/1‑149 (MQ=255)
                                                                        agCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTAt                                                        >  3:62607/1‑149 (MQ=255)
                                                                         gCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTAtt                                                       <  7:23077/149‑1 (MQ=255)
                                                                           gAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACgg                                                                                      <  8:5052/116‑1 (MQ=255)
                                                                                          ttAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACgg                                      <  4:62607/149‑1 (MQ=255)
                                                                                               cccAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAgg                                                                      <  2:63109/112‑1 (MQ=255)
                                                                                               cccAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAgg                                                                      >  1:63109/1‑112 (MQ=255)
                                                                                                         gagaGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAg                                                            >  2:65739/1‑112 (MQ=255)
                                                                                                         gagaGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAg                                                            <  1:65739/112‑1 (MQ=255)
                                                                                                                       gAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGttt                                                                   <  5:153310/91‑1 (MQ=255)
                                                                                                                       gAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGttt                                                                   >  6:153310/1‑91 (MQ=255)
                                                                                                                         gctctccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACa             <  2:31634/136‑1 (MQ=255)
                                                                                                                         gctctccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACa             >  1:31634/8‑143 (MQ=255)
                                                                                                                             tccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgct   >  5:27522/4‑149 (MQ=255)
                                                                                                                              tGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgctg  >  4:50314/1‑149 (MQ=255)
                                                                                                                               caTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgctg  <  6:27522/147‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                        
TAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTG  >  NC_000913/2391279‑2391553

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