Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,315,384	G→A	Y129Y (TAC→TAT)	yciF ←	DUF892 domain‑containing protein YciF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,315,384	0	G	A	95.7%	56.1 / ‑3.4	24	Y129Y (TAC→TAT)	yciF	DUF892 domain‑containing protein YciF
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (10/12);  total (10/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGCTGGCCATGGATTCAAGCATAGATTCGGCTTGCTTTTCCATTGCGTGGGCGTCACGTAACCAGTCATGATAATGTTCAATACGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTC  >  NC_000913/1315134‑1315623                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                ggCTGGCCATGGATTCAAGCATAGATTCGGCTTGCTTTTCTATTGCGTGGGCGTCACGTAACCAGTCATGATAATGTTCAATACGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACgg                                                                                                                                                                                                                                                 >  1:128124/1‑251 (MQ=255)         aTGGATGCAAGCATAGATTCGGCTTGCTTTTCCATTGCGTGGGCGTCACGTAACCAGGCATGATAATGTTCAATACGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGTCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTATGATAACCTaa                                                                                                                                                                                                                                         <  1:390016/251‑1 (MQ=255)                                  gcttTTCCATTGCGTGGGCGTCACGTAACCAGTCATGATAATGTTCAATACGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTc                                                                                                                                                                                                                <  1:91426/251‑1 (MQ=255)                                      ttCCATTGCGTGGGCGTCACGAAACCAGTCATGATCATGTTCAATACGATTCATGACAATGTCATCTTTATACAATGACTAAAAGTTAAAATTCATATTTCCGGCTTTCTTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGACGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTTTTCAGCCCGCGTCGCTAATGTccctt                                                                                                                                                                                                            >  2:208945/1‑249 (MQ=255)                                                     gTCACGTAACCAGTCATGATAATGTTCAATACGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATCACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCAt                                                                                                                                                                                             >  1:137661/1‑251 (MQ=255)                                                                                   tcGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTTATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATc                                                                                                                                                                <  1:168537/249‑1 (MQ=255)                                                                                          aTGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTAGTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTAGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTTCGGCAATCAGTGCGGc                                                                                                                                                        <  2:137661/251‑1 (MQ=255)                                                                                                       ctTTATAAAATGACTAAACGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCg                                                                                                                                           >  2:73662/1‑251 (MQ=255)                                                                                                           ataAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGtttt                                                                                                                                       >  1:381555/1‑251 (MQ=255)                                                                                                                     tAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTGTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTA‑CCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCTCTTCGTTTTTCTCGGTACt                                                                                                                             <  2:128124/250‑1 (MQ=255)                                                                                                                            tAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGAt                                                                                                                       <  1:49067/250‑1 (MQ=255)                                                                                                                                    atatTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCAtt                                                                                                              >  1:39404/1‑251 (MQ=255)                                                                                                                                        ttCAGGCTTTATTTTCGGCTTTCTTATTTAAGATATTAATGGACAGATCAGTCAGTTTAATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGAGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCCTTAGct                                                                                                          >  1:127541/1‑251 (MQ=255)                                                                                                                                          cAGGCTTGATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGAGGGCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGCCCTCATTAGctc                                                                                                         <  2:109043/250‑2 (MQ=255)                                                                                                                                                     ttCGGCTTTCTTATTTACGTTATTAATGGCCAGAGCAGTCAGGTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGTTAACCTAATTGTTCAGCCAGCGTCGCTAATGGCCCGTAACTTGCAATCTCATAATGCTCGATTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGTACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGa                                                                                              <  1:246798/250‑1 (MQ=255)                                                                                                                                                       cGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCt                                                                                           <  1:132124/251‑1 (MQ=255)                                                                                                                                                           tttCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGGGGCCTTTTCTGCTTCCAGGGTTTCGTTCAGAAGCTTCGCTGCTGTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAAGCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCCCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTcc                                                                                        <  2:39404/250‑1 (MQ=255)                                                                                                                                                                     aCGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCCTTAGCTTCTTCAATAAGACCTTCCATTGCcacac                                                                              >  2:369895/1‑250 (MQ=255)                                                                                                                                                                             aaTGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCGTgc                                                                     <  1:73662/251‑1 (MQ=255)                                                                                                                                                                              aTGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATgcg                                                                    <  2:51910/251‑1 (MQ=255)                                                                                                                                                                                              gtTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCAGAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTGCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTc                                                    <  1:390308/251‑1 (MQ=255)                                                                                                                                                                                                                 ttcttcttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACt                                  <  1:225488/250‑1 (MQ=255)                                                                                                                                                                                                                       ttcCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCa                            >  1:322207/1‑250 (MQ=255)                                                                                                                                                                                                                         ccAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAAt                          >  1:223958/1‑250 (MQ=255)                                                                                                                                                                                                                                                gctgctTTACGATAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATACTGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTc  >  2:52276/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                GGCTGGCCATGGATTCAAGCATAGATTCGGCTTGCTTTTCCATTGCGTGGGCGTCACGTAACCAGTCATGATAATGTTCAATACGATTCATGAAAATGTCCTCTTTATAAAATGACTAAAAGTTAAAATTCATATTTCAGGCTTTATTTTCGGCTTTCTTATTTACGTTATTAATGGCCAGATCAGTCAGTTTGATGTCGGTGGCCTTTTCTTCTTCCAGGGTTTCTTTCAGAAGCTTCGCTGCTTTACGGTAACCTAATTGTTCAGCCAGCGTCGCTAATGTCCCGTAACTGGCAATCTCATAATGCTCGACTTTCTGTGCTGCGGCAATCAGTGCGGCATCACGCACTTCGTTTTTCTCGGTACTTTCGATGACCTCATTAGCTTCTTCAATAAGACCTTCCATTGCCACACATTTCATGCGCTTAATTTTCAGGTTCGATTCGGACTCCACAACTTGATCAATACGTTCAATCTGTCCATGAGTTTC  >  NC_000913/1315134‑1315623

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