Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913447,2710CT27.6% 42.8 / 16.1 29A145A (GCG→GCA) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base T (3/5);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATA  >  NC_000913/447126‑447412
                                                                                                                                                 |                                                                                                                                             
gAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGgcggcg                                                                                                                                              >  5:51117/1‑147 (MQ=255)
        gccgccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCgg                                                                                                                                     <  4:96135/148‑1 (MQ=255)
            ccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGc                                                                                                                                <  6:51117/149‑1 (MQ=255)
            ccTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCGgtgatc                                                                                                                                <  2:74726/149‑3 (MQ=255)
                gTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGaa                                                                                                                             <  5:120269/148‑1 (MQ=255)
                      cTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCTCCGGCGCTGCGCCGCCGgtgatc                                                                                                                                >  6:16668/1‑137 (MQ=255)
                      cTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCTCCGGCGCTGCGCCGCCGgtgatc                                                                                                                                <  5:16668/139‑3 (MQ=255)
                       ttAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAc                                                                                                                            <  8:123169/142‑1 (MQ=255)
                       ttAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAAc                                                                                                                            >  7:123169/1‑142 (MQ=255)
                           aGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTg                                                                                                                 <  8:32187/149‑1 (MQ=255)
                                                     tGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTa                                                                                       >  3:36819/1‑149 (MQ=255)
                                                               gCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTAACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCtata                                                                              >  2:22971/1‑148 (MQ=255)
                                                                ccACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCGgtga                                                                                                                                  <  8:37965/95‑1 (MQ=255)
                                                                ccACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGCTGCGCCGCCGgtga                                                                                                                                  >  7:37965/1‑95 (MQ=255)
                                                                  acaGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGCTCTCCGGCGCTGCGCCGCCGgtgatc                                                                                                                                >  3:92954/1‑93 (MQ=37)
                                                                  acaGGCTGAAGAGAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGCTCTCCGGCGCTGCGCCGCCGgtgatc                                                                                                                                <  4:92954/95‑3 (MQ=25)
                                                                           aGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATa                                                                  >  3:40061/1‑148 (MQ=255)
                                                                             aTAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTCACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAAc                                                               <  1:22971/149‑1 (MQ=255)
                                                                                 ccAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCg                                                                                                              >  4:132510/1‑98 (MQ=255)
                                                                                 ccAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCg                                                                                                              <  3:132510/98‑1 (MQ=255)
                                                                                                   cgcTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTc                                                                                            <  6:41168/98‑1 (MQ=255)
                                                                                                            aCTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAgg                                                                                   >  4:106436/1‑98 (MQ=255)
                                                                                                            aCTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAgg                                                                                   <  3:106436/98‑1 (MQ=255)
                                                                                                                 ccGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGccc                                                     <  5:128569/123‑1 (MQ=255)
                                                                                                                 ccGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGccc                                                     >  6:128569/1‑123 (MQ=255)
                                                                                                                  cGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGa                          <  4:36819/149‑1 (MQ=255)
                                                                                                                  cGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGa                          >  1:142927/1‑149 (MQ=255)
                                                                                                                                tAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACca            <  4:40061/149‑1 (MQ=255)
                                                                                                                                    cGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCAcag         >  5:330/1‑148 (MQ=255)
                                                                                                                                           cggcggCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATa  >  4:92607/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GAATGTTTGCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATA  >  NC_000913/447126‑447412

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