Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,117,270 G→A P274L (CCT→CTT)  glpF ← glycerol facilitator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,117,2700GA94.1% 41.1 / NA 18P274L (CCT→CTT) glpFglycerol facilitator
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (8/8);  minor base C (1/0);  total (9/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAAC  >  NC_000913/4117054‑4117495
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
gggTGGTCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGGGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCacaacaa                                                                                                                                                                                                     >  1:368846/1‑246 (MQ=255)
       ccAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAagag                                                                                                                                                                                                                                <  1:28689/213‑1 (MQ=255)
                              ntACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGAGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGttt                                                                                                                                                                   >  1:267926/2‑251 (MQ=255)
                                                         gCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATg                                                                                                                                                                                  <  1:149754/209‑1 (MQ=255)
                                                                        gCGACACGCTAATGATATTGGCATCGTGATCCATGACGACCGCGCGGGAGCTGGTGGTGCTCTGGTCGAGCGCAACGAGATAAGCTTTGTCAGTCATGTTTAATTGGCACGTAGGCATATTACAGCTAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCCGTTCTAACGCAACACGGATACCGCAAGGCAAATTGTGTCTAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAggg                                                                                                                           <  1:208538/249‑1 (MQ=255)
                                                                          gACACGCTAATGATATTGGCATCGTGATCCCTTACGACGGCGCGGGAGCTGGTGGTGCCCTGGCCGTGCGAAACGATATATTTTTTTTCAGTCATGTTTAATTGGCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCTTTTTCTTCca                                                                                                                                                                                                           >  1:170067/1‑167 (MQ=255)
                                                                          gACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCACGGCAAATGGCGACCAATCAGTTTGCGGTAGGCaaa                                                                                                                                                       >  1:250181/1‑219 (MQ=255)
                                                                               gCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAATCTGTTTGTTCTTAAAGAGTTGTGTTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGc                                                                                                                  <  1:295270/251‑1 (MQ=255)
                                                                                                     tCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTgg                                                                                                                                                                                                                          >  1:259931/1‑125 (MQ=255)
                                                                                                           aCGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGATCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTgc                                                                                      <  1:434472/251‑1 (MQ=255)
                                                                                                                    cgGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTACCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTgg                                                                                                                                                                                                                          >  1:386851/1‑110 (MQ=255)
                                                                                                                                                         atTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTGCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAAGGCACCTACAATCGCGCCAACGATAAGGCCGAAAAGCTGCACCAGGAAGTAAGTAATGtc                                                                                           <  1:371981/200‑1 (MQ=255)
                                                                                                                                                                 tCAGTCATGTTTAATTGTCACGTAGTAATATTACAGCGCAGCTTTTTGTTCTGAACGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCACGGCAAATg                                                                                                                                                                                  >  1:162052/1‑105 (MQ=255)
                                                                                                                                                                                tNNNNNGTANTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGagag                                                                                                                             <  2:267926/143‑4 (MQ=255)
                                                                                                                                                                                      gTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCACGGAAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCAc                                                                                                                                                  >  1:190177/1‑116 (MQ=255)
                                                                                                                                                                                            atatTACAGCGAAGCTTTTTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAgt                                                                                                                                                                                        <  1:125075/72‑2 (MQ=255)
                                                                                                                                                                                              atTACNGCGAANCNNTTTGTTCTGACAGAGTTGTNGTTTCCTTTTCTTCCNCAACNCAGATNTCNcn                                                                                                                                                                                           >  2:194552/1‑66 (MQ=255)
                                                                                                                                                                                              atTACAGCGAAGCTTTGTGTTCTGAAAGAGTTGTGGTTTCCTTTTCTTCCACAACACAGAGATCGCAAGGCAAATGGCGACCAAGCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGc                                                                                                                  <  1:194552/140‑1 (MQ=255)
                                                                                                                                                                                                tACAGCGAAGTTTTTTGTTCTGAAAGAGTGGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCTATCAATCAGTTTGCGGGAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATTTCTCTGCCGCCGGTAAAGTCGATATTGCCCTAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAATTCACGCGCTGGGTTCATGGCAAAAc  <  1:294413/250‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                 
GGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTGGTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGCGAAGCTTTTTGTTCTGAAGGAGTTGTGGTTTCCTTTTCTTCCACAACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCGAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCATGGCAAAAC  >  NC_000913/4117054‑4117495

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