Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,065,301 T→C 7.8% D514D (GAT→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,3010TC 7.8% 285.2 / 14.3 102D514D (GAT→GACrpoDRNA polymerase, sigma 70 (sigma D) factor
Reads supporting (aligned to +/- strand):  ref base T (58/36);  new base C (4/4);  total (62/40)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 115 total aligned reads displayed.

AAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACC  >  minE/2065231‑2065370
                                                                      |                                                                     
aaGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGtgatgat                                                                       >  1:506155/1‑71 (MQ=255)
aaGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGtgatgat                                                                       >  1:2890412/1‑71 (MQ=255)
 aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGTAAACGCCGATCGGtgatgatg                                                                      <  1:3166257/71‑1 (MQ=255)
 aGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGtgatgatg                                                                      <  1:2288219/71‑1 (MQ=255)
   aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                    <  1:1880114/71‑1 (MQ=255)
   aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                    <  1:562301/71‑1 (MQ=255)
   aCAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGaa                                                                    <  1:3557579/71‑1 (MQ=255)
     aaGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGtgatgatg                                                                      >  1:1888100/1‑67 (MQ=255)
     aaGATCCGCAAAGTGCTGAAAATCGCCAAAGAGCCAATCTCCAAGGAAACGCCGATCGGtgatgatg                                                                      >  1:3870036/1‑67 (MQ=255)
       gATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGAtt                                                                <  1:624232/71‑1 (MQ=255)
         tcagCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                              >  1:1060041/4‑71 (MQ=255)
         tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                              >  1:2728218/1‑71 (MQ=255)
         tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                              >  1:2656283/1‑71 (MQ=255)
         tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                              >  1:2300954/1‑71 (MQ=255)
         tCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCg                                                              >  1:1088184/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)
                     tGAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                   <  1:2835631/70‑1 (MQ=255)
                      gAAGATTGCGAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGa                                                   <  1:1103652/69‑1 (MQ=255)
                           tCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCt                                                         >  1:188129/1‑58 (MQ=255)
                           tCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCt                                                         >  1:194039/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)
                                aaaGAGCCAATCTCCATGGAAACTCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       <  1:1478647/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)
                                 aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:4015882/1‑70 (MQ=255)
                                 aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:1766298/1‑70 (MQ=255)
                                 aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                      <  1:206824/71‑1 (MQ=255)
                                 aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                      <  1:2413665/71‑1 (MQ=255)
                                 aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                      <  1:3433476/71‑1 (MQ=255)
                                 aaGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGAt                                      <  1:182618/71‑1 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggg                                                      >  1:3158147/1‑53 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTggg                                                      >  1:2621536/1‑53 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                         >  1:1900769/1‑66 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                         >  1:2325880/1‑66 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                         >  1:2764132/1‑66 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAg                                         >  1:2749768/1‑66 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                       >  1:3281185/1‑68 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                       >  1:4238216/1‑68 (MQ=255)
                                   gagCCAATCTCCATGGAAACGCCGATCGGTGACGACGACGATTCGCATCTGGGCGATTTTATCGAGGa                                       >  1:709266/1‑68 (MQ=255)
                                      ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:1928785/1‑65 (MQ=255)
                                      ccAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:2518560/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)
                                            tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                             >  1:1681168/1‑69 (MQ=255)
                                            tcCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                             >  1:1731627/1‑69 (MQ=255)
                                              cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                             >  1:1638111/1‑67 (MQ=255)
                                              cATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTc                             >  1:2465121/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)
                                                 ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:2832599/1‑54 (MQ=255)
                                                 ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:3523336/1‑54 (MQ=255)
                                                 ggAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGa                                       >  1:2218340/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: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)
                                                         cGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGAt              >  1:807155/1‑71 (MQ=255)
                                                            tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                          >  1:936429/1‑40 (MQ=255)
                                                            tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                          >  1:3741131/1‑40 (MQ=255)
                                                            tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                          >  1:3635406/1‑40 (MQ=255)
                                                            tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                          >  1:3348527/1‑40 (MQ=255)
                                                            tCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGa                                          >  1:2779735/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)
                                                                      |                                                                     
AAGACAAGATCCGCAAAGTGCTGAAGATCGCCAAAGAGCCAATCTCCATGGAAACGCCGATCGGTGATGATGAAGATTCGCATCTGGGGGATTTCATCGAGGATACCACCCTCGAGCTGCCGCTGGATTCTGCGACCACC  >  minE/2065231‑2065370

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: