Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 21,619 A→G R71R (AGA→CGG ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091321,6190AG96.7% 84.4 / ‑4.3 31R71R (AGA→CGGribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (17/12);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01

TCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACG  >  NC_000913/21377‑21848
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
tCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCg                                                                                                                                                                                                                                >  1:14219/1‑251 (MQ=255)
 cGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGAccc                                                                                                                                                                                                                                 <  1:120063/249‑1 (MQ=255)
      gcttttCACTGTTTTGAGCCAGACATGAAGCTGATACGAGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGACATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGAACCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAATGCCCCGGCACGTCTGCCCCGCCTGCg                                                                                                                                                                                                                          >  2:185032/3‑251 (MQ=255)
             aCTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGaaaa                                                                                                                                                                                                                    >  1:6304/1‑250 (MQ=255)
                              aTGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGc                                                                                                                                                                                                     >  1:224121/1‑248 (MQ=255)
                                                               cAGGCCCCGCAAGCAGGGTGTGTGCTGACTATTGGTAATTTCGAAGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAACGCCCCGGCCAGGCTGACCCGGCCGCGGGGAAAGCTGCGTTACCTTGCAGAGTGTGGCGTTGAATACATGCTGTGCTGGCGttt                                                                                                                                                                 >  2:42420/1‑251 (MQ=255)
                                                                       gCAAGAAGGGAGTGTGCTGACTATTGGTAATTTCGACTGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTAGTAGGAAGACGGTCGAATGCGCACCATACCGGAGATGGTGATGCCTTTTGAACGTAAAGCCCGGGAACTATTCGCTACCGAGAAAGCACCAGCACTGCTGACCAGGCAGCGGGAAAACCTACGCTACCATGCTGAGTAGGGCGTTGATAACGTgtgttgctgacgtatcgacagt                                                                                                                                                          >  2:206823/1‑229 (MQ=255)
                                                                                   tgtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGt                                                                                                                                            >  1:323738/1‑251 (MQ=255)
                                                                                         gACTATTGGTGATTTCGACGGCGTGCAGCCGGGTCATAGCGCGCTGTTACAGGTGTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAAACTCAACCACTGGAACTGTTTGCTACCGATAAAGTCCCGGCACGGCTGACCCGGCGGCGGGAAAAACTGCGTTACCTTGCAGAGTGGGGCGTTTATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAAcc                                                                                                                                       <  2:82146/250‑1 (MQ=255)
                                                                                         gACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGCTAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAACAACTGCGTTCCCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACcg                                                                                                                                      >  1:268608/1‑251 (MQ=255)
                                                                                               tGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCaa                                                                                                                                 >  1:383935/1‑250 (MQ=255)
                                                                                                   aaatgcGACGGAGTGCATCGCGGTAGTCTCGCGCTGTTTCAGTGCTTGCAGGAAGAAGGTCGCAAGCGCAACTTATCGGTGATGGTGATGCTTGTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCTCCGGCACGGCTGTCCCGGCTGCGGGAAAATCTGCGGTACCTTGCAGATTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGGCAGTTGTTTCGCGGCGTTAACCGCGCAAAAtt                                                                                                                             <  2:335170/245‑1 (MQ=255)
                                                                                                        cGACGGCGGTCATCGCGGTTATCGCGCGCGGTGACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGGGATGGTGATGCTTTTTGTACCTCAACCACTGGAACTGGTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTtcatc                                                                                                                        <  2:224121/250‑1 (MQ=255)
                                                                                                             gCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGAt                                                                                                                  >  1:107041/1‑251 (MQ=255)
                                                                                                                     cgcgGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTggg                                                                                                          >  1:172048/1‑250 (MQ=255)
                                                                                                                       cgGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTg                                                                                                         >  1:136435/1‑250 (MQ=255)
                                                                                                                         gTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCCAAATTTCATCAGCGATCTTCTGGTGa                                                                                                        >  1:378035/1‑249 (MQ=255)
                                                                                                                                 gcgcTGTTAAAGGGCTTGCAGGCAGAAGGACGCAAGCGCAACTTAACGGTGATGGTGATGCTTTTTGCACCTCAACAACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAAATGCGTGAACTTGCAGAGTGAGGCGTTGATTACGGGCTGTGCGGGAGTTTCGAGAGGGGTTTCGGGCGGTGAAACGCGAAAACTTTCCTCAGCGATATTCTTGTGacagttttgcc                                                                                              >  2:271038/1‑241 (MQ=255)
                                                                                                                                      gTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTaa                                                                                         <  1:80751/251‑1 (MQ=255)
                                                                                                                                            gggCTTGCAGGAAGAAGGGCGCAAGCGCACCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTc                                                                                   <  2:323738/251‑1 (MQ=255)
                                                                                                                                                                 caagcgctaATTACCTGTGATTTTGATGATTTTTGCACCTCAACCATTGGAACTGTTAGGTGCCGAAAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCGTGCAGAGTGTGTCGTTGATTACGTGCTGTGCTTGCGTTTCGTCAGTCGGTTCGCGGCGTTAACCGCTCAAAATTTCATCAGCGATCTTCTGGGGAAGCTTTAGCGCGTATAATTTCTGGCCTTAGGTGATGATTTca                                                              <  2:27179/251‑2 (MQ=255)
                                                                                                                                                                 caagcgcaaCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTcc                                                              >  1:326115/1‑251 (MQ=255)
                                                                                                                                                                    gcgtaacgtACCTGTGATGGTGATGTTTGTTGACCCTCGACCACTGGAACTGTTTGCTTCCGATAAAGGCCCGGTACGGCTGACCCGGTTGCTGTAAAAACTGTGGTACCTTGCATAGTGCGTCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCTCGGCTTTAACCGCGTAAAAGTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAAGTTCTTGCCGTAGGTGATGGTGTCCGCt                                                           <  2:233097/243‑1 (MQ=255)
                                                                                                                                                                           ttACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGAGAAGCTTTTCCGCGTAAACTTTCTTGCCGTAGGGGATGATTTCCGCTTTGtcgc                                                    >  2:88362/1‑251 (MQ=255)
                                                                                                                                                                                ggtgatggtgatGCTTTTGGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTTCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGTTGAAGCATTTGCGCGTAAAATTTCTTGCTGTAGGTGATGATTTCCGCTTTGGCGCTGGTc                                               <  2:146145/251‑1 (MQ=255)
                                                                                                                                                                                        tgatgCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGGACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCgcgg                                                                                                                                                <  2:105174/147‑1 (MQ=255)
                                                                                                                                                                                        tgatgCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGGACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCgcgg                                                                                                                                                >  1:105174/1‑147 (MQ=255)
                                                                                                                                                                                                gttGAACCTCAACCCCCGGAACTGTTGGGTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACTTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCt                                <  2:14219/249‑1 (MQ=255)
                                                                                                                                                                                                                 ggAACTGTTTGCTACCGATAAAGCCCCGGCGCGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTg              <  1:118955/251‑1 (MQ=255)
                                                                                                                                                                                                                              aCCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACg  <  1:88362/250‑1 (MQ=255)
                                                                                                                                                                                                                                       gCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGAtt                                                                 >  1:122260/1‑178 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                      
TCGGACGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACG  >  NC_000913/21377‑21848

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