Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,627,454 C→T 32.1% intergenic (+194/‑77) 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,4540CT32.1% 112.7 / 62.9 125intergenic (+194/‑77)rpoH/livJRNA polymerase, sigma 32 (sigma H) factor/leucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (55/29);  new base T (29/11);  total (85/40)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

ATAGCTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGA  >  minE/2627387‑2627522
                                                                   |                                                                    
aTAGCTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTc                                                                      >  1:3237733/1‑68 (MQ=255)
    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:2713451/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2555497/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:1699114/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2919187/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:876819/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2068140/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                      <  1:2287792/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:2019809/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:1107181/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:1723230/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:2753797/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:3310123/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:1081325/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                       >  1:3134148/1‑68 (MQ=255)
                    aaaaTGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:2676513/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:600934/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:1920309/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:1387801/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:2931057/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:3023615/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                                 >  1:221873/1‑69 (MQ=255)
                       aTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGc                                               >  1:2314334/1‑68 (MQ=255)
                              aTGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                       <  1:1109529/69‑1 (MQ=255)
                              aTGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                       <  1:841115/69‑1 (MQ=255)
                               tGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                       >  1:2609688/1‑68 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:2688931/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:3109676/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:859805/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:134756/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:2862916/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:2807280/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:171611/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:604550/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:540370/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:500255/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    <  1:2395650/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                    >  1:2508800/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:2214616/69‑1 (MQ=255)
                                        agagTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCacaa                             <  1:2450874/69‑1 (MQ=255)
                                        agagTATGCTGCTAAAGCACGGGTACTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCacaa                             <  1:792363/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:405143/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:1172861/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:2966794/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:538671/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:840056/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:2846393/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                         <  1:512482/69‑1 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:417262/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:183490/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:3399875/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:2976067/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                        >  1:3022089/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)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACTTCACAacacac                        >  1:2554691/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAataca                         >  1:878634/1‑56 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1130596/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:1777977/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:550690/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:1343352/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:2539198/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:861104/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:370838/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:2642287/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:2659112/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                        >  1:309436/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)
                                                         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:1699310/69‑1 (MQ=255)
                                                                   cATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGa  <  1:2172956/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)
                                                                   |                                                                    
ATAGCTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGGGA  >  minE/2627387‑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: