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	100.0%	53.2 / NA	16	F19F (TTT→TTC)	yjgL	protein YjgL
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (4/12);  total (4/12)

AACCGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAA  >  NC_000913/4475247‑4475731                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               aaCCGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAAc                                                                                                                                                                                                                                            <  1:197720/251‑1 (MQ=255)                 gggcgagacAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGttt                                                                                                                                                                                                                             <  2:93363/242‑1 (MQ=255)                   gcgcggCAGCCCAGGTTCGCAAAGCAGTATTGAAGATGCTGATGGCGTGATGTCTTGTCGCCGTAAATAACGCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTGTTTGGATAGTTTGTGGGGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTTAAGGAGTTCTTCATGAGCAAAATAGCAGATTTAAATGATTCGCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGGGAAGTTTTaa                                                                                                                                                                                                                          <  2:251052/251‑1 (MQ=255)                                    cGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGGTCTTATTCACTCAAGGAGTTCGGCATGAGCAAAATATCAGATTTAAATTATTCGCAACACATTACAGTAGCCGACAATTTCAAACAAAAAAGCGAAGTTGTTAATACCTGGCGTGTTgg                                                                                                                                                                                                          <  2:253994/250‑1 (MQ=255)                                               agtgatggtGCTGATGGCGTGATGACTTTTCGCCGTAAATAACGCAAGGTTTACGGCTAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACTTTCTTATTCACTCAAGGAGTTGGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATGACAGTAGCCGACAATTTCAAACAAAAAAGTGAAGTGTTAAATACCTGGCGTGTTGGAATGAATGAttt                                                                                                                                                                                              <  2:165895/243‑1 (MQ=255)                                                                          tttCGCCGTCAATAACTCAAGGTTTACGACGAGTTTTTGCAATGAGCCATTTTTGATTTTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAAGAGTTTGTCATGAGCAAATAATCAGATTTAAATTATTCTCAACACATTACATTATGCGACAATTTCAAACCAAAAAGAAAATGTTTACATACCTGGCGTGTATGAATGAATGATTTTTCTCGTCTTTCCGGGGGGCAGGATa                                                                                                                                                                    >  2:178923/1‑250 (MQ=255)                                                                               ccGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTCTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTAAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATACCAGa                                                                                                                                                               >  1:17654/1‑250 (MQ=255)                                                                                cGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGTGGCAGGATAACAGAAg                                                                                                                                                             <  1:83924/251‑1 (MQ=255)                                                                                         cTCAAGGTTTACGGCGAGTTTCTGAAAAGAGCGTTTTTTGATCTTTTTTTGTGAGTAAAATTTGTAATGCTTAGGCGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATAttct                                                                                                                                                    >  2:179943/1‑251 (MQ=255)                                                                                                                             tttGATCTTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTCGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCaa                                                                                                                    >  1:160000/1‑247 (MQ=255)                                                                                                                                gATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGatat                                                                                                              <  1:98276/250‑1 (MQ=255)                                                                                                                                                                 aGACGGTCTTATTCACTCAAGGAGGTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATCCCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGTCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGAGTTTAGCAAACGCt                                                                             <  1:69010/250‑1 (MQ=255)                                                                                                                                                                           aTTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGc                                                                   <  1:245645/250‑1 (MQ=255)                                                                                                                                                                                              cATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCt                                                <  1:119101/250‑1 (MQ=255)                                                                                                                                                                                                                                      acattaGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCa        <  1:79690/250‑1 (MQ=255)                                                                                                                                                                                                                                           aGCCGACAATTTCAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGaaa  <  1:131973/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                               AACCGCTGTGCGCAGTGGGCGCGCCAGCCCAGGTTCGCAAAGCAGTAGTGAAGATGCTGATGGCGTGATGACTTTTCGCCGTAAATAACTCAAGGTTTACGGCGAGTTTGTGAAAAGAGCGTTTTTTGATATTTTTTTGTGAGTAAAATTTGTAATGCTTAGACGTTCTTATTCACTCAAGGAGTTCGTCATGAGCAAAATATCAGATTTAAATTATTCTCAACACATTACATTAGCCGACAATTTTAAACAAAAAAGTGAAGTTTTAAATACCTGGCGTGTTGGAATGAATGATTTTGCCCGTATTGCCGGGGGGCAGGATAACAGAAGGAATATTCTTTCTCCTGGAGCATTTTTAGAGTTTTTGGCAAAGATATTTACCCTGGGTTATGTGGATTTTAGCAAACGCTCCAACGAAGCGGGTAGAAATATGATGGCTCATATTAAGTCCTCATCTTATTCTAAAGATACTAATGGCAATGAAA  >  NC_000913/4475247‑4475731

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