Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,138,595 G→C 7.5% intergenic (+103/‑30) ydjM → / → ydjN predicted inner membrane protein regulated by LexA/predicted transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,138,5950GC 7.5% 203.0 / 5.6 67intergenic (+103/‑30)ydjM/ydjNpredicted inner membrane protein regulated by LexA/predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (23/37);  new base C (2/3);  total (25/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTTGAATATAAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCGTGGTGTTCGT  >  minE/1138530‑1138662
                                                                 |                                                                   
tttGAATATAAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGa                                                                <  1:300637/71‑1 (MQ=255)
   gAATATAAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGaaa                                                              >  1:1141495/1‑70 (MQ=255)
     atatAAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAaata                                                            <  1:428306/70‑1 (MQ=255)
      tataAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAaataat                                                          <  1:277985/71‑1 (MQ=255)
      tataAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAaataat                                                          <  1:3479367/71‑1 (MQ=255)
      tataAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAaataa                                                           >  1:1110146/1‑70 (MQ=255)
          agagCCAGCTCAAAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAg                                                       <  1:2799240/70‑1 (MQ=255)
           gagCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGAt                                                     <  1:924899/71‑1 (MQ=255)
              ccAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGAt                                                     <  1:373254/68‑1 (MQ=255)
               cAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAgg                                                 <  1:1928560/71‑1 (MQ=255)
                  cTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGaaa                                                              <  1:312233/55‑1 (MQ=255)
                  cTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGaaa                                                              <  1:1857334/55‑1 (MQ=255)
                  cTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGgag                                               >  1:1118720/1‑70 (MQ=255)
                     acaGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAAc                                            <  1:1460782/70‑1 (MQ=255)
                     acaGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACg                                           <  1:3487009/71‑1 (MQ=255)
                         ttCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAg                                                                 <  1:3077127/45‑1 (MQ=255)
                         ttCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGgag                                               >  1:1046507/1‑63 (MQ=255)
                           cTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGaa                                                               >  1:2003615/1‑45 (MQ=255)
                             gTTAATCTTGCGCCAACACTATGACTGCTACGCAGTTATAGAAATAATAAGATCAGgag                                               <  1:3441730/59‑1 (MQ=255)
                             gTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGgag                                               <  1:621010/59‑1 (MQ=255)
                             gTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGgag                                               <  1:1226451/59‑1 (MQ=255)
                             gTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGgag                                               <  1:1617352/59‑1 (MQ=255)
                              ttAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAaataat                                                          <  1:495581/47‑1 (MQ=255)
                              ttAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAaataat                                                          <  1:3724223/47‑1 (MQ=255)
                              ttAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATc                                                    >  1:2954050/1‑53 (MQ=255)
                              ttAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATc                                                    >  1:3340086/1‑53 (MQ=255)
                                aaTCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACt                                 <  1:1950645/70‑1 (MQ=255)
                                aaTCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACt                                 <  1:51457/70‑1 (MQ=255)
                                      gcgcCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCAt                           >  1:412218/1‑70 (MQ=255)
                                          cAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACg                                           >  1:2052122/1‑50 (MQ=255)
                                          cAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAAtt                      <  1:406686/71‑1 (MQ=255)
                                            acacTATGACTGCTACGCAGTGATAGAAATAATAAGATCCGGAGAACGGGGATGAACTTTCCATTAATTGc                    >  1:1861929/1‑71 (MQ=255)
                                            acacTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGc                    >  1:874033/1‑71 (MQ=255)
                                            acacTATGACTGCTACGCAGTGATAGAAAAAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGc                    >  1:3161277/1‑71 (MQ=255)
                                             cacTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCg                   <  1:2297920/71‑1 (MQ=255)
                                              acTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCg                   >  1:1363814/1‑70 (MQ=255)
                                               cTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 <  1:3235844/71‑1 (MQ=255)
                                               cTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 <  1:269694/71‑1 (MQ=255)
                                               cTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 <  1:2566852/71‑1 (MQ=255)
                                               cTATCACTGCGTAGCAGTCATAGAAATAATAAGATCAGGAGACCGGGGATGAACTTTCCATTAATTGc                    >  1:1526580/1‑68 (MQ=25)
                                               cTATCACTGCGTAGCAGTCATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 >  1:2680169/1‑71 (MQ=37)
                                                tATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 >  1:3070607/1‑70 (MQ=255)
                                                tATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 >  1:296202/1‑70 (MQ=255)
                                                 aTGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACgggg                                        >  1:1017958/1‑46 (MQ=255)
                                                  tGACTGCTACGCAGTTATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACa               <  1:3713607/70‑1 (MQ=255)
                                                  tGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTcc                             <  1:1434354/56‑1 (MQ=255)
                                                  tGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAAc                <  1:1480234/69‑1 (MQ=255)
                                                  tGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAAc                <  1:2387466/69‑1 (MQ=255)
                                                   gACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATc             >  1:1407466/1‑71 (MQ=255)
                                                   gACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATc             >  1:2438544/1‑71 (MQ=255)
                                                   gACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATc             >  1:259485/1‑71 (MQ=255)
                                                    aCTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCg                   <  1:3343315/64‑1 (MQ=255)
                                                      tGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGaa                 >  1:2223523/1‑64 (MQ=255)
                                                        cTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCgtggt        >  1:2875350/1‑71 (MQ=255)
                                                        cTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCgtgg         >  1:3548486/1‑70 (MQ=255)
                                                         tACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCgtggtg       >  1:2456339/1‑71 (MQ=255)
                                                          tagCAGTCATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTc                              <  1:3485450/45‑1 (MQ=255)
                                                          tagCAGTCATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTc                              <  1:86682/45‑1 (MQ=255)
                                                           cGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCgtggtg       <  1:967652/69‑1 (MQ=255)
                                                           cGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCgtggtg       <  1:1357581/69‑1 (MQ=255)
                                                           cGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCgtggtg       <  1:3432979/69‑1 (MQ=255)
                                                             cAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCGTGGTGTTCg   <  1:641395/71‑1 (MQ=255)
                                                             cAGTCATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCGTGGTGTTCg   <  1:405275/71‑1 (MQ=255)
                                                              aGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCGTGGTGTTCGt  <  1:266011/71‑1 (MQ=255)
                                                              aGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCGTGGTGTTCGt  <  1:1667628/71‑1 (MQ=255)
                                                               gTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACt                                 <  1:3192717/39‑1 (MQ=255)
                                                               gTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACt                                 <  1:180999/39‑1 (MQ=255)
                                                                 |                                                                   
TTTGAATATAAGAGCCAGCTCACAGTTCTGTTAATCTTGCGCCAACACTATGACTGCTACGCAGTGATAGAAATAATAAGATCAGGAGAACGGGGATGAACTTTCCATTAATTGCGAACATCGTGGTGTTCGT  >  minE/1138530‑1138662

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: