Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,757,006	G→A	A437T (GCT→ACT)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,757,006	0	G	A	100.0%	121.3 / NA	39	A437T (GCT→ACT)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/20);  total (19/20)

GCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756877‑1757130                                                                                                                                  |                                                                                                                             gCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCaga                                                                                                                     <  2:532735/139‑1 (MQ=255)  ccaccaTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    <  1:115904/139‑1 (MQ=255)       aTCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTc                                                                                                               <  2:53345/139‑1 (MQ=255)                   accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCTGTGTGCCGCAGCTTGTTAAAGAGATCACGTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  2:727601/1‑139 (MQ=255)                   accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    <  2:580551/122‑1 (MQ=255)                   accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCagag                                                                                                                    >  1:580551/1‑122 (MQ=255)                   accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAg                                                                                                   >  1:893329/1‑139 (MQ=255)                   accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGC‑‑AGGCGCTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGca                                                                                                          <  2:650468/130‑1 (MQ=255)                   accaccAACGAAAAAACGGCTCATCAGTTGGTACTGAGC‑‑AGGCGCTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGca                                                                                                          >  1:650468/1‑130 (MQ=255)                        caACTAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACATCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGAGGTCTGGCACACAAAGGTGAc                                                                                              >  2:752920/1‑139 (MQ=255)                          aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             <  1:348196/122‑1 (MQ=255)                          aCGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTg                                                                                                             >  2:348196/1‑122 (MQ=255)                                                ggTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCa                                                                      <  1:752920/139‑1 (MQ=255)                                                       aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTg                                                                                          >  1:398335/1‑112 (MQ=255)                                                       aGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTg                                                                                          <  2:398335/112‑1 (MQ=255)                                                                                ttAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATc                                      <  1:280289/139‑1 (MQ=255)                                                                                      agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  2:1081891/97‑1 (MQ=255)                                                                                      agaTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  1:1081891/1‑97 (MQ=255)                                                                                        aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  2:1208434/1‑139 (MQ=255)                                                                                        aTCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAAGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  2:457281/1‑139 (MQ=255)                                                                                         tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           <  1:1280230/109‑1 (MQ=255)                                                                                         tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAg                                                           >  2:1280230/1‑109 (MQ=255)                                                                                         tCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGt                             >  2:139757/1‑139 (MQ=255)                                                                                         acACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTATAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAAGACCGCATCTGTTCACGt                             >  1:797041/2‑139 (MQ=255)                                                                                          cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCATTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTc                            <  2:797041/139‑1 (MQ=255)                                                                                          cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTc                            <  1:727601/139‑1 (MQ=255)                                                                                          cACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTc                            >  1:334443/1‑139 (MQ=255)                                                                                             ttCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTg                         <  2:893329/139‑1 (MQ=255)                                                                                                     atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGt                                                                                                 <  2:130149/59‑1 (MQ=255)                                                                                                     atgatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGt                                                                                                 >  1:130149/1‑59 (MQ=255)                                                                                                       gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCg                                                                                                                  <  1:983749/40‑1 (MQ=255)                                                                                                       gatTTCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCg                                                                                                                  >  2:983749/1‑40 (MQ=255)                                                                                                           tCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              <  2:381737/120‑1 (MQ=255)                                                                                                           tCTACCGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACg                              >  1:381737/1‑120 (MQ=255)                                                                                                               ccGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCa                                >  2:19116/1‑114 (MQ=255)                                                                                                               ccGTCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCa                                <  1:19116/114‑1 (MQ=255)                                                                                                                  tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          <  1:524216/69‑1 (MQ=255)                                                                                                                  tCTGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTgg                                                                          >  2:524216/1‑69 (MQ=255)                                                                                                                    tGGGTAAAGAACTGACTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGaattaat  <  1:139757/139‑1 (MQ=255)                                                                                                                                  |                                                                                                                             GCCACCATCCTGGCACTGACCACCAACGAAAAAACGGCTCATCAGTTGGTACTGAGCAAAGGCGTTGTGCCGCAGCTTGTTAAAGAGATCACTTCTACTGATGATTTCTACCGTCTGGGTAAAGAACTGGCTCTGCAGAGCGGTCTGGCACACAAAGGTGACGTTGTAGTTATGGTTTCTGGTGCACTGGTACCGAGCGGCACTACTAACACCGCATCTGTTCACGTCCTGTAATATTGCTTTTGTGAATTAAT  >  NC_000913/1756877‑1757130

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