Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	109,960	G→A	R561H (CGT→CAT)	secA →	protein translocation ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	109,960	0	G	A	100.0%	45.9 / NA	16	R561H (CGT→CAT)	secA	protein translocation ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (6/10);  total (6/10)

ATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAG  >  NC_000913/109765‑110185                                                                                                                                                                                                    |                                                                                                                                                                                                                                  aTCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAgggn                                                                                                                                                                            <  1:126395/251‑2 (MQ=255)                                    ggTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAAGTGAAAAAAGTAAAGCCGACTGGCATGTACGTCACGATGCGGTACGGGAAGCAGGTGGCCTGCAGATCATCGGTACCGAGCATCACGAATCCCGTCGTATcn                                                                                                                                                                                                               <  1:49260/181‑2 (MQ=255)                                          ggTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGCCCGCAGAGCAAATTGAAAAAATTAAAGCCGAATGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATgc                                                                                                                                   >  1:354820/1‑251 (MQ=255)                                                         ntCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCACAGAGCAAATTGAAACAATTAAAGCCGACTGGCAGGTACGCCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTCCCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTcc                                                                                                                                                             >  1:44757/2‑210 (MQ=255)                                                                           gcagACGTTGCCGCGCTGGAAAATCAGACCGCCGAGCAAATTGAAAAAATTACAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAg                                                                                                                                       >  1:412986/1‑214 (MQ=255)                                                                              gAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTc                                                                                                                                                                                             <  1:11507/157‑1 (MQ=255)                                                                                                         gCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGGCCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGc                                                                     <  1:260793/250‑1 (MQ=255)                                                                                                                       aaaaaTTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCNNTCGTAt                                                                                                                                                                                                                 >  2:26384/1‑96 (MQ=255)                                                                                                                          aaTTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTgg                                                                                                                                                                                      <  1:169881/120‑1 (MQ=255)                                                                                                                               aaGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCg                                              <  1:26384/251‑1 (MQ=255)                                                                                                                                aGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCAGATCATCGGTACCGATCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCTGTCGTTCTGGTCGTCAGGGGGAGGCTTGTTCTTCCCGGTTCTACCTGTCGATGGAAGATGCGCTGATTCGTATTTGTGCGTCCGACCTAGTGTCCGGCATGATGCTTAAACTGGGTATGAAGCCAGGCGAAGCCAGTGAACACCCGt                                             <  1:284079/251‑1 (MQ=255)                                                                                                                                        ggCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCGGCATATCAGCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGt                                                                                                                                                                                                      <  1:202930/90‑1 (MQ=255)                                                                                                                                            ggTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTaa                                   >  1:451335/1‑249 (MQ=255)                                                                                                                                                cGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATTCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCGGATGCGTATTTTTGCTTCCGACCGAGTATCCGCCATGATGCGTAAACTGGGTATg                                                                         <  1:182429/207‑1 (MQ=255)                                                                                                                                                                             ccTGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTCCGAGGTCGTTCTGGTCGTCCGGGGGACGCTGGTTCTTCCCGTTTCTAcc                                                                                                                                                   >  1:382819/1‑104 (MQ=255)                                                                                                                                                                               tGCATATCATCGGTACCGAGCATCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAACCCATTGAACGCCCGTGTGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAg  <  1:231796/247‑1 (MQ=255)                                                                                                                                                                                                    |                                                                                                                                                                                                                                  ATCCGGCTGCGGTGACTATCGCGACCAATATGGCGGGTCGTGGTACAGATATTGTGCTCGGTGGTAGCTGGCAGGCAGAAGTTGCCGCGCTGGAAAATCCGACCGCAGAGCAAATTGAAAAAATTAAAGCCGACTGGCAGGTACGTCACGATGCGGTACTGGAAGCAGGTGGCCTGCATATCATCGGTACCGAGCGTCACGAATCCCGTCGTATCGATAACCAGTTGCGCGGTCGTTCTGGTCGTCAGGGGGATGCTGGTTCTTCCCGTTTCTACCTGTCGATGGAAGATGCGCTGATGCGTATTTTTGCTTCCGACCGAGTATCCGGCATGATGCGTAAACTGGGTATGAAGCCAGGCGAAGCCATTGAACACCCGTGGGTGACTAAAGCGATTGCCAACGCCCAGCGTAAAGTTGAAAG  >  NC_000913/109765‑110185

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