Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,065,304 A→C 7.6% E515D (GAA→GAC 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,3040AC 7.6% 253.6 / 20.5 106E515D (GAA→GACrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base A (63/35);  new base C (4/4);  total (67/39)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 110 total aligned reads displayed.

ACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCG  >  minE/2065234‑2065371
                                                                      |                                                                   
aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                     <  1:2755709/71‑1 (MQ=255)
aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                     <  1:1880114/71‑1 (MQ=255)
aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                     <  1:3557579/71‑1 (MQ=255)
aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                     <  1:562301/71‑1 (MQ=255)
   aGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAt                                                                  <  1:722895/71‑1 (MQ=255)
      tcagCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                               >  1:1060041/4‑71 (MQ=255)
      tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                               >  1:1088184/1‑71 (MQ=255)
      tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                               >  1:2728218/1‑71 (MQ=255)
      tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                               >  1:2300954/1‑71 (MQ=255)
      tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                               >  1:2656283/1‑71 (MQ=255)
           aaaGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCt                                                          <  1:2307002/71‑1 (MQ=255)
           aaaGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCt                                                          <  1:218585/71‑1 (MQ=255)
               tGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTgggg                                                      <  1:794942/71‑1 (MQ=255)
                 atGAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                    <  1:2734669/70‑1 (MQ=255)
                 atGAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                    <  1:539733/70‑1 (MQ=255)
                   gAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                    <  1:1103652/69‑1 (MQ=255)
                      gATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCa                                               >  1:1906674/1‑71 (MQ=255)
                        tCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCt                                                          >  1:188129/1‑58 (MQ=255)
                         cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                             >  1:4160679/1‑70 (MQ=255)
                         cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                             >  1:1376201/1‑70 (MQ=255)
                         cGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGATTTCATc                                             >  1:2489449/1‑69 (MQ=255)
                          gCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCg                                            <  1:2138166/70‑1 (MQ=255)
                           ccAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          <  1:258792/71‑1 (MQ=255)
                             aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        <  1:4061185/71‑1 (MQ=255)
                             aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        <  1:187460/71‑1 (MQ=255)
                             aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        <  1:597021/71‑1 (MQ=255)
                             aaaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        <  1:2056829/71‑1 (MQ=255)
                              aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTg                                                         <  1:1005285/53‑1 (MQ=255)
                              aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTg                                                         <  1:1843075/53‑1 (MQ=255)
                              aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:1766298/1‑70 (MQ=255)
                              aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:4015882/1‑70 (MQ=255)
                              aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                       <  1:206824/71‑1 (MQ=255)
                              aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                       <  1:2413665/71‑1 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggg                                                       >  1:3158147/1‑53 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          >  1:2325880/1‑66 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          >  1:2749768/1‑66 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          >  1:2764132/1‑66 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          >  1:1900769/1‑66 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                        >  1:1687634/1‑68 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                        >  1:4238216/1‑68 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                        >  1:3281185/1‑68 (MQ=255)
                                gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                        >  1:709266/1‑68 (MQ=255)
                                   ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:2518560/1‑65 (MQ=255)
                                   ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:551090/1‑65 (MQ=255)
                                   ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:1928785/1‑65 (MQ=255)
                                      aTCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATc                                             <  1:2501722/57‑1 (MQ=255)
                                       tctcCATGGAAACGCCGATCGGTGATGATGAAGAt                                                                  <  1:766894/35‑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:775836/1‑57 (MQ=255)
                                         tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          >  1:2152072/1‑57 (MQ=255)
                                         tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                          >  1:2878610/1‑57 (MQ=255)
                                         tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                              >  1:1681168/1‑69 (MQ=255)
                                         tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                              >  1:1731627/1‑69 (MQ=255)
                                           cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                              >  1:2465121/1‑67 (MQ=255)
                                           cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                              >  1:1638111/1‑67 (MQ=255)
                                           cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                              >  1:2182538/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:2922867/71‑1 (MQ=255)
                                            aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCt                         <  1:345118/71‑1 (MQ=255)
                                            aTGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCt                         <  1:3657052/71‑1 (MQ=255)
                                              ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:2539986/1‑54 (MQ=255)
                                              ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:2218340/1‑54 (MQ=255)
                                              ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:2832599/1‑54 (MQ=255)
                                              ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:3523336/1‑54 (MQ=255)
                                              ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:1751771/1‑54 (MQ=255)
                                               gAAACGCCGATCGGTGATGATGAAGATTCGCATCt                                                          >  1:1690168/1‑35 (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:131819/1‑48 (MQ=255)
                                                    gCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:2837606/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:3348527/1‑40 (MQ=255)
                                                         tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                           >  1:3178993/1‑40 (MQ=255)
                                                         tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                           >  1:2779735/1‑40 (MQ=255)
                                                         tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                           >  1:3635406/1‑40 (MQ=255)
                                                         tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                           >  1:936429/1‑40 (MQ=255)
                                                         tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                           >  1:3741131/1‑40 (MQ=255)
                                                          cGGTGATGATGAAGATTCGCTTCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg           >  1:391621/1‑71 (MQ=255)
                                                          cGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTg           >  1:1543789/1‑71 (MQ=255)
                                                             tgatgatgaAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTg                  >  1:4172249/1‑61 (MQ=255)
                                                               atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:3087682/1‑37 (MQ=255)
                                                               atgatgaAGATTCGCATCTGGGGGATTTCATCGAGGa                                        >  1:1852795/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)
                                                                      |                                                                   
ACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACCG  >  minE/2065234‑2065371

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