Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,053,925 A→C 8.3% intergenic (+83/+98) cca → / ← bacA fused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase/undecaprenyl pyrophosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,053,9250AC 8.3% 156.3 / 7.5 60intergenic (+83/+98)cca/bacAfused tRNA nucleotidyl transferase and 2'3'‑cyclic phosphodiesterase/2'nucleotidase/phosphatase/undecaprenyl pyrophosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base C (3/2);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

TGCCGGATGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGAC  >  minE/2053856‑2053993
                                                                      |                                                                    
tGCCGGATGCGGCGTAACCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCa                                                                      <  1:4998899/70‑1 (MQ=255)
 gCCGGATGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAg                                                                    >  1:415947/1‑71 (MQ=255)
 gCCGGATGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAg                                                                    >  1:1407079/1‑71 (MQ=255)
    ggATGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                 <  1:4212813/71‑1 (MQ=255)
    ggATGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                 <  1:5082094/71‑1 (MQ=255)
       tGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                              <  1:1081898/71‑1 (MQ=255)
       tGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGc                                                              <  1:92199/71‑1 (MQ=255)
                aaCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGt                                                     >  1:1842913/1‑71 (MQ=255)
                 aCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGGAAGCGTCGCATCAGGCGtt                                                    >  1:3875910/1‑71 (MQ=255)
                 aCGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGtt                                                    <  1:129242/71‑1 (MQ=255)
                  cGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGtt                                                    >  1:3281464/1‑70 (MQ=255)
                   gCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGa                                                  >  1:328791/1‑71 (MQ=255)
                      ttATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGtt                                                    <  1:1417792/66‑1 (MQ=255)
                      ttATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTg                                               >  1:5198770/1‑71 (MQ=255)
                       tATCCGGCCTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCattgattg                                               >  1:4402974/1‑71 (MQ=2)
                        aTCCGGCCTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGc                                              >  1:3497838/1‑71 (MQ=2)
                        aTCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCa                                                                      <  1:2195553/46‑1 (MQ=255)
                        aTCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCa                                                                      <  1:3992913/46‑1 (MQ=255)
                        aTCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTg                                               <  1:1559422/69‑1 (MQ=255)
                            ggccTACAGGTCGGCAACAGTTGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGa                                          >  1:5129818/4‑71 (MQ=14)
                            ggCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGAt                                         >  1:2698069/1‑71 (MQ=255)
                             gCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                 <  1:5031998/46‑1 (MQ=255)
                             gCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                 <  1:3971231/46‑1 (MQ=255)
                             gCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                 <  1:720948/46‑1 (MQ=255)
                                  cATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCg                                                               >  1:968648/1‑44 (MQ=255)
                                   aTGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGAtt                                                <  1:3024428/58‑1 (MQ=255)
                                   aTGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATgcgg                                     <  1:1681334/69‑1 (MQ=255)
                                    tGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATgcggc                                    <  1:1625650/69‑1 (MQ=255)
                                     gTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAg                                                                    >  1:5472320/1‑36 (MQ=255)
                                     gTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGt                                  >  1:4160627/1‑70 (MQ=255)
                                       cAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGaa                               >  1:4671056/1‑71 (MQ=255)
                                        aaCGACAGTTGTAGGCCTGATAAGACGCGCAAGCGt                                                                 >  1:296514/1‑36 (MQ=255)
                                        aaCGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                              >  1:2545792/1‑71 (MQ=255)
                                          cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAg                                                         >  1:778063/1‑42 (MQ=255)
                                          cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGc                            >  1:688297/1‑71 (MQ=255)
                                          cGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGc                            <  1:483099/71‑1 (MQ=255)
                                              gtttGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCtt                         <  1:3890482/68‑1 (MQ=11)
                                              gtttGTAGGCCTGATAAGACGCGACAGCGTCGCATCAGGCATTGATTGCCGGATGCGGCGTAAACGCCtt                         <  1:4293478/68‑1 (MQ=11)
                                              aGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTa                        <  1:1025895/71‑1 (MQ=255)
                                               tttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGc                                                       <  1:3924670/38‑1 (MQ=37)
                                               gTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCg                                            <  1:4145213/50‑1 (MQ=255)
                                               gTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTAt                       <  1:5202251/71‑1 (MQ=255)
                                                ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCa                                                          >  1:3462187/1‑35 (MQ=37)
                                                ttGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGCAt                                         <  1:4536747/52‑1 (MQ=255)
                                                 tGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAAc                              >  1:2172949/1‑62 (MQ=255)
                                                  gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGTAAACGCCTTATCCg                    >  1:11775/1‑71 (MQ=14)
                                                  gTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGTAAACGCCTTATCCg                    >  1:384209/1‑71 (MQ=14)
                                                     ggCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGcc                 <  1:3706205/71‑1 (MQ=255)
                                                        cTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTAc              >  1:4226647/1‑71 (MQ=255)
                                                         tGATAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGCGGCGta                                 <  1:4217731/51‑2 (MQ=14)
                                                          gATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACg             >  1:5380770/1‑70 (MQ=255)
                                                           aTAAGACGCGCAAGCGTCGCATCAGGCGTTGGTTGCCGAATGAGGCGTAAACGCCTTAt                       >  1:243636/1‑59 (MQ=17)
                                                           aTAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGc            >  1:4815312/1‑70 (MQ=255)
                                                             aaGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAAt         >  1:1385256/1‑71 (MQ=255)
                                                              aGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATc        <  1:5636696/71‑1 (MQ=255)
                                                               gACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCg       >  1:5637974/1‑71 (MQ=255)
                                                                aCGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATgcggc                                    <  1:1660086/41‑1 (MQ=255)
                                                                aCGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGt      >  1:3671525/1‑71 (MQ=255)
                                                                 cgcgCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGtt     >  1:3079450/1‑71 (MQ=255)
                                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGa   >  1:1416237/1‑70 (MQ=255)
                                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGa   >  1:946860/1‑70 (MQ=255)
                                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGAc  <  1:5568931/71‑1 (MQ=255)
                                                                    gcAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGAc  <  1:4357513/71‑1 (MQ=255)
                                                                      |                                                                    
TGCCGGATGCGGCGTAAACGCCTTATCCGGCTACATGTCAACGACAGTTGTAGGCCTGATAAGACGCGCAAGCGTCGCATCAGGCGTTGATTGCCGGATGCGGCGTGAACGCCTTATTCGGCCTACGCAATCGTTGAC  >  minE/2053856‑2053993

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: