Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,627,453 T→C 31.7% intergenic (+193/‑78) 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,4530TC31.7% 109.7 / 66.1 126intergenic (+193/‑78)rpoH/livJRNA polymerase, sigma 32 (sigma H) factor/leucine/isoleucine/valine transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (57/29);  new base C (29/11);  total (86/40)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 140 total aligned reads displayed.

ATAGCTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGG  >  minE/2627387‑2627520
                                                                  |                                                                   
aTAGCTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTc                                                                    >  1:3237733/1‑68 (MQ=255)
    cTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGc                                                                >  1:1828418/1‑68 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:2891451/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:2479585/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:2974663/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:1128474/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:718379/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:2800371/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:394848/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:3018728/1‑58 (MQ=255)
         ccAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGt                                                                     >  1:3413793/1‑58 (MQ=255)
              atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaa                                                      >  1:1317310/1‑68 (MQ=255)
              atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaa                                                      >  1:1287900/1‑68 (MQ=255)
              atCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTGTCATGCATAAAACGaaa                                                     >  1:1224754/1‑68 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:2713451/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:2287792/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:2555497/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:2068140/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:2919187/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:1699114/69‑1 (MQ=255)
               tCATAAAAATGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGCTACGTATAAAACGAAAt                                                    <  1:86059/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)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                     >  1:1723230/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                     >  1:3134148/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:2019809/1‑68 (MQ=255)
               tCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGaaa                                                     >  1:1107181/1‑68 (MQ=255)
                    aaaaTGGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:2676513/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:600934/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:1387801/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:3271402/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:1423567/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:3023615/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:221873/1‑69 (MQ=255)
                    aaaaTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGt                                               >  1:1920309/1‑69 (MQ=255)
                       aTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGc                                             >  1:2934754/1‑68 (MQ=255)
                       aTCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGc                                             >  1:2314334/1‑68 (MQ=255)
                              aTGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                     <  1:841115/69‑1 (MQ=255)
                              aTGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaac                                     <  1:1109529/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:1991927/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  <  1:500255/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  <  1:2395650/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  <  1:1717575/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:538281/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:2862916/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  <  1:2131736/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:540370/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:1253461/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  <  1:1851147/69‑1 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:2688931/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:1124701/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:604550/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:134756/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:817360/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:859805/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaac                                  >  1:2508800/1‑69 (MQ=255)
                                 ttttAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAaacaa                                   <  1:424749/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:405143/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                       <  1:512482/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                       <  1:2846393/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                       <  1:1821713/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                       <  1:1172861/69‑1 (MQ=255)
                                            tATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacaca                       <  1:2481067/69‑1 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                      >  1:417262/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                      >  1:3399875/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                      >  1:3022089/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                      >  1:183490/1‑68 (MQ=255)
                                              tgctgcTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAacacac                      >  1:2976067/1‑68 (MQ=255)
                                                ctgctAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGt                    <  1:2604261/68‑1 (MQ=255)
                                                ctgctAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGt                    >  1:2834637/1‑68 (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:1777977/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:2539198/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:2642287/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:2659112/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:550690/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:3173566/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:1687087/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:1499013/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:40473/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:3034127/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:1343352/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:309436/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:3354798/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacacac                      >  1:3286560/1‑57 (MQ=255)
                                                         cACGGGTAGCTACGTATAAAACGAAATAAAGTGCTGCACAACAACATCACAacaca                       >  1:2934237/1‑56 (MQ=255)
                                                          aCGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAacaaca                                 <  1:2439125/45‑1 (MQ=255)
                                                             ggTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAgaagaa      <  1:2593001/69‑1 (MQ=255)
                                                                 gTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATggg  >  1:2352766/1‑69 (MQ=255)
                                                                  |                                                                   
ATAGCTATTCCAATATCATAAAAATCGGGTATGTTTTAGCAGAGTATGCTGCTAAAGCACGGGTAGTCATGCATAAAACGAAATAAAGTGCTGAAAAACAACATCACAACACACGTAATAACCAGAAGAATGGG  >  minE/2627387‑2627520

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: