Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,627,456 T→C 30.0% intergenic (+196/‑75) rpoH → / → livJ RNA polymerase, sigma 32 (sigma H) factor/leucine/isoleucine/valine transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,627,4560TC30.0% 125.7 / 64.8 133intergenic (+196/‑75)rpoH/livJRNA polymerase, sigma 32 (sigma H) factor/leucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (56/37);  new base C (29/11);  total (85/48)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

Only 100 of 133 total aligned reads displayed.

CTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGA  >  minE/2627391‑2627522
                                                                 |                                                                  
cTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGc                                                                  >  1:1828418/1‑68 (MQ=255)
cTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGc                                                                  >  1:1281061/1‑68 (MQ=255)
cTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGc                                                                  >  1:448190/1‑68 (MQ=255)
          atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaa                                                        >  1:2876972/1‑68 (MQ=255)
          atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaa                                                        >  1:1287900/1‑68 (MQ=255)
          atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaa                                                        >  1:1317310/1‑68 (MQ=255)
          atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTGTCATGCATAAAACGaaa                                                       >  1:1224754/1‑68 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2919187/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2068140/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2287792/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:1699114/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2555497/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2713451/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:876819/69‑1 (MQ=255)
           tCATAAAAATGGGGTATGTTTTAGCAGAGTATCCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:1612487/69‑1 (MQ=255)
           tCATAAAAATGAGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2389189/69‑1 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACg                                                          >  1:3366619/1‑65 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:2458237/1‑68 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:2753797/1‑68 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:3310123/1‑68 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:2019809/1‑68 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:1107181/1‑68 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:1723230/1‑68 (MQ=255)
           tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:1081325/1‑68 (MQ=255)
                aaaaTGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:2676513/1‑69 (MQ=255)
                aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:600934/1‑69 (MQ=255)
                aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:3023615/1‑69 (MQ=255)
                aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:2931057/1‑69 (MQ=255)
                aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:1920309/1‑69 (MQ=255)
                aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:1387801/1‑69 (MQ=255)
                aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:221873/1‑69 (MQ=255)
                   aTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGc                                               >  1:2314334/1‑68 (MQ=255)
                          aTGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                       <  1:841115/69‑1 (MQ=255)
                          aTGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                       <  1:1109529/69‑1 (MQ=255)
                           tGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                       >  1:2609688/1‑68 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:859805/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:2688931/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:2807280/69‑1 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:2862916/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:134756/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:2508800/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:3109676/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:540370/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:2395650/69‑1 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:500255/69‑1 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:171611/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:604550/1‑69 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaa                                     <  1:408182/68‑1 (MQ=255)
                             ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaa                                     <  1:2001540/68‑1 (MQ=255)
                              tttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:1933772/68‑1 (MQ=255)
                                    agagTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCacaa                             <  1:2450874/69‑1 (MQ=255)
                                    agagTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCacaa                             <  1:2214616/69‑1 (MQ=255)
                                    agagTATGCTGCTAAAGCACGGGTACTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCacaa                             <  1:792363/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:840056/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:2966794/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:2846393/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:538671/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:1172861/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:512482/69‑1 (MQ=255)
                                        tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:405143/69‑1 (MQ=255)
                                          tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:3399875/1‑68 (MQ=255)
                                          tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:417262/1‑68 (MQ=255)
                                          tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:2976067/1‑68 (MQ=255)
                                          tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:3022089/1‑68 (MQ=255)
                                          tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:183490/1‑68 (MQ=255)
                                          tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:3033120/1‑68 (MQ=255)
                                            ctgctAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGt                      <  1:1831317/68‑1 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTg                                             >  1:2963701/1‑36 (MQ=16)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGGACAACAACATCACAacacac                        >  1:3160496/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACTTCACAacacac                        >  1:2554691/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAataca                         >  1:878634/1‑56 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:370838/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:550690/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1777977/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1728420/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1582422/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1499013/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:2539198/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1343352/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:861104/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:2642287/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:2659112/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1130596/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:3286560/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:327869/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:3173566/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:309436/1‑57 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacaca                         >  1:2934237/1‑56 (MQ=255)
                                                     cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATAACAacacac                        >  1:3193433/1‑57 (MQ=255)
                                                      aCGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAacaaca                                   <  1:2439125/45‑1 (MQ=255)
                                                         ggTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAgaagaa        <  1:1119225/69‑1 (MQ=255)
                                                         ggTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAgaagaa        <  1:2593001/69‑1 (MQ=255)
                                                             gTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATggg    >  1:2352766/1‑69 (MQ=255)
                                                               cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:1902732/69‑1 (MQ=255)
                                                               cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:2172956/69‑1 (MQ=255)
                                                               cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:1699310/69‑1 (MQ=255)
                                                               cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:2235898/69‑1 (MQ=255)
                                                               cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:258502/69‑1 (MQ=255)
                                                               cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:2774974/69‑1 (MQ=255)
                                                                aTGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:3116938/68‑1 (MQ=255)
                                                                 |                                                                  
CTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGA  >  minE/2627391‑2627522

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: