Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,218,455	C→A	intergenic (+42/‑141)	aceA → / → aceK	isocitrate lyase/isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,218,455	0	C	A	100.0%	180.6 / NA	52	intergenic (+42/‑141)	aceA/aceK	isocitrate lyase/isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (29/23);  total (29/23)

TACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGG  >  NC_000913/4218320‑4218588                                                                                                                                        |                                                                                                                                      tACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAAcgc                                                                                                                                   >  1:46461/1‑140 (MQ=255)           ttCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                                         >  2:215290/1‑140 (MQ=255)              gATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAt                                                                                                                      <  1:109746/140‑1 (MQ=255)              gATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAt                                                                                                                      <  2:559097/140‑1 (MQ=255)                tAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCg                                                                                                                    >  2:480167/1‑140 (MQ=255)                tAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCg                                                                                                                    >  2:507457/1‑140 (MQ=255)                   aGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCgg                                                                                                                   >  2:637297/1‑138 (MQ=255)                   aGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCgg                                                                                                                   <  1:637297/138‑1 (MQ=255)                       acgacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACg                                                                                                             <  1:480167/140‑1 (MQ=255)                         gacTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGAt                                                                                                           >  2:723408/1‑140 (MQ=255)                          acTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAt                                                                                                                      <  1:251065/128‑1 (MQ=255)                          acTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAAt                                                                                                                      >  2:251065/1‑128 (MQ=255)                           cTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGc                                                                                                         <  1:507457/140‑1 (MQ=255)                            tattatTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGcc                                                                                                        >  2:456327/1‑140 (MQ=255)                             attattCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCt                                                                                                       <  2:377670/140‑1 (MQ=255)                               tattCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGa                                                                                                     >  1:486476/1‑140 (MQ=255)                                      ggcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACg                                                                                                             >  2:60741/1‑125 (MQ=255)                                      ggcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACg                                                                                                             <  1:60741/125‑1 (MQ=255)                                      ggcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACg                                                                                              >  2:72067/1‑140 (MQ=255)                                      ggcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACg                                                                                              >  2:504287/1‑140 (MQ=255)                                      ggcggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACg                                                                                              >  2:675387/1‑140 (MQ=255)                                        cggcACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCt                                                                                            >  2:7507/1‑140 (MQ=255)                                           cACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgc                                                                                         >  1:503884/1‑140 (MQ=255)                                           cACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgc                                                                                         >  2:712010/1‑140 (MQ=255)                                           cACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgc                                                                                         >  1:313998/1‑140 (MQ=255)                                              gTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGcgcg                                                                                      <  1:291772/140‑1 (MQ=255)                                                cTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTc                                                                                    <  1:72067/140‑1 (MQ=255)                                                          cgcgCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgc                                                                                         >  2:562161/1‑125 (MQ=255)                                                          cgcgCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTgc                                                                                         <  1:562161/125‑1 (MQ=255)                                                          cgcgCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGcc                                                                          >  2:371684/1‑140 (MQ=255)                                                          cgcgCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGcc                                                                          >  1:720972/1‑140 (MQ=255)                                                           gcgcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGACAATCGGTGCACGAGGCCTGATGCGACGCTTGCGCGTCTTATCAGGCCt                                                                         <  2:367491/140‑1 (MQ=25)                                                            cgcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                        >  2:648726/1‑140 (MQ=255)                                                            cgcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                        <  1:456327/140‑1 (MQ=255)                                                            cgcTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTa                                                                        <  1:215290/140‑1 (MQ=255)                                                                  ccGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCg                                                                  <  2:46461/140‑1 (MQ=255)                                                                      cTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGc                                                              <  1:648726/140‑1 (MQ=255)                                                                           cTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGc                                                              >  2:667899/1‑135 (MQ=255)                                                                           cTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGc                                                              <  1:667899/135‑1 (MQ=255)                                                                                gaaTCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAgg                                                    >  1:223510/1‑140 (MQ=255)                                                                                      cAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                              <  1:371684/140‑1 (MQ=255)                                                                                        gTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATaa                                            <  1:7507/140‑1 (MQ=255)                                                                                                   caacaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTAcgc                                 >  1:166576/1‑140 (MQ=255)                                                                                                     acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                              <  2:537348/125‑1 (MQ=255)                                                                                                     acaacCGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGAt                                              >  1:537348/1‑125 (MQ=255)                                                                                                          cGTTGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATcc                          >  2:529513/1‑140 (MQ=255)                                                                                                            ttGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCgg                        >  1:639289/1‑140 (MQ=255)                                                                                                             tGCTGACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGc                       >  1:415660/1‑140 (MQ=255)                                                                                                                 gACTGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAAtt                   <  2:720972/140‑1 (MQ=255)                                                                                                                              aTAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCTTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgat      <  1:479757/140‑1 (MQ=255)                                                                                                                              aTAAGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgat      >  2:185280/1‑140 (MQ=255)                                                                                                                                aaGGCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCtgatga    <  2:29552/140‑1 (MQ=255)                                                                                                                                  ggCGTTAACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAgg  <  2:313998/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      TACAGGTTACTTCGATAAAGTGACGACTATTATTCAGGGCGGCACGTCTTCAGTCACCGCGCTGACCGGCTCCACTGAAGAATCGCAGTTCTAAGCAACAACAACCGTTGCTGACTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGGTGCACGATGCCTGATGCGACGCTTGCGCGTCTTATCATGCCTACAGCCGTTGCCGAACGTAGGCTGGATAAGGCGTTTACGCCGCATCCGGCAATTCTCTGCTCCTGATGAGG  >  NC_000913/4218320‑4218588

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