Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,507,530 G→A 13.1% P198P (CCG→CCA yicI → predicted alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,507,5300GA13.1% 158.9 / 12.3 61P198P (CCG→CCAyicIpredicted alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base A (4/4);  total (33/28)
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

CAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTGGG  >  minE/2507465‑2507598
                                                                 |                                                                    
cAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTa                                                                 >  1:1830116/1‑71 (MQ=255)
  gACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACa                                                               >  1:19659/1‑71 (MQ=255)
           gaCCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAAcc                                                      <  1:1673012/71‑1 (MQ=255)
            aCCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCg                                                     >  1:3770262/1‑71 (MQ=255)
                ggAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTaa                                                        >  1:1550382/1‑64 (MQ=255)
                ggAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTa                                                         >  1:2749756/1‑63 (MQ=255)
                ggAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGt                                                 >  1:1198779/1‑71 (MQ=255)
                ggAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCAtt                                                                    >  1:550607/1‑52 (MQ=255)
                    ccGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTa                                               <  1:3395604/69‑1 (MQ=255)
                    ccGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATg                                             <  1:771164/71‑1 (MQ=255)
                     cGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATgg                                            >  1:1070326/1‑71 (MQ=255)
                     cGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATgg                                            >  1:2583280/1‑71 (MQ=255)
                     cGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATgg                                            >  1:207777/1‑71 (MQ=255)
                      gggACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTa                                                                 <  1:1656687/49‑1 (MQ=255)
                      gggACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTa                                                                 <  1:67742/49‑1 (MQ=255)
                        gACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGt                                         <  1:3566746/71‑1 (MQ=255)
                         aCGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTa                                        <  1:2859974/71‑1 (MQ=255)
                          cggcggAACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTAc                                       <  1:826651/71‑1 (MQ=255)
                             cggcACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTAc                                       >  1:2652255/1‑68 (MQ=255)
                                 acaAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCaa                                <  1:2808359/71‑1 (MQ=255)
                                 acaAGTACTGAACAGGCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTgg                                    >  1:3543238/1‑67 (MQ=255)
                                   aaGTACTGGACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAat                               >  1:3085507/1‑70 (MQ=255)
                                   aaGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAat                               >  1:385255/1‑70 (MQ=255)
                                   aaGTACTGAACAGGCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGt                                         <  1:1529933/60‑1 (MQ=255)
                                   aaGTACTGAACAGGCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAatc                              <  1:3318776/71‑1 (MQ=255)
                                        cTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATccc                         <  1:1112079/71‑1 (MQ=255)
                                         tGAACAGGCGTATAAAAATATCCCATTCTACATGACTAACCGTGGTTa                                               <  1:3414713/48‑1 (MQ=255)
                                          gAACAGGCGTATAAAAATATCCGGTTCTACATGACTAACCGTGGTTATg                                             <  1:778421/49‑1 (MQ=255)
                                          gAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGt                                   >  1:3678021/1‑59 (MQ=255)
                                           aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTg                                     <  1:3073620/56‑1 (MQ=255)
                                           aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTg                                     <  1:3193039/56‑1 (MQ=255)
                                           aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTg                                     <  1:1525552/56‑1 (MQ=255)
                                           aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                      >  1:1942221/1‑71 (MQ=255)
                                           aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                      >  1:3215584/1‑71 (MQ=255)
                                           aaCAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                      >  1:1447434/1‑71 (MQ=255)
                                             cAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtg                    >  1:38349/1‑71 (MQ=255)
                                             cAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtg                    >  1:3887555/1‑71 (MQ=255)
                                             cAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtg                    >  1:60231/1‑71 (MQ=255)
                                             cAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtg                    >  1:2279987/1‑71 (MQ=255)
                                              aGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAg                      >  1:322360/1‑68 (MQ=255)
                                               ggCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAgtgt                   <  1:1116744/70‑1 (MQ=255)
                                                  gTATAAAAATATCCCATTCTACATGACTAACCGTGGtt                                                >  1:2532668/1‑38 (MQ=255)
                                                   tataAAAATATCCCGTTCTACATGACTAACCGTCGTTATGGGGTACTGGTCAATCATCCCCAGTGTGtctc              <  1:1000515/71‑1 (MQ=255)
                                                    ataAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTctct             >  1:2136915/1‑71 (MQ=255)
                                                    ataAAAATATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTctct             >  1:766891/1‑71 (MQ=255)
                                                     taAAAATATCCCGTTCTACATGACTAACCGTGGTTa                                               >  1:2752612/1‑36 (MQ=255)
                                                     taAAAATATCCCGTTCTACATGACTAACCGTGGTTa                                               >  1:935913/1‑36 (MQ=255)
                                                     taAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTctct             <  1:298671/70‑1 (MQ=255)
                                                      aaaaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAat                               <  1:3877409/51‑1 (MQ=255)
                                                       aaaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCtttt          >  1:3683763/1‑71 (MQ=255)
                                                       aaaaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCttt           <  1:1846123/70‑1 (MQ=255)
                                                       aaaaTATCCCATTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCttt           <  1:3514730/70‑1 (MQ=255)
                                                         aaTATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTg         >  1:1492582/1‑70 (MQ=255)
                                                          atatCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGaa       <  1:961817/71‑1 (MQ=255)
                                                           tatCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAg      <  1:3078695/71‑1 (MQ=255)
                                                             tCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTg    >  1:2075390/1‑71 (MQ=255)
                                                             tCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTg    >  1:3511623/1‑71 (MQ=255)
                                                             tCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCAAGTGTGTCTCTTTTGAAGTg    >  1:2750039/1‑71 (MQ=255)
                                                               ccGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTggg  <  1:3595895/71‑1 (MQ=255)
                                                                cGTTCTACATGACTAACCGTGGTTATGGGGTACTg                                     <  1:1885506/35‑1 (MQ=255)
                                                                cGTTCTACATGACTAACCGTGGTTATGGGGTACTg                                     <  1:3172396/35‑1 (MQ=255)
                                                                 |                                                                    
CAGACGGTAGAGACCTGGAACCGGGACGGCGGCACAAGTACTGAACAGGCGTATAAAAATATCCCGTTCTACATGACTAACCGTGGTTATGGGGTACTGGTCAATCATCCCCAGTGTGTCTCTTTTGAAGTGGG  >  minE/2507465‑2507598

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: