Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,585,151 G→T 11.1% G129G (GGC→GGA)  folD‑II ← bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,585,1510GT11.1% 127.2 / 4.7 45G129G (GGC→GGA) folD‑IIbifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base T (2/3);  total (22/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.51e-01

ATGGCCATCGGGCGGCCAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCATTGAGACGGTCGATCAGCTC  >  NC_002947/2585008‑2585297
                                                                                                                                               |                                                                                                                                                  
aTGGCCATCGGGCGGCCAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGAtg                                                                                                                                              >  2:779151/1‑150 (MQ=255)
     cATCGGGCGGCCAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTa                                                                                                                                         <  1:316601/150‑1 (MQ=255)
         gggCGGCCAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTTCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCGGCAGTGGGATACGCTGGGCCAGGCGGCCGATGTTGTACggg                                                                                                                                     >  1:73829/1‑150 (MQ=255)
             ggcCAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGa                                                                                                                                 >  2:748644/1‑150 (MQ=255)
                cAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACgg                                                                                                                                      >  1:86310/1‑142 (MQ=255)
                cAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACgg                                                                                                                                      <  2:86310/142‑1 (MQ=255)
                    gATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCTGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTc                                                                                                                          <  2:1039376/150‑1 (MQ=255)
                     aTGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTcc                                                                                                                         <  2:223320/150‑1 (MQ=255)
                                   gacgaTCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCgg                                                                                                           >  1:733994/1‑150 (MQ=255)
                                    acgaTCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTCTCCAGCGGGCTGATGATGCCTTTCGGTGTGACAGCTTGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCggg                                                                                                          <  1:519700/150‑1 (MQ=255)
                                             gcgTTCATGCCGTACAGGTCCTGTCCGGTGCTTCCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATgcgc                                                                                                 <  2:252263/150‑1 (MQ=255)
                                             gcgTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATgcgc                                                                                                 >  2:112341/1‑150 (MQ=255)
                                              cgTTCTTGCCGTACAGGACCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCt                                                                                                <  2:340212/150‑1 (MQ=255)
                                                   aTGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGagc                                                                                           >  1:585793/1‑150 (MQ=255)
                                                        gTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCaa                                                                                      <  2:1033771/150‑1 (MQ=255)
                                                           cAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGa                                                                                   <  1:779151/150‑1 (MQ=255)
                                                             ggTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGTCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAgg                                                                                 >  2:995068/1‑150 (MQ=255)
                                                             ggTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGTCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAgg                                                                                 >  2:995036/1‑150 (MQ=255)
                                                              gTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGc                                                                                >  2:666562/1‑150 (MQ=255)
                                                                 cTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTc                                                                             >  1:534912/1‑150 (MQ=255)
                                                                    gCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGTCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAg                                                                          <  1:623913/150‑1 (MQ=255)
                                                                    gCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAg                                                                          <  2:585793/150‑1 (MQ=255)
                                                                    gCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAg                                                                          >  2:623913/1‑150 (MQ=255)
                                                                    gCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAg                                                                          <  2:220513/150‑1 (MQ=255)
                                                                    gCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATACTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAg                                                                          >  1:312945/1‑150 (MQ=255)
                                                                             ttCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTTTGCcgg                                                                 <  2:346704/150‑1 (MQ=255)
                                                                                 agcagGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGTCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCcggcagc                                                             <  1:995036/150‑1 (MQ=255)
                                                                                 agcagGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGTCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCcggcagc                                                             <  1:995068/150‑1 (MQ=255)
                                                                                       gTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGt                                                       >  2:885492/1‑150 (MQ=255)
                                                                                              tgCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGAAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACa                                                <  2:605057/150‑1 (MQ=255)
                                                                                                  tttCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGa                                                                                   <  1:298225/111‑1 (MQ=255)
                                                                                                  tttCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGa                                                                                   >  2:298225/1‑111 (MQ=255)
                                                                                                             aCGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCaa                                                  >  2:538793/1‑133 (MQ=255)
                                                                                                             aCGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCaa                                                  <  1:538793/133‑1 (MQ=255)
                                                                                                             aCGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCg                                        >  2:459476/1‑143 (MQ=255)
                                                                                                             aCGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCg                                        <  1:459476/143‑1 (MQ=255)
                                                                                                                  cgcAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAGGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCggg                            >  1:526726/1‑150 (MQ=255)
                                                                                                                    cagcagTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATACTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGgtc                          <  2:312945/150‑1 (MQ=255)
                                                                                                                          tGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCAtt                    <  1:419734/150‑1 (MQ=255)
                                                                                                                            gAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGagca                                                                                          >  2:288052/1‑78 (MQ=255)
                                                                                                                            gAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGagca                                                                                          <  1:288052/78‑1 (MQ=255)
                                                                                                                            gAATACGCTGGGCCAAGCGGCCGATGTTGTACGGGTGAAAGCCGTCCACATCCTTGTCCGAGCGAATGCTCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCATTga                  >  2:230004/1‑150 (MQ=255)
                                                                                                                                  gCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCATTGAGACGGt            <  2:118921/150‑1 (MQ=255)
                                                                                                                                  atagggCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTACAGGTGTGCCAGCAGCGGCAGTTGAAACAGGAAGACGTCAAGGGCGGGGTCGTGATTGAGACggg            >  1:588200/4‑149 (MQ=255)
                                                                                                                                        ccAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCATTGAGACGGTCGATCa      >  2:811761/1‑150 (MQ=255)
                                                                                                                                           ggcggcCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCATTGAGACGGTCGATCAGCt   <  2:733994/150‑1 (MQ=255)
                                                                                                                                            gcggcCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGAGCTCGAGCAGCAACGCGGCGTCCAGGTGTACCAGCAGCGGCAGTGGTAACATGATGCCGTCGACGGCCGGGTCGTCATTGAGACGGTCGATCAGCTc  >  2:401890/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ATGGCCATCGGGCGGCCAACGATGTTGGACGCGCCGACGATCACCGCGTTCATGCCGTACAGGTCCTGGCCGGTGCTTTCCAGCAGGGTCATGATGCCTTTCGGTGTGCACGGGCGCAGCAGTGGAATACGCTGGGCCAGGCGGCCGATGTTGTACGGGTGGAAGCCGTCCACATCCTTGTCCGGGCGAATGCGCTCGAGCAGCAACGAGGCGTCCAGGTGTGCCGGCAGCGGCAGTTGCAACAGGATGCCGTCGACGGCCGGGTCGTCATTGAGACGGTCGATCAGCTC  >  NC_002947/2585008‑2585297

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