Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,679,484 C→T 9.8% T353M (ACG→ATG)  bioA → adenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,679,4840CT 9.8% 135.3 / 5.1 51T353M (ACG→ATG) bioAadenosylmethionine‑‑8‑amino‑7‑oxononanoate aminotransferase BioA
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (3/2);  total (28/23)
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.62e-01

GGCGTTCTATGACGATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCGTGGCCTGAAGGT  >  NC_002947/5679338‑5679631
                                                                                                                                                  |                                                                                                                                                   
ggCGTTCTATGACGATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCAcgcg                                                                                                                                                  >  1:673859/1‑150 (MQ=255)
ggCGTTCTATGACGATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCAcgcg                                                                                                                                                  >  2:705124/1‑150 (MQ=255)
   gTTCTATGACGATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACACGGCCCTGGCCACGCGCAt                                                                                                                                               <  1:1189526/150‑1 (MQ=255)
      cTATGACGATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGc                                                                                                                                            <  1:705124/150‑1 (MQ=255)
             gATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGcc                                                                                                                                     >  1:389194/1‑150 (MQ=255)
                   ccGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGTGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGca                                                                                                                               >  1:1138259/1‑150 (MQ=255)
                        acTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCCGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACAttt                                                                                                                          <  2:342613/150‑1 (MQ=255)
                             gcgcATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCg                                                                                                                     >  2:303916/1‑150 (MQ=255)
                                         aCTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCGCATg                                                                                                         >  2:477533/1‑150 (MQ=255)
                                             gCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGACGCGCTGGCGACCCTGGCTCCATTCGAGCAGGACAAAGTTTTCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGc                                                                                                     >  1:1053147/1‑150 (MQ=255)
                                             gCACAGCTAAACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGc                                                                                                     >  2:1079665/1‑150 (MQ=255)
                                               acaGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGAGCTAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCg                                                                                                   >  1:81261/1‑150 (MQ=255)
                                               acaGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAACAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTgg                                                                                                                        <  2:631969/129‑1 (MQ=255)
                                               acaGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAACAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTgg                                                                                                                        >  1:631969/1‑129 (MQ=255)
                                                     acacCTGCAACCCGCTGGCGTGCGCTGCAGCGCTGGTGACCCTGGACATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATAc                                                                                             <  2:1053147/150‑1 (MQ=255)
                                                       acCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGc                                                                                           <  2:389194/150‑1 (MQ=255)
                                                        ccGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGcc                                                                                          >  2:61680/1‑150 (MQ=255)
                                                         cGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCa                                                                                         >  2:119456/1‑150 (MQ=255)
                                                          ggCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGc                                                                                                                       <  2:235430/119‑1 (MQ=255)
                                                          ggCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGc                                                                                                                       >  1:235430/1‑119 (MQ=255)
                                                              aCCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTg                                                                                    <  2:673859/150‑1 (MQ=255)
                                                              aCCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTg                                                                                    >  1:224256/1‑150 (MQ=255)
                                                                ccGCTGGCGTGCGCTCCCGCGCTGGCGACCCTGGATAGCTTCGAGCAGGACAACGTGGTCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGc                                                                                  >  2:12622/1‑150 (MQ=255)
                                                                   cTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATg                                                                               <  1:61680/150‑1 (MQ=255)
                                                                      gCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGcc                                                                            <  2:224256/150‑1 (MQ=255)
                                                                        gTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGccag                                                                          <  1:357635/150‑3 (MQ=255)
                                                                                                 gATATCTTCGAGCAGGACAACGTGATCGATCCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAag                                                 <  1:306688/150‑1 (MQ=255)
                                                                                                  atatCTTCGATCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCGTGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCATACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAagg                                                <  1:1152872/150‑1 (MQ=255)
                                                                                                               ggACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGaaa                                                               <  1:313669/122‑1 (MQ=255)
                                                                                                               ggACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGaaa                                                               >  2:313669/1‑122 (MQ=255)
                                                                                                               ggACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGc                                   >  2:233745/1‑150 (MQ=255)
                                                                                                                gACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGc                                                                                           <  2:103789/93‑1 (MQ=255)
                                                                                                                gACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGc                                                                                           >  1:103789/1‑93 (MQ=255)
                                                                                                                       tGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTgg                                                                        <  2:789634/105‑1 (MQ=255)
                                                                                                                       tGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTgg                                                                        >  1:789634/1‑105 (MQ=255)
                                                                                                                       tGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCt                           <  1:303916/150‑1 (MQ=255)
                                                                                                                        gATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGtt                                                          <  1:974621/118‑1 (MQ=255)
                                                                                                                        gATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGtt                                                          >  2:974621/1‑118 (MQ=255)
                                                                                                                              ggCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGa                    >  1:490085/1‑150 (MQ=255)
                                                                                                                              ggCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGa                    <  2:604134/150‑1 (MQ=255)
                                                                                                                                  aacaaGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAgcgc                >  2:562516/1‑150 (MQ=255)
                                                                                                                                    caaGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAgcgcgg              >  1:973470/1‑148 (MQ=255)
                                                                                                                                        gCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATc                                                                   <  2:177725/93‑1 (MQ=255)
                                                                                                                                        gCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATc                                                                   >  1:177725/1‑93 (MQ=255)
                                                                                                                                          tgtGGCCATGCGCGTGGCCAGGGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCGTGGCCt        >  2:28894/3‑150 (MQ=255)
                                                                                                                                          tgtGGCCATGCGCGTGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGc                                                                             >  2:1048412/3‑81 (MQ=255)
                                                                                                                                          tgtGGCCATGCGCGTGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGc                                                                             <  1:1048412/79‑1 (MQ=255)
                                                                                                                                          tgtGGCCATGCGCGTGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGaa                                                                <  2:529888/92‑1 (MQ=255)
                                                                                                                                          tgtGGCCATGCGCGTGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGaa                                                                >  1:529888/3‑94 (MQ=255)
                                                                                                                                           cTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCGTGGCCTg       >  2:252442/1‑150 (MQ=255)
                                                                                                                                             ggCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCGTGGCCTGaa     <  2:1138259/150‑1 (MQ=255)
                                                                                                                                             ggCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCGTGGCCTGaa     >  1:370039/1‑150 (MQ=255)
                                                                                                                                                cACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCTTGGCCTGAAGgt  >  2:375096/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGCGTTCTATGACGATTACCCGACACTGCGCGCATTCCTGCACTCGCACAGCTACACCGGCAACCCGCTGGCGTGCGCTGCAGCGCTGGCGACCCTGGATATCTTCGAGCAGGACAACGTGATCGAGGCCAACAAGGCCCTGGCCACGCGCATGGCCACAGCCACTGCACATTTGGCCGATCATGCTCATGTTGCCGAAATACGCCAGACTGGCATGGCCCTGGCCATCGAAATGGTTCAAGACAAGGCAGGCAAGCTTGCCTACCCCTGGCAGGAGCGCCGTGGCCTGAAGGT  >  NC_002947/5679338‑5679631

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: