Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,433,900	C→T	L71L (CTA→TTA)	def →	peptide deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,433,900	0	C	T	100.0%	53.1 / NA	19	L71L (CTA→TTA)	def	peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/11);  total (8/11)

ACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATC  >  NC_000913/3433664‑3434128                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     aCTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTaatcaat                                                                                                                                                                                                                        <  1:182843/251‑1 (MQ=255)          cacaTCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCtt                                                                                                                                                                                                               >  1:160822/1‑251 (MQ=255)                   ggAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGGTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGGGTTAATCAATCCAGAGCTTTTAGaaaa                                                                                                                                                                                                       >  2:53875/1‑250 (MQ=255)                                            gttacaaaattcgggcGAGCTGCTTCGCAAAGGTGCTAAATCGGTAGAAGAAGTAAATGCAGAAATGCCGCGTATGGTCGATGATATGTTCGAGACGATGTACGCAGAAGAATGTATTGGCCTGGCGGCAACCCAGGTTGATTTCCATCAACGTAACATTGTTAGTTATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAAGCAATCCAGAGCTGTTAGAAAAAAGCGGCGAAACAGGCATTGaagaag                                                                                                                                                                             <  2:60081/236‑1 (MQ=255)                                               aCATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGtt                                                                                                                                                                          <  1:120902/251‑1 (MQ=255)                                                           cGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCATTCAAAGGCTTTTAGAAAAAAGCGGCGAACCAGGCATTgaaaagggttcccttccgttcc                                                                                                                                                              >  2:270256/1‑232 (MQ=255)                                                                         aaGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCTTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCt                                                                                                                                                <  1:293525/251‑1 (MQ=255)                                                                           gTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCttt                                                                                                                                              <  2:216038/251‑1 (MQ=255)                                                                                              aagTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGaacaa                                                                                                                                                      >  1:282701/1‑224 (MQ=255)                                                                                              aagTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGaacaa                                                                                                                                                      <  2:282701/224‑1 (MQ=255)                                                                                                           aaTTCAGCGGATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCAGCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTcgcg                                                                                                              <  2:250818/251‑1 (MQ=255)                                                                                                                            atgatATGGTCGATACAATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATAACCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAAGCAAGCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGt                                                                                             <  2:80610/251‑1 (MQ=255)                                                                                                                                           cGATGTACGCCGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAACACCGTGACGAAAGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGAATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGAGCCGCGCGCAGAGAAAGTTAAAATTCGAGCCCTTGACCGCGACGGTAACCCATTTAAACTg                                                                              >  2:272285/1‑251 (MQ=255)                                                                                                                                                    cAGAAGAAGGTATTGGCCGGGCGGCAACCCAGGATGCTATCCATAAACGTATCATTGTTATTGAAGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCACAACAGGCATTGAAGAAGGCTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGACAGTTAAAATTCGTGCACGTGACCGCGCCGGTAAACCATTTGCACTGGAAGCAgcc                                                                     >  2:106525/1‑249 (MQ=255)                                                                                                                                                                       ggcggcAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCACGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGt                                                   >  1:103393/1‑250 (MQ=255)                                                                                                                                                                        gcggcAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTAt                                                 <  1:206028/251‑1 (MQ=255)                                                                                                                                                                                                  aCGTATCAGTTTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGTCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACACCGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATCCATCACCTGGTcc                       <  2:136893/251‑2 (MQ=255)                                                                                                                                                                                                      aTCATTGTTATTGATGTTTCGGGAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCa                    <  1:270256/250‑1 (MQ=255)                                                                                                                                                                                                                       ttCGGAACACCGTGACGAACGGTTAGTGTTAATCAATCCATAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGCCCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGAGCACCAGGTCGGCAAACTGTTTATGGATTATc  >  2:135079/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                     ACTCCTGAACACATCTCTGGAGATTTATGTCAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATC  >  NC_000913/3433664‑3434128

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