Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 217,996 A→T 18.5% intergenic (‑88/‑18) yafU ← / → yafF predicted inner membrane protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE217,9960AT18.5% 128.6 / 26.1 60intergenic (‑88/‑18)yafU/yafFpredicted inner membrane protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (30/18);  new base T (7/4);  total (37/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.60e-01

CTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTT  >  minE/217928‑218065
                                                                    |                                                                     
ctgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAgg                                                                     >  1:2937070/1‑71 (MQ=255)
 tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                                                    >  1:3919022/1‑71 (MQ=255)
 tgATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGc                                                                    >  1:2751679/1‑71 (MQ=255)
   aTTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGt                                                                  <  1:3954901/71‑1 (MQ=255)
          cGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACg                                                           <  1:2151764/71‑1 (MQ=255)
                    ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                                                 >  1:3326733/1‑71 (MQ=255)
                    ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                                                 >  1:975920/1‑71 (MQ=255)
                    ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                                                 >  1:907351/1‑71 (MQ=255)
                    ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                                                 >  1:4029888/1‑71 (MQ=255)
                    ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGTaatgaa                                                 >  1:3499276/1‑71 (MQ=255)
                    ttCGCCACAGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTGGACGTGGGaatgaa                                                 >  1:3503472/1‑71 (MQ=255)
                      cGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatgaat                                                >  1:4132721/1‑70 (MQ=255)
                        ccACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaatga                                                  <  1:3809680/66‑1 (MQ=255)
                           caGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGc                                                                    >  1:1450720/1‑45 (MQ=255)
                            aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                                               <  1:2076933/49‑1 (MQ=255)
                            aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                                               <  1:2109642/49‑1 (MQ=255)
                            aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                                               <  1:3250859/49‑1 (MQ=255)
                            aGCGATCCGAAATCACTGGCACGTGGAGAATAAGCTGCACTGGCGTCTg                                                               <  1:43725/49‑1 (MQ=255)
                                 tCCGAAACCACTGGCACGTGGAGAATAAGCTGCACCGGCGTCTGGACGTGGTAATGAATGAAGACGACTAc                                    <  1:1041421/71‑1 (MQ=255)
                                  ccGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTAc                                    >  1:1504268/1‑70 (MQ=255)
                                   cGAAACCACTGGCACGTGGAGAATAAGCTTCACAGGCGTCTGGACGTgg                                                        >  1:883595/1‑49 (MQ=255)
                                   cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                                                        >  1:1376534/1‑49 (MQ=255)
                                   cGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTgg                                                        >  1:73805/1‑49 (MQ=255)
                                         cACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAAt                               >  1:541836/1‑68 (MQ=255)
                                           cTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATaagaa                          <  1:1397145/71‑1 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:5087936/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:4980495/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:2602345/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:3158220/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:1258041/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:791585/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:1470093/1‑37 (MQ=255)
                                             ggCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGt                                                          >  1:3970293/1‑37 (MQ=255)
                                              gCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAagag                       >  1:4987817/1‑71 (MQ=255)
                                                aCGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                     <  1:5194263/71‑1 (MQ=255)
                                                 cGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACa                                   <  1:2859138/56‑1 (MQ=255)
                                                 cGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                     <  1:1057529/70‑1 (MQ=255)
                                                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTaa                                                     >  1:3890336/1‑36 (MQ=255)
                                                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaaa                   >  1:250428/1‑70 (MQ=255)
                                                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaaa                   >  1:233039/1‑70 (MQ=255)
                                                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAAc                  <  1:1305037/71‑1 (MQ=255)
                                                   tGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAAc                  <  1:3019170/71‑1 (MQ=255)
                                                      agaATAAGCTGCACAGGCGTCTGGACGTGGTaatg                                                   <  1:2092532/35‑1 (MQ=255)
                                                      agaATAAGCTGCACAGGCGTCTGGACGTGGTaatg                                                   <  1:146182/35‑1 (MQ=255)
                                                      agaATAAGCTGCACAGGCGTCTGGACGTGGTaatg                                                   <  1:567016/35‑1 (MQ=255)
                                                       gaATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag              >  1:2305783/1‑71 (MQ=255)
                                                        aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                     <  1:1058646/63‑1 (MQ=255)
                                                        aaTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag              >  1:955827/1‑70 (MQ=255)
                                                         aTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaa                    <  1:2866950/63‑1 (MQ=255)
                                                         aTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                     <  1:391354/62‑1 (MQ=255)
                                                         aTAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGa                     <  1:444470/62‑1 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGAGGACTACAAAATAAGAAGAGGAAACgcag              >  1:1385863/1‑68 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGaaa                   >  1:2072342/1‑63 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAAc                  >  1:4552892/1‑64 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag              >  1:2820525/1‑68 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag              >  1:314078/1‑68 (MQ=255)
                                                          tAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACgcag              >  1:97438/1‑68 (MQ=255)
                                                            aGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGaa         >  1:3219740/1‑71 (MQ=255)
                                                              cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatt       <  1:506999/71‑1 (MQ=255)
                                                              cTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatt       <  1:4167812/71‑1 (MQ=255)
                                                                 cacaGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAatta      <  1:4206903/69‑1 (MQ=255)
                                                                    aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt  >  1:52085/1‑70 (MQ=255)
                                                                    aGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTAtttt  >  1:8537/1‑70 (MQ=255)
                                                                    |                                                                     
CTGATTTAACCGCAGAGAAGTTCGCCACAGCAATCCGAAACCACTGGCACGTGGAGAATAAGCTGCACAGGCGTCTGGACGTGGTAATGAATGAAGACGACTACAAAATAAGAAGAGGAAACGCAGCAGAATTATTTT  >  minE/217928‑218065

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: