Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,011,880 C→T 5.3% intergenic (+6/‑70) rstB → / → tus sensory histidine kinase in two‑component regulatory system with RstA/inhibitor of replication at Ter, DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,011,8800CT 5.3% 300.4 / 4.2 95intergenic (+6/‑70)rstB/tussensory histidine kinase in two‑component regulatory system with RstA/inhibitor of replication at Ter, DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (57/32);  new base T (3/2);  total (61/34)
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.91e-01

Only 100 of 108 total aligned reads displayed.

TGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAATATTA  >  minE/1011811‑1011950
                                                                     |                                                                      
tggtgCCCGTTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTAcg                                                                       >  1:2570944/1‑71 (MQ=255)
  gtgCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTAcgcg                                                                     <  1:474914/71‑1 (MQ=255)
     cccGTTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCAc                                                                  >  1:2915673/1‑71 (MQ=255)
     cccGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCAc                                                                  >  1:2642375/1‑71 (MQ=255)
       cgtttCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGa                                                                >  1:2726033/4‑71 (MQ=255)
       cgtttCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGa                                                                <  1:1545750/68‑1 (MQ=255)
       cGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGa                                                                >  1:4539383/1‑71 (MQ=255)
         tttCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCAc                                                                  <  1:2031391/66‑1 (MQ=255)
           tctcGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGt                                                            >  1:2777460/1‑71 (MQ=255)
             tcGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTc                                                           >  1:1927605/1‑70 (MQ=255)
               gTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAg                                                        >  1:3602433/1‑71 (MQ=255)
               gTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAg                                                        >  1:3968692/1‑71 (MQ=255)
                tttAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCAc                                                                  <  1:374338/60‑1 (MQ=255)
                tttAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGt                                                       <  1:2225020/71‑1 (MQ=255)
                  tAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTAc                                                                        >  1:2511513/1‑52 (MQ=255)
                  caGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTATGCGCACGGTGGTCaa                                                         >  1:475100/2‑67 (MQ=255)
                    gCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTcac                                                    <  1:811990/70‑1 (MQ=255)
                     cTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTATGCGCACGGTGGTCAAGTcacca                                                  >  1:682009/1‑71 (MQ=255)
                     cTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTATGCGCACGGTGGTCAAGTcacca                                                  >  1:3823099/1‑71 (MQ=255)
                        gCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACg                                                <  1:4203734/70‑1 (MQ=255)
                          cGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACt                                             >  1:768443/1‑71 (MQ=255)
                            ttATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTcac                                                    <  1:1623501/62‑1 (MQ=255)
                             tATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACtgtg                                          >  1:3340326/1‑71 (MQ=255)
                             tATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACtgtg                                          >  1:185240/1‑71 (MQ=255)
                             tATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTAAGCGCACGATGGTCaa                                                         >  1:4319893/1‑56 (MQ=255)
                              aTGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCa                                                                   >  1:1005559/1‑45 (MQ=255)
                              aTGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGa                                                                >  1:1166151/1‑48 (MQ=255)
                              aTGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGAt                                                               >  1:254354/1‑49 (MQ=255)
                              aTGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCaa                                                         >  1:172501/1‑55 (MQ=255)
                              aTGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGc                                         >  1:4569531/1‑71 (MQ=255)
                              aTGGCATAACATCCCGCAATTTAACTCTGCCTGACACTACGCGCAc                                                                  >  1:4000592/1‑46 (MQ=255)
                                ggCATAACATCCCGCAATTTACCTCTGCCTGACACTAc                                                                        <  1:2695230/38‑1 (MQ=255)
                                  cATAACATCCCGCAATTTACCTCTGCCTGACACTATGCGCACGGTGGTCAAGTCACCACCACTGTGc                                         <  1:2216810/67‑1 (MQ=255)
                                  cATAACATCCCGCAATTTACCTCTGCCTGACACTATGCGCACGGTGGTCAAGTCACCACCACTGTGc                                         <  1:2375065/67‑1 (MQ=255)
                                   aTAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCtata                                     >  1:2611694/1‑70 (MQ=255)
                                    tAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGt                                                       <  1:3506851/51‑1 (MQ=255)
                                    tAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaa                                    >  1:996436/1‑70 (MQ=255)
                                     aaCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  >  1:3135099/1‑71 (MQ=255)
                                     aaCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  >  1:620/1‑71 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCAc                                                                  <  1:1339489/38‑1 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGAt                                                               >  1:4666634/1‑41 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTca                                                     >  1:4464569/1‑51 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  >  1:4063803/1‑70 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  >  1:369222/1‑70 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  <  1:1380512/70‑1 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  >  1:1241094/1‑70 (MQ=255)
                                      aCATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAAt                                 <  1:262992/71‑1 (MQ=255)
                                      aCATCCCGCAATTTACCTATGCCTGACACTACGCGCAc                                                                  <  1:41805/38‑1 (MQ=255)
                                       cATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCaa                                                         <  1:937920/46‑1 (MQ=255)
                                       cATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCaa                                                         <  1:347264/46‑1 (MQ=255)
                                        aTCCCGCAATTTACCTCTGCCTGACACTACGCGCACGa                                                                >  1:2409716/1‑38 (MQ=255)
                                          cccGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAg                              >  1:2184092/1‑70 (MQ=255)
                                           ccGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTa                            >  1:1797735/1‑71 (MQ=255)
                                           ccGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTa                            >  1:2852909/1‑71 (MQ=255)
                                            cGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGa                                               >  1:527718/1‑51 (MQ=255)
                                            cGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGa                                               >  1:2637493/1‑51 (MQ=255)
                                            cGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTAt                           <  1:1181111/71‑1 (MQ=255)
                                              cAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTa                            >  1:955793/1‑68 (MQ=255)
                                              cAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTAtgt                         >  1:4674842/1‑71 (MQ=255)
                                                aTTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTa                            <  1:978595/66‑1 (MQ=255)
                                                 tttACCTCTGCCTGACACTACGCGCACGATGGTCAAGTc                                                      <  1:677232/39‑1 (MQ=255)
                                                 tttACCTCTGCCTGACACTACGCGCACGATGGTCAAGTc                                                      <  1:288422/39‑1 (MQ=255)
                                                  ttACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTAtgttgt                      >  1:629502/1‑70 (MQ=255)
                                                   tACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCtat                                      >  1:3262630/1‑53 (MQ=255)
                                                   tACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCt                                        >  1:4304755/1‑51 (MQ=255)
                                                   tACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTaa                    >  1:1052132/1‑71 (MQ=255)
                                                   tACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTa                     >  1:1989124/1‑70 (MQ=255)
                                                    aCCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAAc                   >  1:2203901/1‑71 (MQ=255)
                                                        ctGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTAtgtt                        >  1:1996724/1‑62 (MQ=255)
                                                        ctGCCTGACACTACGCCCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTa                 >  1:1840981/1‑69 (MQ=255)
                                                          gCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTaa                <  1:640885/68‑1 (MQ=255)
                                                           ccTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGt                             >  1:1742874/1‑54 (MQ=255)
                                                           ccTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGt                             >  1:3081929/1‑54 (MQ=255)
                                                           ccTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTaaa               <  1:1552788/68‑1 (MQ=255)
                                                            cTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATaaaa                                  <  1:1386062/48‑1 (MQ=255)
                                                            cTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTa                            >  1:3865802/1‑54 (MQ=255)
                                                             tGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAg                              >  1:2452472/1‑51 (MQ=255)
                                                             tGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAg                              >  1:1038228/1‑51 (MQ=255)
                                                               acacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGtt         <  1:3765521/70‑1 (MQ=255)
                                                               acacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGt          <  1:3353922/69‑1 (MQ=255)
                                                               acacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTa        <  1:2664716/71‑1 (MQ=255)
                                                               acacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTa        <  1:2057627/71‑1 (MQ=255)
                                                               acacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAACTATGTTGTAACTAAAGTGGt          <  1:282468/69‑1 (MQ=255)
                                                                cacTACGCGCACGATGGTCAAGTCACCACGACTGTGCta                                       >  1:2033077/1‑39 (MQ=255)
                                                                cacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTaa       <  1:1047095/71‑1 (MQ=255)
                                                                cacTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTa        <  1:2732933/70‑1 (MQ=255)
                                                                   tACGCGCAGGATGGTCAAGTCACCACGACTGTGCtat                                      <  1:3493641/37‑1 (MQ=255)
                                                                   tACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAata     >  1:1860552/1‑70 (MQ=255)
                                                                    aCGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTAtgttg                       <  1:3373481/51‑1 (MQ=255)
                                                                    aCGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTaa                    >  1:3958108/1‑54 (MQ=255)
                                                                    aCGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAatat    >  1:2506489/1‑70 (MQ=255)
                                                                    aCGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAatat    >  1:1396465/1‑70 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAg                              >  1:1068406/1‑43 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTaa                    >  1:3693185/1‑53 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTaa                    >  1:2934904/1‑53 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAatat    <  1:1770062/69‑1 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAatat    <  1:2048514/69‑1 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAat      >  1:1536969/1‑67 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAAtatta  >  1:1350408/1‑71 (MQ=255)
                                                                     cgcgCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAAtatta  >  1:3350105/1‑71 (MQ=255)
                                                                     |                                                                      
TGGTGCCCGCTTCTCGTTTAGCTGGCCGTTATGGCATAACATCCCGCAATTTACCTCTGCCTGACACTACGCGCACGATGGTCAAGTCACCACGACTGTGCTATAAAATAAGTATGTTGTAACTAAAGTGGTTAATATTA  >  minE/1011811‑1011950

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: