Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,712,841 T→C 8.0% T755A (ACG→GCG)  metH ← cobalamin‑dependent methionine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,712,8410TC 8.0% 138.7 / 2.3 50T755A (ACG→GCG) metHcobalamin‑dependent methionine synthase
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base C (2/2);  total (26/24)
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.99e-01

TTCTGCTCGCGGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTACCCGCGCCGAA  >  NC_002947/2712708‑2712979
                                                                                                                                     |                                                                                                                                          
ttCTGCTCGCGGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCtt                                                                                                                             >  1:182643/1‑149 (MQ=255)
  cTGCTCGCGGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATc                                                                                                                               >  4:99544/1‑145 (MQ=255)
  cTGCTCGCGGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATc                                                                                                                               <  3:99544/145‑1 (MQ=255)
        gcgTGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCTCGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTgg                                                                                                                     <  1:261371/149‑1 (MQ=255)
        gcgGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTTCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTg                                                                                                                      <  6:43689/148‑1 (MQ=255)
        gcgGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCtt                                                                                                                       <  8:110928/147‑1 (MQ=255)
        gcgGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTg                                                                                                                      <  6:38725/148‑1 (MQ=255)
                    cAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGt                                                                                                         >  1:183218/1‑149 (MQ=255)
                     aGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTc                                                                                                        >  2:212638/1‑149 (MQ=255)
                       gATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACACCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGc                                                                                                      >  2:136635/1‑149 (MQ=255)
                         tcttctCCGCCGGTTCCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGCGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCt                                                                                                    >  8:83542/1‑149 (MQ=255)
                         tcttctCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGCGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCt                                                                                                    >  2:167170/1‑149 (MQ=255)
                                cgccgGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTc                                                                                               <  7:137047/147‑1 (MQ=255)
                                cgccgGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCgg                                                                                             >  3:88121/1‑149 (MQ=255)
                                      tACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGa                                                                                       >  3:268764/1‑149 (MQ=255)
                                        ccATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGAt                                                                                      >  4:8890/1‑148 (MQ=255)
                                          aTCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGaa                                                                                   >  3:9993/1‑149 (MQ=255)
                                             aCGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACgg                                                                                <  2:183218/149‑1 (MQ=255)
                                                ccAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCCGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGAt                                                                             <  6:332522/149‑1 (MQ=255)
                                                    ggTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAgg                                                                         >  5:190087/1‑149 (MQ=255)
                                                     gTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGt                                                                        <  4:268764/149‑1 (MQ=255)
                                                      tCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGt                                                                        >  8:242388/1‑148 (MQ=255)
                                                        gacgaTGTCGTATCCGTTACAGCCCAGCACCACACCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGa                                                                      <  1:136635/148‑1 (MQ=255)
                                                                        ttACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGCGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCAt                                                     <  7:83542/149‑1 (MQ=255)
                                                                        ttACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGCGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCAt                                                     <  1:167170/149‑1 (MQ=255)
                                                                         tACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCAt                                                     <  4:88121/148‑1 (MQ=255)
                                                                           cAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGAc                                                  <  2:182643/149‑1 (MQ=255)
                                                                             gCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGAcgc                                                >  2:147092/1‑149 (MQ=255)
                                                                                 aGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTg                                                                                                                            >  5:3152/1‑69 (MQ=255)
                                                                                 aGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTg                                                                                                                            <  6:3152/69‑1 (MQ=255)
                                                                                     ccacGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGAc                                        <  7:242388/149‑1 (MQ=255)
                                                                                          ccGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGac                                   >  5:141869/1‑149 (MQ=255)
                                                                                           cGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCAt                                                     <  2:68692/130‑1 (MQ=255)
                                                                                           cGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCAt                                                     >  1:68692/1‑130 (MQ=255)
                                                                                           cGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGacc                                  >  4:81284/1‑149 (MQ=255)
                                                                                           cGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGacc                                  <  3:8890/149‑1 (MQ=255)
                                                                                           cGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGac                                   >  4:304617/1‑148 (MQ=255)
                                                                                                      ttGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCg                           <  6:152251/145‑1 (MQ=255)
                                                                                                      ttGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCg                           >  5:152251/1‑145 (MQ=255)
                                                                                                        gCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAgg                      >  5:286538/1‑148 (MQ=255)
                                                                                                            aTGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACAt                 <  6:190087/149‑1 (MQ=255)
                                                                                                                 gTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTa             <  3:304617/148‑1 (MQ=255)
                                                                                                                 gTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTa             <  1:147092/148‑1 (MQ=255)
                                                                                                                 gTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTAc            >  4:304004/1‑149 (MQ=255)
                                                                                                                    cACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTAccc          <  1:212638/148‑1 (MQ=255)
                                                                                                                     aCGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTACCcgc        >  2:134722/1‑149 (MQ=255)
                                                                                                                       gTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTACCCgcgc      >  6:293056/1‑149 (MQ=255)
                                                                                                                       gTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTACCCgcgc      >  6:298097/1‑149 (MQ=255)
                                                                                                                          gCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCgg                                              <  6:245296/106‑1 (MQ=255)
                                                                                                                             tttCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTACCCGCGCCGaa  >  3:111638/1‑147 (MQ=255)
                                                                                                                                     |                                                                                                                                          
TTCTGCTCGCGGGCCGTTTGCAGGATCTTCTCCGCCGGTACCATCACGCCAAGGTCGACGATGTCGTATCCGTTACAGCCCAGCACCACGCCGACAATGTTCTTGCCGATGTCGTGCACGTCGCCTTTCACCGTGGCCATCAGGATCTTGCCCTTGGCTTCCGGCTTGTCGCCTTTTTCGGCTTCGATGAACGGGATCAGGTGAGCAACGGCCTGTTTCATGACGCGGGCCGACTTGACCACCTGCGGCAGGAACATTTTACCCGCGCCGAA  >  NC_002947/2712708‑2712979

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