Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,045,864 T→A 6.6% intergenic (+626/+169) yqiK → / ← rfaE conserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,045,8640TA 6.6% 240.5 / 6.7 88intergenic (+626/+169)yqiK/rfaEconserved hypothetical protein/fused heptose 7‑phosphate kinase and heptose 1‑phosphate adenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (31/51);  new base A (2/4);  total (33/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAGGGGG  >  minE/2045794‑2045930
                                                                       |                                                                  
gcaAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCt                                                                    >  1:3826505/1‑71 (MQ=255)
  aaGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                                  >  1:1376259/1‑71 (MQ=255)
     acgCTGAAACGGGAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGt                                                                 >  1:1708430/3‑70 (MQ=0)
      ggCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGtcat                                                              <  1:1936945/71‑1 (MQ=255)
         tGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCg                                                                  <  1:748236/64‑1 (MQ=255)
           aaaCGTGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                                           <  1:3550215/69‑1 (MQ=255)
           aaaCGGGAAAGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcat                                                           <  1:3128594/69‑1 (MQ=11)
           aaaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                                           <  1:2923315/69‑1 (MQ=255)
           aaaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                                                         >  1:673150/1‑71 (MQ=255)
            aaCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATAc                                                         <  1:2633942/70‑1 (MQ=255)
             aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                                                   >  1:2261239/1‑59 (MQ=255)
             aCGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTc                                                                   >  1:3282549/1‑59 (MQ=255)
               gggAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTgg                                                      <  1:1049747/71‑1 (MQ=11)
               gggAAAGCCTCTCCCGGAGAAGAGGGCTTTTAATAAGGAAAGGGTTATGATGAAGCACGTCATCATACTgg                                                      <  1:3376571/71‑1 (MQ=11)
                 gAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                                                   >  1:2125140/1‑71 (MQ=255)
                  aaaGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgt                                                  <  1:484967/71‑1 (MQ=255)
                    aGCCCCTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatc                                                          <  1:1334196/61‑2 (MQ=9)
                       cccTCCCGAGGAAGGGGCCATAAATAAGGAAAGGGTCATGATGAAGCTACTcatcatcg                                                         <  1:994508/59‑3 (MQ=255)
                       cccTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt                                             <  1:3639705/71‑1 (MQ=255)
                         cTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgt                                                >  1:1881117/1‑66 (MQ=255)
                         cTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTa                                           <  1:3464314/71‑1 (MQ=255)
                         cTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTa                                           >  1:3654444/1‑71 (MQ=255)
                         cTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTa                                           >  1:1424878/1‑71 (MQ=255)
                           cccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACt                                         >  1:897779/1‑71 (MQ=255)
                            ccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt                                                                 >  1:11342/1‑46 (MQ=255)
                            ccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGt                                                                 >  1:1171079/1‑46 (MQ=255)
                            ccGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTg                                        <  1:1935560/71‑1 (MQ=255)
                             cGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGt                                       <  1:2794176/71‑1 (MQ=255)
                              gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                                           <  1:2934175/50‑1 (MQ=255)
                              gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTcatcat                                                           <  1:2451304/50‑1 (MQ=255)
                              gaagaaGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      <  1:321771/71‑1 (MQ=255)
                               aagaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      >  1:3731076/1‑70 (MQ=255)
                               aagaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      >  1:2619352/1‑70 (MQ=255)
                                agaagGGGCCTTGTATAAGGAAAGGTTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGt                                       >  1:2305512/1‑68 (MQ=255)
                                 gaagGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      <  1:1890929/68‑1 (MQ=255)
                                    ggggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgt                                                >  1:2122189/1‑55 (MQ=255)
                                    ggggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgt                                                >  1:2667839/1‑55 (MQ=255)
                                    ggggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                                                 <  1:3131370/54‑1 (MQ=255)
                                    ggggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtg                                                 <  1:1989837/54‑1 (MQ=255)
                                    ggggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGttt                                <  1:881928/71‑1 (MQ=255)
                                       ggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt                                            <  1:3154011/57‑1 (MQ=255)
                                       ggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTgttgtt                                            <  1:2941019/57‑1 (MQ=255)
                                       ggCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTg                                        <  1:3865645/61‑1 (MQ=255)
                                        gCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTccc                             <  1:2440058/71‑1 (MQ=255)
                                         ccTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTTCTGTTAAGTTTccc                             >  1:3193328/1‑70 (MQ=255)
                                         ccTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCg                            <  1:2919349/71‑1 (MQ=255)
                                         ccTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGATTccc                             >  1:595478/1‑70 (MQ=255)
                                          cTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttgtgttgt                                             <  1:3029819/53‑1 (MQ=255)
                                          cTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGa                           >  1:1380385/1‑71 (MQ=255)
                                          cTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGa                           <  1:1786667/71‑1 (MQ=255)
                                             gTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGa                           <  1:1747711/68‑1 (MQ=255)
                                              tataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGtt                                 <  1:1645400/61‑1 (MQ=255)
                                              aataAGGAAAGGGTTATGATGAAGCACGTCATCATAc                                                         >  1:3496380/2‑37 (MQ=11)
                                              aataAGGAAAGGGTTATGATGAAGCACGTCATCATACTg                                                       <  1:1770919/38‑1 (MQ=11)
                                               ataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc                          <  1:1473262/67‑1 (MQ=255)
                                               ataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc                          <  1:3658789/67‑1 (MQ=255)
                                               ataAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGAc                          <  1:2119217/67‑1 (MQ=255)
                                                taAGGAAAGGGTTATGATGAAGCTCGTCATCATACTg                                                       <  1:2250128/37‑1 (MQ=33)
                                                taAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGa                           >  1:3412248/1‑65 (MQ=255)
                                                  aGGAAAGGGTTATGATGAAGCTCGTCATCATACTGgttg                                                   <  1:594639/39‑1 (MQ=40)
                                                  aGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTc                               >  1:2682945/1‑59 (MQ=255)
                                                   ggAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCg                            <  1:1113711/61‑1 (MQ=255)
                                                   ggAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCg                            <  1:3240690/61‑1 (MQ=255)
                                                    gAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTTTGTTGTTACTGTTAAGTTTCCCGACTTACTaaca                 <  1:1091101/71‑1 (MQ=255)
                                                     aaaGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTg                                        >  1:954642/1‑47 (MQ=40)
                                                      taGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      >  1:1094889/2‑48 (MQ=255)
                                                      aaGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      >  1:3733424/1‑48 (MQ=255)
                                                           ttATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTa                                     <  1:935167/44‑1 (MQ=255)
                                                            tATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGtt                                      <  1:153998/42‑1 (MQ=255)
                                                             atgatgAAGCTCGTCATCATACTGGTTGTgttgtt                                            <  1:2637366/35‑1 (MQ=40)
                                                             atgatgAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTa                       >  1:3399207/1‑56 (MQ=255)
                                                              tgatgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCag        <  1:3218432/70‑1 (MQ=255)
                                                               gatgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGttt                                <  1:3270620/45‑1 (MQ=255)
                                                                atgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGt                                  <  1:26256/42‑1 (MQ=255)
                                                                atgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGt                                  <  1:991015/42‑1 (MQ=255)
                                                                atgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaa                   >  1:2594083/1‑57 (MQ=255)
                                                                atgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaa                   >  1:3045629/1‑57 (MQ=255)
                                                                atgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAgg     <  1:364460/71‑1 (MQ=255)
                                                                 tgaAGCTCGTCATCATACTGGTTGTGTTGTTACTTTTAAGTTTCCCGACTTACTAACAACTCATCAGAggg    <  1:732493/71‑1 (MQ=255)
                                                                 tgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTaa                                    >  1:3811719/1‑39 (MQ=255)
                                                                 tgaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAggg    <  1:3064798/71‑1 (MQ=255)
                                                                  gaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaacaa                >  1:393412/1‑58 (MQ=255)
                                                                  gaAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTaacaa                >  1:221587/1‑58 (MQ=255)
                                                                    aGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTc                               >  1:1955551/1‑41 (MQ=255)
                                                                    aGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAggggg  <  1:3009946/70‑1 (MQ=255)
                                                                    aGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAggggg  <  1:3932634/70‑1 (MQ=255)
                                                                    aGCTCGTCATCATACTGGTTGAGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAggggg  <  1:458582/70‑1 (MQ=255)
                                                                    aGCACGTCATCATACTGGTGATACTCTTAGTGATTAGCTTCCAGGCTTACTaag                  <  1:1602202/54‑2 (MQ=255)
                                                                      cTCGTCATCATACTGGGTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACt              <  1:3433658/56‑1 (MQ=255)
                                                                       |                                                                  
GCAAGGGGCTGAAACGGGAAAGCCCCTCCCGAAGAAGGGGCCTTGTATAAGGAAAGGGTTATGATGAAGCTCGTCATCATACTGGTTGTGTTGTTACTGTTAAGTTTCCCGACTTACTAACAACTCATCAGAGGGGG  >  minE/2045794‑2045930

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: