Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	945,259	A→G	V100A (GTC→GCC)	ycaC ←	putative hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	945,259	0	A	G	100.0%	86.6 / NA	25	V100A (GTC→GCC)	ycaC	putative hydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (8/17);  total (8/17)

CATTAAAGGTCCCTGACGCGTCAGTCACGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGACAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAACTAATGGGCCGTTAGGACCGGTTTCAAAACTGGT  >  NC_000913/945125‑945386                                                                                                                                       |                                                                                                                                cATTAAAGGTCCCTGACGCGTCAGTCACGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCt                                                                                                                            <  2:356418/140‑1 (MQ=255)             cTGACGCGTCAGTCACGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCt                                                                                                                >  1:312614/1‑140 (MQ=255)                aCGCGTCAGTCACGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCg                                                                                                             <  1:213723/140‑1 (MQ=255)                    gtcagtcaCGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTAGCAGCTTTTACAAAATCTTCGTTAt                                                                                                         >  2:793252/1‑140 (MQ=255)                        gtcaCGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCa                                                                                                     <  2:213169/140‑1 (MQ=255)                                     cGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATAttt                                                                                        >  1:508352/1‑140 (MQ=255)                                        cAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTccc                                                                                     <  2:312614/140‑1 (MQ=255)                                              cttcttCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGa                                                                               <  2:83599/140‑1 (MQ=255)                                              cttcttCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGa                                                                               <  2:226498/140‑1 (MQ=255)                                              cttcttCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGa                                                                               <  1:283494/140‑1 (MQ=255)                                                 cttcAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCa                                                                            <  1:111008/140‑1 (MQ=255)                                                 cttcAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCa                                                                            <  2:169252/140‑1 (MQ=255)                                                 cttcAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCa                                                                            >  2:197848/1‑140 (MQ=255)                                                    cAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATg                                                                         <  2:147127/140‑1 (MQ=255)                                                       tGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTaa                                                                      >  2:329343/1‑140 (MQ=255)                                                            gAAAGCCCCGGGCATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTCTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTGATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGgtgc                                                                 >  2:102654/1‑139 (MQ=255)                                                                  gcCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCggg                                                           >  1:594706/1‑140 (MQ=255)                                                                              aCGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCttt                                               <  2:332412/140‑1 (MQ=255)                                                                                    aCTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTc                                         <  1:318162/140‑1 (MQ=255)                                                                                          gTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAAc                                   >  1:560431/1‑140 (MQ=255)                                                                                                acacCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCttt                                               >  2:589500/1‑122 (MQ=255)                                                                                                acacCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCttt                                               <  1:589500/122‑1 (MQ=255)                                                                                                            atTAACTGTTTTTTACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAACTAATGGGCCGTTAGGAcc                 <  2:594706/140‑1 (MQ=255)                                                                                                                      ttttACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAACTAATGGGCCGTTAGGACCGGTTTCaaaa       <  1:436846/140‑1 (MQ=255)                                                                                                                         tACCTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAACTAATGGGCCGTTAGGACCGGTTTCAAAACtg    <  2:560431/140‑1 (MQ=255)                                                                                                                           ccTGTCGCTTTGGCAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAACTAATGGGCCGTTAGGACCGGTTTCAAAACtggt  <  2:604154/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                CATTAAAGGTCCCTGACGCGTCAGTCACGACAAAGACGTCAAACCCTTCTTCAATGGCTGAAAGCGCCGGGAATGCCACGCAAACTTCGGTTACCACACCGGCAATAATTAACTGTTTTTTACCTGTCGCTTTGACAGCTTTTACAAAATCTTCGTTATCCCAGGCGTTAATATTTCCCGGGCGAGCAATGTAAGGTGTATCGGGGAATTGTGCTTTCAGTTCTGGAACTAATGGGCCGTTAGGACCGGTTTCAAAACTGGT  >  NC_000913/945125‑945386

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