Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	5,201	A→G	intergenic (+181/‑33)	thrC → / → yaaX	threonine synthase/DUF2502 domain‑containing protein YaaX

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	5,201	0	A	G	100.0%	69.9 / NA	21	intergenic (+181/‑33)	thrC/yaaX	threonine synthase/DUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (10/11);  total (10/11)

TATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGT  >  NC_000913/5076‑5324                                                                                                                               |                                                                                                                            tataAAAAATTTATGGAAAAAAAGGACAGGA‑AAAAAGGAAAAATTCTAAAAAAATGCGGTAATTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATaa                                                                                                                      <  2:42987/133‑1 (MQ=255)     aaaaaaTTAGGAAAAAAAAGGACAGGA‑AAAAAAAAAAAATTAACAATAATGGCGGTAACTTAGAAATTAGGATGGGGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACggg                                                                                                                  <  2:47618/130‑1 (MQ=255)     aaaaaaTTAAGGAAAAAAAAGAAAGGA‑AAAAAGGAAAAATTCTAAATAAATGGGGTAACTAAGAGATTAGGATGGGGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACggg                                                                                                                  <  2:399930/130‑1 (MQ=255)        aaaaTAAGGAAAAAAAAGAAAGGGAAAAAAGAAAAAATCCTAAAAAAATGGTAAAATTAAAAGATAAGGATGGGGGAAAAAAAAAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCa                                                                                                                <  2:5531/133‑1 (MQ=255)            tagggAAAAACAGGATAGGA‑AAAAAGAATAATATCTAAATAGATGGGGTAATTAAGAGATTAGGATGGCGGAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGAt                                                                                                           <  2:205124/130‑1 (MQ=255)              tGGAGAAAAAGGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCACCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATaa                                                                                                         <  2:213933/133‑1 (MQ=255)              ggaaaaaaaaaGACATAA‑AAAAAGAAAATATCCTAAAAAATCGCGTTAACTAAGAGATTAGGATTGTGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATaa                                                                                                         <  2:175352/128‑1 (MQ=255)                    aaaaGGAAAGGA‑AAAAAGGAAAAATTCTAAATAAATGCGGTAACTAAGAGATTAGGATTGCGAAAAAAAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAg                                                                                                   <  2:318301/133‑1 (MQ=255)                                                                       aTTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCt                                                 >  1:573136/1‑133 (MQ=255)                                                                         tAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                                               >  2:311278/1‑133 (MQ=255)                                                                           ggtttGCGAAAAATAAAATCCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                                             <  1:343185/130‑1 (MQ=255)                                                                           ggATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAATAGATGCAAGCTAGCGTCCTCGCACTGTCCCTGGTTCTGGTc                                             >  2:177696/1‑133 (MQ=255)                                                                           ggATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                                             >  1:372084/1‑133 (MQ=255)                                                                                 cGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTCCTCGCACGTTCACTAGTTCTGCTCGCTccc                                       >  1:299243/1‑133 (MQ=255)                                                                                    aaaaaaacaacCGCCGTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATg                                    <  2:399523/128‑1 (MQ=255)                                                                                                gccgTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt                        >  1:562037/1‑133 (MQ=255)                                                                                                gccgTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt                        >  1:477596/1‑133 (MQ=255)                                                                                                gccgTTCTCATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCt                        >  2:88375/1‑133 (MQ=255)                                                                                                        cgtCGAGTATTCTCCCGATATCGCGCAAAAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAAt                <  2:299243/131‑1 (MQ=255)                                                                                                        cATCGAGTAATCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAAt                <  1:311278/133‑1 (MQ=255)                                                                                                                 aTCTCCGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGt       >  2:350133/1‑133 (MQ=255)                                                                                                                      cGGATATCGGCCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGt  >  2:479989/1‑133 (MQ=255)                                                                                                                               |                                                                                                                            TATGAAGAAATTATGGAGAAAAATGACAGGG‑AAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGGAAATTACGTTAGTCCCGT  >  NC_000913/5076‑5324

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