Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 22,373 T→C intergenic (+25/‑18) ribF → / → ileS bifunctional riboflavin kinase/FMN adenylyltransferase/isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091322,3730TC88.6% 81.8 / 1.2 35intergenic (+25/‑18)ribF/ileSbifunctional riboflavin kinase/FMN adenylyltransferase/isoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (14/17);  total (16/19)
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.70e-01

GGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAG  >  NC_000913/22126‑22606
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
ggCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACg                                                                                                                                                                                                                                        <  1:417646/251‑1 (MQ=255)
               ccGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCtga                                                                                                                                                                                                                         <  1:50679/250‑1 (MQ=255)
                         gTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACtt                                                                                                                                                                                                               >  2:168394/1‑250 (MQ=255)
                            gCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATaa                                                                                                                                                                                                            <  1:353924/251‑1 (MQ=255)
                                gTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATc                                                                                                                                                                                                         >  1:31876/1‑250 (MQ=255)
                                  aTTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAAc                                                                                                                                                                                                      <  1:100152/251‑1 (MQ=255)
                                      gCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATTGACCTTTACGGTCGCCATATTCAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGAGTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGGGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTg                                                                                                                                                                                                  <  2:339586/251‑1 (MQ=255)
                                               tGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCaa                                                                                                                                                                                                       >  1:169440/1‑237 (MQ=255)
                                               tGGAAGTGCATTTGTTAGATGTTGCAAGGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTATCCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCaa                                                                                                                                                                                                       <  2:169440/237‑1 (MQ=255)
                                                            gTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTc                                                                                                                                                                             <  2:62753/250‑1 (MQ=255)
                                                               aGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCg                                                                                                                                                                          >  2:242434/1‑250 (MQ=255)
                                                                   gTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGc                                                                                                                                                                      <  1:242434/250‑1 (MQ=255)
                                                                      gCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATTAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAAGTTGCCGGAAACAGGGTTCCCGATGCGTg                                                                                                                                                                   <  2:368928/250‑1 (MQ=255)
                                                                      gCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTg                                                                                                                                                                   <  1:259910/250‑1 (MQ=255)
                                                                                         gCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCTTGGCGATCTCGCCAAGCGCGa                                                                                                                                                >  2:64962/1‑250 (MQ=255)
                                                                                         ccccacccgacgcccttccgaaccaaaaaaaaCGCAATGAGCAGCGAATTGCGGCGCTGGAAGAACTGAAAGCGCAGAATGCGCGTGATGAAATAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGGTATGTAAACAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGaa                                                                                                                                               <  2:55895/227‑1 (MQ=255)
                                                                                            atatACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAATCGCGAAccg                                                                                                                                            >  1:202698/1‑250 (MQ=255)
                                                                                            atatACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAAccg                                                                                                                                            >  1:384277/1‑250 (MQ=255)
                                                                                                           tgcGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGt                                                                                                                             >  1:428699/1‑251 (MQ=255)
                                                                                                                 aaaaaaTACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGaa                                                                                                                                               >  1:55895/1‑227 (MQ=255)
                                                                                                                            aaTGAGCAGCTATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCt                                                                                                            <  2:428699/251‑1 (MQ=255)
                                                                                                                                  cagcGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACgg                                                                                                      >  1:72750/1‑251 (MQ=255)
                                                                                                                                             gTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTg                                                                                           >  1:465726/1‑251 (MQ=255)
                                                                                                                                                                    aGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAAcc                                                                    >  1:422941/1‑251 (MQ=255)
                                                                                                                                                                     gatttcgcttgagtaaTTAACCGCCCGCGAATTTTTTGAGCTAACAAAACCGGCTTAAGCCTGTGTTGTAATCAAACCGAAACACTGAACCGAGAATCTGATGAGTGACTATAAATCAACCCGGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGGTGGACTGATGATTATCTGTTCGGCATCATCCGTGCGTCTAAAAGAGGCAAAATAACCt                                                                   <  2:461781/237‑1 (MQ=255)
                                                                                                                                                                     gATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCt                                                                   <  1:31845/251‑1 (MQ=255)
                                                                                                                                                                      aTTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCtt                                                                  >  2:74131/1‑251 (MQ=255)
                                                                                                                                                                                         ccgcccgcGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCcct                                               <  2:132147/251‑1 (MQ=255)
                                                                                                                                                                                                  aaTTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTCAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGACCTCCTTATGCg                                      >  2:136881/1‑251 (MQ=255)
                                                                                                                                                                                                               aaCAAAACCGGCTTAAGCCTGTGTTGTAATCAAACCGAAACTCGTAACCGAGAAACTGATGATTGACTAGAATTCAACCCTTATTTTGCCGGAAACAGGGTTCCCGATGCGTGGCGCTCTCGCCAAGCGCGACCCGGGAATGCTGGCGCGTTTGACTGATTCGTATCTGTACGGCAGCATCCTTGCGGCTAAAAAAGGCAAAAAAACCCTCATTCGGCATGATTGCCCTCCTTCTGCGAATGTCAGCAtgc                         <  2:302410/251‑3 (MQ=255)
                                                                                                                                                                                                                aCAAAACCAGCTTAAGCCTGTTATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAAGCAACCCTGAATGTGCCGGAAACAGGGTGCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTGATGCGAATGGCAGCATTCa                        <  2:258524/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tATGTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGACACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGCTCTGTACGGCATCCTCCGAGCGGCTAAAAACGGCAAAAAAACATTCATTCTGCATGATGGCCCTCCTTATCAGAATGGCAACAGTCatttgggtccataggttcaca    >  2:305248/1‑231 (MQ=255)
                                                                                                                                                                                                                                      aagTAATCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATGCATATTGGTCACTCGGTTAACAAg  <  2:482113/249‑1 (MQ=255)
                                                                                                                                                                                                                                           aTCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCAt                                                               >  2:63529/1‑185 (MQ=255)
                                                                                                                                                                                                                                           aTCAAACCGAAACACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCGGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCAt                                                               <  1:63529/185‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         
GGCAAACATCGGAACACGCCCAACGGTTGCCGGTATTCGCCAGCAGCTGGAAGTGCATTTGTTAGATGTTGCAATGGACCTTTACGGTCGCCATATACAAGTAGTGCTGCGTAAAAAAATACGCAATGAGCAGCGATTTGCGTCGCTGGACGAACTGAAAGCGCAGATTGCGCGTGATGAATTAACCGCCCGCGAATTTTTTGGGCTAACAAAACCGGCTTAAGCCTGTTATGTAATCAAACCGAAATACGGAACCGAGAATCTGATGAGTGACTATAAATCAACCCTGAATTTGCCGGAAACAGGGTTCCCGATGCGTGGCGATCTCGCCAAGCGCGAACCCGGAATGCTGGCGCGTTGGACTGATGATGATCTGTACGGCATCATCCGTGCGGCTAAAAAAGGCAAAAAAACCTTCATTCTGCATGATGGCCCTCCTTATGCGAATGGCAGCATTCATATTGGTCACTCGGTTAACAAG  >  NC_000913/22126‑22606

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