Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,642,784 T→A 10.0% P74P (CCA→CCT ypfJ ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,642,7840TA10.0% 136.5 / 7.8 50P74P (CCA→CCTypfJconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/21);  new base A (3/2);  total (27/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 = 8.52e-01

GTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAA  >  minE/1642720‑1642851
                                                                |                                                                   
gtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCt                                                                <  1:89460/70‑1 (MQ=255)
 tCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTa                                                               <  1:2121330/70‑1 (MQ=255)
      aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAc                                                         >  1:881544/1‑71 (MQ=255)
      aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAc                                                         >  1:2226175/1‑71 (MQ=255)
         tgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                      <  1:943473/71‑1 (MQ=255)
         tgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                      <  1:859839/71‑1 (MQ=255)
         tgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                      <  1:252180/71‑1 (MQ=255)
          gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                      <  1:1612845/70‑1 (MQ=255)
             ttCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGtt                                                   >  1:444033/1‑70 (MQ=255)
                    gTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAt                                                >  1:961851/1‑66 (MQ=255)
                    gTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAt                                                >  1:1425229/1‑66 (MQ=255)
                    gTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAt                                                >  1:3324329/1‑66 (MQ=255)
                     ttGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgt                                          <  1:2376230/71‑1 (MQ=255)
                      tGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGAcc                                                        >  1:3216675/1‑56 (MQ=255)
                      tGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                 >  1:2827736/1‑63 (MQ=255)
                      tGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgt                                          >  1:1070948/1‑70 (MQ=255)
                       gCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTaa                                                              >  1:3745499/1‑49 (MQ=255)
                       gCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                        >  1:2448164/1‑71 (MQ=255)
                       gCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgtt                                         <  1:1704458/70‑1 (MQ=255)
                        ccAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCt                                                                <  1:2254743/46‑1 (MQ=255)
                        ccAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCt                                                                <  1:2119749/46‑1 (MQ=255)
                        ccAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGAttgttgtt                                         <  1:855793/69‑1 (MQ=255)
                         cAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTgg                                       <  1:3752530/70‑1 (MQ=255)
                         cAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTgg                                       <  1:824030/70‑1 (MQ=255)
                           gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTGGGa                                     >  1:351067/1‑70 (MQ=255)
                           gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTGGGa                                     >  1:3646630/1‑70 (MQ=255)
                           gAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTAGGGCTAATTGACCGCGTTGATTGCTGTTGGGa                                     >  1:2135943/1‑70 (MQ=255)
                                aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACcgc                                                      <  1:1282692/48‑1 (MQ=255)
                                aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCgg                               >  1:853570/1‑71 (MQ=255)
                                 ccGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGCCCGCGTTGAttgttgt                                          >  1:2566813/1‑59 (MQ=255)
                                     ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGcc                          <  1:1482181/71‑1 (MQ=255)
                                      gTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGccc                         <  1:2079174/71‑1 (MQ=255)
                                        gAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCg                        >  1:224524/1‑70 (MQ=255)
                                             ttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                   >  1:3421468/1‑70 (MQ=255)
                                             ttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGtcat                   >  1:2473393/1‑70 (MQ=255)
                                              tGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGt                      >  1:847699/1‑66 (MQ=255)
                                                  gCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGa                                                 <  1:1593337/35‑1 (MQ=255)
                                                    ttCGTCTTCATTTGGGCTAATTGACCGCGTTGAtt                                               >  1:573742/1‑35 (MQ=255)
                                                     tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCg            <  1:1826230/69‑1 (MQ=255)
                                                     tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGt          >  1:3695849/1‑71 (MQ=255)
                                                          ttCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                        >  1:3652310/1‑36 (MQ=255)
                                                          ttCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                        >  1:1251970/1‑36 (MQ=255)
                                                          ttCATTTGGGCTAATTGACCGCGTTGAttgttgttg                                        >  1:3326141/1‑36 (MQ=255)
                                                          ttCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCg                                <  1:3158454/44‑1 (MQ=255)
                                                           tCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAAt     >  1:1432580/1‑70 (MQ=255)
                                                            cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGc                              >  1:779590/1‑44 (MQ=255)
                                                            cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc    <  1:2681094/70‑1 (MQ=255)
                                                            cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc    <  1:3492781/70‑1 (MQ=255)
                                                            cATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATc    <  1:1098881/70‑1 (MQ=255)
                                                             aTTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCaa  <  1:922170/71‑1 (MQ=255)
                                                              tttGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTg                            <  1:2821461/44‑1 (MQ=255)
                                                                |                                                                   
GTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAA  >  minE/1642720‑1642851

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: