Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 556,997 A→G 10.8% T16A (ACG→GCG)  rpsN → 30S ribosomal protein S14

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947556,9970AG10.8% 115.5 / 4.9 46T16A (ACG→GCG) rpsN30S ribosomal protein S14
Reads supporting (aligned to +/- strand):  ref base A (21/20);  new base G (3/2);  total (24/22)
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.53e-01

GACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTGC  >  NC_002947/556855‑557132
                                                                                                                                              |                                                                                                                                         
gACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGACGCTTCAGCTCGCgg                                                                                                                                      >  2:318079/1‑146 (MQ=255)
gACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGACGCTTCAGCTCGCgg                                                                                                                                      <  1:318079/146‑1 (MQ=255)
            accaccACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAGAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGc                                                                                                                         <  2:37275/147‑1 (MQ=255)
                ccacTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTaagaa                                                                                                                   <  1:159944/149‑1 (MQ=255)
                 cacTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTaagaag                                                                                                                  >  8:51698/1‑149 (MQ=255)
                 cacTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGACCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTaag                                                                                                                     >  5:51657/1‑146 (MQ=255)
                      cTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGc                                                                                                                         >  3:243052/1‑137 (MQ=255)
                      cTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGc                                                                                                                         <  4:243052/137‑1 (MQ=255)
                         gTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGTGTAGCAAAGTTCGCTAAGAAGCGTgc                                                                                                             <  1:171654/147‑1 (MQ=255)
                          ttCGGACGACGAAGGGCGCGCTCTGCTGCTTTCTTTCAAGTTCCCGTTCCGAAAGTTATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTgctga                                                                                                          <  7:114600/148‑1 (MQ=255)
                             ggACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTgctgagct                                                                                                       >  8:69855/1‑148 (MQ=255)
                             ggACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTACCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTCGCGGtttttca                                                                                                                               >  5:183073/1‑118 (MQ=38)
                             ggACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTACCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGTGAGCTGACGCTTCAGCTCGCGGtttttca                                                                                                                               <  6:183073/124‑7 (MQ=38)
                                cgacgaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGaaa                                                                                                   <  1:36336/149‑1 (MQ=255)
                                    gaAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACAGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGa                                                                                               >  1:50430/1‑149 (MQ=255)
                                      aGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGAc                                                                                              >  2:110523/1‑148 (MQ=255)
                                                 tgctgcGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAAc                                                                                    <  8:110697/147‑1 (MQ=255)
                                                    ggcGTGCTTTCAAGTTCCCCTTCCCCAACTGATTGCAGTAGGAAAATGGCCAAGACGACCAGGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGaa                                                                               <  8:103253/148‑1 (MQ=255)
                                                       gTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGc                                                                            >  7:226090/1‑149 (MQ=255)
                                                           tttCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCaaa                                                                                                                               >  3:29062/1‑94 (MQ=255)
                                                           tttCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCaaa                                                                                                                               <  4:29062/94‑1 (MQ=255)
                                                               aaGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCtctc                                                                       >  1:138026/1‑146 (MQ=255)
                                                                 gTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCagaa                                                                  <  7:51698/149‑1 (MQ=255)
                                                                   tCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCagaa                                                                  >  8:175107/1‑147 (MQ=255)
                                                                   tCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCagaa                                                                  >  1:107894/1‑147 (MQ=255)
                                                                        ttCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGtt                                                           >  6:103854/1‑149 (MQ=255)
                                                                            gCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGc                                                       >  1:317059/1‑149 (MQ=255)
                                                                                    ttGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTg                                               <  7:69855/149‑1 (MQ=255)
                                                                                       gAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCg                                            <  5:279851/149‑1 (MQ=255)
                                                                                          tAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGctc                                         >  1:148401/1‑149 (MQ=255)
                                                                                              aaaaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCa                                     <  5:103854/149‑1 (MQ=255)
                                                                                               aaaTGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCa                                     <  2:148401/148‑1 (MQ=255)
                                                                                                    gCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAag                                 >  7:243785/1‑147 (MQ=255)
                                                                                                      cAATAAGAGCATGAAAAACCGCTAGCTGAAGCGTCAGCTCACCACCAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAagcag                              <  1:110523/149‑1 (MQ=255)
                                                                                                       aagaagAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGcc                            <  2:50430/149‑1 (MQ=255)
                                                                                                         gaagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCAc                          <  2:317059/149‑1 (MQ=255)
                                                                                                         gaagaGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCAc                          >  2:263779/1‑149 (MQ=255)
                                                                                                                  gAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGc                                      <  2:56323/128‑1 (MQ=255)
                                                                                                                  gAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGc                                      >  1:56323/1‑128 (MQ=255)
                                                                                                                  gAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGc                   >  4:124932/1‑147 (MQ=255)
                                                                                                                  gAAAAACCGCGAGCTGAAGCGACAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCAGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGACCGTTTCGCTGCCGTTGTCGCTCTTCAGAAGCAGCCACGCGACGcaa                 >  7:73509/1‑147 (MQ=255)
                                                                                                                   aaaaaCCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTa                 <  2:107894/148‑1 (MQ=255)
                                                                                                                    aaaaCCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAgc               <  2:138026/149‑1 (MQ=255)
                                                                                                                     aaaCCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAgcg              >  1:210143/1‑149 (MQ=255)
                                                                                                                         cgtgAGCTGACGCTTCAGCTCGCGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTc          >  1:282142/4‑149 (MQ=255)
                                                                                                                                 gAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTgc  >  1:195314/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                         
GACATCACCCTGACCACCACTGCTCGTTCGGACGACGAAGGCCGCGCTCTGCTGCGTGCTTTCAAGTTCCCGTTCCGCAACTGATTGGAGTAGGAAAATGGCCAAGAAGAGCATGAAAAACCGCGAGCTGAAGCGTCAGCTCACGGTAGCAAAGTTCGCTAAGAAGCGTGCTGAGCTGAAAGCGACCATCGTCAACCTGAACGCCTCTCCAGAAGAGCGTTTCGCTGCCGTTGTCGCTCTGCAGAAGCAGCCACGCGACGCTAGCGCCTCGCGCCTGC  >  NC_002947/556855‑557132

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