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_000913	21,619	0	A	G	96.0%	69.7 / NA	25	R71R (AGA→CGG)	ribF	bifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (16/8);  minor base T (0/1);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

CGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGAT  >  NC_000913/21382‑21856                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              cGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGc                                                                                                                                                                                                                                  <  1:54466/251‑1 (MQ=255)            gTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACt                                                                                                                                                                                                                        <  2:53253/250‑1 (MQ=255)                        aCATGAAGCTGATACGCGGCATACACAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTAGTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGTCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCa                                                                                                                                                                                                           <  2:189828/251‑1 (MQ=255)                                     aCGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTg                                                                                                                                                                                              >  1:75815/1‑251 (MQ=255)                                          gCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTAc                                                                                                                                                                                         >  2:65633/1‑251 (MQ=255)                                                         gCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGTCGGCGTGCATCGCGGTCATAGCGCGCTGTTACAGGGCTTGCAGGAAGCAGGGCGAAAGAGCACCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACAGCTCAAGCCCCGGCCCGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGGGGCCTTGCTTACGTGCTGTgcgtgcgt                                                                                                                                                                          >  1:262449/1‑251 (MQ=255)                                                                         agGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTcg                                                                                                                                                          >  2:60191/1‑251 (MQ=255)                                                                               tgtgCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACATTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGTCAGGCGTTTCGCGGCGt                                                                                                                                                    >  2:106081/1‑251 (MQ=255)                                                                                               aaTTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAAtt                                                                                                                                     <  1:156521/250‑1 (MQ=255)                                                                                                    cGACGGCGTGCATCGCGGTAATCGAGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCACCCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAACCTGCGTTACCTTGCAGAGTGTGTCGTTGATTACGTGCTGGGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAAcagctccaactttcttc                                                                                                                                >  2:164217/1‑234 (MQ=255)                                                                                                                          cgcgcgCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTTCCGGTGATGGTGATTCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGAGTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCtttt                                                                                                         <  2:264254/251‑1 (MQ=255)                                                                                                                                              gCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCCCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACAGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCg                                                                                     >  1:255403/1‑251 (MQ=255)                                                                                                                                                              aagcgcaaCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCCAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCg                                                                     >  1:159927/1‑251 (MQ=255)                                                                                                                                                                     aCTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGgc                                                              >  1:76100/1‑251 (MQ=255)                                                                                                                                                                           cGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGt                                                        <  2:76100/251‑1 (MQ=255)                                                                                                                                                                            ggtgatggtgatGCTTTTTAAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGTCAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGGGATGATTTCCGCTTTGGCGCTGGt                                                        >  2:220399/1‑250 (MQ=255)                                                                                                                                                                                  ggtgatgCTTTTTGAACCTAAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGGTGCGGGAACAACTGCGTTAACTTGCAGAGTGTAGCATTCATTACGTGATGTGCGTGCGTTTCGCCAGGCGCTTCGCGGCGTTCACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGAGTAACATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTggacggga                                                   >  1:238369/1‑243 (MQ=255)                                                                                                                                                                                           ttttGAACCCCAACCACTGGAACTGTTTGCTACCGATAAAGCGCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTGGCAGAGTGTGGCGTTGATTACGGGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTTCCGTAGGTGATGATTTCCGCTTTGGCGCTAGTCGTGAAGGCGAttt                                          <  2:255403/249‑1 (MQ=255)                                                                                                                                                                                              tGAACCTCAACCACTGGAACTGTTTGCTACCGTTCAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGt                                     >  2:96279/1‑251 (MQ=255)                                                                                                                                                                                                              gAACTGTTTGCTACCGATAAAGCCCCGGCACGGCTGAACCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGAGTTTCGACCGGCGTTTCGCGGCGTTAACCGCGAAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTATTGTTATTAAAGAAAGCTgg                     >  2:194363/1‑251 (MQ=255)                                                                                                                                                                                                                  tGTTTGCTACCGATAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAACCTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGTCATg                 >  2:148764/1‑251 (MQ=255)                                                                                                                                                                                                                        cTACCGATAAAGCCACGGCACGGCTGACCCGGCTGCGGGAACAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGGGATGCTTTCCGCTTTGGCGCTGGTCGTGCAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATAc           >  2:251398/1‑251 (MQ=255)                                                                                                                                                                                                                            cGATAACGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCt       <  2:140326/251‑1 (MQ=255)                                                                                                                                                                                                                              aTAAAGCCCCGGCACGGCTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTc     <  1:65633/251‑1 (MQ=255)                                                                                                                                                                                                                                  aGCCCCGGCACGGCTGACCCGGCTGCGGGAAACACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCTTTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGat  >  1:224874/1‑250 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              CGATTTTCACTGTTTTGAGCCAGACATGAAGCTGATACGCGGCATACATAATCTCAGCCAGGCCCCGCAAGAAGGGTGTGTGCTGACTATTGGTAATTTCGACGGCGTGCATCGCGGTCATCGCGCGCTGTTACAGGGCTTGCAGGAAGAAGGGCGCAAGCGCAACTTACCGGTGATGGTGATGCTTTTTGAACCTCAACCACTGGAACTGTTTGCTACCGATAAAGCCCCGGCAAGACTGACCCGGCTGCGGGAAAAACTGCGTTACCTTGCAGAGTGTGGCGTTGATTACGTGCTGTGCGTGCGTTTCGACAGGCGTTTCGCGGCGTTAACCGCGCAAAATTTCATCAGCGATCTTCTGGTGAAGCATTTGCGCGTAAAATTTCTTGCCGTAGGTGATGATTTCCGCTTTGGCGCTGGTCGTGAAGGCGATTTCTTGTTATTACAGAAAGCTGGCATGGAATACGGCTTCGAT  >  NC_000913/21382‑21856

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