Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,384 T→A 14.3% intergenic (+152/‑1627) yibH → / → ysaB hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,3840TA14.3% 174.5 / 21.2 70intergenic (+152/‑1627)yibH/ysaBhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/40);  new base A (4/6);  total (24/46)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

CTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGG  >  minE/2571317‑2571453
                                                                    |                                                                     
cTTTAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTACCCGTCTTGATc                                                                     <  1:681628/71‑1 (MQ=255)
atttAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATc                                                                     >  1:1204721/2‑71 (MQ=255)
  ttAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTg                                                                   >  1:1385251/1‑71 (MQ=255)
   tAATATCCAGCATGTTAAGAAGTTCAGGGATAGCTGTGATCTCATTGGATTTCTTATCCGTCTTGATCTgt                                                                  >  1:3859165/1‑70 (MQ=255)
      tatTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTAACCGTCTTGTTCTGCCGc                                                               <  1:1658843/70‑1 (MQ=255)
           aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                          <  1:2701468/70‑1 (MQ=255)
           aTTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGAt                                                          <  1:3677991/70‑1 (MQ=255)
            ttATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                        <  1:2009661/71‑1 (MQ=255)
             tATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTc                                                                     >  1:779794/1‑57 (MQ=255)
             tATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGa                                                           >  1:2418661/1‑67 (MQ=255)
              aTGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATgag                                                       <  1:3209317/70‑1 (MQ=255)
               tGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGAc                                                     >  1:355609/1‑71 (MQ=255)
                 ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                                    <  1:3269259/70‑1 (MQ=255)
                 ttAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                                    <  1:209001/70‑1 (MQ=255)
                   aataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtat                                                  <  1:257301/70‑1 (MQ=255)
                    ataaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtata                                                 <  1:2815788/70‑1 (MQ=255)
                     taaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACa                                               >  1:13632/1‑71 (MQ=255)
                      aaGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAt                                              >  1:3522116/1‑71 (MQ=255)
                        gTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcg                                            <  1:1944225/71‑1 (MQ=255)
                             tGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGa                                                           >  1:506791/1‑52 (MQ=255)
                               gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                        <  1:3675681/53‑1 (MQ=255)
                               gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACtat                                                  >  1:3271316/1‑59 (MQ=255)
                               gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCgag                                       <  1:2084695/70‑1 (MQ=255)
                               gATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCga                                        <  1:3808151/69‑1 (MQ=255)
                                  gAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                   >  1:869814/1‑71 (MQ=255)
                                     tgtgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATga                                                        <  1:2270/47‑1 (MQ=255)
                                       tgCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAat                               >  1:3041879/1‑70 (MQ=255)
                                         ctctCATTGGATTTTTTATCCGTCTTGTTCTGCCGCa                                                              <  1:3611129/37‑1 (MQ=255)
                                          tctcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGa                                                           >  1:3002213/1‑39 (MQ=255)
                                            tcATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATgaga                                                      >  1:2672705/1‑42 (MQ=255)
                                             cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                                    <  1:1372972/43‑1 (MQ=255)
                                             cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACt                                                    <  1:254577/43‑1 (MQ=255)
                                             cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCaa                                >  1:3461047/1‑63 (MQ=255)
                                             cATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                        >  1:29333/1‑71 (MQ=255)
                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAAtgatga                          <  1:2517521/68‑1 (MQ=255)
                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                       <  1:3512756/71‑1 (MQ=255)
                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                       <  1:3235391/71‑1 (MQ=255)
                                              aTTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGata                        <  1:1621088/70‑1 (MQ=255)
                                               ttGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATgatgat                         <  1:2860854/68‑1 (MQ=255)
                                               ttGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGacat                       >  1:2080036/1‑70 (MQ=255)
                                               ttGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATg                      <  1:1646094/71‑1 (MQ=255)
                                               ttGGATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATg                      <  1:260330/71‑1 (MQ=255)
                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                          <  1:3469100/50‑1 (MQ=255)
                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                          <  1:3769408/50‑1 (MQ=255)
                                                tGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGatat                       >  1:1556753/1‑69 (MQ=255)
                                                 ggATTTCTTATCCGTCTTGATCTGTCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                     <  1:1184069/70‑1 (MQ=255)
                                                 ggATTTCTTATCCGTCTTGATCTGCCCGATGACGAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                     <  1:2467650/70‑1 (MQ=255)
                                                 ggATTTCTTATCCGTCTTGATCTGCCCGATGACCAGACTGTGCATTGTTGAGAACGCACTAATGACATGa                     <  1:1156558/70‑1 (MQ=255)
                                                  gATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgt                   <  1:1932645/71‑1 (MQ=255)
                                                  gATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAg                    <  1:1563958/70‑1 (MQ=255)
                                                    ttttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttg                 <  1:670240/71‑1 (MQ=255)
                                                    ttttttATCCGTCTTGTTCTGCCGCAGGATCAGACTATACATCGTCGAGATCGCAATGATGATAAGAgttg                 >  1:2667978/1‑71 (MQ=255)
                                                     tttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttg                 <  1:3012631/70‑1 (MQ=255)
                                                      ttttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAgttgtt               >  1:2647521/1‑71 (MQ=255)
                                                        ttATCCGTCTTGTTCTGCCGCAGGATGAGACTATACAtcgtc                                          >  1:1731894/1‑42 (MQ=255)
                                                         tATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcc             <  1:234069/70‑1 (MQ=255)
                                                          aTCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtc              >  1:332396/1‑68 (MQ=255)
                                                           tCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTtcctcc          <  1:3341236/71‑1 (MQ=255)
                                                            ccGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                            <  1:2739097/52‑1 (MQ=255)
                                                            ccGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgat                            <  1:3576466/52‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCTGGATGAGACTATACATCGTCGAGATCGCAatga                             <  1:3886979/48‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGc                                  <  1:592404/43‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                           <  1:2090128/50‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatgatg                           <  1:3506040/50‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                             <  1:1473688/48‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                             <  1:453366/48‑1 (MQ=255)
                                                               tCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGc       >  1:419092/1‑70 (MQ=255)
                                                                cTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCg                                   <  1:14261/41‑1 (MQ=255)
                                                                 ttGTTCTGCCGCAGGATGAGACTATACATCGTCga                                        >  1:3343464/1‑35 (MQ=255)
                                                                   gTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAatga                             <  1:424019/44‑1 (MQ=255)
                                                                    ttCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACgg  <  1:2401273/70‑1 (MQ=255)
                                                                    |                                                                     
CTTTAATATTCATTATGTTAATAAGTTCTGGATGAGTGTGCTCTCATTGGATTTTTTATCCGTCTTGTTCTGCCGCAGGATGAGACTATACATCGTCGAGATCGCAATGATGATATGAGTTGTTCCTCCTGCGACGG  >  minE/2571317‑2571453

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: