Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,995,269 C→A 69.9% A170E (GCG→GAG)  dapF → diaminopimelate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,995,2690CA69.9% 38.0 / 32.7 46A170E (GCG→GAG) dapFdiaminopimelate epimerase
Reads supporting (aligned to +/- strand):  ref base C (8/6);  new base A (16/16);  total (24/22)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.84e-01

CTTCGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGC  >  NC_000913/3995126‑3995411
                                                                                                                                               |                                                                                                                                              
cttggcaCCTTCCGCCGTGCCGTTTCGGGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGAAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGGCGATGATGTCGATACCGaggcgg                                                                                                                                           <  2:323948/143‑1 (MQ=255)
      aCCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGaaa                                                                                                                                     <  1:338987/149‑1 (MQ=255)
       ccTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGGCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAAc                                                                                                                                    <  1:162720/149‑1 (MQ=255)
            cgccgTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTg                                                                                                                               <  1:79608/149‑1 (MQ=255)
               cgtgccgtTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTc                                                                                                                            <  1:218844/149‑1 (MQ=255)
               cgtgccgtTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTc                                                                                                                            <  1:308849/149‑1 (MQ=255)
                      tttCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCt                                                                                                                     >  1:252012/1‑149 (MQ=255)
                      ttccgcgCTAACAAAGCGGAAAAGCCCTATAGTATGCGCGCCCCCGAGCAGACAACCTTATTCGGCCTGTTTTCCACAGCAAATCCTCATTGCACTATTCACGTCGATGATGTCGATACCCAGGCGGTAGAAACGCTTGGTCCTGTTCt                                                                                                                     <  2:354711/146‑1 (MQ=255)
                       tccgcgCTACTAGAGCGGAAAAGACCTATAGTAGGCGCGCCGCCGAGCAGACAAGCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTg                                                                                                                    <  1:375587/147‑1 (MQ=255)
                          gcgcTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCGGTTCAGGaa                                                                                                                 >  2:533029/1‑149 (MQ=255)
                          gcgcTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGCCGATATTGTTGAAGCCGAGGCGGTAGAAACTCTTGGTCCTGTTCTGGaa                                                                                                                 >  1:681615/1‑149 (MQ=255)
                                  aaaGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGa                                                                                                         >  1:396175/1‑149 (MQ=255)
                                  aaaGCGGAAAAGACCTATATTATGCGCGCCGCCGAGAAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATGGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGa                                                                                                         >  2:611349/1‑149 (MQ=255)
                                   aaGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAg                                                                                                        >  1:667202/1‑149 (MQ=255)
                                              aCCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGa                                                                                             >  1:145955/1‑149 (MQ=255)
                                               ccTATATTATGCGGGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAg                                                                                            <  1:686778/149‑1 (MQ=255)
                                                       acgctcgccgcCCAGCAGACCATCTTATGCGGCGTGGTGTCGATGGGAAATCCTCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCaa                                                                                    <  1:256928/144‑1 (MQ=255)
                                                           gcgcCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAAAGCCTGGTCCTGTTCTGGAAACCCACGAGCGTGTTCCGGGGCGCGCCCATATc                                                                                >  2:101518/1‑149 (MQ=255)
                                                             gccgccGCGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCgg                                                                              <  2:561402/149‑1 (MQ=255)
                                                                       aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCaa                                                                    >  1:692482/1‑149 (MQ=255)
                                                                       aGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGCTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTCGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCTATATCGGTTTTATGCaa                                                                    <  2:541512/149‑1 (MQ=255)
                                                                        gACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAg                                                                   >  2:636675/1‑149 (MQ=255)
                                                                           actCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGGGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTgg                                                                <  2:63014/147‑1 (MQ=255)
                                                                           aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTgg                                                                >  1:94477/1‑149 (MQ=255)
                                                                           aaTCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAACATCGGTTTTATGCAAgggg                                                                >  2:441531/1‑149 (MQ=255)
                                                                           aaTCTTATGCGGCGTGGTGTCGATGGGAAATCAGCCACCCTGCATGCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGCAAGCCACGAGCGTTTTCCGGAGCGCGCCAAAGTCGGTTTTATGCCAGTgg                                                                >  2:486379/1‑149 (MQ=255)
                                                                               ttATTCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                                                            <  1:710265/149‑1 (MQ=255)
                                                                               ttATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTaa                                                            <  1:393250/149‑1 (MQ=255)
                                                                                 aTGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTGTATGCAAGTGGTTAAgc                                                          >  2:195154/1‑149 (MQ=255)
                                                                                     ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa                                                      <  1:378071/149‑1 (MQ=255)
                                                                                     ggcgTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGa                                                      >  1:13679/1‑149 (MQ=255)
                                                                                               cGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGt                                            <  2:426226/149‑1 (MQ=255)
                                                                                                      aaaTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAAGGCATGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTAcgcg                                     >  1:143741/1‑149 (MQ=255)
                                                                                                        aTCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGtt                                   >  2:693746/1‑149 (MQ=255)
                                                                                                         tCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGttt                                  >  1:64886/1‑149 (MQ=255)
                                                                                                         tCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGGTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGttt                                  >  2:457760/1‑149 (MQ=255)
                                                                                                           cGCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGGGTTTAt                                >  1:553099/1‑149 (MQ=255)
                                                                                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATg                               >  1:22406/1‑149 (MQ=255)
                                                                                                            gCATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAACATCGGGTTTATGCAAGTGGTTAAGCGCGACCATACTCGTTTAAGCGTTTATg                               >  2:164002/1‑149 (MQ=255)
                                                                                                             cATTGCGTGATTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCCAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGa                              >  2:215245/1‑149 (MQ=255)
                                                                                                                 gCGTGAGTCAGGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGt                          <  2:341645/149‑1 (MQ=255)
                                                                                                                  cGTGATTCATGTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTg                         <  2:323886/149‑1 (MQ=255)
                                                                                                                     gATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATTCATTTGTTTAAGCGAGAGCAGATTCGGTTACGCGTTTATGGGCGTgggg                      >  2:700151/1‑149 (MQ=255)
                                                                                                                       tcccggTCCATGCTGGCGGTACTGAGGCGGTAGAAACGCCTGGGCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGGGCGCGCCAATATCGGTTTTCTGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCa                    <  2:608480/145‑1 (MQ=255)
                                                                                                                            gTCGATGATGTCGATACCGAGGCGGTAGAAACGCTTGGTCCTGTTCGGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGgaga               <  2:711568/149‑1 (MQ=255)
                                                                                                                                   atgTCGATCACGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAg        <  1:486379/149‑1 (MQ=255)
                                                                                                                                         aTACCGAGGCGGTAGAAACGCGTGGTCCTGTTCTGGACAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGc  <  2:116817/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CTTCGAACCTTCCGCCGTGCCGTTTCGCGCTAACAAAGCGGAAAAGACCTATATTATGCGCGCCGCCGAGCAGACAATCTTATGCGGCGTGGTGTCGATGGGAAATCCGCATTGCGTGATTCAGGTCGATGATGTCGATACCGCGGCGGTAGAAACGCTTGGTCCTGTTCTGGAAAGCCACGAGCGTTTTCCGGAGCGCGCCAATATCGGTTTTATGCAAGTGGTTAAGCGCGAGCATATTCGTTTACGCGTTTATGAGCGTGGGGCAGGAGAAACCCAGGCCTGC  >  NC_000913/3995126‑3995411

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