Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,171,893 Δ1 bp 8.2% coding (465/1836 nt) glmS ← L‑glutamine/D‑fructose‑6‑phosphate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,171,8930G. 8.2% 151.6 / 8.0 49coding (465/1836 nt)glmSL‑glutamine/D‑fructose‑6‑phosphate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (20/25);  new base . (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTC  >  NC_002947/6171750‑6172039
                                                                                                                                               |                                                                                                                                                   
tCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                              <  1:36065/150‑1 (MQ=255)
tCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAAGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                              <  1:167608/150‑1 (MQ=255)
         tGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCAcc                                                                                                                                     <  1:732275/150‑1 (MQ=255)
                   ccAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTATCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTaa                                                                                                                           <  1:607830/150‑1 (MQ=255)
                     aGCTGGTCAGAAGAAAGGTAGTGCTGCCCAGGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAgg                                                                                                                         <  2:569000/150‑1 (MQ=255)
                           tCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCg                                                                                                                   >  2:180341/1‑150 (MQ=255)
                            cGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGcggcgg                                                                                                                  <  1:756812/150‑1 (MQ=255)
                                  ccAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTGAGCGCCGTGCAGGCGCTTCACCGACGCCCTTAAGGCGGCGGtgaagc                                                                                                            <  1:136895/150‑6 (MQ=255)
                                   cAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGCTGGGGAGCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGTTTCACAGACGAAATAAGGGCGTCGGTCAGGTc                                                                                                           >  1:772881/1‑150 (MQ=255)
                                           tctcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCTGCACGGCGCttac                                                                                                                                       >  2:478107/1‑112 (MQ=255)
                                           tctcGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCTGCACGGCGCttac                                                                                                                                       <  1:478107/114‑3 (MQ=255)
                                              cGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGctct                                                                                                >  1:10721/1‑150 (MQ=255)
                                                     gCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                              <  2:716718/97‑1 (MQ=255)
                                                     gCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAgg                                                                                                                                              >  1:716718/1‑97 (MQ=255)
                                                              gATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAggg                                                                                          <  2:395318/140‑1 (MQ=255)
                                                              gATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAggg                                                                                          >  1:395318/1‑140 (MQ=255)
                                                                  aCCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAgg                                                                            <  1:426963/150‑1 (MQ=255)
                                                                    cAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTg                                                                          >  2:455775/1‑150 (MQ=255)
                                                                          gCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGAt                                                                    <  2:10721/150‑1 (MQ=255)
                                                                            tACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGa                                                                  >  1:79255/1‑150 (MQ=255)
                                                                            tACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGa                                                                  <  1:180341/150‑1 (MQ=255)
                                                                            tACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTATTCAGGGTGTGGTGGATCAGGTGGACGATGa                                                                  <  2:78064/150‑1 (MQ=255)
                                                                               cACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGt                                                                                                                                                    >  1:95742/1‑66 (MQ=255)
                                                                               cACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGt                                                                                                                                                    <  2:95742/66‑1 (MQ=255)
                                                                               cACTGCATGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCCGCGCCACCCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCCTAAGGGCGTCGGTCAGATCTGGAATGCGCTTCAGGGAGTGGTGGATCAGGTGGACGATGAcat                                                               >  2:215915/1‑148 (MQ=255)
                                                                                aCTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTc                                                              <  2:79255/150‑1 (MQ=255)
                                                                                 ttGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                             >  2:399601/2‑150 (MQ=255)
                                                                                  tGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGgcgc                                                                                                                                           >  1:613272/1‑71 (MQ=255)
                                                                                  tGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGgcgc                                                                                                                                           <  2:613272/71‑1 (MQ=255)
                                                                                     gTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTAt                                                         <  2:655742/150‑1 (MQ=255)
                                                                                           ccaccaGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCt                                                   >  1:700248/1‑150 (MQ=255)
                                                                                                     gTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGCCGATGACTTCGGTATTGGTCGGCGAAGTCaa                                         <  2:772881/150‑1 (MQ=255)
                                                                                                     gTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTAAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGaa                                         >  1:647367/1‑150 (MQ=255)
                                                                                                      tCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAg                                        >  1:505651/1‑150 (MQ=255)
                                                                                                            tgctTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAg                                  <  1:399601/150‑1 (MQ=255)
                                                                                                               ttGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGccc                               <  1:108868/150‑1 (MQ=255)
                                                                                                                     cTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATa                                   <  1:255635/140‑1 (MQ=255)
                                                                                                                     cTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATa                                   >  2:255635/1‑140 (MQ=255)
                                                                                                                        aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtgg                                                                                     >  1:286292/1‑87 (MQ=255)
                                                                                                                        aTCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTgtgg                                                                                     <  2:286292/87‑1 (MQ=255)
                                                                                                                         tCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATATCCCAGGCCCTTCa                     >  2:221848/1‑150 (MQ=255)
                                                                                                                                cAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                             <  2:429502/103‑1 (MQ=255)
                                                                                                                                cAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCg                                                             >  1:429502/1‑103 (MQ=255)
                                                                                                                                    ggtgaaGCGCCTGCACGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCgcgc          <  2:442230/146‑1 (MQ=255)
                                                                                                                                    ggtgaaGCGCCTGCACGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCgcgc          >  1:442230/5‑150 (MQ=255)
                                                                                                                                     cGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATATCCCAGGCCCTTCAGCTCTTCGCGCa         <  1:221848/150‑1 (MQ=255)
                                                                                                                                       tAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCCAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGc       >  2:543905/1‑150 (MQ=255)
                                                                                                                                           cgcCGTGCAGGAGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGACGTAAAAGCGTAGCGCAGGCCCTTCAGCTCTTCGCGCAGctcct   >  2:662382/1‑150 (MQ=255)
                                                                                                                                            gcCGTGCAGGCGCTTCACCGCCGCCTGAATGGCGTCGGTCAGGTCTGGATTGCTCTTCAGGGTCTGCCGATTTAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGctcctc  <  1:215915/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
TCAGTAACCTGACGCAGCGCCAGCTGGTCGGACGCCAGGAAGTTCTCGCCAAGGCCCAGGCCGATAACCAGCGGGCTACCACTGCGTGCTGCCACCAGGCGGTCCGGCTGCTTGGCACTGATCAGCGCCAGGCCGTAAGCGCCGTGCAGGCGCTTCACCGCCGCCTTAAGGGCGTCGGTCAGGTCTGGAATGCTCTTCAGGGTGTGGTGGATCAGGTGGACGATGACTTCGGTATCGGTCTGCGAAGTGAAGACATAGCCCAGGCCCTTCAGCTCTTCGCGCAGCTCCTC  >  NC_002947/6171750‑6172039

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