Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,072,048 A→T 15.9% intergenic (‑88/‑131) ygjH ← / → ebgR conserved hypothetical protein/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,072,0480AT15.9% 186.9 / 25.7 75intergenic (‑88/‑131)ygjH/ebgRconserved hypothetical protein/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base T (7/5);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAGGT  >  minE/2071978‑2072117
                                                                      |                                                                     
tAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGa                                                                       >  1:2473560/1‑71 (MQ=255)
 aaaaTTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:2909390/1‑71 (MQ=255)
  aaaTTTTACTAAAGTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:3327308/1‑70 (MQ=255)
  aaaTTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:588453/1‑70 (MQ=255)
  aaaTTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:2612748/1‑70 (MQ=255)
  aaaTTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:4073224/1‑70 (MQ=255)
  aaaTTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:134517/1‑70 (MQ=255)
  aaaTTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATc                                                                     <  1:3975109/71‑1 (MQ=255)
      tttACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACAt                                                                 >  1:3016862/1‑71 (MQ=255)
        tACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGTTCACAtt                                                                <  1:3486635/70‑1 (MQ=255)
          cTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGt                                                             <  1:501291/71‑1 (MQ=255)
          cTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGt                                                             <  1:3588682/71‑1 (MQ=255)
             aaCTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAg                                                          >  1:3332902/1‑71 (MQ=255)
             aaCTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAg                                                          >  1:1127344/1‑71 (MQ=255)
             aaCTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAg                                                          >  1:2716116/1‑71 (MQ=255)
             aaCTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAg                                                          >  1:3794556/1‑71 (MQ=255)
              aCTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      >  1:3943840/1‑58 (MQ=255)
               cTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGcc                                                        >  1:2501550/1‑71 (MQ=255)
                  taGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCt                                                     <  1:2419774/71‑1 (MQ=255)
                  taGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCt                                                     >  1:174119/1‑71 (MQ=255)
                   aGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGAt                                                                      <  1:3365210/53‑1 (MQ=255)
                    gCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGTTCACATTCGTGAGCCGCtgt                                                   <  1:4126042/71‑1 (MQ=255)
                    gCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACAt                                                                 <  1:1593311/57‑1 (MQ=255)
                     cAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATc                                                                     <  1:3692542/52‑1 (MQ=255)
                     cAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATc                                                                     <  1:446386/52‑1 (MQ=255)
                       aaaGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATcac                                                                   >  1:2689647/1‑52 (MQ=255)
                       aaaGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATcac                                                                   >  1:2876583/1‑52 (MQ=255)
                       aaaGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATcac                                                                   >  1:1267557/1‑52 (MQ=255)
                        aaGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCtg                                                    <  1:3688487/66‑1 (MQ=255)
                         aGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGTTCACATTCGTGAGCCGCTGTGTTAc                                              >  1:308740/1‑71 (MQ=255)
                         aGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGc                                                      <  1:1308207/63‑1 (MQ=255)
                         aGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTAc                                              <  1:1385280/71‑1 (MQ=255)
                         aGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTAc                                              >  1:952593/1‑71 (MQ=255)
                          gTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTg                                                            >  1:1025259/1‑56 (MQ=255)
                          gTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTAcc                                             <  1:593190/71‑1 (MQ=255)
                                  cAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGc                                     >  1:937328/1‑71 (MQ=255)
                                     tCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATca                                                                    <  1:3630228/37‑1 (MQ=255)
                                     tCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCg                                                              <  1:3157155/43‑1 (MQ=255)
                                     tCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCg                                                              <  1:837431/43‑1 (MQ=255)
                                       atGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCa                                 >  1:3514114/2‑70 (MQ=255)
                                          tAGGCTAAAAATGGAGAATGCAGGCGTGATCACAt                                                                 <  1:3129576/35‑1 (MQ=255)
                                          tAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGa                             <  1:2839159/71‑1 (MQ=255)
                                           aGGCTAAAAATGGAGAATGCAGGCGTGTTCACATTCGTGAGCCGCtgt                                                   >  1:2844746/1‑48 (MQ=255)
                                            ggCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGaaa                           >  1:641826/1‑71 (MQ=255)
                                            ggCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGaa                            <  1:2896277/70‑1 (MQ=255)
                                               tAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGaa                            >  1:3627955/1‑67 (MQ=255)
                                                 aaaaTGGAGAATGCAGGCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACgcgc                      <  1:411437/71‑1 (MQ=255)
                                                  aaaTGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGtt                                                <  1:1695928/44‑1 (MQ=255)
                                                   aaTGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCtt                    <  1:402377/71‑1 (MQ=255)
                                                   aaTGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCtt                    <  1:651200/71‑1 (MQ=255)
                                                     tGGAGAATGCAGGCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCttta                  >  1:2601346/1‑71 (MQ=255)
                                                     tGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCttta                  >  1:723056/1‑71 (MQ=255)
                                                     tGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCttta                  >  1:674126/1‑71 (MQ=255)
                                                     tGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCttta                  <  1:1308135/71‑1 (MQ=255)
                                                      ggAGAATGCAGGCGTGTTCACATTCGTGAGCCGCTGTGTTAc                                              >  1:331566/1‑42 (MQ=255)
                                                       gagaATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGaaa                           >  1:14286/1‑60 (MQ=255)
                                                        agaATGCAGGCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCtttattt               <  1:3136967/71‑1 (MQ=255)
                                                           aTGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACt            >  1:279844/1‑71 (MQ=255)
                                                           aTGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGCTACAGCGTCAAAGAAACGCGCTTTATTTACt            >  1:2577925/1‑71 (MQ=255)
                                                            tGCAGGCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACt            >  1:2597699/1‑70 (MQ=255)
                                                            tGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGa                             <  1:3454742/53‑1 (MQ=255)
                                                            tGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACt            >  1:2042714/1‑70 (MQ=255)
                                                             gCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTAc             >  1:2384943/1‑68 (MQ=255)
                                                              catgCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGaa         <  1:1022873/68‑1 (MQ=255)
                                                              cAGGCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGaa         <  1:794984/71‑1 (MQ=255)
                                                              cAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGaa         <  1:3909210/71‑1 (MQ=255)
                                                              cAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGa          <  1:3883158/70‑1 (MQ=255)
                                                              cAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGa          <  1:655804/70‑1 (MQ=255)
                                                               aGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGaaa        <  1:2327213/71‑1 (MQ=255)
                                                                 gCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGaaa                           >  1:1375469/1‑50 (MQ=255)
                                                                 gCGTGTTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGaaa                           >  1:3713981/1‑50 (MQ=255)
                                                                  cGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCttta                  <  1:1607485/58‑1 (MQ=255)
                                                                  cGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGaaaa       >  1:1535650/1‑69 (MQ=255)
                                                                    tGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAg                              <  1:729342/44‑1 (MQ=255)
                                                                    tGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCttta                  <  1:824979/56‑1 (MQ=255)
                                                                    tGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAg    <  1:3059523/70‑1 (MQ=255)
                                                                      aTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGaaaa       >  1:1474843/1‑65 (MQ=255)
                                                                      aTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAg    >  1:2132493/1‑68 (MQ=255)
                                                                      aTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAg    >  1:2858049/1‑68 (MQ=255)
                                                                      aTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAGGt  >  1:1982877/1‑70 (MQ=255)
                                                                      aTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAGGt  >  1:900696/1‑70 (MQ=255)
                                                                      aTCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAGGt  >  1:390915/1‑70 (MQ=255)
                                                                      |                                                                     
TAAAATTTTACTAAACTATAGCAAAAGTTTTTCTCAATCCTGTAGGCTAAAAATGGAGAATGCAGGCGTGATCACATTCGTGAGCCGCTGTGTTACCGTTACAGCGTCAAAGAAACGCGCTTTATTTACTGAAAACAGGT  >  minE/2071978‑2072117

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