Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913447,2750GC27.6% 50.2 / 21.8 29A144G (GCT→GGT) cyoEprotoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base G (10/11);  new base C (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.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATA  >  NC_000913/447134‑447412
                                                                                                                                             |                                                                                                                                         
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                          >  1:142927/1‑149 (MQ=255)
                                                                                                          cGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGa                          <  4:36819/149‑1 (MQ=255)
                                                                                                                        tAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACca            <  4:40061/149‑1 (MQ=255)
                                                                                                                            cGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCAcag         >  5:330/1‑148 (MQ=255)
                                                                                                                                   cggcggCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATa  >  4:92607/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                         
GCCGCCTGGTAATCCTTAAAGCGGAAAATGGCGATGGCATAGGAGTGAGGCATCTGCCACAGGCTGAAGATAGCCAGCAGGATCGCTGCGCCGCTATCGAACTCACCGGTTACCGCACAGTAGCCGATCACCGGCGGCGCAGCGCCGGAGAGCGAACCAATCAACGTGCCGTAGACAGAGTGGCGTTTCATGTACAGGCTATAAACGCCGACATAAACCACAAAGCCCATCACCCCCAGCCAGCAGGCCAGCGGATTCGCGCCAAACCACAGCAGCATA  >  NC_000913/447134‑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: