Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	466,148	G→A	P389S (CCA→TCA)	ybaE ←	uncharacterized protein YbaE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	466,148	0	G	A	100.0%	57.6 / NA	18	P389S (CCA→TCA)	ybaE	uncharacterized protein YbaE
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/10);  total (8/10)

AATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGA  >  NC_000913/465921‑466365                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          aaTTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCAGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGAATGAAGTTCTATCGGTAGGAGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTgg                                                                                                                                                                                                    >  1:273093/1‑250 (MQ=255)               ccATCATGAGGTCTGCGGGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCg                                                                                                                                                                                      <  1:285416/251‑1 (MQ=255)                           cTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGtt                                                                                                                                                                          >  1:273782/1‑251 (MQ=255)                                   gCCTGTAGGTTCGTGTCGTCCCTGTTTTTTGCGTTATGAAAAATAATTGTTAGTGCACAGCCTTCCGCTGCCAGGGTCGCCTGTAGGCGTTCTGCCATGGTATTAAGTTCTATCGGTAGGTGATAACCCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGtaataa                                                                                                                                                                  <  2:258746/251‑1 (MQ=255)                                         aGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCgg                                                                                                                                                             <  1:43682/250‑1 (MQ=255)                                            gtcgtgtcgtCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATtg                                                                                                                                                         <  2:29010/251‑1 (MQ=255)                                               gtgtcgtCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCAGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGCTCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAAACGGATtgatg                                                                                                                                                      >  2:98431/1‑251 (MQ=255)                                                 gtcgtcCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGaa                                                                                                                                                    <  1:211991/251‑1 (MQ=255)                                                        cAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAGGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGGTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGga                                                                                                                                             <  2:131410/251‑1 (MQ=255)                                                               ttGCTTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCAc                                                                                                                                       <  2:69666/250‑1 (MQ=255)                                                                                       gTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCCGAGGGAGAGt                                                                                                              >  1:50095/1‑251 (MQ=255)                                                                                                 ttCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGCTAACCCAGCGTCCAGGTTTTCGTTCGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTTCTCCGGTTTTTTCCGCCTTCTAACGTTTGTTATAAACCCGCTTTATTAATCCTGGCTATCCCTTTTCTCGCCTGCCCGATGGAGAATCGGGGACttt                                                                                                    >  2:108095/1‑251 (MQ=255)                                                                                                         ccAGTGTCGCCTGTAGGCGTTCTGCCATGTTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAGGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAAc                                                                                            <  2:141376/251‑1 (MQ=255)                                                                                                          cAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTGGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATTAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGGGAGTCGGGGGCTTTTGCGCAACg                                                                                           >  2:151050/1‑251 (MQ=255)                                                                                                                                                                              ttttCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCAAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTTATGAATAATGGAGATCACTTTTCGCGCCTGCCAAAGGGAGAGTCGGGGACTTTTGCGCCAAGTCAAATAACAAAAACCTAAACTGATGCCGCTAATGAACCGGCTGACCCGCTGCCTCTCCTCCGGtttt                       >  2:38461/1‑250 (MQ=255)                                                                                                                                                                                    gTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCGGGCTGACCCGTTGCAGCGCCTCCGGGTTGCCGATg                 <  1:108095/251‑1 (MQ=255)                                                                                                                                                                                                cATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACg     <  1:11179/251‑1 (MQ=255)                                                                                                                                                                                                   ccGGTACCTGCCAATGCGGAATAGTCCAGCCTGACAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTTTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACCTTTTCGCAACGGCAAATAGAAAAAACCTAAACTGATGCCGCCACTGACCTGGCCGACCCGTTGCAGACCCTCCGGGTTGCCGATGGTGATTTGCCCgggg  >  1:140784/1‑250 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                          AATTAATCTGTCGCCCATCATGAGGTCTGCGTGTGCCTGTAGGGTCGTGTCGTCCCAGTTTTTTGCGTTATGAAAAATAATTGTGAGTTCACAGCCTTCCGCTGCCAGTGTCGCCTGTAGGCGTTCTGCCATGGTATGAAGTTCTATCGGTAGGTGATAAACCAGCGTCAAGGTTTTCGGTAGTTTGACTTCATCCGGTACCTGCCAATGCGGAATAGTCCAGCCTGGCAGTAATGCATGACTGGCGGTGATCAGGTTTTCTCCGACTTCTAACGTTTGTAATAAACCGGATTGATGAATAATGGAGATCACTTTTCGCGCCTGCCAGAGGGAGAGTCGGGGACTTTTGCGCAACGTCAAATAGCAAAAACCTAAACTGATGCCGCTACTGACCTGGCTGACCCGTTGCAGCTCCTCCGGTTTGCCGATGGTGATTTGCACGGGA  >  NC_000913/465921‑466365

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