Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,335,568	G→A	19.4%	L526L (CTG→TTG)	fumB ←	anaerobic class I fumarate hydratase (fumarase B)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,335,568	0	G	A	19.4%	37.8 / 4.4	21	L526L (CTG→TTG)	fumB	anaerobic class I fumarate hydratase (fumarase B)
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base A (2/2);  total (10/11)
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.75e-01

GTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGAT  >  CP009273/4335323‑4335798                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       gtgaGCGTCAGTGCGTTGGTGTTAAGAATGAACGGCTTGGGTGTTCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTTACTATTCACCAGCAAGCGTTAAATTACAATCGCAAAAGAGGTTGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCg                                                                                                                                                                                                            <  1:229854/271‑1 (MQ=255)                                aCGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGATTCCATACCAAACTCCg                                                                                                                                                                             >  1:56844/1‑274 (MQ=255)                                                      ggTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGAGCGGGGCCGAAGAGGTTACTTAGTGCAGTGCGAGCACTGATTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTGGTCATCGACCAGGATAAACGACGGGAAAACTTCTAATCTGATTTTCCCAATAGCTTCCATACCCAGCCACGTATACACGACCTACTCCAAAtta                                                                                                                                                      >  2:401432/1‑273 (MQ=255)                                                                             aTACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                                                                                                                                                                                                   >  1:188672/1‑207 (MQ=255)                                                                                 tGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAAGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                                                                                                                                               <  2:188665/207‑1 (MQ=255)                                                                                             aGCAAGCGTGAAAGTACAATCGCAAAAGAGGTAGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGAccgccg                                                                                                               <  1:54393/275‑1 (MQ=255)                                                                                                                    caaaGGGGTGTCTAGGTGGTTCACCTGGCAGCATTCTGCTAGGCGCGGGGCCGAAGAGGTTACTTAGTGGAGTTAGCTCACTGTTTGTTGANGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCTGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAgccgcc                                                                                        <  1:401447/274‑1 (MQ=38)                                                                                                                                        tCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGc                                                                                                                                                                                              >  2:270484/1‑153 (MQ=255)                                                                                                                                            cTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTcc                                                                                                                                                                                          <  1:270495/153‑1 (MQ=255)                                                                                                                                                   caTGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATa                                                                                                                                                                                                                                    >  2:30093/1‑104 (MQ=255)                                                                                                                                                   caTGCTGCCAGGCGATGGGCCGAAGAGGTTACTTAGTGCCGTTCGCGCACTGGTAGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGATAGAAGCAGCCGTGTTTATGACAGGCGCCGGTAACccgc                                                            >  1:108396/1‑272 (MQ=255)                                                                                                                                                       ctgcCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACg                                                                                                                                                                                                                                <  1:30093/104‑1 (MQ=255)                                                                                                                                                         gcCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTAcc                                                                                                       <  1:216642/223‑1 (MQ=255)                                                                                                                                                                    ggCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGATTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGATTGATGcg                                                                                                                                              >  2:6857/1‑172 (MQ=255)                                                                                                                                                                        gAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTctgc                                                                                                                                          <  1:6857/173‑1 (MQ=255)                                                                                                                                                                         aaGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGCCGGTAAACTGCTGACTGCGGTTACCTTTCGCCAGc                                    >  1:78168/1‑274 (MQ=255)                                                                                                                                                                                            ttCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAg                                                                                                                                                                                                   >  2:135344/1‑96 (MQ=255)                                                                                                                                                                                                cgcACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAg                                                                                                                                                                                               <  1:135346/96‑1 (MQ=255)                                                                                                                                                                                                           gtgccGATTCNCTGGAAGAAGTCGTTACCGTTGTCACCGACCAAGTTAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGAt  <  2:56844/270‑1 (MQ=255)                                                                                                                                                                                                                                          tGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGa                                                                              >  1:238081/1‑167 (MQ=255)                                                                                                                                                                                                                                              aTCGACCATGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACAcg                                                                          <  2:238072/167‑1 (MQ=255)                                                                                                                                                                                                                                                      ggATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCTGATAAGc                                                                                                                                                                      >  1:28580/1‑67 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                       GTTAGCGTCAGTGCGCTGGTGTTAAGAATTAACGGCTTGAGTGTGCGGGTCATGGTGAGTTCGCAAATCAGGGAAAATACTGACTATTCACCAGCAAGCGTGAAATTACAATCGCAAAAGAGGTGGCCAGGGGGATCACCTGGCAGCATGCTGCCAGGCGCTGGGCCGAAGAGGTTACTTAGTGCAGTTCGCGCACTGTTTGTTGACGATTTGCTGGAAGAAGTCGTTACCTTTGTCATCGACCAGGATAAACGCCGGGAAATCTTCTACTTCGATTTTCCAGATAGCTTCCATACCCAGCTCCGGATAAGCGACGCACTCCAGATGCTTGATGCTCTGCTGCGCCAGTACCGCCGCCGGACCGCCGATGCTACCGAGGTAGAAGCCGCCGTGTTTATGACACGCGTCGGTAACCTGCTGACTGCGGTTACCTTTCGCCAGCATGATCATGCTGCCGCCGTGGGATTGCAGCAGAT  >  CP009273/4335323‑4335798

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