Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,080,879 C→G 8.5% P236P (CCG→CCC lysC ← aspartokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,080,8790CG 8.5% 139.6 / 4.7 47P236P (CCG→CCClysCaspartokinase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base G (2/2);  total (25/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

TCGGGCCGAGGATCTTGAAGGCCACGCCCGGGGTATCCGGCACGCCGCGAATGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGATCTTCTCCAGGCGACGGGCCTGTG  >  NC_002947/5080737‑5081027
                                                                                                                                              |                                                                                                                                                    
tCGGGCCGAGGATCTTGAAGGCCACGCCCGGGGTATCCGGCACGCCGCGAATGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGaccct                                                                                                                                               <  2:718426‑M2/150‑6 (MQ=255)
  gggCCGAGGATCTTGAAGGCCACGCCCGGGGTATCCGGCACGCCGCGAATGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACcctcc                                                                                                                                             <  1:937300/150‑1 (MQ=255)
                     ccACGCCCGGGGTATCCGGCACGCCGCGAATGGTCAGCTTTGCTTCATCCGGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCtt                                                                                                                                           >  1:811942/1‑133 (MQ=255)
                     ccACGCCCGGGGTATCCGGCACGCCGCGAATGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCtt                                                                                                                                           <  2:811942/133‑1 (MQ=255)
                                      ggCACGCCGCGAATGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTAc                                                                                                         <  2:859948/150‑1 (MQ=255)
                                                  aTGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCAGGCCCAGCGGAACGTTGTACTTGCCGGCGAAc                                                                                             >  1:315240/1‑150 (MQ=255)
                                                      tCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTAATCAATGGTGTTGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCa                                                                                         <  1:883881/150‑1 (MQ=255)
                                                       cAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCAc                                                                                        >  2:1056422/1‑150 (MQ=255)
                                                                ttCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGAt                                                                               >  1:851598/1‑150 (MQ=255)
                                                                ttCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGAt                                                                               >  1:285298/1‑150 (MQ=255)
                                                                   atcaCGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTg                                                                            >  1:102951/1‑150 (MQ=255)
                                                                    tcaCGGTTGAAGGCGATACCGGAAATGATCGGCTGTTGCAGGGATCCCTCTTCATCACTGGTAATGAGGGTACTCGGACCGTCCTCTACGCTGTGTAGCACGCCCGGCGGAATGCTGGAATTCCCAACGCACTCCACCGAACGAATCTgc                                                                           >  2:586527/1‑150 (MQ=255)
                                                                    tcaCGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTgc                                                                           <  2:285298/150‑1 (MQ=255)
                                                                    tcaCGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTATTAATAAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTCTACTTGACGGCGAACTCCACCGAACAGATCTgc                                                                           >  2:630922/1‑150 (MQ=255)
                                                                            gAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCtt                                                                   <  2:851598/150‑1 (MQ=255)
                                                                               ggCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGa                                                                <  1:1056422/150‑1 (MQ=255)
                                                                               ggCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTCCGGGtaccctcattaccattgatgaaga                                                                                                                            <  2:1349812‑M2/89‑25 (MQ=255)
                                                                               ggCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTCCGGGtaccctcattaccattgatgaaga                                                                                                                            >  1:1349812‑M2/1‑65 (MQ=255)
                                                                                gCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTCCGGGTaccctcattaccat                                                                                                                                      <  1:673068‑M2/79‑15 (MQ=255)
                                                                                gCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTCCGGGTaccctcattaccat                                                                                                                                      >  2:673068‑M2/1‑65 (MQ=255)
                                                                                       cGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTCTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCaa                                                           <  1:815354/147‑1 (MQ=255)
                                                                                       cGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTCTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCaa                                                           >  2:815354/1‑147 (MQ=255)
                                                                                                tCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCAttt                                               >  2:1390522/1‑150 (MQ=255)
                                                                                                tCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCAttt                                               >  1:243321/1‑150 (MQ=255)
                                                                                                   gCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCa                                            <  2:102951/150‑1 (MQ=255)
                                                                                                   cgatctCCATGGATTCCTCTTCATCAATGGTAATGGCGCAATCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCa                                            <  1:74236/145‑1 (MQ=255)
                                                                                                    cTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAg                                           <  1:718426/150‑1 (MQ=255)
                                                                                                      gTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTg                                                      <  1:1134427/137‑1 (MQ=255)
                                                                                                      gTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTg                                                      >  2:1134427/1‑137 (MQ=255)
                                                                                                            tGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCtt                                   <  1:1390522/150‑1 (MQ=255)
                                                                                                                 tCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAg                              >  2:1237732/1‑150 (MQ=255)
                                                                                                                   ctctTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGaa                               >  1:572798/1‑147 (MQ=255)
                                                                                                                   ctctTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGaa                               <  2:572798/147‑1 (MQ=255)
                                                                                                                   ctctTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGt                            >  1:1224844/1‑150 (MQ=255)
                                                                                                                    tctTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCa                                            <  2:196461/133‑1 (MQ=255)
                                                                                                                    tctTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCa                                            >  1:196461/1‑133 (MQ=255)
                                                                                                                          tcaATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCCACCGGAGCGATAAGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGAtcttc                     >  2:1416676/1‑150 (MQ=255)
                                                                                                                          tcaATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGAtcttc                     >  1:950625/1‑150 (MQ=255)
                                                                                                                                 tAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTc                                              <  1:1075745/118‑1 (MQ=255)
                                                                                                                                 tAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTc                                              <  1:1075927/118‑1 (MQ=255)
                                                                                                                                 tAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTc                                              >  2:1075927/1‑118 (MQ=255)
                                                                                                                                 tAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTc                                              >  2:1075745/1‑118 (MQ=255)
                                                                                                                                 tAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCCTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGATCTTCTCCAGGc              >  2:310667/1‑150 (MQ=255)
                                                                                                                                   aTGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGATCTTCTCCAGGCGa            >  2:957397/1‑150 (MQ=255)
                                                                                                                                     gAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGATCTTCTCCAGGCGACg          <  1:630922/150‑1 (MQ=255)
                                                                                                                                      aGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTgc                                                                           <  1:397400/84‑1 (MQ=255)
                                                                                                                                      aGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTgc                                                                           >  2:397400/1‑84 (MQ=255)
                                                                                                                                             ccGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGATCTTCTCCAGGCGACGGGCCtgtg  >  2:753740/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCGGGCCGAGGATCTTGAAGGCCACGCCCGGGGTATCCGGCACGCCGCGAATGGTCAGCTTTGCTTCATCACGGTTGAAGGCGATACCGGAAATGATCGGCTGTTCCATGGATTCCTCTTCATCAATGGTAATGAGGGTACCCGGACCCTCCTTGAAGCTGTGCAGCACGCGCAGCGGAACGTTGTACTTGCCGGCGAACTCCACCGAACGGATCTGCAGCACCTTGGAACCAAGGCTGGCCATTTCCAGCATCTCTTCGAAGGTGATCTTCTCCAGGCGACGGGCCTGTG  >  NC_002947/5080737‑5081027

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