Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,216,329	C→T	I284I (ATC→ATT)	aceB →	malate synthase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,216,329	0	C	T	90.0%	46.9 / ‑1.4	21	I284I (ATC→ATT)	aceB	malate synthase A
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (10/8);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

AGGGCAGTGGTCCCTATTTCTATCT‑GCCGAAAACCCAGTCCTGGCAGGAAGCGGCCTGGTGGAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATCAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAG‑CACAATAACCAGGTGCTCAACAAAGTAAAAGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCAC  >  NC_000913/4216091‑4216568                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  aGGGCAGTGTTCCCTATTTCTTTCT‑GCCGAAAACCCAGTCTTGGCAGGAAGCGGCCTGTGGGATCGAAGTCTTCAGCTAGGCAGAAGATCGCTTTAATCTTCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCTCTGCGTGACCATATTGTTTGTCTGAACTGCGGTCGTTGGGAGTACATCTTCAGCTATATTAAAACGTTGaaa                                                                                                                                                                                                                                       <  1:287547/251‑1 (MQ=255)    gcatgggTCCATATTTGTATCTCGCC‑AAAACCCAGTCCGGGCAGGAAGCGGCCTGGTGGAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGGCGTGTCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGTATTACATCTTCAGCTATATCAAAACGTTGAAAAAc                                                                                                                                                                                                                                    <  2:142174/246‑1 (MQ=255)                aTTTCTATCT‑GCCGAAAACCCAGTCCTGGCAGGAAGCGGCCTGGTGGAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATc                                                                                                                                                                                                                          <  1:448600/249‑1 (MQ=255)                     tATCT‑GCCGAAAACCCAGTCCTGGCAGGAAGCGGCCTGGTGGAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCt                                                                                                                                                                                                                   >  2:201662/1‑251 (MQ=255)                      aTCT‑GTCGAAACCCCAGTCCTGCAAGGACGCGGCCTGGTGGAGCGAAGTCTCCAGCTATGCAGAAGACCGCTTTAATCTGCCACGCGGCGCCATCAAGGCGACGTTGCTGATGGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCt                                                                                                                                                                                                                   <  1:449401/250‑1 (MQ=255)                             cGAAAACCCAGTCCTGGCAGGAAGCGGCCTGGTGGAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCGGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCcagac                                                                                                                                                                                                            <  1:149128/251‑1 (MQ=255)                                                           ggtggAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTc                                                                                                                                                                              <  1:201662/251‑1 (MQ=255)                                                                              ttATGCAGGAGATCGTTTTAATCTGCCGCGCGGCACCATCCAGACGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGCTGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTATGAACTGCGGTCGATGAGATTACAACTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTACTGACAGACAGACAGGCAGTAACGATGGATAAAcctttacagactgattacacactcc                                                                                                                                                           >  2:91681/2‑228 (MQ=255)                                                                               tATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATAAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGACTGCTTACTCACGCCt                                                                                                                                                          >  1:270449/1‑251 (MQ=255)                                                                                  gCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGtt                                                                                                                                                       <  1:5541/251‑1 (MQ=255)                                                                                        gaTCGCTTTAATCTGCCGCGAGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCACGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATCTCAAAACGTTGAACAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCCTTCCTGAATGCTTACTCACACCTGTTGATTa                                                                                                                                                  >  2:211936/1‑250 (MQ=255)                                                                                                              gCACCATCAAGGCGCCGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGGCTGAACTGCGGTCGTTGGAATTACATCTTCAGCTCTATTAAAACGTTGAAAAAATATCCCGATCGCGTCCTGCCAGACAGCCAGGCAGTGCCGCTGGATAAACCATTACTGAATGCTTACTCACGCCTGGTGATTAAAACCTGCCATAAC‑AGCGGTGCt                                                                                                                          >  1:66106/1‑251 (MQ=255)                                                                                                               cACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCGGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAA‑CGCGGTGCtt                                                                                                                         >  1:81854/1‑251 (MQ=255)                                                                                                                       aGGCGACGTTGCTGATTGAAAGGCTGCCCGCAGCGTGCCAGATGGAGGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAAATGCGGGAGATGGGATTGCATCTTACGCTATATTAAAACGTTGACAAACTATCCCTATCGCGTCATGCCAGACAGGAAGGCAGTGACGATGGATAATCGATTACTGCATGCTTACTACCGCCAGTTAATTAAAACCTGC‑ATAAAACCCGGTGCTATTGCGATg                                                                                                                 >  2:234803/1‑251 (MQ=255)                                                                                                                                                               tcagagtgtgctctaccgatctcTGCGTTACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACGCACGCCTGTTGATTAAACCCTGCCATAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAagag                                                                          <  1:38035/228‑1 (MQ=255)                                                                                                                                                                       taatgcTTCACGTGCTGCGTGACCGTATTGTTTGTCTGAACAGAGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGCAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAGAACCTGCCATAAC‑CGTGGTGCTTTTGCGCGGGGCGTCATGGCGGAGTTTATTCCAAGCAAAGATGAAGAGACACaaaac                                                                 >  2:94296/6‑248 (MQ=255)                                                                                                                                                                                     cTGCGTGACCATATTGTTGGACTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGACTGCTTACTCACGCCTGTTGATTAAAACCTGCCCTAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGGGGCGTTTCTTCCGAGCAAAGATGAagag                                                                          >  2:38035/1‑228 (MQ=255)                                                                                                                                                                                                    gttgGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAG‑CACAATAACCAGGTGCTCAACAACGTAAAAGCGGATa                                    >  1:37344/1‑250 (MQ=255)                                                                                                                                                                                                                             aTTACACCTTCAGATATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCGGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAG‑CACAATAACCAGGTGCTCAACAAAGTAAAAGCGGATAAATCGCTGGAAGCCAATAACGGTCAc          <  1:445714/251‑1 (MQ=255)                                                                                                                                                                                                                                  aTCTTCAGCAATATTAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCAAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTACAACCTGCCATAAA‑CGCGACGCTTTTGCGCTGGGCGGCATGGCGACGTTTATTCCGAGCAAAGATGAAGAG‑CACAATCACCAGGTGCTACACACTGAGAAAGCGGATAAATCGCTGGAAGCACATAACGGTCACAATgg     >  2:265588/1‑251 (MQ=255)                                                                                                                                                                                                                                      tcaccTATATAAAAACGTTGAAAAACTATCCCGATCGCTTCCTGCCAGACGGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACGCACGCCTGTTGATTAAAACCTGCCATAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGCGTCGTTTATTCCGAGCAAAGATGAAGAG‑CACAATAACCAGGTGCTCAACAAAGTAAAAGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGcac  <  1:416079/246‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  AGGGCAGTGGTCCCTATTTCTATCT‑GCCGAAAACCCAGTCCTGGCAGGAAGCGGCCTGGTGGAGCGAAGTCTTCAGCTATGCAGAAGATCGCTTTAATCTGCCGCGCGGCACCATCAAGGCGACGTTGCTGATTGAAACGCTGCCCGCCGTGTTCCAGATGGATGAAATCCTTCACGCGCTGCGTGACCATATTGTTGGTCTGAACTGCGGTCGTTGGGATTACATCTTCAGCTATATCAAAACGTTGAAAAACTATCCCGATCGCGTCCTGCCAGACAGACAGGCAGTGACGATGGATAAACCATTCCTGAATGCTTACTCACGCCTGTTGATTAAAACCTGCCATAAA‑CGCGGTGCTTTTGCGATGGGCGGCATGGCGGCGTTTATTCCGAGCAAAGATGAAGAG‑CACAATAACCAGGTGCTCAACAAAGTAAAAGCGGATAAATCGCTGGAAGCCAATAACGGTCACGATGGCAC  >  NC_000913/4216091‑4216568

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