Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,967,689 C→G 16.7% G115A (GGC→GCC) ‡ yciK ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,967,6890CG16.7% 66.6 / 9.7 30G115A (GGC→GCC) ‡yciKputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (12/13);  new base G (2/3);  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 = 5.20e-01

GCCCCAGTAAGCCCGGCCCTTGCGCCCCACGCTGCTGCTGGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGTTTCCAGGTTGAACGGGATGACCAC  >  NC_002947/1967541‑1967837
                                                                                                                                                    |                                                                                                                                                    
gcCCCAGTAAGCCCGGCCCTTGCGCCCCACGCTGCTGCTGGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGcc                                                                                                                                                     >  1:555314/1‑150 (MQ=255)
      gTAAGCCCGGCCCTTGCGCCCCACGCTGCTGCTGGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTTCATCACGCGCATGAAGTTGTCGCCCGACaa                                                                                                                                               <  1:128222/150‑1 (MQ=255)
            ccGGCCCTTGCGCCCCACGCTGCTGCTGGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGGGCGACAACttcat                                                                                                                                         <  1:597138/150‑5 (MQ=255)
                                      tgGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCGCGACAACTGc                                                                                                                                           >  2:63487/1‑122 (MQ=255)
                                      tgGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGc                                                                                                                                           <  1:63487/122‑1 (MQ=255)
                                              acgaccGAAGCGTCTTCCGACAGCTACAGCAGCGGCAGCAGAGTGCTGGTGATCATGAGGGTAGCATGGACGTGGAGGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCGCCAGCGGCGTGCACGGGCCGATGATAGACGCAttg                                                                                                       >  2:51126/4‑150 (MQ=255)
                                              accaccGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCAttg                                                                                                       >  1:242105/1‑150 (MQ=255)
                                                caccGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCAttgtt                                                                                                     <  2:125931/150‑1 (MQ=255)
                                                  ccGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCGGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGa                                                                                                   <  1:757974/150‑1 (MQ=255)
                                                     aaGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGCCAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCa                                                                                                <  2:312226/150‑1 (MQ=255)
                                                              ccGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTa                                                                                       <  2:555314/150‑1 (MQ=255)
                                                                          gcagcGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAggc                                                                                    <  1:267859/141‑1 (MQ=255)
                                                                          gcagcGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAggc                                                                                    >  2:267859/1‑141 (MQ=255)
                                                                           cagcGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATGACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCACTGTTGAGCAGGCCATCTAGGCGGCCGAACTg                                                                          <  1:557022/150‑1 (MQ=255)
                                                                                  agcagGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGACGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCg                                                                   >  1:699109/1‑150 (MQ=255)
                                                                                                     tGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACt                                                                                                                                             >  2:137562/1‑57 (MQ=255)
                                                                                                     tGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACt                                                                                                                                             <  1:137562/57‑1 (MQ=255)
                                                                                                             gCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTgg                                        >  1:745347/1‑150 (MQ=255)
                                                                                                               aTTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGa                                                                  <  1:92123/122‑1 (MQ=255)
                                                                                                               aTTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGa                                                                  >  2:92123/1‑122 (MQ=255)
                                                                                                                     gTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCtcgt                                                >  1:530909/1‑134 (MQ=255)
                                                                                                                     gTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCtcgt                                                <  2:530909/134‑1 (MQ=255)
                                                                                                                                        gagcaGTTGTCGGGCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGtt                         >  2:326616/5‑138 (MQ=255)
                                                                                                                                        gagcaGTTGTCGGGCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGtt                         <  1:326616/134‑1 (MQ=255)
                                                                                                                                        aTGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGc                             >  1:547211/1‑134 (MQ=255)
                                                                                                                                        aTGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGc                             <  2:547211/134‑1 (MQ=255)
                                                                                                                                         tGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCGTTCGAGCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGTTTCCAGGTTGAACg            <  2:353199/150‑1 (MQ=255)
                                                                                                                                          gAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGTTTCCAGGTTGAACgg           >  1:565073/1‑150 (MQ=255)
                                                                                                                                             gTTGTCGGGCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGc                                  >  1:62819/1‑124 (MQ=255)
                                                                                                                                             gTTGTCGGGCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGc                                  <  2:62819/124‑1 (MQ=255)
                                                                                                                                                   gCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGTTTCCAGGTTGAACGGGATGaccac  >  1:195181/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCCCAGTAAGCCCGGCCCTTGCGCCCCACGCTGCTGCTGGTGAACACCACCGAAGCGTCTTCCGACAGCTTCAGCAGCGGCAGCAGGGTGCTGGTGAGCATGAAGGTGGCATTGACGTTGATGTGCATCACGCGCATGAAGTTGTCGCCCGACAACTGCTCCAGCGGCGTGCGCGGGCCGATGATCGACGCATTGTTGAGCAGGCCATCTAGGCGGCCGAACTGGCCTTCGATCATGGCCGCCAGCTCGTCGTACTGGTGTGGCAGGGCGGTTTCCAGGTTGAACGGGATGACCAC  >  NC_002947/1967541‑1967837

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