Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,1930AG28.6% 34.0 / 19.8 28F105L (TTC→CTC) nuoNNADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base A (11/9);  new base G (2/6);  total (13/15)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGAAAGCGTAACCGACCAGGCCAAACAGCGGCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATA  >  NC_000913/2391062‑2391326
                                                                                                                                   |                                                                                                                                     
gCGGAAAGCGTAACCGACCAGGCCAAACAGCGGCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATcc                                                                                                                                >  1:183678/1‑139 (MQ=255)
    aaaGCGTAACCGACCAGGCCAAACAGCGGCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCttgt                                                                                                                            <  1:73741/139‑1 (MQ=255)
                               ggCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGt                                                                                                 >  2:491771/1‑139 (MQ=255)
                               ggCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGt                                                                                                 >  2:202425/1‑139 (MQ=255)
                                gCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCttgttg                                                                                                                          <  1:237509/113‑1 (MQ=255)
                                gCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCttgttg                                                                                                                          >  2:237509/1‑113 (MQ=255)
                                gCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTa                                                                                                <  2:154259/139‑1 (MQ=255)
                                gCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTa                                                                                                <  2:551932/139‑1 (MQ=255)
                                                gATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAgg                                                                                              <  1:132506/125‑1 (MQ=255)
                                                gATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAgg                                                                                              >  2:132506/1‑125 (MQ=255)
                                                    ccGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGa                                                                            <  2:522003/139‑1 (MQ=255)
                                                       aGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGc                                                                         >  2:154103/1‑139 (MQ=255)
                                                       aGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGc                                                                         >  2:364539/1‑139 (MQ=255)
                                                            cAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGcc                                                                    <  2:86486/139‑1 (MQ=255)
                                                            cAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGcc                                                                    <  1:491771/139‑1 (MQ=255)
                                                              gagaCGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCcaa                                                                  >  2:471774/1‑139 (MQ=255)
                                                                   gCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATa                                                             <  1:364539/139‑1 (MQ=255)
                                                                    ccAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATAc                                                            >  1:467194/1‑139 (MQ=255)
                                                                    ccAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATAc                                                            <  2:45787/139‑1 (MQ=255)
                                                                        aTGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGc                                                        >  1:238023/1‑139 (MQ=255)
                                                                         tGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGcc                                                       <  1:153038/139‑1 (MQ=255)
                                                                         tGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGcc                                                       <  2:183678/139‑1 (MQ=255)
                                                                             tGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGgt                                                   >  1:514864/1‑139 (MQ=255)
                                                                              ggCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGgtg                                                  >  2:262515/1‑139 (MQ=255)
                                                                                                                aTTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAgcg                <  2:238023/139‑1 (MQ=255)
                                                                                                                  tAACACCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAgcggc              <  2:83556/139‑1 (MQ=255)
                                                                                                                    acacCAACAGGTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGt            <  2:467194/139‑1 (MQ=255)
                                                                                                                              gTAGAGCTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATa  >  1:558002/1‑139 (MQ=255)
                                                                                                                                   |                                                                                                                                     
GCGGAAAGCGTAACCGACCAGGCCAAACAGCGGCAAAGAGATCAGTTCGATACCGAGGAACAGAGACGCCAGATGGTTGGCATTCGCCAGCAGGATCCCGCCCAGCGCGGCAATTAACACCAACAGGTAGAACTCATCCTTGTTGTCGTTATAGCCTTCAAGCCACGGGTAGGCGAAAGTACAGGTGGCGAGGCTCGCCAACAATACCAGCCCGGTGTAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATA  >  NC_000913/2391062‑2391326

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