Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 370,911 C→T 6.9% P131L (CCT→CTT)  dgcA → dimethylglycine dehydrogenase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947370,9110CT 6.9% 168.4 / 2.2 58P131L (CCT→CTT) dgcAdimethylglycine dehydrogenase subunit
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base T (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

CATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAC  >  NC_002947/370768‑371053
                                                                                                                                               |                                                                                                                                              
catcCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCCACGCTGATGTCGCCATCGGGCATTCGCGAACCtgtg                                                                                                                                            >  8:61130/1‑148 (MQ=255)
     cGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCCGGCATTCGCGAACCTGTGCACcgc                                                                                                                                      >  4:50700/1‑149 (MQ=255)
      gCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACcgcg                                                                                                                                     >  8:135411/1‑149 (MQ=255)
            ccAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCt                                                                                                                               >  8:90130/1‑149 (MQ=255)
             cAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTg                                                                                                                              >  6:317782/1‑149 (MQ=255)
                 aCCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCaaa                                                                                                                          <  3:209308/149‑1 (MQ=255)
                   ctggCCGACGCCCTGCACAAGCTTGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAAc                                                                                                                        <  5:275933/146‑1 (MQ=255)
                       ccGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGATCCGGTGCACCGCGCCACCTGCAAAACCATc                                                                                                                    >  6:235139/1‑149 (MQ=255)
                       ccGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATc                                                                                                                    <  8:121499/149‑1 (MQ=255)
                        cgccGCCATGCCCAAGCATGGCGCCAAGATCATGATCCAGCTTACCCACACCGCCCGTCCCGCGCGCTGGGACGGCTTCAACTGGCCGACGCCGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCg                                                                                                                   <  8:12956/146‑1 (MQ=255)
                         gACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCa                                                                                                                                  <  3:225356/133‑1 (MQ=255)
                         gACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCa                                                                                                                                  >  4:225356/1‑133 (MQ=255)
                         gACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGa                                                                                                                  <  3:7898/149‑1 (MQ=255)
                          aCGCCATGCACAAGCATGGCGCCAAGATCATGGTCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAg                                                                                                                 <  7:260119/149‑1 (MQ=255)
                            gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg                                                                                                                >  4:255088/1‑148 (MQ=255)
                            gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg                                                                                                                >  4:250772/1‑148 (MQ=255)
                            gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg                                                                                                                <  3:255088/148‑1 (MQ=255)
                            gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACTTGTGCACCGCGCCACCTGCAAAACCATCGAgg                                                                                                                <  3:250772/148‑1 (MQ=255)
                            gCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGt                                                                                                               >  5:159495/1‑149 (MQ=255)
                                  cacaAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAaga                                                                                                         <  7:61130/149‑1 (MQ=255)
                                       gCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCt                                                                                                    >  2:83063/1‑149 (MQ=255)
                                        cATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTg                                                                                                   >  3:260985/1‑149 (MQ=255)
                                            gcgcCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCgg                                                                                               <  6:57069/149‑1 (MQ=255)
                                             cgcCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCggg                                                                                              <  5:317782/149‑1 (MQ=255)
                                                    aTCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCgg                                                                                       <  6:238058/149‑1 (MQ=255)
                                                           tCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTAc                                                                                <  4:136269/149‑1 (MQ=255)
                                                                     ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACgcgc                                                                            <  2:42094/143‑1 (MQ=255)
                                                                     ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACgcgc                                                                            >  1:42094/1‑143 (MQ=255)
                                                                     ccACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAAACATAGAGGTGGAAGAGATCTGGCGGCTAATCGGCAACTAAGCGCAGCCTg                                                                      >  4:42281/1‑149 (MQ=255)
                                                                       acaTGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCa                                                                                                                      >  2:278185/1‑99 (MQ=255)
                                                                            ggCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTc                                                                <  7:90130/148‑1 (MQ=255)
                                                                                  cgctcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGACGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCaa                                                         <  3:9425/149‑1 (MQ=255)
                                                                                   gctcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCa                                                          >  5:198631/1‑147 (MQ=255)
                                                                                     tcgcGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAACACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCGaaa                                                        >  4:277532/1‑147 (MQ=255)
                                                                                       gcgcTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAgg                                                    >  6:265698/1‑149 (MQ=255)
                                                                                          cTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGggcg                                                 <  7:135411/149‑1 (MQ=255)
                                                                                                        aCTGGCCGACGCTTATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                   <  3:42281/149‑1 (MQ=255)
                                                                                                        aCTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                   <  4:260985/149‑1 (MQ=255)
                                                                                                        aCTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                   >  3:78102/1‑149 (MQ=255)
                                                                                                         cTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                   >  5:272737/1‑148 (MQ=255)
                                                                                                            gCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAt                               >  8:35170/1‑149 (MQ=255)
                                                                                                             ccGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTATGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAAtt                              <  1:83063/149‑1 (MQ=255)
                                                                                                              cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAg                                                                          >  3:272830/1‑104 (MQ=255)
                                                                                                              cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAg                                                                          <  4:272830/104‑1 (MQ=255)
                                                                                                              cGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTg                             >  3:228517/1‑149 (MQ=255)
                                                                                                               gACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCGGGACGGCGTGGAATTGt                            <  7:303745/149‑1 (MQ=255)
                                                                                                                  gCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCgg                         <  7:35170/149‑1 (MQ=255)
                                                                                                                       tGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTg                    <  3:44895/149‑1 (MQ=255)
                                                                                                                            ccATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCa               <  6:198631/149‑1 (MQ=255)
                                                                                                                              aTCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGc             <  2:299102/149‑1 (MQ=255)
                                                                                                                               tCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGca            <  4:228517/149‑1 (MQ=255)
                                                                                                                                 gggCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGcaca          >  1:59611/1‑149 (MQ=255)
                                                                                                                                       tCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCg    >  4:132279/1‑149 (MQ=255)
                                                                                                                                         gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                   <  1:63449/116‑1 (MQ=255)
                                                                                                                                         gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTg                                   >  2:63449/1‑116 (MQ=255)
                                                                                                                                         gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa                       >  2:189273/1‑128 (MQ=255)
                                                                                                                                         gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCa                       <  1:189273/128‑1 (MQ=255)
                                                                                                                                         gcgAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAc  >  1:308603/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CATCCCGCACTTCCAGAACCTTGCCGACGCCATGCACAAGCATGGCGCCAAGATCATGATCCAGATTACCCACATGGGCCGTCGCTCGCGCTGGGACGGCTTCAACTGGCCGACGCTGATGTCGCCATCGGGCATTCGCGAACCTGTGCACCGCGCCACCTGCAAAACCATCGAGGTGGAAGAGATCTGGCGGGTAATCGGCAACTACGCGCAGGCTGCGCGTCGCGCCAAAGAGGGCGGCCTGGACGGCGTGGAATTGTCGGCAGTGCACCAGCACATGATCGAC  >  NC_002947/370768‑371053

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