Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,256,859 C→T 5.8% R60H (CGT→CAT)  rpoC ← RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,256,8590CT 5.8% 231.3 / 2.1 68R60H (CGT→CAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (31/33);  new base T (2/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTAACTTC  >  minE/2256794‑2256929
                                                                 |                                                                      
cACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCTTAATCTTTTACCGGCCCAAAGATACGGgcg                                                                   <  1:1673729/71‑1 (MQ=255)
cACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGgcg                                                                   <  1:1521080/71‑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)
            ggCGCTTGTACTTGCCGCACAGGCACTCGTAGTCCTTCACCGGGCCGAAGATCTTGGCGCAGAACAGGCCg                                                       >  1:547731/1‑71 (MQ=255)
            ggCGCTTGTACTTACCGCACAGGCTCTCGTAATCTTTTACCGGCCCAAAGATAc                                                                        >  1:298424/1‑54 (MQ=255)
                  tGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGt                                                 >  1:1622404/1‑71 (MQ=255)
                   gTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCa                                                    >  1:258386/1‑67 (MQ=255)
                   gTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGtt                                                >  1:739846/1‑71 (MQ=255)
                   gTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGt                                                 >  1:3443924/1‑70 (MQ=255)
                     aCTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCAc                                                   >  1:1136026/1‑66 (MQ=255)
                     aCTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCAc                                                   >  1:3716180/1‑66 (MQ=255)
                       ttACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTgg                                            >  1:3562054/1‑71 (MQ=255)
                        tACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGt                                           >  1:2436464/1‑71 (MQ=255)
                          ccGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACTTTCTGGtt                                          <  1:2982129/70‑1 (MQ=255)
                          ccGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGtt                                          >  1:1429869/1‑70 (MQ=255)
                            gCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATAc                                                                        >  1:4226562/1‑38 (MQ=255)
                            gCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGtt                                                >  1:265228/1‑62 (MQ=255)
                            gCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGtt                                          >  1:2468014/1‑68 (MQ=255)
                             cacaGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGt                                           <  1:1637287/66‑1 (MQ=255)
                             cacaGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGt                                           <  1:445381/66‑1 (MQ=255)
                             cacaGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGaa                                      <  1:123186/71‑1 (MQ=255)
                                 ggCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTa                                  <  1:395793/71‑1 (MQ=255)
                                 ggCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTa                                  <  1:1505071/71‑1 (MQ=255)
                                  gCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTAc                                 <  1:3306271/71‑1 (MQ=255)
                                   cACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGcc                                                        <  1:1704308/47‑1 (MQ=255)
                                    aCTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACgg                               <  1:1487578/71‑1 (MQ=255)
                                    aCTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACg                                <  1:4181078/70‑1 (MQ=255)
                                      tCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTa                                  <  1:2777307/66‑1 (MQ=255)
                                       cGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGaaaa                                                            >  1:743604/1‑39 (MQ=255)
                                        gTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAg                            <  1:3275453/70‑1 (MQ=255)
                                        gTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAg                            <  1:1440229/70‑1 (MQ=255)
                                        gTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGt                           <  1:2014580/71‑1 (MQ=255)
                                         tAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACg                                <  1:3996170/65‑1 (MQ=255)
                                         tAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGtt                          <  1:2354691/71‑1 (MQ=255)
                                         tAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGtt                          <  1:3327321/71‑1 (MQ=255)
                                          aaTCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAg                            >  1:3482977/1‑68 (MQ=255)
                                           aTCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGcc                                                        >  1:1659029/1‑39 (MQ=255)
                                           aTCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGa                        >  1:1257734/1‑71 (MQ=255)
                                              ttttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGaa                                      >  1:2106230/1‑54 (MQ=255)
                                              ttttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGaa                                      >  1:4421492/1‑54 (MQ=255)
                                              ttttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATg                      >  1:936954/1‑70 (MQ=255)
                                               tttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCAc                                                   >  1:2845458/1‑40 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTGTGGTTTGAACGTACGGTa                             >  1:2938719/1‑61 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGaa                                      >  1:1694250/1‑52 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTa                             >  1:3225976/1‑61 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTa                             >  1:3060998/1‑61 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTa                             >  1:2447908/1‑61 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGAt                       <  1:1960950/67‑1 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGAt                       <  1:3284622/67‑1 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGAt                       <  1:4133918/67‑1 (MQ=255)
                                                ttACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGtt                   >  1:1802535/1‑71 (MQ=255)
                                                 tACCGGCCCAAAGATATGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGttt                  >  1:284008/1‑71 (MQ=255)
                                                 tACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGttt                  >  1:2535651/1‑71 (MQ=255)
                                                 tACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGttt                  >  1:2960344/1‑71 (MQ=255)
                                                 tACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGtt                   >  1:737328/1‑70 (MQ=255)
                                                  aCCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGtt                   <  1:320276/69‑1 (MQ=255)
                                                    cGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATg                      >  1:928874/1‑64 (MQ=255)
                                                         cAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCttt          <  1:2863439/71‑1 (MQ=255)
                                                           aaGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTa                                  <  1:4282311/45‑1 (MQ=255)
                                                           aaGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCttttt        <  1:4385838/71‑1 (MQ=255)
                                                            aGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCttttt        >  1:2844118/1‑70 (MQ=255)
                                                            aGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTa       <  1:4342694/71‑1 (MQ=255)
                                                             gATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTaa      <  1:38990/71‑1 (MQ=255)
                                                              aTACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTAAc     >  1:1584575/1‑71 (MQ=255)
                                                               tACTGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTAACt    <  1:146928/71‑1 (MQ=255)
                                                                 cGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACg                                    >  1:1839761/1‑37 (MQ=255)
                                                                 cGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTAACTTc  <  1:3720016/71‑1 (MQ=255)
                                                                 cGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTATCCGGCTTTTTAACTTc  <  1:299158/71‑1 (MQ=255)
                                                                 |                                                                      
CACGGTGTTTCAGGCGCTTGTACTTACCGCACAGGCACTCGTAATCTTTTACCGGCCCAAAGATACGGGCGCAGAAAAGGCCGTCACGTTCTGGTTTGAACGTACGGTAGTTGATGGTTTCCGGCTTTTTAACTTC  >  minE/2256794‑2256929

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: