Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,013,255 A→T 17.2% intergenic (+83/+48) prfC → / ← potF‑I peptide chain release factor 3/putrescine‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,013,2550AT17.2% 73.7 / 9.7 35intergenic (+83/+48)prfC/potF‑Ipeptide chain release factor 3/putrescine‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (22/7);  new base T (3/3);  total (25/10)
Fisher's exact test for biased strand distribution p-value = 3.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

GCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCTCGTCGC  >  NC_002947/1013113‑1013403
                                                                                                                                              |                                                                                                                                                    
gcgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCCGCCTAGCGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCtg                                                                                                                                                  <  8:36447/147‑1 (MQ=18)
 cgcTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTgg                                                                                                                                               <  2:180958/149‑1 (MQ=25)
      aTGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTACCAAAGGAATAACCAGACCTGGAGGGCATGTTGTACCGGTGGGAGCg                                                                                                                                          >  2:209678/1‑149 (MQ=38)
         gAAGAGCGCTGGCCGGACATCAAGTTACGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACCCCCCCTACCACAGGAATCACCAGACCTGGAGGGGAGGTGATAGCGGTGGGCGCGGGc                                                                                                                                       >  4:77828/1‑149 (MQ=38)
                                                   cACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCAGTGCGAGCGGGCACGCCcgctcccac                                                                                                                         >  1:267132/1‑115 (MQ=37)
                                                             gCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATAc                                                                                   >  7:9551/1‑149 (MQ=255)
                                                                                                 cACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCg                                               >  2:28762/1‑149 (MQ=255)
                                                                                                  aCGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCttaga                                              >  4:200917/1‑146 (MQ=255)
                                                                                                    gCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCg                                                      >  2:277248/1‑139 (MQ=255)
                                                                                                    gCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCg                                                      <  1:277248/139‑1 (MQ=255)
                                                                                                    gCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCg                                            >  5:120709/1‑149 (MQ=255)
                                                                                                             cacaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGcc                                    >  8:276253/1‑148 (MQ=255)
                                                                                                              acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGaa                                                              >  2:354826/1‑121 (MQ=18)
                                                                                                              acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGaa                                                              <  1:354826/121‑1 (MQ=18)
                                                                                                              acaGGTATCACATGACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                                  >  8:290649/1‑149 (MQ=25)
                                                                                                              acaGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                                  >  3:366124/1‑149 (MQ=255)
                                                                                                                aGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTc                                  >  5:77685/1‑147 (MQ=255)
                                                                                                                  gtaTCACATGACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAgc                              >  6:145798/3‑149 (MQ=37)
                                                                                                                        cctgACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTc                                                <  3:79022/122‑1 (MQ=255)
                                                                                                                        ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa                         >  4:117870/1‑148 (MQ=255)
                                                                                                                        ccAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGAGCCCGCGAACCGCCGCTAGGCGGCCCCATTTCAATTACCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGACCATCTTc                                                >  4:79022/1‑125 (MQ=255)
                                                                                                                        catgACCTCCAGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa                         >  4:62290/4‑148 (MQ=38)
                                                                                                                         cAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGcc                                    >  2:158642/1‑136 (MQ=255)
                                                                                                                         cAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGcc                                    <  1:158642/136‑1 (MQ=255)
                                                                                                                         cAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa                         >  3:230551/1‑147 (MQ=255)
                                                                                                                         atgACCTCCAGGTCTGGTGATTCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATa                         >  3:163091/3‑147 (MQ=37)
                                                                                                                           gACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGc                     <  6:77685/149‑1 (MQ=255)
                                                                                                                           gACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGc                     >  3:169244/1‑149 (MQ=255)
                                                                                                                                 gAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGc                 <  8:9551/147‑1 (MQ=255)
                                                                                                                                 caGGTCTGGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGc                 <  3:62290/146‑1 (MQ=255)
                                                                                                                                      cATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGcc                <  5:339897/143‑1 (MQ=255)
                                                                                                                                      cATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGcc                >  6:339897/1‑143 (MQ=255)
                                                                                                                                         gtgATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCtc       >  5:213036/1‑149 (MQ=255)
                                                                                                                                            aTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCtcgtc    >  2:285859/1‑149 (MQ=255)
                                                                                                                                            aTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCtcgtc    >  7:346888/1‑149 (MQ=255)
                                                                                                                                              aCCTGTGGGAGCGGGCGTGCCCGCGACGGGCCGCTAGGCGGCCCCAGTGCAATTTCCGGGTTTGATACTGGTACAGCTTCGCGTACGAATACGGTCGCTCTTCGCCGGCCGCGCCTCCAGCGCATACAGCTGGCGCATCCCCTCGgcgc  >  4:261638/1‑149 (MQ=37)
                                                                                                                                              aCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCtcgtc    <  4:71186/147‑1 (MQ=255)
                                                                                                                                              aCCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCtcgtc    >  3:71186/1‑147 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GCGCTGATGGAAGAGCGCTGGCCGGACATCAAGTTCCGCGCTACCCGCGAGCACCATTAAGGCCTGAAATCGGGGCTGCTTTGCAGCCCTTCGCGGGCACGCCCGCTCCCACAGGAATCACCAGACCTGGAGGTCATGTGATACCTGTGGGAGCGGGCGTGCCCGCGAAGGGCCGCTAGGCGGCCCCATTTCAATTTCCGGCTTTGATACTGGTCCAGATTCGCGTACGAATACGGTCGATCTTCGCCGGCATCGCCTCCAGCGCATACAGCTTGCCCATCACCTCGTCGC  >  NC_002947/1013113‑1013403

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: