Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,759,985 G→A 9.3% L206L (CTG→CTA) ‡ folEB → GTP cyclohydrolase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,759,9850GA 9.3% 118.2 / 3.0 43L206L (CTG→CTA) ‡folEBGTP cyclohydrolase I
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (2/2);  total (21/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

GGTCAGGATCGGCTACTCCTCGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACGGGCT  >  NC_002947/3759838‑3760126
                                                                                                                                                   |                                                                                                                                             
ggTCAGGATCGGCTACTCCTCGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTgg                                                                                                                                              <  5:53290/149‑1 (MQ=255)
      gATCGGCTACTCCTCGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGc                                                                                                                                        <  3:86748/149‑1 (MQ=255)
        tCGGCTACTCCTCGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCca                                                                                                                                      <  2:144558/149‑1 (MQ=255)
               ctcctcGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGc                                                                                                                                >  4:194980/1‑148 (MQ=255)
                tcctcGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGcc                                                                                                                               >  8:198244/1‑148 (MQ=255)
                tcctcGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGcc                                                                                                                               >  4:16480/1‑148 (MQ=255)
                  ctcGACCTGTCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAAc                                                                                                                            <  7:120076/149‑1 (MQ=255)
                  ctcGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCCAGCGACGcgttgc                                                                                                                                   <  8:201588/142‑6 (MQ=255)
                  ctcGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCCAGCGACGcgttgc                                                                                                                                   >  7:201588/1‑137 (MQ=255)
                        cTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGc                                                                                                                       >  8:203815/1‑148 (MQ=255)
                          gCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCg                                                                                                                      >  8:14160/1‑147 (MQ=255)
                                gCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAGCCGTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCCGCAACGCGTCGCTGGCGCCGCCCCATAGCCAACGCAGCGCAGCTc                                                                                                                >  8:140151/1‑147 (MQ=255)
                                 cTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGa                                                                                                              <  8:891/148‑1 (MQ=255)
                                 cTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGa                                                                                                              <  7:140151/148‑1 (MQ=255)
                                              tcCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCAc                                                                                                   <  6:232155/146‑1 (MQ=255)
                                              tcCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACg                                                                                                  <  7:14160/147‑1 (MQ=255)
                                                  ggCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCt                                                                                            >  8:152721/1‑149 (MQ=255)
                                                                gACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGa                                                                              >  3:77454/1‑149 (MQ=255)
                                                                   ggCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGAtct                                                                           <  7:203815/149‑1 (MQ=255)
                                                                               aCCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGc                                                               >  1:103533/1‑149 (MQ=255)
                                                                                  ttCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCt                                                            >  7:36178/1‑149 (MQ=255)
                                                                                     aaGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGa                                                          <  7:198244/148‑1 (MQ=255)
                                                                                             gCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACg                                                 >  6:167035/1‑149 (MQ=255)
                                                                                              ccACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGccc                                                             <  5:277507/136‑1 (MQ=255)
                                                                                              ccACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGccc                                                             >  6:277507/1‑136 (MQ=255)
                                                                                              ccACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGt                                                >  5:223483/1‑149 (MQ=255)
                                                                                              ccACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGt                                                >  8:21333/1‑149 (MQ=255)
                                                                                                    cTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGaa                                            <  7:152721/147‑1 (MQ=255)
                                                                                                    cTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAgg                                          >  5:948/1‑149 (MQ=255)
                                                                                                         aGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCt                                     <  6:223483/149‑1 (MQ=255)
                                                                                                           tCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTg                                   >  1:87774/1‑149 (MQ=255)
                                                                                                               cgcCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCa                               >  8:104918/1‑149 (MQ=255)
                                                                                                                ggcGTCGGGCTGTGGCTGCGCAGCAACGCGTCGCTGGCGACGGCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCAc                              <  2:103533/147‑1 (MQ=255)
                                                                                                                gcCGTGGCGCAGTGGCTGCGCAGGAACGCGTCTCTGGCGACGCCAGATAGCCAACGCAGGGAAGCTCGAGTGAGCGTCACGGTACCTGCCGATGCGCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGGGCCTTGCGCAc                              <  4:77454/149‑1 (MQ=255)
                                                                                                                                      gcaACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTaaaa        <  3:16480/149‑1 (MQ=255)
                                                                                                                                       caACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAAc       >  8:190981/1‑149 (MQ=255)
                                                                                                                                        gtggcgTCGCCAGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACg      >  8:86278/4‑149 (MQ=255)
                                                                                                                                        aaCGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACg      <  5:139218/149‑1 (MQ=255)
                                                                                                                                         aCGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACgg     <  3:194980/149‑1 (MQ=255)
                                                                                                                                         aCGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACgg     <  5:167035/149‑1 (MQ=255)
                                                                                                                                           gcgTCGCCAGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACGGGc   <  7:86278/149‑1 (MQ=255)
                                                                                                                                            cgTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACGGGCt  <  1:211745/149‑1 (MQ=255)
                                                                                                                                              tCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGa                                                      >  3:111168/1‑95 (MQ=255)
                                                                                                                                              tCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGa                                                      <  4:111168/95‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
GGTCAGGATCGGCTACTCCTCGACCTGCCCGTGCTCGGCCGCTCTCTCCCGGCAGTTGATTGAGGACGGCTTCCTGCGAACCTTCAAGGACCAGCCACTGCTCGAAGTCAGCGCCGTCGCGCAGTGGCTGCGCAGCAACGCGTCGCTGGCGACGCCACATAGCCAACGCAGCGAAGCTCGAGTGCGCGTCACGGTACCTGCCGATGCCCCGGATCTCGGGTTGCTGGCCCTGATCGACCACGTCGAAGGCGCCTTGCGCACGCCGGTTCAGACCGCCGTAAAACGGGCT  >  NC_002947/3759838‑3760126

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