Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,075,552 G→C 22.6% A210G (GCG→GGG)  yggF ← fructose 1,6 bisphosphatase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,075,5520GC22.6% 60.9 / 3.3 38A210G (GCG→GGG) yggFfructose 1,6 bisphosphatase isozyme
Reads supporting (aligned to +/- strand):  ref base G (21/7);  new base C (2/7);  total (23/15)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

GATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTGCGCAGGTTGTCCGTCAGAGGCAAA  >  NC_000913/3075350‑3075776
                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
gATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATTGTGTacatcacatc                                                                                                                                                                                                         >  1:520158/1‑228 (MQ=255)
  tatCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGTCGCACCGCCGATGGCGTACATCACATCGt                                                                                                                                                                                                       >  1:205419/1‑228 (MQ=255)
                    gTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGCCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCATTGGGGTACATCAAATCGTAAGGATTATACTgaccat                                                                                                                                                                                     >  1:662220/1‑223 (MQ=255)
                          cgagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATTCTTTAACGGCACAGGCAGACATCCCGCCTTCTGGCGCACCGCCGATGTTGTacatca                                                                                                                                                                                                             >  1:501251/1‑198 (MQ=255)
                          cgagcgagTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCCTCCGCACCGCCCATGTGGTAACTCCCATCGAAAGGATTATCCTGccccccgcccag                                                                                                                                                                              >  2:578915/1‑219 (MQ=255)
                                       acGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCCAGCCTTCTGGCGCACCGCCGTTGGAGTACATCACATCGTAAGGATTATCCTGCCAGCAGGACAGCACGCTGgcggcg                                                                                                                                                                 >  2:28769/1‑229 (MQ=255)
                                         ggTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGCCGAc                                                                                                                                                               >  1:132041/1‑229 (MQ=255)
                                                   aaCGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAACAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGCAATCACGCCTTCTGGCGCACCGCCGATGGGGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCCCCAccc                                                                                                                                                     >  1:330152/1‑229 (MQ=255)
                                                      cctctttctctgccAACGCTTGCGCCCCTCCTCCGCAATCCGCGTTTTTCCCGTGTAATCCCCTTTAGCCTGGAAAAATCCAAACCGCCCGCCCTGCATGCCCCCGCTAAATGGTTAAACGGAAAGGGCAAAAACCACGCTTTTGGGCCCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGc                                                                                                                                                  <  2:520158/216‑1 (MQ=255)
                                                           ttGCCTTGCAACGCTTCCGCTCCTGCGCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTCCATATTCCCGCCTAATGCTTTAAAGGCCCAGGACGAAAACCCGCCTTCTGGCCCACCGCCCGTTCTGTATATCACATCAGAAGGCATCACCAGCCCACAGGTCAGCACGCTggcgggcgcaacggaaaccgacaggag                                                                                                                                             >  2:377517/1‑207 (MQ=255)
                                                                     ccgttggcgcCCCTCCTCGGAAATGCGCCTATTTCCCGTATACTCCCTTTTAACGGGAAAAAGGCAAATCAGCCGGCCCTGCATGCCCCCGCCTAATGTGTAACGGTCACAGGCAAAAACCCCCGCTTCTGGCCCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTc                                                                                                                                   <  2:19364/223‑1 (MQ=255)
                                                                        ctgccgcTCCGGCTCTGCAATCTGCCAATTTTCCTGTTAACCCCCTTTAGCCGGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACg                                                                                                                                <  1:405742/225‑1 (MQ=255)
                                                                        cTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGCTGGTGTACATCACATCGTAAGGACTATCCTGCCAGCAGGTCAGCACGCTGCCGGCGACATTGCCGTCCGGCAGGGCAAAAACCTTCACg                                                                                                                                >  1:398166/1‑229 (MQ=255)
                                                                        cTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTc                                                                                                                                   >  1:19364/1‑226 (MQ=255)
                                                                        cTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCCCACCGCCGATAGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCCGGGCAAAAACCCTCACg                                                                                                                                >  1:187768/1‑229 (MQ=255)
                                                                          tGCCCTCCTGCCCGGCAATCTGCCGATTTCCCTTTTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTGGGCCTGCATCTCCCCGCCTAATCCTTTAACGGCACAGGAAGAAATCACGCCTTTGGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGcc                                                                                                                              <  1:468079/229‑1 (MQ=255)
                                                                             ccccctgctaccaaactcgCCGTTTTTCCGTGTAACCCCTTTTAGGCTGGAAAAAGTCAAACAGCCGGCCCTGCATGTCCCCCGCTAATGCTTTAACGGCACGGGCAAAAATCACGCCTTCTGGCCCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAg                                                                                                                           <  2:223163/211‑1 (MQ=255)
                                                                                                                       tAGCCTGGAAAAAGTAAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTACCGGCACAGGAAGAAATCACGCCGTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATc                                                                                 <  1:647785/229‑1 (MQ=255)
                                                                                                                        tgtgctggaaaatctcaccaggCTCGGCTTGTATGTCCCCGCGTAGTGTTTTAACGGCACGGCCAGAAACCCCGCTTTCTGCCGCCCCGCCGATGGTGTAAATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATc                                                                                 <  2:657706/208‑1 (MQ=255)
                                                                                                                              gCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGTCCAGCACGCTGGCGGCGACACCCCCATCCGTCTGGGCAAAAACCTTCAGCCCCAGCTCGTTTGATTCCCCAATCCCGGCACACAgcaccctctgtacaccggta                                                                          >  2:653039/1‑210 (MQ=255)
                                                                                                                                                cGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACCCAGGCGCGGTTTATCCAGCGTAACCATGCGCACATTAt                                                          >  1:469945/1‑228 (MQ=255)
                                                                                                                                                  gCCCGCATACCCCCCCTTAAGGCTTAACCTGCACAGGCAGAAATCACGCCTTGTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCa                                                       <  1:578916/229‑1 (MQ=255)
                                                                                                                                                        atatCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCCCCCCCAGCCCGGGTGCTTCCTCCATCGCGGCCCCCAGGCCCGGTTTTTCCCGCGTAACCATgccgcgctttacacgcggt                                                  >  2:218649/1‑211 (MQ=255)
                                                                                                                                                                      tGCTTTAATGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCACCTGCGTTGCTTCCTCAATCGCGGCACCCAGGCGCGGTTTACCCCCCGTAACCCTCCGCAGCTTATCCCGCGGTTTGCCTAACGCTTTc                                   >  2:296851/1‑229 (MQ=255)
                                                                                                                                                                      tGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCACCGGTTTGCCTAACGCttt                                    >  1:373018/1‑228 (MQ=255)
                                                                                                                                                                      tGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGGGTTGCTTCCTCAATCGCGGCACCCCGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCAAACGCctt                                    >  1:90509/1‑228 (MQ=255)
                                                                                                                                                                           ttcgcgctacaGGCAAAATTCCCGCTTCCGGCCCACCCGCCGAGGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCa                               <  2:373018/221‑1 (MQ=255)
                                                                                                                                                                           tAACGGCACAGGCAGAAATCCCGCCTTCTGGCCCCCCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGgc                                                                                           <  1:538454/168‑1 (MQ=255)
                                                                                                                                                                              gcgacacaGGCAAAAATCCCGCAGTCTGGCTGACCGCCGATGGTGTATATCACAACGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGt                            <  1:148040/225‑1 (MQ=255)
                                                                                                                                                                                  ccgcagCACAAATCCCGCCTTCTGGCCACCCGCCGGTGGTGTACATCACATCGTAAGGATTACCCTGCCAGGAGGTCAGCCCGCTGGCGGCGACATCGCCATCCGGGAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTcgcg                        <  2:621652/224‑1 (MQ=255)
                                                                                                                                                                                  cacaGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACCCAGCGCCGGTTTATCCAGCGTAACCCTGCGCCGCTTATCCCGCGGTTTCCCTAACGCTTTCGCCCCCTTgcg                        >  2:687533/1‑229 (MQ=255)
                                                                                                                                                                                  caacggCAGAAACCCCGCCTGCGGGCCACCCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTgcg                        <  1:611868/225‑1 (MQ=255)
                                                                                                                                                                                           aaaTCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTGCGCAGGTTGTc               >  1:555679/1‑229 (MQ=255)
                                                                                                                                                                                                ccGCCGTCTGGCGCCCGCGCGATTGGTTACGTCCCATCAGAATGTCTATCCGGCCAGCAGGTCAGCACGCTGGCGGCGAAATCGCCCGCCGGCAGTGCAAAAACCTTGCCGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTGCGCAGGTTGTCCGTCa          <  2:469945/228‑1 (MQ=255)
                                                                                                                                                                                                 ccttgtcgggCGCCCCGCCCAGGGGTTACAACACATCGGAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTGCGCAGGTTGTCCGTCa          <  2:555678/220‑1 (MQ=255)
                                                                                                                                                                                                 cGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTACCCATGCGAAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTGCGCAGGTTGTCCGTCag         >  1:442728/1‑229 (MQ=255)
                                                                                                                                                                                                 cGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACCCTTTCCCCACCTTACCCAGGGTTTCCGCCag         >  2:152247/1‑229 (MQ=255)
                                                                                                                                                                                                        tGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGCTTTGCCTAACGCTTTCGCAACGTTccgcctgtctcttctacccacctcc  >  1:612203/1‑208 (MQ=255)
                                                                                                                                                                                                            |                                                                                                                                                                                                                                  
GATATCGTTCCCCCTCACCAGTTCATCGAGCGAGTACACACGGTTAACGTCAACGCCCATTGCCTTGCAACGCTTGCGCTCCTGCTCTGCAATCTGCCGATTTTCCGTGTAATCCCCTTTAGCCTGGCAAAAGTCAATCAGCTCGGCCTGCATATCCCCGCCTAATGCTTTAACGGCACAGGCAGAAATCACGCCTTCTGGCGCACCGCCGATGGTGTACATCACATCGTAAGGATTATCCTGCCAGCAGGTCAGCACGCTGGCGGCGACATCGCCATCCGGCAGGGCAAAAACCTTCACGCCCAGCTGTGTTGCTTCCTCAATCGCGGCACTCAGGCGCGGTTTATCCAGCGTAACCATGCGCAGCTTATCCAGCGGTTTGCCTAACGCTTTCGCCACGTTGCGCAGGTTGTCCGTCAGAGGCAAA  >  NC_000913/3075350‑3075776

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