Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,679,581 G→A G55R (GGA→AGA)  kdgK → 2‑dehydro‑3‑deoxygluconokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,679,5810GA90.9% 48.6 / NA 24G55R (GGA→AGA) kdgK2‑dehydro‑3‑deoxygluconokinase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  major base A (12/8);  minor base C (2/0);  total (15/9)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AAACCACCACAGCGCAAACTAACGCTAATTTTTTACAGATCAGGTTCACGACTATGTCCAAAAAGATTGCCGTGATTGGCGAATGCATGATTGAGCTTTCCGAGAAAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAGCCGCCGCCAAATTCTGGCTGGAGAGTGAGCAGTCTGCGGCGATTTGCGAAGAGCTGGCGAATTTCGATTATCTCTACCTGAGCGGG  >  NC_000913/3679366‑3679823
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
aaaCCACCACAGCGCAAACTAACGCTAATTTTTTACAGATCAGGTTCACGACTATGTCCAAAAAGATTGCCGTGATTGGCGAATGCATGATTGAGCTTTCCGGGAAAGGCGCGGACGTTAAGCGCGTTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTATCCGGTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGcc                                                                                                                                                                                                                 <  1:351773/251‑1 (MQ=255)
aaaCCACCACAGCGCAAACTAACGCTAATTTTTTACAGATCAGGTTCACGACTATGTCCAAAAAGATTGCCGTGATTGGCGAATGCATGATTGAGCTTTCCGGGAAAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGcc                                                                                                                                                                                                                 <  1:351761/251‑1 (MQ=255)
        acaGCGCAAACTAACGCTAATTTTTTACAGATCAGGTTCACGACTATGTCCACAAAGATTGCCGTGATTGGCGAATGCATGATTGAGCTTTCCGAGAAAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGATATCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGc                                                                                                                                                                                                             >  1:162131/1‑247 (MQ=255)
                                              cACGACTATGTCCAAAAAGATTGCAGTGATTGGCGAATGCATGATTTAGATTTCCGAGCCAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGCTACCCTGAACACTTCCGTCTATCTCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGCGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATg                                                                                                                                                                       >  1:319827/1‑247 (MQ=255)
                                                                           ttGGCGAATGCATGATTGAGCTTTCCGAGAAAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGCGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGaaa                                                                                                                                                                   >  1:4111/1‑222 (MQ=255)
                                                                                                   ccGAGAAAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACg                                                                                                              >  1:3830/1‑251 (MQ=255)
                                                                                                                ggACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAg                                                                                                                              >  1:190015/1‑222 (MQ=255)
                                                                                                                 naCGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGTTCCTGCGGCATTACACGTTCATTACGCAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGCTGCTGGACGCCTGGCACGGCGAGCACGTTGATACTTCCCTGa                                                                                                                                                                                  >  1:64875/2‑169 (MQ=255)
                                                                                                                 gACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCGGTCTATATCGCAGGTCAGGTCGATCCTGCGGCATTAAGAGTTCATTACGTACCGGCGCTGAGAACGGACAGTTTTAGCCAGAAGATGCTGGACGCCTGGCACGGAGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACTGCCAGCACCGGCGAGCGTACGTTc                                                                                                           >  1:255588/1‑240 (MQ=255)
                                                                                                                         gcgcGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTACCGGCGCTGAGAACGGACAGTTTTCGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGCACGa                                                                                          >  1:349848/1‑249 (MQ=255)
                                                                                                                          cgcgGTTTCGGCGGCGATACCCTGAACGCTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAg                                                                                        <  1:5944/250‑1 (MQ=255)
                                                                                                                                        cgATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTGACCGTTCATTACGTAACGGCGCTGAGAACGGACNNTTTTAgnn                                                                                                                                                                                                                                  >  2:19457/1‑96 (MQ=255)
                                                                                                                                         gATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGc                                                                                                                                                                                                        >  1:7468/1‑123 (MQ=255)
                                                                                                                                          aTACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCGGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAGCCGCCGCCAAATTctgg                                                                       <  1:294417/251‑1 (MQ=255)
                                                                                                                                                                 anNGCCCGTNNGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGc                                                                                                                                                                                                        <  2:7468/99‑1 (MQ=255)
                                                                                                                                                                         tCAGGGCGAGCCTGCGGCATTAACCGTTCAGTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAGCCGCCGCCAAATACTGGCTGGAGAGTGAGCAGTCTGCGGCGATTTGc                                         <  1:271467/250‑1 (MQ=255)
                                                                                                                                                                            ngTCGATCCTGCGGCCTTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAGCCGCCGCCAAATTCTGGCTGGAGAGTGAGCAGTCTGCGGCGATTTGCGa                                       >  1:91000/2‑249 (MQ=255)
                                                                                                                                                                            ggTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAgccgccgc                                                                                 >  1:77211/1‑207 (MQ=255)
                                                                                                                                                                             gTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGCCGCCTGGCACGGCGNNAACGTTNNt                                                                                                                                                                                            >  2:11590/1‑99 (MQ=255)
                                                                                                                                                                             gTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTGTTAGCCAGCAGATGCTGGACGCCTGGCCCGGCGAGAACGTTGATACTTCCCTGAc                                                                                                                                                                                 <  1:11590/110‑1 (MQ=255)
                                                                                                                                                                                           aanaacngnncnnaccNTAACGGCGCTGAGAACGGACANNTTTAGCCAGNAGAGGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGNTGGAAAACCNNCGNCCGGGNCNNNCNNNCANTGAAACCGACAGCACCGGCNAGNNTACGTTNTACTACTGGCGTATCGAAgccgccgc                                                                                 <  2:77211/177‑1 (MQ=255)
                                                                                                                                                                                             tAACCGTTCATTACGTAACGGCGCTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGCAAACCGTCTGCCGGGCCTTTACTACATNGCANNNNNCannnnnnnnnn                                                                                                                     >  2:223331/1‑144 (MQ=255)
                                                                                                                                                                                               aNCNTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCNAGCTNANNCTNGNCGNCTGGc                                                                                                                                                                                                             <  2:162131/64‑1 (MQ=255)
                                                                                                                                                                                                                ngcgcTGAGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGTGGAACGAAGCCGCCGCCAAATTCTGGCTGGAGAGTGAGCAGTCTGCGGCGATTTGCGAAGAGCTGGCGAATTTCGATTATCTCTACCTGAGCggg  >  1:86301/2‑250 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
AAACCACCACAGCGCAAACTAACGCTAATTTTTTACAGATCAGGTTCACGACTATGTCCAAAAAGATTGCCGTGATTGGCGAATGCATGATTGAGCTTTCCGAGAAAGGCGCGGACGTTAAGCGCGGTTTCGGCGGCGATACCCTGAACACTTCCGTCTATATCGCCCGTCAGGTCGATCCTGCGGCATTAACCGTTCATTACGTAACGGCGCTGGGAACGGACAGTTTTAGCCAGCAGATGCTGGACGCCTGGCACGGCGAGAACGTTGATACTTCCCTGACCCAACGGATGGAAAACCGTCTGCCGGGCCTTTACTACATTGAAACCGACAGCACCGGCGAGCGTACGTTCTACTACTGGCGGAACGAAGCCGCCGCCAAATTCTGGCTGGAGAGTGAGCAGTCTGCGGCGATTTGCGAAGAGCTGGCGAATTTCGATTATCTCTACCTGAGCGGG  >  NC_000913/3679366‑3679823

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: