Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,974,131 G→T 12.5% P127P (CCC→CCA)  ubiG ← ubiquinone biosynthesis O‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,974,1310GT12.5% 102.0 / 2.8 32P127P (CCC→CCA) ubiGubiquinone biosynthesis O‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/16);  new base T (2/2);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

CTTCAGGATGTACTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATCTGCCGGTATTCCACCTGCACGCCGG  >  NC_002947/1974002‑1974237
                                                                                                                                 |                                                                                                          
cTTCAGGATGTACTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACg                                                                                                      >  2:163980/1‑136 (MQ=255)
cTTCAGGATGTACTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACg                                                                                                      >  1:355503/1‑136 (MQ=255)
     tgttGTTCTCGTCGCAGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGTAGAACACCTGGCCGCCCGGCTTGACCATGCTGTAGCAGGCGCGAATCTCCGAAGACGGGTCGGGCACGTGCTc                                                                                                 <  2:666878/133‑1 (MQ=255)
     ggATGTACTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTc                                                                                                 <  1:163980/136‑1 (MQ=255)
     ggATGTACTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTc                                                                                                 <  2:316975/136‑1 (MQ=255)
            cTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCAtt                                                                                          >  2:480278/1‑136 (MQ=255)
            cTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCAtt                                                                                          >  2:247384/1‑136 (MQ=255)
            cTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCAtt                                                                                          >  2:460797/1‑136 (MQ=255)
            cTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCAtt                                                                                          >  2:483725/1‑136 (MQ=255)
              cGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTc                                                                                        >  2:178243/1‑136 (MQ=255)
                      acgaTGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCcaggcagg                                                                                <  2:267143/136‑1 (MQ=255)
                        gaTGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTg                                                                              <  2:312826/136‑1 (MQ=255)
                           ggCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGAcc                                                                           <  1:178243/136‑1 (MQ=255)
                                gcagGTTGGCGTTGGGGTTGCGGTTGATGGTGGAGACGAACAGCTGGCGGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  2:626445/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  1:123026/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  2:557947/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  2:355503/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  2:582023/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  1:483725/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  1:480278/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  1:460797/136‑1 (MQ=255)
                                gcagGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTc                                                                      <  1:311465/136‑1 (MQ=255)
                                                gTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCTGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCAtt                                                      >  2:145817/1‑136 (MQ=255)
                                                gTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCTGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCAtt                                                      >  2:132009/1‑136 (MQ=255)
                                                gTTGCGGTTGATGGTGGAGAAGAACACATGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCAtt                                                      >  1:412102/1‑136 (MQ=255)
                                                      gTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCg                                                <  2:16139/136‑1 (MQ=255)
                                                                           cTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATc                           >  2:519236/1‑136 (MQ=255)
                                                                           cTGGCCGCCCCGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATc                           >  1:625452/1‑136 (MQ=255)
                                                                              gccgccCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCTGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATCTGc                        <  1:132009/136‑1 (MQ=255)
                                                                              gccgccCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCTGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATCTGc                        <  1:145817/136‑1 (MQ=255)
                                                                                                   gTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATCTGCCGGTATTCCACCTGCACGCCg   >  2:576740/1‑136 (MQ=255)
                                                                                                    tAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATCTGCCGGTATTCCACCTGCACGCCgg  >  2:75108/1‑136 (MQ=255)
                                                                                                                                 |                                                                                                          
CTTCAGGATGTACTCGGCGCCGACGATGGCCAGCAGGTAGGCCTTGGGGTTGCGGTTGATGGTGGAGAAGAACACCTGGCCGCCCGGCTTGACCATGCGGTAGCAGGCGCGAATCACCGAAGACGGGTCGGGCACGTGCTCGAGCATTTCCAGGCAGGTGACCACGTCGAACTGTTCAGGCATTTCCTCGGCCAGGGCTTCGGCGGTGATCTGCCGGTATTCCACCTGCACGCCGG  >  NC_002947/1974002‑1974237

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