Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,073,952 G→A A72T (GCA→ACA)  yihV → 6‑deoxy‑6‑sulfofructose kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,073,9520GA97.1% 105.0 / NA 35A72T (GCA→ACA) yihV6‑deoxy‑6‑sulfofructose kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/13);  minor base C (1/0);  total (22/13)
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

CAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCC  >  NC_000913/4073718‑4074183
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
cAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTg                                                                                                                                                                                                                         <  1:75964/251‑1 (MQ=255)
 aaaaTGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTgg                                                                                                                                                                                                                        >  1:157279/1‑251 (MQ=255)
              aaTGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGGGATGACGACACGGGCAACAGACTGCTGACAGAGCTGGAATCCTGGGGCGTTAACAccc                                                                                                                                                                                                           >  1:3385/1‑251 (MQ=255)
                 gaacATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCACCAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGt                                                                                                                                                                                                         <  2:201515/247‑1 (MQ=255)
                    cATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCTCATCTATTACGTGGAAGGGTTACCGACTGAGAGCGGTAAATTCGTGGCGAGAAATTATACGGAAGTTTGTGGCGTGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTTGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCGGGAATCCTGGGGCGTTAACACCCGTTaca                                                                                                                                                                                                     <  2:130183/251‑1 (MQ=255)
                         ttCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGATACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTTTTGGGCGCGGAGGGGATGACGACACGGGCAACAGCCTGCTGACAGAGCTAGAATCCTGGGGCGTTAACACCCGTTACACCaa                                                                                                                                                                                                 >  1:200127/1‑250 (MQ=255)
                                         aGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACAACCAACGGTATAACCAGGCg                                                                                                                                                                                 >  1:152809/1‑250 (MQ=255)
                                                     gATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAAt                                                                                                                                                                    <  2:154485/251‑1 (MQ=255)
                                                      aTGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGCAGTTGGGGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATCACCAGGCGAAATCTTCGCAATc                                                                                                                                                                   >  2:11563/1‑251 (MQ=255)
                                                       tGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATcc                                                                                                                                                                  >  1:282353/1‑251 (MQ=255)
                                                        ggATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGGCGATTTTATTGGTCGCGTAGGAGATGACGACCCGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAAAACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCCATccc                                                                                                                                                                 >  2:306930/1‑250 (MQ=255)
                                                            cGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGACGCGTAGGTGATGACGACTCGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTGCACCAAACGGTATACCCAGGCGAAATCTTCGCAATCCGCcat                                                                                                                                                             <  1:339601/251‑1 (MQ=255)
                                                                            aaGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGGAGGTGATGTCGACACGGGCAACAGCCTGCTGACAGAGGTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATTGTGGATAccca                                                                                                                                              >  2:341902/1‑248 (MQ=255)
                                                                                        cGGAGAGCGGCAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGGGATGACGACACGGGCAACAGCCTGCTGCCAGAACTGGAATCCTGGGGGGGTAACACCCGTTAAACCACACGGGATAACCAGGCGAAATCTTCGCAATCCGCCATCCTGGTGGATACCAAAGGCGAGCGGAt                                                                                                                                  >  2:247898/1‑250 (MQ=255)
                                                                                            gagCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCa                                                                                                                             <  1:212145/251‑1 (MQ=255)
                                                                                              gCGGTAAATACGTGGCGAGAAATTATACGGAAGTGGGTGGCGGGCCAGCGGCGACGGCTGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCGGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAACGGCGAGCGGATAATCAtt                                                                                                                           <  2:15286/251‑1 (MQ=255)
                                                                                               cGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGGGCGGGTAATCAtt                                                                                                                           >  2:269409/1‑250 (MQ=255)
                                                                                                     aTACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAac                                                                                                                                                                                                               <  2:42687/160‑1 (MQ=255)
                                                                                                     aTACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAac                                                                                                                                                                                                               >  1:42687/1‑160 (MQ=255)
                                                                                                     aTACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTAcc                                                                                                                    <  1:269409/251‑1 (MQ=255)
                                                                                                       aCGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTAcccc                                                                                                                  <  1:218404/251‑1 (MQ=255)
                                                                                                             cGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGGATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTAACCCAGCCCg                                                                                                            >  2:289518/1‑251 (MQ=255)
                                                                                                                             atttGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACg                                                                                            <  2:201130/249‑1 (MQ=255)
                                                                                                                                       gggCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGtt                                                                                  >  1:152888/1‑251 (MQ=255)
                                                                                                                                          ccAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTgg                                                                                >  1:103330/1‑250 (MQ=255)
                                                                                                                                                    cGGCAGCGGTTGCGGCGGCAAGGCTGGCGGCGCCGGTCGTTTTTATTGGCCGCGTAGGTGATGACAACACGGGCAACAGCCTGCTGACAGAGCTGCGATCCTGGGGCGTTAACACCCCTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACACCAGCCCGGCTCTGCTGCCTGACGCAGCGTGGTTGGAGGCAATTGAt                                                                     >  2:128710/1‑251 (MQ=255)
                                                                                                                                                                        aGGCTGGGGGCGCCGGTCGATTTTATTGGTCGCGTAGGTGCTGACGACCCGGGCACCAGCCTGCTGACAGAGCTGGCATCCTGGGGCGTTACCACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGAtgttgt                                                 >  2:183287/1‑251 (MQ=255)
                                                                                                                                                                             gggggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTgg                                            <  2:44385/251‑1 (MQ=255)
                                                                                                                                                                               gggCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGc                                           >  2:96382/1‑250 (MQ=255)
                                                                                                                                                                                             tttATTGGTCGCGTAGGTGATGACGACACGGGAAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGAGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGACCGGATCTGCTGCATGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTAAGTGGGATGTTGTGCTGGCAGATGTACGCTGGc                             >  1:15768/1‑250 (MQ=255)
                                                                                                                                                                                                   ggccgcgTAGGTGATGCCGACACGGCCAAAAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTCACACCCGTTACACCAAACGGTATAACCAAGCGTAACCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACGACCCTCGCCCGGCTCTGCGGCCTGACGCCGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGGGGTTGTGGTGGCAGATGTACGCTGGCACGACgg                      >  2:207953/4‑251 (MQ=255)
                                                                                                                                                                                                   ggTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGAACGACgg                      >  1:61061/1‑251 (MQ=255)
                                                                                                                                                                                                           aGGTGATGAAGACACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTaaaa              <  2:250778/251‑1 (MQ=255)
                                                                                                                                                                                                                     gacACGGGCAACAGCCTGCTGACAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCAc    >  1:147218/1‑251 (MQ=255)
                                                                                                                                                                                                                       cacGGGCAACAGCCTGCTGACAGAGCTGGAATCCCGGGGCGTTAACACCCACTACCCCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGGACGCTGGCACGACGGCGCTAAAAAAGCCTTCAccc  >  2:142498/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                       
CAAAATGTTTAAGGAATGACCATGATTCGTGTTGCTTGTGTAGGTATAACCGTGATGGATCGCATCTATTACGTGGAAGGGTTACCGACGGAGAGCGGTAAATACGTGGCGAGAAATTATACGGAAGTTGGTGGCGGGCCAGCGGCGACGGCAGCGGTTGCGGCGGCAAGGCTGGGGGCGCAGGTCGATTTTATTGGTCGCGTAGGTGATGACGACACGGGCAACAGCCTGCTGGCAGAGCTGGAATCCTGGGGCGTTAACACCCGTTACACCAAACGGTATAACCAGGCGAAATCTTCGCAATCCGCCATCATGGTGGATACCAAAGGCGAGCGGATAATCATTAACTACCCCAGCCCGGATCTGCTGCCTGACGCAGAGTGGTTGGAGGAAATTGATTTCTCTCAGTGGGATGTTGTGCTGGCAGATGTACGCTGGCACGACGGCGCTAAAAAAGCCTTCACCC  >  NC_000913/4073718‑4074183

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