Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,256,801 T→C 5.4% K79K (AAA→AAG rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,256,8010TC 5.4% 212.7 / 3.8 74K79K (AAA→AAGrpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (40/30);  new base C (2/2);  total (42/32)
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.70e-01

CATACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAA  >  minE/2256732‑2256871
                                                                     |                                                                      
cATACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGtt                                                                       >  1:1164304/1‑71 (MQ=255)
cATACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGtt                                                                       >  1:2926021/1‑71 (MQ=255)
  tACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGttt                                                                      >  1:2054233/1‑70 (MQ=255)
  tACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGttt                                                                      >  1:3580205/1‑70 (MQ=255)
  tACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGttt                                                                      >  1:4248753/1‑70 (MQ=255)
     gCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGttt                                                                      <  1:1548305/67‑1 (MQ=255)
       tCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGgc                                                                 >  1:4181056/1‑70 (MQ=255)
         aCGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGgcgc                                                               <  1:4425974/70‑1 (MQ=255)
         aCGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCt                                                              >  1:2324267/1‑71 (MQ=255)
         aCGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCt                                                              >  1:3452429/1‑71 (MQ=255)
         aCGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCt                                                              >  1:518503/1‑71 (MQ=255)
              gTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTAc                                                         <  1:2703991/71‑1 (MQ=255)
               tACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTg                                                            <  1:1798711/67‑1 (MQ=255)
                 cTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                      >  1:956842/1‑71 (MQ=255)
                  tttAGTCTGGGTCACTTTAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCa                                                                    >  1:4294430/1‑56 (MQ=255)
                  tttAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCa                                                                    >  1:3968336/1‑56 (MQ=255)
                  tttAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCa                                                                    >  1:2420958/1‑56 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTttacttac                                                     <  1:2455670/69‑1 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                      >  1:1893311/1‑68 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                      >  1:3504351/1‑68 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTa                                                      >  1:2355965/1‑68 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTAc                                                     <  1:2859699/69‑1 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCg                                                   >  1:2964687/1‑71 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCg                                                   >  1:833099/1‑71 (MQ=255)
                    tAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACATACCg                                                   >  1:1261062/1‑71 (MQ=255)
                      gTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                 <  1:4183438/71‑1 (MQ=255)
                       tCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGcac                                                <  1:2688547/71‑1 (MQ=255)
                         tGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAg                                              >  1:3574906/1‑71 (MQ=255)
                           ggTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAg                                              <  1:1342881/69‑1 (MQ=255)
                              cACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACt                                         >  1:465325/1‑71 (MQ=255)
                              cACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACt                                         >  1:2773609/1‑71 (MQ=255)
                               aCTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTc                                        >  1:4043579/1‑71 (MQ=255)
                                cTTCGACGCCGCACTTTTCGCAGATCACGCCGCGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCg                                       >  1:4149514/1‑71 (MQ=255)
                                cTTCGACGCCGCACTTTTCGCAGATCACGCCGCGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCg                                       >  1:3602722/1‑71 (MQ=255)
                                 ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACt                                                        <  1:3063959/53‑1 (MQ=255)
                                 ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                 <  1:2852610/60‑1 (MQ=255)
                                 ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                 <  1:4205750/60‑1 (MQ=255)
                                 ttCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGca                                                 <  1:292678/60‑1 (MQ=255)
                                   cAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACt                                                        >  1:3805605/1‑51 (MQ=255)
                                   cAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACt                                                        >  1:2822008/1‑51 (MQ=255)
                                   cAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACt                                                        >  1:2084539/1‑51 (MQ=255)
                                   aaaCGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCTCTTGTACt                                                        >  1:4313246/2‑51 (MQ=255)
                                       gccgcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCtt                                <  1:3641170/71‑1 (MQ=255)
                                          gcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             <  1:267855/71‑1 (MQ=255)
                                          gcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             <  1:2852369/71‑1 (MQ=255)
                                          gcACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             <  1:721655/71‑1 (MQ=255)
                                           cACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTAc                            <  1:155797/71‑1 (MQ=255)
                                            aCTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTAc                                                         >  1:3581393/1‑41 (MQ=255)
                                            aCTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTAc                                                         >  1:4218799/1‑41 (MQ=255)
                                            aCTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTAc                                                         >  1:3431304/1‑41 (MQ=255)
                                            aCTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTAc                                                         >  1:1675775/1‑41 (MQ=255)
                                            aCTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTAcc                           <  1:3956610/71‑1 (MQ=255)
                                             cTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCttt                               <  1:4324115/66‑1 (MQ=255)
                                              ttCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGcac                                                <  1:2212390/48‑1 (MQ=255)
                                               tctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGc                                                  >  1:1482223/1‑45 (MQ=255)
                                               tctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             >  1:3535335/1‑66 (MQ=255)
                                               tctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             >  1:634881/1‑66 (MQ=255)
                                                ctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCgga                        >  1:2215632/1‑69 (MQ=255)
                                                ctcACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCgg                         <  1:3844695/69‑1 (MQ=255)
                                                   acaGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTAc                            >  1:2744608/1‑63 (MQ=255)
                                                   acaGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                     >  1:1594818/1‑70 (MQ=255)
                                                    caGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGccc                      <  1:235584/68‑1 (MQ=255)
                                                    caGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCaaa                   <  1:4395260/71‑1 (MQ=255)
                                                     aGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             >  1:2945820/1‑60 (MQ=255)
                                                     aGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             >  1:240790/1‑60 (MQ=255)
                                                     aGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCaaa                   <  1:1589896/70‑1 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTAc                                                     >  1:3654770/1‑35 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                     >  1:39560/1‑67 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                     >  1:2587188/1‑67 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCa                     >  1:2564452/1‑67 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                 <  1:1145780/71‑1 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                 <  1:2046605/71‑1 (MQ=255)
                                                      gATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGa                 <  1:2041196/71‑1 (MQ=255)
                                                       aTGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             >  1:963957/1‑58 (MQ=255)
                                                       aTGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTa                             >  1:2376058/1‑58 (MQ=255)
                                                              cACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCTTAATCTTTTACCGGCCCAAAGATACGGgcg         <  1:1673729/71‑1 (MQ=255)
                                                              cACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGgcg         <  1:1521080/71‑1 (MQ=255)
                                                               aCGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTa                                     <  1:3575386/42‑1 (MQ=255)
                                                                cGGTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCGTAGTCCTTGATCGGGCCAAAGATCTTGGCGCa       <  1:1482452/71‑1 (MQ=255)
                                                                 ggTGCTTCAGGCGCTTGTACTTGCCGCACAGGCATTCGTAGTCCTTGATCGGGCCAAAGATCTTGGCGCa       <  1:3895856/70‑1 (MQ=255)
                                                                   tgtTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaa    <  1:891174/71‑1 (MQ=255)
                                                                   tgtTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaa    <  1:3884487/71‑1 (MQ=255)
                                                                     tttCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaaaa  <  1:2673864/71‑1 (MQ=255)
                                                                     tttCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaaa   <  1:4409524/70‑1 (MQ=255)
                                                                     |                                                                      
CATACGCTCACGGCGTACTTTAGTCTGGGTCACTTCAACGCCGCACTTCTCACAGATGACGCCACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAA  >  minE/2256732‑2256871

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: