Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	554,631	T→A	T7T (ACT→ACA)	cysS →	cysteine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	554,631	0	T	A	100.0%	80.7 / NA	23	T7T (ACT→ACA)	cysS	cysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (13/10);  total (13/10)

GAAATATGGGTATTATACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACTCTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGTTT  >  NC_000913/554554‑554707                                                                              |                                                                             gAAATATGGGTATTATACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATAcact                                                                            >  1:121041/1‑80 (MQ=255)  aaaTATGGGTATTATACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGa                                                                          >  2:121024/1‑81 (MQ=255)       tGGGTATTATACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACg                                                                      <  1:120997/80‑1 (MQ=255)         ggTATTATACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCa                                                                   <  2:120998/81‑1 (MQ=255)           tattatACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCaa                                                                  >  1:121026/1‑80 (MQ=255)             ttatACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCaaaaa                                                               <  1:120999/81‑1 (MQ=255)                taCGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAagag                                                            <  2:121000/81‑1 (MQ=255)                   gCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGaa                                                         <  2:121001/81‑1 (MQ=255)                    cAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAAt                                                        >  1:121040/1‑81 (MQ=255)                       cTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAAttt                                                      >  2:121025/1‑80 (MQ=255)                       cTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAAttt                                                      >  1:121034/1‑80 (MQ=255)                          aaatACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGc                                                  <  1:121002/78‑1 (MQ=255)                                ccACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTc                                            <  1:121003/81‑1 (MQ=255)                                   cacaTGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACg                                         <  1:121004/81‑1 (MQ=255)                                       tGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCgg                                     >  2:121032/1‑81 (MQ=255)                                         tCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCggg                                    <  1:121005/80‑1 (MQ=255)                                           tAAACGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCgggg                                   >  2:121038/1‑79 (MQ=255)                                            aaaCGGAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAg                                >  1:121036/1‑81 (MQ=255)                                                 gAATCTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCgg                            >  2:121044/1‑80 (MQ=255)                                                     cTTCGATGCTAAAAATCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATg                        >  2:121035/1‑80 (MQ=255)                                                                 aaaTCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATGTACGTGTGTGGaa           >  2:121030/1‑81 (MQ=255)                                                                 aaaTCTTCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATGTACGTGTGTGGaa           >  1:121033/1‑81 (MQ=255)                                                                      ttCAATACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCAcc      >  1:121055/1‑81 (MQ=255)                                                                          aTACACTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGttt  <  2:121006/81‑1 (MQ=255)                                                                              |                                                                             GAAATATGGGTATTATACGCAACTCAATTACCCACACATGTCTAAACGGAATCTTCGATGCTAAAAATCTTCAATACTCTGACACGCCAAAAAGAGGAATTTAAGCCTATTCACGCCGGGGAAGTCGGCATGTACGTGTGTGGAATCACCGTTT  >  NC_000913/554554‑554707

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