Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,219,615	C→T	D340D (GAC→GAT)	aceK →	isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,219,615	0	C	T	100.0%	161.9 / NA	46	D340D (GAC→GAT)	aceK	isocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/23);  total (23/23)

GCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCAT  >  NC_000913/4219485‑4219747                                                                                                                                   |                                                                                                                                     gCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTcgcg                                                                                                                             <  2:544188/140‑1 (MQ=255)   tACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGcgcc                                                                                                                           <  1:387454/140‑1 (MQ=255)     cctcgAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGcgccgc                                                                                                                         >  2:380351/4‑140 (MQ=255)        cgAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGa                                                                                                                      >  2:705590/1‑140 (MQ=255)           atatCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAag                                                                                                                   <  1:421406/140‑1 (MQ=255)                  gTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTc                                                                                                            <  1:125427/140‑1 (MQ=255)                    ttATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTg                                                                                                          <  2:768477/140‑1 (MQ=255)                    ttATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTg                                                                                                          <  2:119571/140‑1 (MQ=255)                            aGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCAc                                                                                                  >  2:390747/1‑140 (MQ=255)                              ggCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGt                                                                                                >  2:246537/1‑140 (MQ=255)                                cTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTg                                                                                                          >  2:366968/1‑128 (MQ=255)                                cTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTg                                                                                                          <  1:366968/128‑1 (MQ=255)                                 tGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCTAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCg                                                                                             >  2:62366/1‑140 (MQ=255)                                 tGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCg                                                                                             >  2:469119/1‑140 (MQ=255)                                 tGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCg                                                                                             >  1:764394/1‑140 (MQ=255)                                 tGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCg                                                                                             >  2:752500/1‑140 (MQ=255)                                 tGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCg                                                                                             <  1:246537/140‑1 (MQ=255)                                   tAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTg                                                                                           >  2:287265/1‑140 (MQ=255)                                    aaTGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGc                                                                                          <  1:156868/140‑1 (MQ=255)                                    aaTGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGc                                                                                          <  1:287265/140‑1 (MQ=255)                                       gAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTg                                                                                       <  1:705590/140‑1 (MQ=255)                                        aGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGc                                                                                      >  2:27121/1‑140 (MQ=255)                                         gCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCt                                                                                     >  2:169085/1‑140 (MQ=255)                                          cAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTa                                                                                    >  2:605135/1‑140 (MQ=255)                                           aGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTAt                                                                                   <  1:169085/140‑1 (MQ=255)                                              tCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCaa                                                                                <  2:132410/140‑1 (MQ=255)                                               cATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAAc                                                                               >  1:744888/1‑140 (MQ=255)                                               cATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAAc                                                                               <  2:764727/140‑1 (MQ=255)                                               cATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAAc                                                                               <  2:516344/140‑1 (MQ=255)                                                     aGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCTAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGa                                                                         <  1:62366/140‑1 (MQ=255)                                                     aGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGa                                                                         <  1:752500/140‑1 (MQ=255)                                                       cgcCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGaaa                                                                       >  2:724538/1‑140 (MQ=255)                                                       cgcCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGaaa                                                                       >  2:41264/1‑140 (MQ=255)                                                         ccGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAaga                                                                     <  1:142489/140‑1 (MQ=255)                                                                        aTGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTggg                                                      <  2:764394/140‑1 (MQ=255)                                                                          ggTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGcc                                                    >  2:535286/1‑140 (MQ=255)                                                                           gTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCg                                                   >  2:171356/1‑140 (MQ=255)                                                                                   ggTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCgg                                           <  1:535286/140‑1 (MQ=255)                                                                                    gggtTTTCGCTGCCGGGCTTTGATCGGCTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTGCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGa                                          <  2:218826/138‑1 (MQ=255)                                                                                          aCGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCa                                    >  2:315650/1‑140 (MQ=255)                                                                                          aCGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCa                                    >  2:238861/1‑140 (MQ=255)                                                                                            gCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAgg                                  <  1:171356/140‑1 (MQ=255)                                                                                                                      aaGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAg        <  1:315650/140‑1 (MQ=255)                                                                                                                       aGTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGc       >  2:57644/1‑140 (MQ=255)                                                                                                                        gTCATCAAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCg      <  1:469119/140‑1 (MQ=255)                                                                                                                            tcaAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCat  >  2:504523/1‑140 (MQ=255)                                                                                                                            tcaAAGATAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGGTTGAAAACTTTGTGCTGGAGAAGCGGCat  >  2:128785/1‑140 (MQ=255)                                                                                                                                   |                                                                                                                                     GCTACCGCGAATATCTCGTTTATCTACAGGGCTGTAATGAGCAGTTCATTGAAGCGCCGGGTATTCGTGGAATGGTGATGTTGGTGTTTACGCTGCCGGGCTTTGATCGGGTATTCAAAGTCATCAAAGACAGGTTCGCGCCGCAGAAAGAGATGTCTGCCGCTCACGTTCGTGCCTGCTATCAACTGGTGAAAGAGCACGATCGCGTGGGCCGAATGGCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCAT  >  NC_000913/4219485‑4219747

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: ‑