Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,934,687 G→T 13.3% P236T (CCG→ACG)  visC ← ubiquinone biosynthesis oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,934,6870GT13.3% 90.6 / 3.1 30P236T (CCG→ACG) visCubiquinone biosynthesis oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (12/14);  new base T (2/2);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GTGTCGCTGGCGCAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCAACGCACGCTGG  >  NC_002947/5934553‑5934818
                                                                                                                                      |                                                                                                                                   
gtgtCGCTGGCGCAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCggg                                                                                                                                    >  2:68448/1‑136 (MQ=255)
gtgtCGCTGGCGCAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCggg                                                                                                                                    >  2:615090/1‑136 (MQ=255)
   gcGCTGGCGCAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGggtg                                                                                                                                 <  1:384456/135‑1 (MQ=255)
          cgcAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGAcc                                                                                                                          <  1:162686/136‑1 (MQ=255)
          cgcAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGAcc                                                                                                                          <  1:296775/136‑1 (MQ=255)
          cgcAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGAcc                                                                                                                          <  1:545570/136‑1 (MQ=255)
                ggCCCGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCAACAACTGCTCGGCAAGCTCCGGGGTGGTCGACCAGACGa                                                                                                                    <  1:531476/136‑1 (MQ=255)
                ggCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGa                                                                                                                    <  1:335326/136‑1 (MQ=255)
                                ggTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCctg                                                                                                    <  2:115593/136‑1 (MQ=255)
                                        tGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCAt                                                                                            >  2:524716/1‑136 (MQ=255)
                                        tGAAGCACCTCGCCCACACGCCCTTCGAAGACCCGCTCAATGCCCTGCAAGAACGCCGCTACAACCATAGCCATAAACTGCACGCCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGGCAt                                                                                            >  1:81002/1‑136 (MQ=255)
                                             cACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCgg                                                                                       <  2:234812/136‑1 (MQ=255)
                                                   gCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGTGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAAc                                                                                 >  1:40265/1‑136 (MQ=255)
                                                   gCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGTGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAAc                                                                                 >  1:330542/1‑136 (MQ=255)
                                                                tCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGccc                                                                    >  1:373092/1‑135 (MQ=255)
                                                                 cGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAg                                                                   >  2:28019/1‑136 (MQ=255)
                                                                              cAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCa                                                      >  1:100170/1‑136 (MQ=255)
                                                                              cAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCa                                                      >  1:235157/1‑136 (MQ=255)
                                                                                        aGGAAGGCCGCTTCATCCATCGCCATCTACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGATCACCAGTCCTGCTGGCCATCGCGGGTCAACGGGAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCt                                            <  2:64863/136‑1 (MQ=255)
                                                                                         ggAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTg                                           >  2:653178/1‑136 (MQ=255)
                                                                                         ggAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTg                                           >  2:331575/1‑136 (MQ=255)
                                                                                         ggAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGACAGCGGACCTTCGTCAGTGAAGCGCTg                                           >  1:604441/1‑136 (MQ=255)
                                                                                                  cTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGt                                  >  2:608204/1‑136 (MQ=255)
                                                                                                    tCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCCGGCCGTgg                                >  2:623534/1‑136 (MQ=255)
                                                                                                        ccatcGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTg                            <  2:370749/136‑1 (MQ=255)
                                                                                                          atcGCCATCAACTGCTCGGCATGCTCCGTGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGt                          <  2:330542/136‑1 (MQ=255)
                                                                                                          atcGCCATCAACTGCTCGGCATGCTCCGTGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGt                          <  2:40265/136‑1 (MQ=255)
                                                                                                                        cTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCaa            >  1:144571/1‑136 (MQ=255)
                                                                                                                          cGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCAacg          <  1:68448/136‑1 (MQ=255)
                                                                                                                           ggCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCAacgc         <  1:89524/136‑1 (MQ=255)
                                                                                                                           ggCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCAacgc         <  1:567102/136‑1 (MQ=255)
                                                                                                                                  tCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCAACGCACGCTgg  <  1:431789/136‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GTGTCGCTGGCGCAGCGGCACGCACACCCGTGGGTCGGCCTGAAGCACCTCGCCCAGACGCCCTTCGAAGGCCCGCTCCAGGGCCTGCAGGAAGGCCGCTTCATCCATCGCCATCAACTGCTCGGCATGCTCCGGGGTGGTCGACCAGACGATCGAGCACCAGTCCTGCTGGCCATCGCGGGTCAACGGCAGGAACGCCAGCGGACCTTCGTCAGTGAAGCGCTGCCAGGCCGTGGCCTGGTGCCCGGCGCTGCAACGCACGCTGG  >  NC_002947/5934553‑5934818

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: ‑