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	95.7%	65.0 / ‑2.6	23	coding (211/510 nt)	def	peptide deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/10);  minor base . (0/1);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.13e-01

CAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTGAAA  >  NC_000913/3433694‑3434142                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   cAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACa                                                                                                                                                                                                        <  1:188000/251‑1 (MQ=255)                     atTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAATCCGTGACGAACGGTTAGTGTTAATCAAGCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGcc                                                                                                                                                                                     <  1:239992/250‑1 (MQ=255)                       tCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTg                                                                                                                                                                                   <  1:115751/250‑1 (MQ=255)                         cGGACGAGCGGCTTCGCAAAGTTGCTAAACAGGTAGAAGAAGCGATTACAGCAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGACTGGCGGCAACCCAGGCTGATATCAATCAACGTATCATTGTTATTGATGTTTAGGACAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAACGCGGCGAAACAGGCATTGAAGAAGGttcaatgtcg                                                                                                                                                                                >  1:109127/1‑243 (MQ=255)                             cGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATcc                                                                                                                                                                            >  1:203993/1‑251 (MQ=255)                                 cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTg                                                                                                                                                                         >  1:205491/1‑250 (MQ=255)                                 cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTg                                                                                                                                                                         >  1:48479/1‑250 (MQ=255)                                 cGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGGTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGCAGGATGCCTGTCGATCCCTg                                                                                                                                                                         >  2:51910/1‑250 (MQ=255)                                       cGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACg                                                                                                                                                                  >  1:101915/1‑251 (MQ=255)                                          aaaGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGc                                                                                                                                                               >  1:29702/1‑251 (MQ=255)                                                                                tCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGcc                                                                                                                          <  2:101915/250‑1 (MQ=255)                                                                                tCAGCGTATCGTCGATGATATGGTCGAGACGATGTACGCAGAATAAGTTATGGGCCTGGCGGCAACCCAGGTGTATATCCATCAACGTATCATTGTTATTGATGTGTCGGAAATCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGcc                                                                                                                          <  1:143549/250‑1 (MQ=255)                                                                                    cGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGACAACAGTGACGAACGGTTAGTGATAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTg                                                                                                                      >  1:266888/1‑250 (MQ=255)                                                                                                     gTTCGAGACGAGGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGATGCCTGTCGATCCCTGAACAACGTGCTTTAGGGCCGCGCTCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCAt                                                                                                    <  2:267838/251‑1 (MQ=255)                                                                                                               aTGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTgg                                                                                           >  1:35570/1‑250 (MQ=255)                                                                                                                            aaggaattgacttGGCGGCAACCCAGGTGGATATCCATCAACGAATCATTGTTATTGATGCTTAGGAAAACCGTGACGAACGGTTAGAGTTAATCAAACCAAAGCTTTTAGAAAAAAGCGGCGAAAAAGGCATTGAAGAAGGTTGCCTGTAGATCAATGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAAATGGAAGCAGACGGTCTg                                                                             >  2:157421/13‑251 (MQ=255)                                                                                                                                    ggCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCa                                                                      <  2:205491/250‑1 (MQ=255)                                                                                                                                    ggCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTGAGAAAAAAGCGGCGAAACAGGCAGTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCGTTAGTGCCGCGCGCAGAGAAAGTTAAAATGCGCGCGCTTGACCGCGACGGTAAGCCATTTGAACTGGAAGCAGACGGTCTGTTAGCCa                                                                      <  2:48479/250‑1 (MQ=255)                                                                                                                                       cTGGCGGCAACCCAGGATGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGAGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTg                                                                  >  1:221705/1‑251 (MQ=255)                                                                                                                                          gcggcAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTGCGGAAAACCGTGACGAACGG‑TAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTAtt                                                              <  2:87158/251‑1 (MQ=255)                                                                                                                                                                                                aaCCGTGACGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCAcc         <  1:226497/251‑1 (MQ=255)                                                                                                                                                                                                   agTGACGTACGGTTAGTGTTAATCAATCCAGAGCTTTTAGACAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCt      >  2:16918/2‑251 (MQ=255)                                                                                                                                                                                                       aCGAACGGTTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACTGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTGaaa  <  1:26638/251‑1 (MQ=255)                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   CAGTTTTGCAAGTGTTACATATTCCGGACGAGCGGCTTCGCAAAGTTGCTAAACCGGTAGAAGAAGTGAATGCAGAAATTCAGCGTATCGTCGATGATATGTTCGAGACGATGTACGCAGAAGAAGGTATTGGCCTGGCGGCAACCCAGGTTGATATCCATCAACGTATCATTGTTATTGATGTTTCGGAAAACCGTGACGAACGGCTAGTGTTAATCAATCCAGAGCTTTTAGAAAAAAGCGGCGAAACAGGCATTGAAGAAGGTTGCCTGTCGATCCCTGAACAACGTGCTTTAGTGCCGCGCGCAGAGAAAGTTAAAATTCGCGCCCTTGACCGCGACGGTAAACCATTTGAACTGGAAGCAGACGGTCTGTTAGCCATCTGTATTCAGCATGAGATGGATCACCTGGTCGGCAAACTGTTTATGGATTATCTGTCACCGCTGAAA  >  NC_000913/3433694‑3434142

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