Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,577,937	T→C	E147E (GAA→GAG)	maeB ←	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,577,937	0	T	C	82.8%	55.9 / 2.3	29	E147E (GAA→GAG)	maeB	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (14/10);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAG  >  NC_000913/2577693‑2578166                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                      ggCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAata                                                                                                                                                                                                                                 <  1:93337/251‑1 (MQ=255)     cTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAatat                                                                                                                                                                                                                                >  1:129811/1‑248 (MQ=255)           gataACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCCGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACAtt                                                                                                                                                                                                                       >  1:144659/1‑251 (MQ=255)                       tGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTa                                                                                                                                                                                                            >  1:180/1‑250 (MQ=255)                                aaCCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCg                                                                                                                                                                                                   >  1:272029/1‑250 (MQ=255)                                   cACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACAACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGt                                                                                                                                                                                               >  1:251446/1‑251 (MQ=255)                                                              cgcTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAgcg                                                                                                                                                                    >  1:64458/1‑251 (MQ=255)                                                                                         tgctgcGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCg                                                                                                                                         >  1:311202/1‑251 (MQ=255)                                                                                             gcgtacacccctccggacatcaccaACCGCACGTTGGAAATGTTTTTCCCCAACTCGCGCAAGCCGTTGACGAGGGAGGCAGTGATGATAATGGCCGTGCCGTGCTGATCGTTGTGGAATACATGAATTTTCATCCGTTTGCGCAGTTTCTGCTATATATAGAAACTTTCTGGCGCTTTAATGTCTTCGAGGTGGATTCCGCCGAAGGTTGGTTCGAGCGCGGCGACATCGTCTCTAAATTTGTCCgtggc                                                                                                                                     <  2:78347/231‑5 (MQ=255)                                                                                              cGGCCCCCGCGCCGGAAACCACCATCCGCAAGCCGGAGATGTTTTTCACCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGGCTTCGAGGTTGATGCCGCCGAAGGTGGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCTGGGTCg                                                                                                                                    <  2:64458/251‑1 (MQ=255)                                                                                                    ccggcccGGAAACCCCCATACGCACGGCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATCCCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGACACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGATCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCg                                                                                                                              <  2:91188/246‑1 (MQ=255)                                                                                                          cGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATTATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGACTATTCATCCGCTCGCGCAGTTTCTGCTCCATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATCAATTTGTCCGGGTCGAGGTCGTCAACt                                                                                                                        >  2:257090/1‑251 (MQ=255)                                                                                                           ggAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACtt                                                                                                                       >  1:129065/1‑251 (MQ=255)                                                                                                                       ccGCACGACGGAGATGTCTTTCTCCACCACGCGGAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTTCCGTGCTGATCGTCGTGGAATTCCGGAATATTCATCCGCTCGCGCATTTTCTGCCCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCTCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATa                                                                                                           <  1:288223/251‑1 (MQ=255)                                                                                                                            cGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCa                                                                                                      <  1:312326/251‑1 (MQ=255)                                                                                                                                        ttttCTCCACCTCGCACAAGCCGTTTTGGATGGCGTCAGTGCTGATAAGTGCCGTGTCGTGCTGATCGGCGTTGCATACCTGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATAGAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATTCCGTCTAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGaa                                                                                           <  2:372178/250‑1 (MQ=255)                                                                                                                                             tcCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCAAGTTCGTCa                                                                                                                           <  1:206278/213‑1 (MQ=255)                                                                                                                                             tcCACCACGCGAAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGCATATTAATCCGCTCGCGCAGATACTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGGTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGGCAAGTTCGTCa                                                                                                                           >  2:206278/1‑213 (MQ=255)                                                                                                                                                   acGCGCAAGCCGTTGAGGATGGCGGCAGCGATGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACa                                                                                >  1:94307/1‑250 (MQ=255)                                                                                                                                                                     aTGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCAt                                                             <  1:248465/251‑1 (MQ=255)                                                                                                                                                                       ggcggcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATc                                                            >  1:61544/1‑250 (MQ=255)                                                                                                                                                                              gTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGttt                                                    <  1:313884/251‑1 (MQ=255)                                                                                                                                                                              gTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGttt                                                    >  1:287574/1‑251 (MQ=255)                                                                                                                                                                                                    gtagATCTTCGGGGAATACCGGAATATTCATCCGCGCGCGCTGTTTCGGTTCATTATAGAAACATGCTTGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACATCTGCAATAAATTTGTCCGGGTCGAGTTCTTCAACTTCTATGTCAAATACATCAAGCCCGGCGTATTCCTTAAACAGAACGCCCTTGCCTTCCATCCCCGGTTTGCCTGCCAGCGCGCCATTGtt                               <  2:272029/247‑1 (MQ=255)                                                                                                                                                                                                      tGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGcc                            >  1:238206/1‑251 (MQ=255)                                                                                                                                                                                                                     aCCGGAATATTCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATGCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCAccg              <  1:83623/250‑1 (MQ=255)                                                                                                                                                                                                                               tCATCCGCTCGCGCAGTTTCTGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGtt    >  1:115678/1‑250 (MQ=255)                                                                                                                                                                                                                                cATCCGCTCGCGCAGTTTATGCTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCAAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCCCCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCCGCACCGCCGTCCCGTTa   >  2:128039/1‑250 (MQ=255)                                                                                                                                                                                                                                   ccGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGGATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTag  >  1:77659/1‑248 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                      GGCCCTGATAGATAACGCCTTTTGAATCGCAAACCACGATGTTATGTTTTTGCAGACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAACCCCAGCACCGCCGTACCGTTAG  >  NC_000913/2577693‑2578166

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