Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,219,936 G→A Q447Q (CAG→CAA)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,219,9360GA100.0% 59.5 / NA 24Q447Q (CAG→CAA) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (13/11);  total (13/11)

GCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGAC  >  NC_000913/4219703‑4220170
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
gCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACAt                                                                                                                                                                                                                           >  1:19653/1‑251 (MQ=255)
   gACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCCTTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACAtttt                                                                                                                                                                                                                        >  1:156449/1‑251 (MQ=255)
            naGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCg                                                                                                                                                                                                                >  1:202706/2‑250 (MQ=255)
                                         gCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTGTCCCTGGCGACATGCGGTTTAAAAACTTCTGTGTCAccc                                                                                                                                                                                  <  1:258849/251‑1 (MQ=255)
                                                          tAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTCTATTGAGCGGCGGATGGTGCCGCTCACTATATGGCTGGCACAAGTGGAAGGTCAGCAGTTGCGCGACCCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCgct                                                                                                                                                                                                                                   >  1:128406/1‑185 (MQ=255)
                                                                                 gcagcGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTGTTATGATGACTATGAAATTTg                                                                                                                                          <  1:87410/251‑1 (MQ=255)
                                                                                       gAAAAAATCACCGATCTTGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGGGGAAGGTCAGCAGTTGCGCGACGCCAGTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTg                                                                                                                                                                                                                   <  1:184592/172‑1 (MQ=255)
                                                                                               cACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGATCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGaa                                                                                                                             >  1:210454/1‑250 (MQ=255)
                                                                                                           cGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCcgcn                                                                                                                <  1:30044/251‑2 (MQ=255)
                                                                                                                            nccATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTT‑CCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCt                                                                                              >  1:204329/2‑251 (MQ=255)
                                                                                                                            gCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTCACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTCCAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCgcgc                                                                                               >  1:329158/1‑251 (MQ=255)
                                                                                                                                            gcggcgGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGAc                                                                                   <  1:112068/247‑1 (MQ=255)
                                                                                                                                                                      ggctggAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTTCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACGTCGGTGTCACCCTTCACGGGCGTTTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGGTAATTTCCTCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCg                                                       <  1:53672/249‑1 (MQ=255)
                                                                                                                                                                           gAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCCCCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGcc                                                <  1:22498/251‑1 (MQ=255)
                                                                                                                                                                                ngnagNAGGTCNGCAGNTGCGCGACGCCATTGNAGAATACGGTAACGCTATTNNCNAACTTNCCgct                                                                                                                                                                                                                                   <  2:128406/63‑1 (MQ=255)
                                                                                                                                                                                        gTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATg                                                                                                                                                 <  1:309391/141‑1 (MQ=255)
                                                                                                                                                                                                     cgACGCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACCTTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGATTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGACGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGGGGTAAAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTc                      >  1:302734/1‑251 (MQ=255)
                                                                                                                                                                                                         gCCATTGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACCTGCTGTTTAAAAACTTCGGTGTCCCCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCtccctgac                                                                                                                                 >  1:180801/1‑133 (MQ=255)
                                                                                                                                                                                                            aTTGAAGAATACGGTAACGCTATTCGCCAACTTGCAGCTGCTAACATTTTCACTGGCGACATACTGATTCAAAACTTCGGTGTCACCCGACACGGGCGTGTGGTTTTTTATGATTACGCTGAAAtt                                                                                                                                            >  1:135774/1‑126 (MQ=255)
                                                                                                                                                                                                              tGAAGNATACGNTNNCGCTATTCGCCAACTTGCCNCTGCTAACATTTTCCNTGGCNACATgntgntn                                                                                                                                                                                                     >  2:131222/1‑66 (MQ=255)
                                                                                                                                                                                                              tGAAGAATACGGTAACGCTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTcc                                  <  1:131222/230‑1 (MQ=255)
                                                                                                                                                                                                                         nnnnnnnnnntCNNNNNCTTNCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGTCATCccgccgcc                                                                                                   <  2:202706/144‑1 (MQ=255)
                                                                                                                                                                                                                         gTAACGCTAGTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTGAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACTATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCTGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGGTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGAc  <  1:291995/251‑1 (MQ=255)
                                                                                                                                                                                                                              ncTATTCGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTAc                                                                                                                                      >  1:79896/2‑114 (MQ=255)
                                                                                                                                                                                                                                    cGCCAACTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCgg                                                                                                                                                                                          >  1:172809/1‑56 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                          
GCGGACACCCAGGAGTTTGAAAACTTTGTGCTGGAGAAGCGGCATATTTCCCCGGCATTAATGGAATTACTGCTTCAGGAAGCAGCGGAAAAAATCACCGATCTCGGCGAACAAATTGTGATTCGCCATCTTTATATTGAGCGGCGGATGGTGCCGCTCAATATCTGGCTGGAACAAGTGGAAGGTCAGCAGTTGCGCGACGCCATTGAAGAATACGGTAACGCTATTCGCCAGCTTGCCGCTGCTAACATTTTCCCTGGCGACATGCTGTTTAAAAACTTCGGTGTCACCCGTCACGGGCGTGTGGTTTTTTATGATTACGATGAAATTTGCTACATGACGGAAGTGAATTTCCGCGACATCCCGCCGCCGCGCTATCCGGAAGACGAACTTGCCAGCGAACCGTGGTACAGCGTCTCGCCGGGCGATGTTTTCCCGGAAGAGTTTCGCCACTGGCTATGCGCCGAC  >  NC_000913/4219703‑4220170

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