Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,117,109 T→G 100% Q38P (CAA→CCA)  glpK ← glycerol kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,117,1090TG94.5% 95.5 / ‑1.9 37Q38P (CAA→CCA) glpKglycerol kinase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base G (17/18);  total (19/18)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGC  >  NC_000913/4116963‑4117250
                                                                                                                                                  |                                                                                                                                              
tGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGTCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGGTTGGc                                                                                                                                              >  1:74505/1‑149 (MQ=255)
              ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTg                                                                                                                                >  1:286643/1‑149 (MQ=255)
              ttCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTg                                                                                                                                <  1:708347/149‑1 (MQ=255)
               tCGTAATTCCGATAGCGGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGc                                                                                                                               <  2:120770/149‑1 (MQ=255)
                    aTACCGATAGCTGCAATTTGATAGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTTGGGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGacac                                                                                                                          >  1:570577/1‑149 (MQ=255)
                      aCCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTTGAGCTTTGGGTGGCCCAGATTTCCTTTGGGTCGTGTTCTACCCAACCTGGCTTTGGGTTGATTGGCTCAAATTCGCGCTGCGACACGc                                                                                                                        <  1:475216/149‑1 (MQ=255)
                         gcttgccgCAATTTGATCGGAAATGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTaa                                                                                                                     <  2:152765/142‑1 (MQ=255)
                          agaGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAAt                                                                                                                    <  2:53290/147‑1 (MQ=255)
                               tGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCCCAAATTCGCGCTGCGACACGCTAATGatat                                                                                                               >  1:674411/1‑149 (MQ=255)
                                gCAATTCGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCCCAAATTCGCGCTGCGACACGCCAAAGAAAtt                                                                                                              >  2:232550/1‑149 (MQ=255)
                                 cAATTGGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTg                                                                                                             <  1:76730/149‑1 (MQ=255)
                                        aTCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGt                                                                                                      <  1:133792/149‑1 (MQ=255)
                                          aggAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCGCAGATTTAGATTGGGTCGTGTTTTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATcagga                                                                                                    >  1:376692/2‑145 (MQ=255)
                                             aaCTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATcc                                                                                                 <  1:563336/149‑1 (MQ=255)
                                                tGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCAtt                                                                                              <  1:523547/149‑1 (MQ=255)
                                                     tCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGTGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACgcc                                                                                         >  1:59377/1‑147 (MQ=255)
                                                           ttCGCCAGCACCTCTACCCGCCGGGGGCTTTGGGTGGCCCAGGTTTCCATTGGGTCGTGTTCTTCCCATCCTGGTTTTGGGTAGAATGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACcgcgcg                                                                                   <  2:121266/149‑1 (MQ=255)
                                                               gcAGCACTTCGACCAGCGTGGAGCTTTGCGGTGCCCACATTTCCATTGTGTCTTGTTCTACCCAACCTTGTTTTGGGTAGATTGGCTCAAATTCTCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                               <  2:441043/148‑1 (MQ=255)
                                                               ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCATATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACTCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAg                                                                               <  2:441018/149‑1 (MQ=255)
                                                               ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAACTTCGCGCTGCGACACGCTACTGATATTGGCATCGTGATCCATTACGACCGGCCGGGAg                                                                               >  1:423094/1‑149 (MQ=255)
                                                               ccAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATGCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGTGCGGGAg                                                                               >  1:423028/1‑149 (MQ=255)
                                                                                 tGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGGCCTGGTCg                                                            >  1:414333/1‑149 (MQ=255)
                                                                                   gAGCTTTGGGTGGCCCAGATTTCCAGTGGGTCGTGTTCTACCCAACCAGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGTGCTGGTG‑GTCGCCTGGTCGac                                                          >  1:409837/1‑148 (MQ=255)
                                                                                        ttggttggCCCAGCTTTCCATTGGTCCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCCCAAAGTCGCGCTGCGACACCCTAATGATATTGGCTTCGTGATCCAGTACGACCGCGCGGGAGCTGGTG‑GTGGCCTGGTCGAGCGCaa                                                     <  2:94929/149‑1 (MQ=255)
                                                                                             tGGCCCAGATTTTCATTGGGTCGTGTTCTACCCAAACTGGTTTTGGGTAGCTTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCAGCGTGAACCATTACGACAGCGCGGGAGCTGGGT‑GTGACCTGGTCTAGCGCAACGagt                                                >  2:9803/1‑147 (MQ=255)
                                                                                                         cgctggggTCGTGTTCTGCCCAACCTGGTTTTGGGTAGTTTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCa                                    <  2:168120/144‑1 (MQ=255)
                                                                                                            ttGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACCACCGAG‑GGGGTCTACTGTGGGAACTGGTAGATCGCAACGATATATTTTTTTTCAGTc                                 >  1:530404/1‑149 (MQ=255)
                                                                                                                 tCGTGTTCTACCCAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGGCATgat                            >  2:620005/1‑147 (MQ=255)
                                                                                                                            ccACCCTGGTTTTGGGGAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCg                 <  2:202177/149‑1 (MQ=255)
                                                                                                                            ccAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTTGCATCGTGATCGATTACGACCGCGCGGGAGCTAGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCCGGCATGTGTAATAGTccca                 >  1:135433/1‑148 (MQ=255)
                                                                                                                            ccAACCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGACCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCg                 <  1:237688/149‑1 (MQ=255)
                                                                                                                              aaCCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTGAGGTTTAATTGTCACGTa               >  1:502439/1‑149 (MQ=255)
                                                                                                                              aaCCTGGTTTTGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTa               <  1:231307/149‑1 (MQ=255)
                                                                                                                                      tttGGGTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATAtt       >  2:712312/1‑149 (MQ=255)
                                                                                                                                         gggTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTTTCCCGTAGTCATACTACa    >  1:336170/1‑149 (MQ=255)
                                                                                                                                          ggTAGATTGGCTCAAATTCGCGCTGCGACACGCTAGTGATATTGGCATCGTGATCCATTACGACCGCGCGGGGGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGCCATGTTTAATTGTCCCGTAGTCATATTACAg   <  2:535965/149‑1 (MQ=255)
                                                                                                                                           gTAGATTGGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTTATTGTCCGGTAGTCATATTACAGc  >  2:585916/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TGGTTTCACGCTGGTTCGTAATACCGATAGCTGCAATTTGATCGGAACTGATATCGGCTTTCGCCAGCACTTCTACCAGCGTGGAGCTTTGGGTGGCCCAGATTTCCATTGGGTCGTGTTCTACCCAACCTGGTTTTGGGTAGATTTGCTCAAATTCGCGCTGCGACACGCTAATGATATTGGCATCGTGATCCATTACGACCGCGCGGGAGCTGGTG‑GTGCCCTGGTCGAGCGCAACGATATATTTTTTTTCAGTCATGTTTAATTGTCCCGTAGTCATATTACAGC  >  NC_000913/4116963‑4117250

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