Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,599,224 A→G V73A (GTC→GCC)  opgB ← phosphoglycerol transferase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,599,2240AG96.3% 82.7 / ‑6.4 27V73A (GTC→GCC) opgBphosphoglycerol transferase I
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/10);  minor base T (1/0);  total (17/10)
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.75e-01

TTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATCCTCACTAAACCGGGCCG  >  NC_000913/4598990‑4599467
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
ttATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTc                                                                                                                                                                                                                                       <  1:218022/249‑1 (MQ=255)
                                                                        gCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCTCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGGCAGCCCCAGCACAATGCCTATACCCGGCAAAAGGGAGTTGCTTGCGCCAGCACCTGTCAGGCGGGTGGTTAAGGTATAGAGATCCGCGGCGTGAATACCATCGCCAGTAAAAGAGTCGCTGGcttatagggt                                                                                                                                                             <  2:158609/251‑10 (MQ=255)
                                                                           gggCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTACCCAAAATGGTGCGGATGATGGAGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAAAGTATTTGCTGGCGCCAGCCCCGGTCAGGCTGTTGGGTAAGGTATAGAGAACCGCGTCGTTAATAACCTAGCCAGTAAAATAGTCGCttacaaaccaggaga                                                                                                                                                           >  2:244871/1‑236 (MQ=255)
                                                                                              ccAGCGCCAGTAAGAGCGCCAGAAGGCTGTAACCAAAATGGCGCGGATGCTGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGATAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCCGTAAAATAGTCGCTGGCAAACAGGGAGATATTTAAAACGACAACTAGc                                                                                                                                       >  1:216789/1‑251 (MQ=255)
                                                                                                     cAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGACCACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAACATGTATTTGCTGTCGCCCGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTCAAATAGTCGATTGCAAACAGGGTGATATTTAAACCGACAAATAGCCCCatc                                                                                                                                 >  2:119010/1‑248 (MQ=255)
                                                                                                      aGTAAGAGCGCCAGCAGGCTGTAACCAAAAGGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGGATGTGCTGGCGCCAGCACCGGTCAGGCTGTTGTTTAAGGTATATAGAACCGCTTCGTTAATACCATCGCCAGTAAAATAGTCTCTGGCaaa                                                                                                                                                                   <  1:34845/215‑1 (MQ=255)
                                                                                                                          gTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCAACGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCGGCAAacc                                                                                                           >  1:105119/1‑251 (MQ=255)
                                                                                                                              ccAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAaccacc                                                                                                        >  1:274287/1‑250 (MQ=255)
                                                                                                                                      gTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTAc                                                                                               >  2:68518/1‑251 (MQ=255)
                                                                                                                                         cGGATGATGGAGACGACGGCGCAGGATCCAGCACAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGACGATACCCAGCAAAATGTATTAGCTGGCGCCAGCACCGGTCAGGCTGATGGTTAAGGTATAGAGAACCACGTCGTTAATACCATCGCCAGTAAAATCGACGATGGCAAACAGGGTGACATTTAAAACGACCAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTc                                                                                            >  2:36374/1‑251 (MQ=255)
                                                                                                                                          ggATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAAcc                                                                                                                                                                                                            <  2:62679/138‑1 (MQ=255)
                                                                                                                                          ggATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAAcc                                                                                                                                                                                                            >  1:62679/1‑138 (MQ=255)
                                                                                                                                                    gacgacGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCAt                                                                                 <  1:249222/251‑1 (MQ=255)
                                                                                                                                                      cgacgGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGc                                                                               >  1:171702/1‑251 (MQ=255)
                                                                                                                                                                              gcACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGCGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGaa                                                        >  1:181965/1‑250 (MQ=255)
                                                                                                                                                                                                    ccatcacaCTGCCGATTCCCGGCAAAATGTTTTTGCTGGCGCCAGCACCGTTCAGGCTTTTGGTAAAGGTATAGAGACCAGCGTCGTTAATACCATCGCCATGAAAATAGGCTCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGCCCCGCTTTCCAGGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGATAGTAGTTCt                                 <  2:7975/247‑1 (MQ=255)
                                                                                                                                                                                                              gCCGATACCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGATGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCAAGTAAAATAGGCGCTGGCAAACAGGGTGATCTTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGAGTTCCGTACCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGAAGTTCTGACAACTGGc                       >  1:54808/1‑251 (MQ=255)
                                                                                                                                                                                                                    aCCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATcc                  <  2:274287/250‑1 (MQ=255)
                                                                                                                                                                                                                    aCCCGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCGTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATCCt                 <  1:115955/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAg                                                                      >  2:226807/1‑196 (MQ=255)
                                                                                                                                                                                                                      ccGGCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAg                                                                      <  1:226807/196‑1 (MQ=255)
                                                                                                                                                                                                                         gCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATCTTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCAt                                                                                 >  1:308336/1‑182 (MQ=255)
                                                                                                                                                                                                                         gCAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGGTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAATCGACAAATAGCCCCAGCACCGTTAACGGGGCTGCAAATCACCAGGTGTTACGTCCCGTTTTCCAt                                                                                 <  2:308336/182‑1 (MQ=255)
                                                                                                                                                                                                                          cAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAACAGGGCGAAAGAGCGTAGTTCTGACAACGGGCGATCCTCACTa            >  2:9862/1‑250 (MQ=255)
                                                                                                                                                                                                                          cAAAATGTATTTGCTGGCGCCAGCACCGGTCAGGCTGATGGATAGGGTCTAGAGAACCGCGTCGTTCATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGAGATATTTAAAACGACAACTAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGCAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATCCTCACTa            >  1:299380/1‑250 (MQ=255)
                                                                                                                                                                                                                              aTGTATTTGCTGGCGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCCTCGCCAGTAAAATAGTCGATGGCAAACAGGGAGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCAACAGGAGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATCCTCACTAACCCg       >  1:270210/1‑251 (MQ=255)
                                                                                                                                                                                                                                   tttGCTGGCGCCAGCACCGGTCCGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATCCTCACTAAACCGGGCCg  >  2:27486/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   
TTATAATAAGCCGCAAAGTCCGGGTCGCCGTCGCGTGACTGGGATTTCACCAGTTCCGTTATCTGACGAAATGCCGGGCTGGCGTCCACCGAACCCAGCGCCAGTAAGAGCGCCAGCAGGCTGTAACCAAAATGGTGCGGATGATGGCGACGACGGCGCAGGATCCAGCCCAGCGCACCGAACACCGCTGTCAGCCCCAGCACAATGCCGATACCCGGCAAAATGTATTTGCTGACGCCAGCACCGGTCAGGCTGTTGGTTAAGGTATAGAGAACCGCGTCGTTAATACCATCGCCAGTAAAATAGTCGCTGGCAAACAGGGTGATATTTAAAACGACAAATAGCCCCAGCACCGTTAACGTGGCTGCAAACCACCAGGTGTTACGTCCCGCTTTCCATGCGTAAATCAGCACAGAGGCGAGAAAAAGGGCGAAAGAGAGTAGTTCTGACAACGGGCGATCCTCACTAAACCGGGCCG  >  NC_000913/4598990‑4599467

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