Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,577,987	C→A	G131C (GGC→TGC)	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,987	0	C	A	92.0%	65.0 / ‑5.5	25	G131C (GGC→TGC)	maeB	malate dehydrogenase (oxaloacetate‑decarboxylating) (NADP(+))
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (12/11);  total (13/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 = 9.39e-01

GACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGT  >  NC_000913/2577747‑2578220                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             gACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGATGTTGATGCCGCAGAAggttggt                                                                                                                                                                                                                                    <  1:368074/251‑1 (MQ=255)            aCCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCTCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGcgcg                                                                                                                                                                                                                           <  2:387502/249‑1 (MQ=255)                  agGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAAc                                                                                                                                                                                                                   <  1:393133/251‑1 (MQ=255)                      tCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTc                                                                                                                                                                                                                <  1:340083/250‑1 (MQ=255)                           cAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATaaa                                                                                                                                                                                                          <  1:259935/251‑1 (MQ=255)                                   tgctgcGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTcc                                                                                                                                                                                                   <  1:369774/250‑1 (MQ=255)                                    gcggcGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTGTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGAGCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGGTCGAGCGCGGCGACAACTTCAATAAAGTTGTCCg                                                                                                                                                                                                  <  1:57587/247‑1 (MQ=255)                                                    cGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACt                                                                                                                                                                                 >  1:252979/1‑251 (MQ=255)                                                              caTCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCaa                                                                                                                                                                       <  2:236306/251‑1 (MQ=255)                                                                    cACGTCGGAGATGTTTTTCTCCACCACAAGCCAGCCGGTGAGGATGGCGGCAGCGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCACTATAGAAACATTCTGGCGCTTTATTGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCACGGCGTCAACTTCAATACATTTGTCCGGGTCGCGTTCGTAACCTTCAATGTCCAAT‑aact                                                                                                                                                                >  2:443523/1‑248 (MQ=255)                                                                                       tcCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGAGCTGATCGTAGCGGAATACCGGAATATTCATCCGCCCGCGCAGTTTCTGTTCAATATAGAACCATTCTGGCGCTTTAATGTCCTCGAGGTTGATGCCGCAGAAGGTTGGTTCGCGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGCGTTCGTCAACTTCCATGTCAACT‑ACATCAATCCCGGCGAATTTCtt                                                                                                                                             >  2:390802/1‑251 (MQ=255)                                                                                           ccacgatcaatccGTTGAGGATGGTGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGCATACCGGAATATTCATCCGCTCGCCCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGACGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATCAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGACAAAT‑ACATCAATCCTGGCGAATTTCTTAAAc                                                                                                                                         >  2:202511/12‑251 (MQ=255)                                                                                                        gTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGc                                                                                                                            >  1:229579/1‑251 (MQ=255)                                                                                                                 ggcggcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGTGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAC‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCa                                                                                                                   <  2:324182/251‑1 (MQ=255)                                                                                                                     gcAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCgg                                                                                                               >  2:171347/1‑251 (MQ=255)                                                                                                                                 aTTGCAGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATAAGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTt‑gcctgcc                                                                                                    >  1:300172/1‑249 (MQ=255)                                                                                                                                       gtgccgtgcTGATCGTCGTGGAATTCCGGAATATTCATCCGCTTGCGCAGTTTCTGTTCAATATAGAAACATGCTGGCGCTTTAATGGCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGATTTCGTCAACGGCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCGGCCAgcgcgca                                                                                            <  2:66154/251‑2 (MQ=255)                                                                                                                                                     gtcgtGGAATACCGGAATATTCCTCCGCTCGCGCAGTTTCTGTTCAATATAGAACCATTATGGCGCTTTAATGTCTTCGCGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAAcc                                                                              >  1:164751/1‑251 (MQ=255)                                                                                                                                                                 cGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAACCCC‑AGCAccgccg                                                                 >  1:342137/1‑251 (MQ=255)                                                                                                                                                                  ggAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAACCCC‑AGCACCGCCGt                                                                >  1:393032/1‑251 (MQ=255)                                                                                                                                                                            tCCGCTCGCGCAGGTTCTGTTCAATATAGAAACCTTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATAAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAACCCC‑AGCACAGCCGTACCGTTag                                                        >  2:9544/1‑249 (MQ=255)                                                                                                                                                                                  cgcgTAGTTTCTGTTAAAAAAAAAAACTTTCGGGCTCTTTATGGCCTTCGGGGTTGAGGCCGCAGAGGGTTGGTTCGGGGGCGGCAACAACTTCAAAAATTTTGCCCGGGCCGAGTTCGTCAATTACAATGTCAAAT‑ACTTCAATCCCGGCGAATTTCTTAAACAAAACGCCCTTGCCTTCCACCCCCGTTTT‑TCCTGCAAGCGCGCCAGGGTTGCCTACCCCC‑AGCACCGCCGTACCGTTAGAGATCAcc                                                <  1:249541/251‑1 (MQ=255)                                                                                                                                                                                        gTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGCGGATGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCAAGTACGTCACCTTCTATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTACCTTCCATCACCGTTTT‑GCCTGCCAGCTCGCCAATGTTGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGccacc                                          >  1:71143/1‑251 (MQ=255)                                                                                                                                                                                                                cTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCAGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg                  >  1:267910/1‑251 (MQ=255)                                                                                                                                                                                                                       tttAATGTCTTCGGGGTTGAGGCCGAAGAAGTTTGTTGCGAGCGGGGCGACAACTTAAATAAATTTGCCGGGGTCGAGTTCGTCAACTCAAATGTAAAAT‑ACATCAACCCCGGGAAATTTCTAAAAAAGAACCCCCTTGCCTTCCACCCCGGGTTT‑GCCTGCCAGCGCGCAAATGTTGCCAACCCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt            <  1:171347/250‑1 (MQ=255)                                                                                                                                                                                                                                 tgggggtggatcccgcagAAGGTTGGTCCGAGCGGGGCAACATCTTAAAAAAATTTGCCGGGGTCGATTTCGCAAACTCAAATG‑CAAATTACATCAACCCCGGCAATTTTCTTAAACAGACCGCCTTTCCCTCCCA‑CCCCGGTTTTGCTGGCCAGCGCCCAAATGT‑GCCCAACCCCCAGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTGACGGGt  <  1:126652/238‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             GACCCAGCGCTACCAGCAGGTTCATACAGGCGATTGCTGCGGCACCCGCGCCGGAAACCACCATCCGCACGTCGGAGATGTTTTTCTCCACCACGCGCAAGCCGTTGAGGATGGCGGCAGTGCTGATAATTGCCGTGCCGTGCTGATCGTCGTGGAATACCGGAATATTCATCCGCTCGCGCAGTTTCTGTTCAATATAGAAACATTCTGGCGCTTTAATGTCTTCGAGGTTGATGCCGCCGAAGGTTGGTTCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAAT‑ACATCAATCCCGGCGAATTTCTTAAACAGAACGCCCTTGCCTTCCATCACCGGTTT‑GCCTGCCAGCGCGCCAATGTTGCCTAACCCC‑AGCACCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGT  >  NC_000913/2577747‑2578220

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