Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,126,899 G→C 13.3% A190G (GCC→GGC)  PP_2744 ← ribose‑phosphate pyrophosphokinase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,126,8990GC13.3% 85.6 / 6.1 30A190G (GCC→GGC) PP_2744ribose‑phosphate pyrophosphokinase family protein
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base C (2/2);  total (16/14)
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.90e-01

GCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAG  >  NC_002947/3126750‑3127036
                                                                                                                                                     |                                                                                                                                         
gcGTCTGCCCCGTGCTGACCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGGTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAggg                                                                                                                                           >  2:275096/1‑150 (MQ=255)
         ccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGCCCTg                                                                                                                                       <  2:406224/145‑1 (MQ=255)
         ccGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGCCCTg                                                                                                                                       >  1:406224/1‑145 (MQ=255)
                  tCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCgg                                                                                                                         >  1:642001/1‑150 (MQ=255)
                         gttgAATACGATGACGGTGCTGGCTGCCACGTTACCGACCAGTTGTTCGCCGGTCAGGCCTGGCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGttt                                                                                                                  >  2:592849/4‑150 (MQ=255)
                           cgAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTTCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGa                                                                                                                <  1:325930/150‑1 (MQ=255)
                            gAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGAc                                                                                                               <  1:185029/150‑1 (MQ=255)
                                                 tgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCggc                                                                                               <  2:515018/145‑2 (MQ=255)
                                                 tgccACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCggc                                                                                               >  1:515018/1‑144 (MQ=255)
                                                     aCGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGAcc                                                                                      <  2:642001/150‑1 (MQ=255)
                                                           ccGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGAcc                                                                                <  1:498597/150‑1 (MQ=255)
                                                                                    ctgcctgcTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCa                                                       >  1:181061/1‑150 (MQ=255)
                                                                                       cctgcTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGc                                                    <  1:361767/150‑1 (MQ=255)
                                                                                        ctgctgGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGAGATAGGGCCTGGAGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGGGACACGACAGTGAGAGTGTGGTGACGGAGTAGCGGGGCGAgggg                                                   >  1:321629/1‑148 (MQ=255)
                                                                                           ctgGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTg                                                <  1:667931/150‑1 (MQ=255)
                                                                                           ctgGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTg                                                <  1:275096/150‑1 (MQ=255)
                                                                                                         ccATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTggg                                  >  1:736017/1‑150 (MQ=255)
                                                                                                         ccATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTggg                                  >  2:537946/1‑150 (MQ=255)
                                                                                                         ccATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTggg                                  >  2:537898/1‑150 (MQ=255)
                                                                                                                      cGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGtt                     >  2:537300/1‑150 (MQ=255)
                                                                                                                       gagaCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTc                    >  2:490731/1‑150 (MQ=255)
                                                                                                                          aCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCg                                                                                                 <  2:313845/70‑1 (MQ=255)
                                                                                                                          aCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCg                                                                                                 >  1:313845/1‑70 (MQ=255)
                                                                                                                          aCCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAg                 <  2:362497/150‑1 (MQ=255)
                                                                                                                              gTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCGCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCg             >  2:434839/1‑150 (MQ=255)
                                                                                                                                    ccAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCAt                             <  2:368843/128‑1 (MQ=255)
                                                                                                                                    ccAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCAt                             >  1:368843/1‑128 (MQ=255)
                                                                                                                                    ccAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTc       >  2:263104/1‑150 (MQ=255)
                                                                                                                                         agGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAg  <  2:181061/150‑1 (MQ=255)
                                                                                                                                                ccAGGCCCTGGCGCACCTGCT‑GGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTa   >  2:224490/1‑141 (MQ=255)
                                                                                                                                                ccAGGCCCTGGCGCACCTGCT‑GGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTa   <  1:224490/141‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                         
GCGTCTGCCCCGTGCTGATCAGGTCGTCGAATACGATGACGGTGCTGCCTGCCACGTTACCGACCAGGTGTTCGCCGCTCAGGCCTGCCTGCTGGCGGTGTTTCTCCATCATGGCCACCGAGACCGGTCGCGCCAGCAGGTGCGCCAGGGCCTGGCGAAACTGCTCGGCGCGTTTGACGCCACCGATATCCGGTGACACGACCGTGACCGCCTGCTCACCGAGTAGCGGGGCCAGGCGTTGGGCGAATAGCTGGGTGCATTGCAGGTTCCAGCTCGGGATTCGGTAG  >  NC_002947/3126750‑3127036

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