Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 862,786 C→T E242E (GAG→GAA ybiY ← putative pyruvate formate‑lyase activating enzyme YbiY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913862,7860CT96.7% 101.1 / NA 30E242E (GAG→GAAybiYputative pyruvate formate‑lyase activating enzyme YbiY
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/17);  minor base G (0/1);  total (12/18)
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.98e-01

GCACCAGCGCATTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACG  >  NC_000913/862543‑863021
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
gCACCAGCGCATTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTcgg                                                                                                                                                                                                                                      >  1:404633/1‑251 (MQ=255)
          aTTTTTGTGCGCTTTAATGCGGTCGATGACAGTATCCCGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGAGGTTAACCCTTTCCGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGATTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATc                                                                                                                                                                                                                            >  2:192923/1‑251 (MQ=255)
            ttttGTGCGCTTTAATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCg                                                                                                                                                                                                                           <  1:124414/250‑1 (MQ=255)
             tttGTGCGCTTTAATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCgg                                                                                                                                                                                                                          <  2:175774/250‑1 (MQ=255)
                                 cGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGttt                                                                                                                                                                                                     <  1:123698/251‑1 (MQ=255)
                                 cGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGtt                                                                                                                                                                                                      <  1:248231/250‑1 (MQ=255)
                                     gAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTtcgtc                                                                                                                                                                                                 <  1:390171/251‑1 (MQ=255)
                                     gAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTtcgt                                                                                                                                                                                                  <  1:124179/250‑1 (MQ=255)
                                                    cAGTGTGGTGATGGTTGGTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCGGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTCATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCt                                                                                                                                                                                   <  2:399414/250‑1 (MQ=255)
                                                     aGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCGGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGa                                                                                                                                                                                 <  2:114776/251‑1 (MQ=255)
                                                          ggTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAg                                                                                                                                                                            <  2:376080/251‑1 (MQ=255)
                                                            tCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCg                                                                                                                                                                          >  2:321416/1‑251 (MQ=255)
                                                              aTGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGc                                                                                                                                                                        <  1:405267/251‑1 (MQ=255)
                                                              aTGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGc                                                                                                                                                                        >  2:405267/1‑251 (MQ=255)
                                                                       aTCCGCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGGTCGTCGGCGGCAAAATCGGTAATGGCTGTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGAt                                                                                                                                                                <  1:29242/250‑1 (MQ=255)
                                                                            cGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGAt                                                                                                                                                          >  1:11101/1‑251 (MQ=255)
                                                                                    cgcgGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGc                                                                                                                                                  >  2:143752/1‑251 (MQ=255)
                                                                                     gcgGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGcc                                                                                                                                                 >  1:326076/1‑251 (MQ=255)
                                                                                      cgGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGTGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGTTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTGTTGcca                                                                                                                                                <  2:55072/251‑2 (MQ=255)
                                                                                                                           gTCGAGCAGTTCTGGCGCATCAAGCGGTTGTTCCGGTGCGTCATAGGGCAGATTAAGTAAGTGATAGTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGGGCAGTTCGTCGGCGGCAAAAGCGGTAATGGCTTTTACAGAGTTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTgg                                                                                                           <  1:190216/251‑1 (MQ=255)
                                                                                                                                      cTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGAAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGCCGGCGGCAAAATCGGCAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTActg                                                                                                >  1:26809/1‑249 (MQ=255)
                                                                                                                                                aaGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGTCGGCGTTACCGTCGGTCCACt                                                                                      <  2:21942/251‑1 (MQ=255)
                                                                                                                                                  gCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTGCGTCGGCATTAAAGCGCTGAATCAGCGGCACGCGGATAATGATTTTGTTGCCCGCGGCGGCGAGTTTGTTCAGGTGATCCAGGACTCTGGCGGCGTTACCGTCGGTCCACTGt                                                                                    <  1:397152/251‑1 (MQ=255)
                                                                                                                                                                          gCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACgtgt                                                            <  1:26676/251‑1 (MQ=255)
                                                                                                                                                                                                       cccAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCCGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTACACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCTatat                                >  2:197785/1‑250 (MQ=255)
                                                                                                                                                                                                         cAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAgg                             >  1:386351/1‑251 (MQ=255)
                                                                                                                                                                                                         cAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAgg                             >  1:222154/1‑251 (MQ=255)
                                                                                                                                                                                                                 gggaGGGCAGAAAATTAATTTCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGaa                     <  2:288881/248‑1 (MQ=255)
                                                                                                                                                                                                                 ggTAGGGCAGAAAATGAATTTCGCCAACGTGCAGTTCGTCTGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTGCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTTGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGGCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAgaga                      <  2:65738/250‑1 (MQ=255)
                                                                                                                                                                                                                                    ttCGCCAACGTGCAGTTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACg  >  2:3306/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                           
GCACCAGCGCATTTTTGTGCGCTTTAATGCGGTCGCTGAGCGTGTCCAGTTTCAGTGTGGTCATGGTTGTTATCCTCGTAAGGTCGCGGTTAACCCTTTCTGGCAGGCATACTGCTGGGCAAAGTCGAGCAGTTCTGGCGCATCAAGCGGTTTTTCCGGGGCGTCATAGGGCAGATTAAGTAAGTGATATTTGTTGATGCCCAGCGTGTGGTAGGGCAGAAAATGAATTTCGCCAACGTGCAGCTCGTCGGCGGCAAAATCGGTAATGGCTTTTACAGAGGTTTCGTCGGCATTAAAGCCCTGAATCAGCGGCACGCGGATAATGATTTTTTTGCCCGCTGCGGCGAGTTTTTTCAGGTTATCCAGCACTCTGGCGGCGTTACCGTCGGTCCACTGTTTAAACGGCGCGTCGGCAACGTGTTTTAAATCGGCAAGAAACAGATCGATATAGGGCAGAGAAGGGGCGATATATTTCCACG  >  NC_000913/862543‑863021

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