Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,669,965 C→G 6.9% L148V (CTG→GTG)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,669,9650CG 6.9% 175.0 / 2.2 55L148V (CTG→GTG) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base C (22/29);  new base G (2/2);  total (24/31)
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

GCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTC  >  NC_002947/5669821‑5670111
                                                                                                                                                |                                                                                                                                                  
gCCGAGTACAAGGCTGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTgg                                                                                                                                                 <  5:261122/148‑1 (MQ=255)
gCCGAGTACAAGACAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTgg                                                                                                                                                 >  7:272493/1‑148 (MQ=255)
   gAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGt                                                                                                                                             <  7:87082/149‑1 (MQ=255)
      tACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGc                                                                                                                                          <  7:9700/149‑1 (MQ=255)
         aaGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGc                                                                                                                                       >  8:111573/1‑149 (MQ=255)
         aaGACAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGc                                                                                                                                       >  6:124737/1‑149 (MQ=255)
                    cATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGaaa                                                                                                                            <  2:14147/149‑1 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCTGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGc                                                                                                                         <  4:88165/148‑1 (MQ=255)
                        aaaCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCa                                                                                                                        >  2:300842/1‑149 (MQ=255)
                          aCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTCATCACTTCGACCTGGAAGTGGCTGCCTACCCGGAAATGCAcc                                                                                                                      >  1:67205/1‑149 (MQ=255)
                          aCGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATAGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCAcc                                                                                                                      >  4:96983/1‑149 (MQ=255)
                              aTTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTgg                                                                                                                                           <  4:180587/124‑1 (MQ=255)
                              aTTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTgg                                                                                                                                           >  3:180587/1‑124 (MQ=255)
                                              gtgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGa                                                                                                   <  5:206522/148‑1 (MQ=255)
                                              gtgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAg                                                                                                  >  1:195984/1‑149 (MQ=255)
                                              gtgACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAg                                                                                                  <  7:111573/149‑1 (MQ=255)
                                                   cTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTAATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGAACACGCAGGCGCGCAACTACGAGGCGGa                                                                                             >  5:337035/1‑149 (MQ=255)
                                                      ccGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGAt                                                                                            <  8:272493/147‑1 (MQ=255)
                                                       cGTCGGGCATGGTCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTg                                                                                         <  7:181668/149‑1 (MQ=255)
                                                         tCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTTACCTGGTCGAGTTCATCCGCGAGGTGAACGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCGGTAAATGCACTGTCAGGCGCGCAACTTCGTGGCGGATCTGGc                                                                                       <  2:263881/149‑1 (MQ=255)
                                                         tCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGc                                                                                       <  5:124737/149‑1 (MQ=255)
                                                                  ggcatggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTtgt                                                                              >  3:147735/1‑149 (MQ=255)
                                                                     atggCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGca                                                                           <  2:170301/149‑1 (MQ=255)
                                                                        gCCAGTGGTGAACTGCGCTATGCCAGTGACGTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGcac                                                                          >  1:31510/1‑147 (MQ=255)
                                                                        gCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACaa                                                                        <  1:43174/149‑1 (MQ=255)
                                                                         ccAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAg                                                                       <  3:253113/149‑1 (MQ=255)
                                                                          cAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAgg                                                                      >  1:325478/1‑149 (MQ=255)
                                                                                    cTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCgg                                                            >  8:68004/1‑149 (MQ=255)
                                                                                         cTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCg                                                       <  4:147735/149‑1 (MQ=255)
                                                                                           aTGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGAGCGCTGATGACTTCCACGTGGAAGTGGCTGCCTACCCGGATATGCACCCGCAGGCGCGCAACTTCGTGGCGGATCTGGGGAACTTTGTGGACAAGGTCAAGGCCGGt                                                           <  5:60887/143‑1 (MQ=255)
                                                                                           aTGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGc                                                               <  8:302471/139‑1 (MQ=255)
                                                                                           aTGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGc                                                               >  7:302471/1‑139 (MQ=255)
                                                                                           aTGCCAGTGACCTCGTCCAGATAATCCCCGACGAGACCGCTCATCACATCCACCACCAACACGCTCCCTACCCGCAAAACCACCCGCAGGCGCGCAACTTCGAGGCGCATCTGGCCAACTTTGTGCACAAGGACAAGGCCgga                                                           >  6:60887/1‑142 (MQ=255)
                                                                                               cAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAgcg                                                 >  7:63293/1‑149 (MQ=255)
                                                                                                    aCCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCAt                                             >  7:229488/1‑148 (MQ=255)
                                                                                                     ccTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCa                                           >  1:165408/1‑149 (MQ=255)
                                                                                                          tCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAg                                      <  3:96983/149‑1 (MQ=255)
                                                                                                            gAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGt                                     <  8:229488/148‑1 (MQ=255)
                                                                                                                tCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActt                                 <  2:106130/148‑1 (MQ=255)
                                                                                                                 cATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTActtc                                <  3:265878/148‑1 (MQ=255)
                                                                                                                    ccGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGcc                                                        <  8:169771/121‑1 (MQ=255)
                                                                                                                     cgcgAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCaa                           >  3:58019/1‑149 (MQ=255)
                                                                                                                          ggagACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCAcc                                         <  6:178957/130‑1 (MQ=255)
                                                                                                                           gagaCCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGTGCCATCACCCAGTACTTCTTCAACGCCGa                     <  1:300842/149‑1 (MQ=255)
                                                                                                                           gagaCCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGa                     <  2:31510/149‑1 (MQ=255)
                                                                                                                                cGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACGGCGCCATCACCCAGTACTTCTTCAACGCCGACAGc                 >  3:88215/1‑148 (MQ=255)
                                                                                                                                    agccACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctactt            >  4:273107/4‑149 (MQ=255)
                                                                                                                                      ccACTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACGAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttc           >  4:122531/2‑148 (MQ=255)
                                                                                                                                        aCTTCCAGGTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttctac        <  3:273107/149‑1 (MQ=255)
                                                                                                                                        aCTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttcta         <  2:165408/148‑1 (MQ=255)
                                                                                                                                         cTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAg                                                                       <  2:151164/85‑1 (MQ=255)
                                                                                                                                         cTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAg                                                                       >  1:151164/1‑85 (MQ=255)
                                                                                                                                         cTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCCCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGctacttctact       >  8:310096/1‑149 (MQ=255)
                                                                                                                                             cACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTtcgt   <  8:63293/149‑1 (MQ=255)
                                                                                                                                              aCCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTTCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTtcgtc  <  2:67205/149‑1 (MQ=255)
                                                                                                                                                cTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTtcgtc  >  1:49233/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GCCGAGTACAAGGCAGCCGGCATCAAACGCATTGTCGCGCTGCGTGGTGACCTGCCGTCGGGCATGGGCATGGCCAGTGGTGAACTGCGCTATGCCAGTGACCTGGTCGAGTTCATCCGCGAGGAGACCGCTGATCACTTCCACCTGGAAGTGGCTGCCTACCCGGAAATGCACCCGCAGGCGCGCAACTTCGAGGCGGATCTGGCCAACTTTGTGCACAAGGTCAAGGCCGGTGCCGACAGCGCCATCACCCAGTACTTCTTCAACGCCGACAGCTACTTCTACTTCGTC  >  NC_002947/5669821‑5670111

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: ‑