Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,065,319 G→C 8.2% G520G (GGG→GGC rpoD → RNA polymerase, sigma 70 (sigma D) factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,065,3190GC 8.2% 281.5 / 12.1 97G520G (GGG→GGCrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base G (58/31);  new base C (4/4);  total (62/35)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 102 total aligned reads displayed.

TGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGC  >  minE/2065249‑2065387
                                                                      |                                                                    
tGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTgggg                                                                      <  1:794942/71‑1 (MQ=255)
  atGAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                                    <  1:2734669/70‑1 (MQ=255)
  atGAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                                    <  1:539733/70‑1 (MQ=255)
   tGAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                                    <  1:2835631/70‑1 (MQ=255)
    gAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                                    <  1:1103652/69‑1 (MQ=255)
       gATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCa                                                               >  1:1906674/1‑71 (MQ=255)
          cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                                             >  1:1376201/1‑70 (MQ=255)
          cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                                             >  1:4160679/1‑70 (MQ=255)
          cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGATTTCATc                                                             >  1:2489449/1‑69 (MQ=255)
           gCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGGAGATTCGCATCTGGGGGATTTCATCg                                                            <  1:2342018/70‑1 (MQ=255)
           gCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCg                                                            <  1:2138166/70‑1 (MQ=255)
            ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          <  1:258792/71‑1 (MQ=255)
              aaaGAGCCAATCTCCATGGAAACTCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        <  1:1478647/71‑1 (MQ=255)
              aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        <  1:597021/71‑1 (MQ=255)
              aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        <  1:4061185/71‑1 (MQ=255)
              aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        <  1:2056829/71‑1 (MQ=255)
              aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        <  1:187460/71‑1 (MQ=255)
               aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:4015882/1‑70 (MQ=255)
               aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:1766298/1‑70 (MQ=255)
               aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                                       <  1:2413665/71‑1 (MQ=255)
               aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                                       <  1:3433476/71‑1 (MQ=255)
               aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                                       <  1:206824/71‑1 (MQ=255)
               aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                                       <  1:182618/71‑1 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          >  1:1900769/1‑66 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          >  1:2749768/1‑66 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          >  1:2325880/1‑66 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                                        >  1:709266/1‑68 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                                        >  1:3281185/1‑68 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                                        >  1:4238216/1‑68 (MQ=255)
                 gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                                        >  1:1687634/1‑68 (MQ=255)
                    ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:551090/1‑65 (MQ=255)
                    ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:2518560/1‑65 (MQ=255)
                    ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:1928785/1‑65 (MQ=255)
                       aTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                                             <  1:2501722/57‑1 (MQ=255)
                        tctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggggg                                                                     <  1:606987/48‑1 (MQ=255)
                        tctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:1569285/1‑58 (MQ=255)
                         ctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCAccc                                                >  1:1045689/1‑68 (MQ=255)
                         ctcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                                              >  1:1165973/1‑70 (MQ=255)
                          tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTc                                                                <  1:940158/51‑1 (MQ=255)
                          tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          >  1:2878610/1‑57 (MQ=255)
                          tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          >  1:775836/1‑57 (MQ=255)
                          tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          >  1:2152072/1‑57 (MQ=255)
                          tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                                              >  1:1731627/1‑69 (MQ=255)
                          tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                                              >  1:1681168/1‑69 (MQ=255)
                            cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                                              >  1:2465121/1‑67 (MQ=255)
                            cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                                              >  1:2182538/1‑67 (MQ=255)
                            cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                                              >  1:1638111/1‑67 (MQ=255)
                            cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAg                                           <  1:2795790/70‑1 (MQ=255)
                            cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGc                                          >  1:687364/1‑71 (MQ=255)
                             aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCg                                                            <  1:3701678/52‑1 (MQ=255)
                             aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCt                                         <  1:3657052/71‑1 (MQ=255)
                             aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCt                                         <  1:2922867/71‑1 (MQ=255)
                             aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCt                                         <  1:345118/71‑1 (MQ=255)
                               ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:2832599/1‑54 (MQ=255)
                               ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:2539986/1‑54 (MQ=255)
                               ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:2218340/1‑54 (MQ=255)
                               ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:3523336/1‑54 (MQ=255)
                               ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:1751771/1‑54 (MQ=255)
                                gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGcc                                      >  1:2476499/1‑71 (MQ=255)
                                gAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGc                                       >  1:1920338/1‑70 (MQ=255)
                                  aaCGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGc                                    <  1:2930853/71‑1 (MQ=255)
                                   aCGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCAccc                                                <  1:1770446/58‑1 (MQ=255)
                                    cGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTg                                  <  1:692985/71‑1 (MQ=255)
                                     gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                                          <  1:2488736/46‑1 (MQ=255)
                                     gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:2837606/1‑48 (MQ=255)
                                     gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:131819/1‑48 (MQ=255)
                                       cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:187037/1‑46 (MQ=255)
                                       cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:2610645/1‑46 (MQ=255)
                                       cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATaccac                                                  >  1:1242852/1‑52 (MQ=255)
                                       cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGAt                               >  1:807155/1‑71 (MQ=255)
                                          tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:936429/1‑40 (MQ=255)
                                          tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:3741131/1‑40 (MQ=255)
                                          tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:2779735/1‑40 (MQ=255)
                                          tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:3348527/1‑40 (MQ=255)
                                          tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:3178993/1‑40 (MQ=255)
                                          tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                                           >  1:3635406/1‑40 (MQ=255)
                                           cGGTGATGATGAAGATTCGCTTCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg                           >  1:391621/1‑71 (MQ=255)
                                           cGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg                           >  1:1543789/1‑71 (MQ=255)
                                           cGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg                           >  1:4055775/1‑71 (MQ=255)
                                              tgatgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTg                                  >  1:4172249/1‑61 (MQ=255)
                                                atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:1852795/1‑37 (MQ=255)
                                                atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGa                                                        >  1:3087682/1‑37 (MQ=255)
                                                atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATacc                                                    <  1:3098027/41‑1 (MQ=255)
                                                atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATacc                                                    <  1:3759643/41‑1 (MQ=255)
                                                atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGc                                    <  1:4213970/57‑1 (MQ=255)
                                                   atgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCAACCTCGAGCTGCCGCTGGATTATGCGcccacc                   >  1:1085443/1‑71 (MQ=255)
                                                     gaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCg                  >  1:2760360/1‑70 (MQ=255)
                                                        taTTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaaa               <  1:3013004/69‑1 (MQ=255)
                                                        gATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAg              <  1:3175188/71‑1 (MQ=255)
                                                        gATTCACATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAg              <  1:821671/71‑1 (MQ=255)
                                                            cGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtg          >  1:3078644/1‑71 (MQ=255)
                                                                 cTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgc     >  1:4130535/1‑71 (MQ=255)
                                                                  tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaaa               <  1:2886222/60‑1 (MQ=255)
                                                                  tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGaaa               <  1:1083180/60‑1 (MQ=255)
                                                                  tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgcg    >  1:2504290/1‑71 (MQ=255)
                                                                  tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgcg    >  1:1795159/1‑71 (MQ=255)
                                                                  tGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCtgcgtgcg    >  1:3295723/1‑71 (MQ=255)
                                                                    ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc             >  1:3133233/1‑60 (MQ=255)
                                                                    ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGc             >  1:173167/1‑60 (MQ=255)
                                                                    ggggATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGc  <  1:4214128/71‑1 (MQ=255)
                                                                      |                                                                    
TGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCGAAAGCCTGCGTGCGGC  >  minE/2065249‑2065387

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