Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,756,374	G→A	G226D (GGC→GAC)	pykF →	pyruvate kinase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,756,374	0	G	A	98.2%	154.2 / NA	54	G226D (GGC→GAC)	pykF	pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (30/23);  total (30/24)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAA  >  NC_000913/1756235‑1756513                                                                                                                                            |                                                                                                                                            tttGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTAAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCt                                                                                                                                          >  2:238725/1‑143 (MQ=255)            aCAAGGCGTAGACGTTGTTGCTGCTTCCGTGATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAGAATCGAAAACCAGGAAGACCTCAACAACTTcg                                                                                                                               <  2:487275/143‑1 (MQ=255)                 gCGTAGACTTTGTTGCTGCTTCCTGTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGaa                                                                                                                          <  2:338130/143‑1 (MQ=255)                  cGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGaaa                                                                                                                         >  2:638680/1‑143 (MQ=255)                       aCTTGGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTc                                                                                                                    <  2:845018/143‑1 (MQ=255)                          tttttgCTGCTTCCTTTCTTCGTAAGCGCTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGCAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGaa                                                                                                                 <  2:382291/140‑1 (MQ=255)                          ttgttgCTGCTTCCTTTATTTGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGaa                                                                                                                 <  1:736527/143‑1 (MQ=255)                               gctgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGGGAAAACATCGACATCATCTCCAAAATCGAAACCCAGGAAGACCTCAACAAATGCGACGAAATCCTCGAAGCctc                                                                                                            >  2:502520/1‑143 (MQ=255)                                ctgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCGCATCATCTCCAAAATCGAAAACCAGGAAGACCTAGAGAACTGCGACGAAGTCCCAGAAGccacc                                                                                                           >  2:718735/1‑140 (MQ=255)                                 tgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAACCATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTACGACGCAATCCTCGAAGCCTCTg                                                                                                          >  1:591092/1‑143 (MQ=255)                                 tgctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTGCGACGAAATCCTCGGCGCCTCTg                                                                                                          >  1:138169/1‑143 (MQ=255)                                  gctTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGa                                                                                                         <  1:28251/143‑1 (MQ=255)                                              cGTAAGCGGGCTGACGCTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAACACCATGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGt                                                                                             <  2:206708/143‑1 (MQ=255)                                               gTAAGCGTTCTGACGTTATCGAAATACGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGtt                                                                                            <  1:153964/143‑1 (MQ=255)                                                  aGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATCGTTgcg                                                                                         >  2:542926/1‑143 (MQ=255)                                                   gCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACACCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTgcgc                                                                                        >  1:192708/1‑143 (MQ=255)                                                     gTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGt                                                                                      <  1:253599/143‑1 (MQ=255)                                                       tCTGACGTTATCGAAAGCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTgg                                                                                    <  2:428605/143‑1 (MQ=255)                                                         tGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCATCTGACGGCGTCATGCCTGCGAGtatcc                                                                                  >  1:921402/1‑139 (MQ=255)                                                             gTTATCGAAATCCGTGAGCAACAGAAAGCGCACGGCAGCGAAAACATACAAATCAACCACAAAATCGAAAACCAGGAAGACCTCAACAACTCAGACGACATCATACAAGCCTCTGACGGCAACATAGTAGCACGTgaacacgt                                                                              >  2:397597/1‑136 (MQ=255)                                                             gTTATCGAAAGCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCt                                                                              <  2:187540/143‑1 (MQ=255)                                                              ttATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCta                                                                             >  1:484776/1‑142 (MQ=255)                                                                 tCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGgt                                                                          <  1:501658/143‑1 (MQ=255)                                                                    aaaTCCGTGCGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCCACAACTTCGACGAAATCCTCGCAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTa                                                                       >  1:9375/1‑143 (MQ=255)                                                                    aaaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTa                                                                       >  1:443262/1‑143 (MQ=255)                                                                     aaTCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAg                                                                      <  1:418662/143‑1 (MQ=255)                                                                       tCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATGGAAAACCAGGAAGACCTCAACAACTTCGCCGAAAGCCTCGAAGCCTCTGACGGCCCCACGGCTGCGCGTGGCGCCCTGCGTCTTGaa                                                                    >  1:155561/1‑143 (MQ=255)                                                                               cACCTGAAAGCGCACGGCGGCGAACACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCCCAACAACCGCGACGCAATCCTCGCAGCCTCTGACGGCATCACGCTTGCGCGTGGCCACCCGGGCGCAGGACTCCCGGt                                                            >  2:393265/1‑143 (MQ=255)                                                                                          gcACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGCTTAGAAATCCCGGTAGAAGCAGttc                                                 >  1:823121/1‑142 (MQ=255)                                                                                          gcACGGCGGCGAAAACATCCACATCAGCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTa                                                 <  2:922769/143‑1 (MQ=255)                                                                                            aCGGCGGCGAAAACCTCCACACCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATc                                               <  2:363743/143‑1 (MQ=255)                                                                                                cggcgAAAACATCCACATCAGCTCCAAACTCGAAAACCAGGAAGACCTCAACCACTTCGACGAAATCCTCGAAGCCTCGGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCg                                           <  2:402230/143‑1 (MQ=255)                                                                                                cggccAACACATCCACACCCTCTCCAAACTCTAAAACCAGGAACACCTCAGCAACTTCGACGAAATCCTCTCAGCCTCTGACGCCATCATGGTTGCGCGTGTCGACCCGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCg                                           <  2:389192/143‑1 (MQ=255)                                                                                                    gAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAACTCCCGGTAGAAGAAGTTATGTTGGCCCa                                       >  2:410841/1‑143 (MQ=255)                                                                                                      aaaCATCCACATCATCTCCAAAATCGAAAACCCAGAACACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAACTCCCGGTAGAAGAAGTTATCACCGACCaga                                     >  2:200655/1‑143 (MQ=255)                                                                                                         cATCCACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGCGTAGAAATCCCGCTCGAAGAAGTTATCTTCGCCCaggaga                                  >  1:323996/1‑143 (MQ=255)                                                                                                            ccACATCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatga                               <  1:238725/143‑1 (MQ=255)                                                                                                              acaTCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGAGGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAGGTTATCTTCGCCCGGAAgaggatg                             >  1:501151/1‑143 (MQ=255)                                                                                                              acaTCATCTCCAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatg                             >  1:300051/1‑143 (MQ=255)                                                                                                                agcatcTCCAAAATCGAAAACCAGGCAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGGAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAgatgatgat                           <  2:871057/141‑1 (MQ=255)                                                                                                                        cagaaTCGAAAACCAGGAAGACCTCAACAACTCCGACGAAATCCGCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt                   <  1:45518/140‑1 (MQ=255)                                                                                                                        cAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCGCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt                   >  1:74225/1‑143 (MQ=255)                                                                                                                        cAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGCTCGAACAAt                   >  2:153849/1‑143 (MQ=255)                                                                                                                        cAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt                   >  1:558749/1‑143 (MQ=255)                                                                                                                        cAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATATTCGCCCAGAAGATGATGATCGAAAAAt                   >  1:934946/1‑143 (MQ=255)                                                                                                                        cAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGGAGAAGAAGTTATCTTCGCCCAGAAGATGATGCTCGAACAAt                   >  2:852002/1‑143 (MQ=255)                                                                                                                        cAAAATCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATACTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGGGTAGCAATCCCGGTAGCAGAAGGTATCTTCGCCCGTAAGATGATGATCGaaaaag                   >  1:30154/1‑142 (MQ=255)                                                                                                                        cAAAACCGAAAACCAGGAAGATCTCAACAACTTCGACGCAATGCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAAt                   <  1:161585/143‑1 (MQ=255)                                                                                                                             tCGAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATc              >  1:37910/1‑143 (MQ=255)                                                                                                                               gAAAACCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGCTGATCGAAAAATGTATCCg            >  1:706734/1‑143 (MQ=255)                                                                                                                                  aaCCGGGAGGGCCGCAACAGCGTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGc         <  2:481993/143‑1 (MQ=255)                                                                                                                                  aaCCAGGAAGACCTCAACAACTTCGACGAAACCCTCGAAGCCTCTGACGGCAACATGGTTGCGCGTGACGACGTGGGTGTAGCAATCCCGGTAGAAGAGGTTATCTTCGGCCAGAGGACGCGGACCGCAAAAAGTAACCGTGc         >  2:274840/1‑143 (MQ=255)                                                                                                                                   aCCAGGAAGACCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCa        <  1:812075/143‑1 (MQ=255)                                                                                                                                         aaGACCTCAACAACTGCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTaaa  <  2:251770/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            TTTGGTTGCGAACAAGGCGTAGACTTTGTTGCTGCTTCCTTTATTCGTAAGCGTTCTGACGTTATCGAAATCCGTGAGCACCTGAAAGCGCACGGCGGCGAAAACATCCACATCATCTCCAAAATCGAAAACCAGGAAGGCCTCAACAACTTCGACGAAATCCTCGAAGCCTCTGACGGCATCATGGTTGCGCGTGGCGACCTGGGTGTAGAAATCCCGGTAGAAGAAGTTATCTTCGCCCAGAAGATGATGATCGAAAAATGTATCCGTGCACGTAAA  >  NC_000913/1756235‑1756513

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