Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,337,394	A→G	S522G (AGC→GGC)	acnA →	aconitate hydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,337,394	0	A	G	96.5%	71.5 / ‑5.7	30	S522G (AGC→GGC)	acnA	aconitate hydratase 1
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (11/17);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

GATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCT  >  NC_000913/1337162‑1337634                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 gATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTcggtg                                                                                                                                                                                                                                >  1:423894/1‑251 (MQ=255)             cTTGCTGCCGAACAAACCCGTAACTCTGGGCCTCATGCGGCAACCATGGGTAAAAGCGTCGCTGCAACCGGGTCCGCCAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGCTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTCACTCTTGGCCGCTGCCCGATCCTATCGCACCGGCACTCAAAAAAGGCGCTTTAACCGTCGGAGCGGTGCTGTCCg                                                                                                                                                                                                                    >  2:206215/1‑251 (MQ=255)               tGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCgg                                                                                                                                                                                                                   <  1:470174/250‑1 (MQ=255)                gCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTCTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCCCCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGc                                                                                                                                                                                                                  <  1:350460/250‑1 (MQ=255)                  tGGCTAAAAAAGCCGTAACTCAGGGACAGAATCGGAAACCCAGGGACAAAGCATCGATGGCACAGGGTTCGAAAGTCGGTACTGATTCTCTGGCAAAAGCGAAACTGACACCGTATCACGACGAACTGGGGTTTACAATTGTGGGATAAGATTGTACCACCTGTATTGTTAACTCTGGGCAGCTGCCCACGCCTATCGAATCAGAAATCAACAAAGTTGATTTAACCGTcgcagcggagcggtcaaggaa                                                                                                                                                                                                                >  2:338170/1‑231 (MQ=255)                        aaaaaGCGGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGGCCGGCAACCGTAAc                                                                                                                                                                                                         <  2:445104/251‑1 (MQ=255)                            agcctgaaCTCTGGTCCTCAAGCGGCAGCCATGGGGCAAGGTGTCGCTGGCACCGGGTTCTAAAGTCGGGTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCGCTACGAACTGGGGTTTATCCTTGTGGGATACAGTTGTACCACCTGTATGGGTAACTCTGGGCCGCTGCCCGATTCTATCTTTACGGCAAGCAAAAAAGGTGATTTAACTGTCTGTGCTGTGGTGTCCTGCAACCGTATCTTTg                                                                                                                                                                                                     <  2:530388/245‑1 (MQ=255)                                  aaCTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGcc                                                                                                                                                                                               <  1:348968/251‑1 (MQ=255)                                   aCTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTGTGAAGGCCg                                                                                                                                                                                              <  2:211652/251‑1 (MQ=255)                                        gggCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatcc                                                                                                                                                                                         <  2:141012/251‑1 (MQ=255)                                          gCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTatccat                                                                                                                                                                                       <  2:449454/251‑1 (MQ=255)                                               aaGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCt                                                                                                                                                                                  >  1:117439/1‑251 (MQ=255)                                                agcagCAACCATGAGTCAAAGCGTCGCTGGCACCGGGTACGAAAGTCGTTTCTGATTATATGGCAAAAGCGAAACTGACACCGTATCTCGACGATCTGGGGATTAACCTTGTGAGATACGATTGTACACCCTGTAGTGGTAACTCTGAGCCGCTGCCCGATCCTATAGAAACGGCAATCACAAAACGCGATTTAACCGTCGGTGCGGTGCTGTCCAGCAAACGACACTTTACAGGCAGTATCCATCCGc                                                                                                                                                                                   >  2:338101/5‑249 (MQ=255)                                                                  aaGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAACCTGACACCGTATCTCGGCGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCAGCCGCTGGTTAAAACTAActggctg                                                                                                                                                               <  2:342738/251‑1 (MQ=255)                                                                        cGCTGGCACCGGGTTCGTAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACGGGGGTTTAACCTTGTGGGATCCGGTTGTCCCCCCTGTATTGGTAACTCTGGGCCGCTACCCGATCCTAGCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCGTCCGCTGGTTAAAACTAACTGGATGGCCTc                                                                                                                                                          <  2:427471/250‑1 (MQ=255)                                                                          cgggctccGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACCCCGTATCGCGACGAACTGGTGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTTGCCTcgcc                                                                                                                                                       <  2:423894/245‑1 (MQ=255)                                                                                ccGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCt                                                                                                                                                 >  1:368293/1‑251 (MQ=255)                                                                                               gTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATgc                                                                                                                                  <  1:155962/251‑1 (MQ=255)                                                                                                     gATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGc                                                                                                                            >  2:146113/1‑251 (MQ=255)                                                                                                             ggCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGAGTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTTAAGGCCTTATCGATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGaa                                                                                                                       <  1:545146/248‑1 (MQ=255)                                                                                                                   aGCGAAACTGACACCGGATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCGCCTGTATTGGTAACTCTGGGCCGCGGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTTGCGGGAAATATGAata                                                                                                              <  2:117439/251‑1 (MQ=255)                                                                                                                     cGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGAGGGAAATATGAATATc                                                                                                            <  2:466479/251‑1 (MQ=255)                                                                                                                                                tGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCAt                                                                                 >  1:62913/1‑251 (MQ=255)                                                                                                                                                                  gATACGGTTGTACAGCCTGTATTAGTAACTCTGGGCCGCGGCCGGACCCTACCGAAACGGAAATCAAAAAAGGCGATTTAACCGGCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCg                                                               <  2:158618/251‑1 (MQ=255)                                                                                                                                                                   aTACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCgg                                                              >  2:291780/1‑251 (MQ=255)                                                                                                                                                                            tACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCCACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTg                                                      >  1:378381/1‑250 (MQ=255)                                                                                                                                                                                                  gtgCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAGGCTATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTTAAGGCCGTAGCCAGCCGCTGGTGAAAACTAACTGGCGGGCCTCGCCGCCGCTGTTGGTTGCCTATTCGCTGGCGGGAAATATGAATATCAACCGGGGTTCGTATCCTTTCGGCCATGATCTCAAAGGCTATCCGGTGTAGCTGAAAGATATTTAGCCTTCGGctc                                <  2:579082/248‑3 (MQ=255)                                                                                                                                                                                                         tGCCCGATCCTACCGACACGGCAACCAAAAAAGGCGATTTAACCGTAGGTGCGGTGCTGTCCGGCAACCGTACATTAGAAGGCCGAATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGCAATATGGATATCAACCTGGCTTCTGAGCCTCTCGGCCACGCTCGCAAAGGCGAACCGGTTTATTTGCAAGATATCttgacctcggccaaaagaatt                        >  2:109537/1‑231 (MQ=255)                                                                                                                                                                                                           cccGATCCTATAGAAACGGAAATCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAAAATGAATATCAAACTGGCTTCTGAGCCTATCGGCCATGATCGCAACGGCGATCCGGTTTATCTGAAAGATATATGGCACTCGGCAAAAGAAATTGc                      >  2:27505/1‑251 (MQ=255)                                                                                                                                                                                                                               aaTCAAAAAAGGCGATTTAACCGTCGGTGCGGTGCGGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTAGCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATGGCCCGTGCGGTAGAACAAGtct  <  2:547370/251‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                 GATGGCCGCAGGCTTGCTGGCGAAAAAAGCCGTAACTCTGGGCCTCAAGCGGCAACCATGGGTCAAAGCGTCGCTGGCACCGGGTTCGAAAGTCGTTTCTGATTATCTGGCAAAAGCGAAACTGACACCGTATCTCGACGAACTGGGGTTTAACCTTGTGGGATACGGTTGTACCACCTGTATTGGTAACTCTGGGCCGCTGCCCGATCCTATCGAAACGGCAATCAAAAAAAGCGATTTAACCGTCGGTGCGGTGCTGTCCGGCAACCGTAACTTTGAAGGCCGTATCCATCCGCTGGTTAAAACTAACTGGCTGGCCTCGCCGCCGCTGGTGGTTGCCTATGCGCTGGCGGGAAATATGAATATCAACCTGGCTTCTGAGCCTATCGGCCATGATCGCAAAGGCGATCCGGTTTATCTGAAAGATATCTGGCCATCGGCACAAGAAATTGCCCGTGCGGTAGAACAAGTCT  >  NC_000913/1337162‑1337634

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