Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,579,023	C→A	14.1%	Q120K (CAA→AAA)	talA →	transaldolase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,023	0	C	A	14.1%	56.9 / 4.5	29	Q120K (CAA→AAA)	talA	transaldolase A
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base A (2/2);  total (12/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.04e-01

ACAAACTGGCGGTCAATTTCGGTGCTGAAATCCTCAAAATCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACA  >  NC_000913/2578886‑2579150                                                                                                                                          |                                                                                                                                acgaaCGGGCGGTAGATTTCGGTGCTGATAGCCTCAAGATCGTGCCCGGGCGCGTGTCAACAGAAGTGGATGCACGCCTCTCTGTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGcaacaa                                                                                                                            <  2:464473/140‑1 (MQ=255)          cGGTCAATTTCGGTGCTGAAATCCTCAAAATCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTTAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTgag                                                                                                                   <  1:595994/143‑1 (MQ=255)                        gCTGAAATCCTCAAAATCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCt                                                                                                     >  1:336372/1‑143 (MQ=255)                                  tCAAAATCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGCAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTg                                                                                           >  1:402/1‑143 (MQ=255)                                      aaTCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCtt                                                                                       <  2:319721/143‑1 (MQ=255)                                       ttCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTc                                                                                      <  1:501163/142‑1 (MQ=255)                                             cccGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAACAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTg                                                                                <  1:452995/143‑1 (MQ=255)                                                 gTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCGTCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGTCTTCGACCTGggaa                                                                            <  2:159051/143‑1 (MQ=255)                                                  tCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCTAGCTTGCTTCGACATAggaat                                                                           >  1:268204/1‑142 (MQ=255)                                                                  gTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATAAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCaga                                                           >  2:618246/1‑143 (MQ=255)                                                                   ttGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCggaa                                                          >  1:112186/1‑143 (MQ=255)                                                                       tGCACGGCTCTCTTTTGAGAAAGAAAGGAGGATTGAGAGGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGc                                                      <  2:183267/143‑1 (MQ=255)                                                                                ctctTTTGATAAAGAAAAGAGGATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAaag                                             <  1:74145/143‑1 (MQ=255)                                                                                  ctttgcaTAAAGAAAAGAGTATGGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAaagaa                                           <  1:264155/137‑1 (MQ=255)                                                                                  ctTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAaagaa                                           <  1:474350/143‑1 (MQ=255)                                                                                     ttGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGt                                        >  1:515093/1‑143 (MQ=255)                                                                                        aTAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAGGAGCTGGAAAAAGAAGGTAtt                                     >  2:584463/1‑143 (MQ=255)                                                                                        aTAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTAtt                                     >  2:167685/1‑143 (MQ=255)                                                                                         tAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGGGGACTTGTATCAGCAACAAGGCGTGGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTa                                    <  2:156897/143‑1 (MQ=255)                                                                                          aaagaaaagaGCATTGAGAAGGCGCGCCATCTGGTGGACTTGTCTCAGTAGCAATGCGTTGCGAAATCACGCATGCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTaa                                   <  1:425012/143‑1 (MQ=255)                                                                                                aagaGTATTGAGAAGGCTTGCCATCTGGTGGACTTCCACCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCaa                             <  1:314415/143‑1 (MQ=255)                                                                                                 cgagTATTGAGAAGGCGAGCCATCTGGTGGACTTGTATCAGCGACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAAc                            <  2:355385/142‑1 (MQ=255)                                                                                                 agagTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAAc                            >  1:591748/1‑143 (MQ=255)                                                                                                                   gcCGTCTGGTGGAATCGGATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGAGCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGGTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCt          <  1:506074/143‑1 (MQ=255)                                                                                                                           gtggACTTGTATCAGCAACAAGGCGTTGAGACATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGcaca  <  2:285770/143‑1 (MQ=255)                                                                                                                           gtggACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGcaca  >  1:12997/1‑143 (MQ=255)                                                                                                                             ggACTTGTATCAGAAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGctgctg                >  1:269721/1‑127 (MQ=255)                                                                                                                             ggACTTGTATCAGAAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGctgctg                >  1:605446/1‑127 (MQ=255)                                                                                                                             ggACTTGTATCAGAAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGctgctg                <  2:269721/127‑1 (MQ=255)                                                                                                                             ggACTTGTATCAGAAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGctgctg                <  2:605446/127‑1 (MQ=255)                                                                                                                                          |                                                                                                                                ACAAACTGGCGGTCAATTTCGGTGCTGAAATCCTCAAAATCGTACCCGGTCGCGTGTCAACAGAAGTTGATGCACGCCTCTCTTTTGATAAAGAAAAGAGTATTGAGAAGGCGCGCCATCTGGTGGACTTGTATCAGCAACAAGGCGTTGAGAAATCACGCATTCTGATCAAGCTGGCTTCGACCTGGGAAGGAATTCGCGCGGCAGAAGAGCTGGAAAAAGAAGGTATTAACTGCAACCTGACGCTGCTGTTTTCTTTTGCACA  >  NC_000913/2578886‑2579150

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