Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,311,509 T→C 15.1% intergenic (+20/‑47) hslU → / → menA molecular chaperone and ATPase component of HslUV protease/1,4‑dihydroxy‑2‑naphthoate octaprenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,311,5090TC15.1% 125.4 / 14.8 53intergenic (+20/‑47)hslU/menAmolecular chaperone and ATPase component of HslUV protease/1,4‑dihydroxy‑2‑naphthoate octaprenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (17/28);  new base C (6/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 6.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCTGGATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTA  >  minE/2311439‑2311574
                                                                      |                                                                 
cATCTGGATGCGGTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTtcat                                                                   <  1:805469/71‑1 (MQ=255)
 aTCTGGATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTtcatc                                                                  <  1:224553/71‑1 (MQ=255)
      gATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTtcatca                                                                 <  1:3066828/67‑1 (MQ=255)
      gATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTg                                                              <  1:1784148/70‑1 (MQ=255)
      gATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGt                                                             <  1:2686092/71‑1 (MQ=255)
        tGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGttt                                                           <  1:63271/71‑1 (MQ=255)
                   cAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTg                                                 >  1:2825287/1‑70 (MQ=255)
                   cAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTGa                                                <  1:3560332/71‑1 (MQ=255)
                   cAGATGAAGATATGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATggg                                                     >  1:59232/1‑66 (MQ=255)
                     gATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              <  1:2805898/71‑1 (MQ=255)
                      aTGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAg                                             <  1:3277467/71‑1 (MQ=255)
                        gAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAg                                             <  1:3270973/69‑1 (MQ=255)
                         aaGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGc                                           <  1:2919860/70‑1 (MQ=255)
                          aGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              >  1:196118/1‑66 (MQ=255)
                           gATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGcccc                                        <  1:181756/71‑1 (MQ=255)
                            aTCTTAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCa                                       >  1:3466071/1‑71 (MQ=255)
                            aTCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAgg                                            >  1:1391033/1‑66 (MQ=255)
                             tCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCAt                                      >  1:568882/1‑71 (MQ=255)
                              cTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGAt                                                        >  1:1147588/1‑52 (MQ=255)
                                 aGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTg                                                 <  1:2018413/56‑1 (MQ=255)
                                 aGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTg                                                 <  1:1679339/56‑1 (MQ=255)
                                  gCCGTTTTATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTa                                 >  1:1895459/1‑71 (MQ=255)
                                     gTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              >  1:335520/1‑55 (MQ=255)
                                     gTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTg                              >  1:864911/1‑71 (MQ=255)
                                      ttttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              <  1:2417407/54‑1 (MQ=255)
                                      ttttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              <  1:39275/54‑1 (MQ=255)
                                      ttttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              <  1:3926996/54‑1 (MQ=255)
                                      ttttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGa                                                <  1:2951851/52‑1 (MQ=255)
                                      ttttATCCTATAATCGCGTTCAATCAATTTCACCATTGTTTGATGGGGCTGAAAGGCCCCAttttt                                  <  1:2895614/66‑1 (MQ=255)
                                        ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                         >  1:2143447/1‑41 (MQ=255)
                                        ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGa                                                         >  1:3227448/1‑41 (MQ=255)
                                        ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGaaa                                              >  1:2199427/1‑52 (MQ=255)
                                        ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgcg                           >  1:3497955/1‑71 (MQ=255)
                                        ttATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgc                            <  1:1683236/70‑1 (MQ=255)
                                        ttATCCTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGgcg                           >  1:3361168/1‑71 (MQ=255)
                                         tATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCAtt                                     <  1:3376771/60‑1 (MQ=255)
                                           tCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTg                                                 >  1:3227655/1‑46 (MQ=255)
                                           tCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCAtt                                     <  1:2479327/58‑1 (MQ=255)
                                            ccTATAATCGCGTTCAATCATTTTCACCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGta                       >  1:862161/1‑71 (MQ=255)
                                             cTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGtat                      <  1:2236212/71‑1 (MQ=255)
                                             cTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGtat                      >  1:1204373/1‑71 (MQ=255)
                                             cTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGta                       <  1:1580596/70‑1 (MQ=255)
                                                 aaTCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGtattat                   >  1:3733435/1‑70 (MQ=255)
                                                  aTCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGa                 <  1:1887511/71‑1 (MQ=255)
                                                  aTCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGa                 <  1:264342/71‑1 (MQ=255)
                                                     gcgTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACt               >  1:924377/1‑70 (MQ=255)
                                                        ttCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCAt                                      <  1:386155/44‑1 (MQ=255)
                                                          cAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa         >  1:3488445/1‑71 (MQ=255)
                                                           aaTCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa         <  1:325881/70‑1 (MQ=255)
                                                           aaTCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGaacaa         <  1:2231868/70‑1 (MQ=255)
                                                            aTCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAacaaca       <  1:2189522/71‑1 (MQ=255)
                                                             tCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAacaaca       >  1:536490/1‑70 (MQ=255)
                                                              cATTTTCACCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACaaa     >  1:1175224/1‑71 (MQ=255)
                                                                ttttCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGtattat                   <  1:995156/55‑1 (MQ=255)
                                                                ttttCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAAtt   <  1:1745574/71‑1 (MQ=255)
                                                                  ttCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTa  <  1:1293219/70‑1 (MQ=255)
                                                                   tcatcaTTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGt                        >  1:761408/1‑47 (MQ=255)
                                                                      |                                                                 
CATCTGGATGCGTTGGTGGCAGATGAAGATCTGAGCCGTTTTATCCTATAATCGCGTTCAATCATTTTCATCATTGTTTGATGGGGCTGAAAGGCCCCATTTTTATTGGCGCGTATTATGACTGAACAACAAATTA  >  minE/2311439‑2311574

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: