Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,475,493	T→C	F19F (TTT→TTC)	yjgL →	protein YjgL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,475,493	0	T	C	96.7%	94.7 / NA	30	F19F (TTT→TTC)	yjgL	protein YjgL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (13/16);  minor base A (1/0);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTA  >  NC_000913/4475254‑4475710                                                                                                                                                                                                                                                |                                                                                                                                                                                                                          gtgCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACaaaaa                                                                                                                                                                                                                  <  2:5491/249‑1 (MQ=255)          gggCGAGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGCTATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTAACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGtttt                                                                                                                                                                                                       >  2:328386/1‑251 (MQ=255)                          gTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGtt                                                                                                                                                                                       <  1:526419/251‑1 (MQ=255)                                aaaGCAGTAGTGAAGCTGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGCTTTAAATTATTCTCAAAACATTACATTAGCAGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGAGTGTTGgaatg                                                                                                                                                                                 >  2:164917/1‑251 (MQ=255)                                aaaGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGgaatg                                                                                                                                                                                 >  2:31806/1‑251 (MQ=255)                                  aGCAGTAGTGAAGATGCTGATGGTGTTATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACGTGGCGTTTTtgaatgaa                                                                                                                                                                               <  2:103125/251‑1 (MQ=255)                                                                      cGCCGTAAATAACTCAAGGTTTACGGCGAGGTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATAGCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTTGAATGAATGATTTGTCCCGTAGTTCCGGGGGTCATTATAACAg                                                                                                                                           <  2:360196/251‑1 (MQ=255)                                                                           tAAATAACTCAAGGATTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGCGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAAAAGAAgg                                                                                                                                       >  1:375427/1‑250 (MQ=255)                                                                                 aCTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATAttc                                                                                                                                <  1:163556/251‑1 (MQ=255)                                                                                    cAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTGTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATGCACTCAAGGAGTTCGTCAGGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTTAAGTTTTAAATGCCTGGCGTGTTGGATTGATTGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATAttctt                                                                                                                              <  1:37383/250‑1 (MQ=255)                                                                                               gCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCa                                                                                                                  <  1:164232/251‑1 (MQ=255)                                                                                                cgattttGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAt                                                                                                                 <  1:320485/247‑1 (MQ=255)                                                                                                   gTTTGTGAAAATAGCGATTTTAGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGAAGTTCTTATTAACTAAAGCAGTTCGACATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCAGACAATTTAAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGTGGGCAGGATAACAGAAGGAATATTCTTTCTCATGGAGAAtttt                                                                                                              >  2:527530/1‑251 (MQ=255)                                                                                                     ttGTGAAAAGAGCGTTTTTTGATATTTTTTAGTGAGTAAAATTTGTAATGCTTAGACGATCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATAAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAACAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCAGGGGGGCAGTATAACAGCAGGAATATTATTTCTCCTGGAGCAttttt                                                                                                             >  2:80284/1‑250 (MQ=255)                                                                                                        tgAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTag                                                                                                           >  1:9514/1‑249 (MQ=255)                                                                                                                                   tgtgAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCTTGAGCAAAATATCAGATTTAAATTATTCTAAACACATTACATTCGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTgggg                                                                              >  2:220291/1‑250 (MQ=255)                                                                                                                                    gtgAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGTCAGGATAACAGAAgg                                                                                                                                       <  2:399685/193‑1 (MQ=255)                                                                                                                                    gtgAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAgg                                                                                                                                       >  1:399685/1‑193 (MQ=255)                                                                                                                                                            aCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCa                                                     <  1:31806/251‑1 (MQ=255)                                                                                                                                                                           cAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGa                                       <  1:451235/250‑1 (MQ=255)                                                                                                                                                                                    cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg                             >  1:363439/1‑251 (MQ=255)                                                                                                                                                                                    cGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATgg                             >  1:113878/1‑251 (MQ=255)                                                                                                                                                                                       cATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCt                           <  1:544814/250‑1 (MQ=255)                                                                                                                                                                                           aGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTTGCCGACAATTTCAAACAAAAAAGGGAAGTTGTAAATACCTGGCGTGTTGGAATTAATGATTTTGCCCGTAGTGCCGGGTCGCATGATTACAGAAGAAATATTCTTTCTCCTGGAGCATTTTTAGAGTATTTGGCAAATATAGTGACCCTGGCTTTTGTTGATTTTAGCAAACGCTCCAATGAAGCGGGTATAAATATCATGGCTCatat                      <  2:555995/251‑1 (MQ=255)                                                                                                                                                                                               aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                                                >  1:78796/1‑157 (MQ=255)                                                                                                                                                                                               aaaTATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGGGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCAtt                                                                                                                <  2:78796/157‑1 (MQ=255)                                                                                                                                                                                                 atatCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTc                <  1:432188/251‑1 (MQ=255)                                                                                                                                                                                                   atCAGATTTCAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTTAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGAACGTATTGCCGGGGTGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAACGATATTTACCATGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGTCTCTTATTAAGTCCt              >  2:55974/1‑251 (MQ=255)                                                                                                                                                                                                   atCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCAGATTAAGTCCt              <  1:467252/251‑1 (MQ=255)                                                                                                                                                                                                                tattCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGCGAAGTTTTAAATACCTGGCGTGACGGAATGAATGATTTTGCCCGTATTGACGGGGGGCAGGATAACAGAAGGAATATTATTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTa  >  1:77966/1‑250 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                          GTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTA  >  NC_000913/4475254‑4475710

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: ‑