Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 335,162 G→A 78.3% G130G (GGC→GGT ybaK ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE335,1620GA78.3% 116.7 / 32.0 60G130G (GGC→GGTybaKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (5/8);  new base A (17/30);  total (22/38)
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.82e-01

AAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGT  >  minE/335096‑335222
                                                                  |                                                            
aaaTTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGccaccg                                                           >  1:2195544/1‑70 (MQ=255)
       gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCTATATCCAGTCCGCGCTTGCCACCGGAAACAt                                                    <  1:2608908/70‑1 (MQ=255)
       gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa                                                   <  1:2166512/71‑1 (MQ=255)
       gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACAt                                                    <  1:1219997/70‑1 (MQ=255)
       gCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACAt                                                    <  1:354985/70‑1 (MQ=255)
        cATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa                                                   >  1:4053717/1‑70 (MQ=255)
        cATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa                                                   >  1:735941/1‑70 (MQ=255)
          tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt                                                <  1:2848182/71‑1 (MQ=255)
          tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt                                                <  1:2924940/71‑1 (MQ=255)
          tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt                                                <  1:2959134/71‑1 (MQ=255)
          tCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt                                                <  1:3297988/71‑1 (MQ=255)
          tCGAGGATCTTTGCCAGATCGCCTGACGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAAt                                                <  1:1138807/71‑1 (MQ=255)
                aTCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGaa                                                        <  1:61388/57‑1 (MQ=255)
                     tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt                                      <  1:2532597/70‑1 (MQ=255)
                     tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt                                      <  1:3610426/70‑1 (MQ=255)
                     tGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt                                      <  1:3178271/70‑1 (MQ=255)
                      gCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAtt                                     <  1:1378834/70‑1 (MQ=255)
                      gCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAtt                                     <  1:3481334/70‑1 (MQ=255)
                       ccAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCt                                   <  1:1803517/71‑1 (MQ=255)
                       ccAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCt                                   <  1:305597/71‑1 (MQ=255)
                        cAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCtt                                  <  1:166438/71‑1 (MQ=255)
                        cAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCtt                                  <  1:180308/71‑1 (MQ=255)
                          gATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTtgt                                <  1:1450490/71‑1 (MQ=255)
                             cGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:3282782/1‑56 (MQ=255)
                             cGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:156035/1‑56 (MQ=255)
                              gCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCgg                            >  1:3481094/1‑71 (MQ=255)
                              gCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCgg                            >  1:1205851/1‑71 (MQ=255)
                                  gccgccAGTTCGATATCCATTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc                                    <  1:1421727/59‑1 (MQ=255)
                                  gccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc                                    <  1:442779/59‑1 (MQ=255)
                                  gccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc                                    <  1:3427281/59‑1 (MQ=255)
                                   ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAg                                              >  1:87097/1‑48 (MQ=255)
                                   ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAg                                              >  1:3260895/1‑48 (MQ=255)
                                   ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAg                                              >  1:843015/1‑48 (MQ=255)
                                   ccgccAGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc                                    <  1:3218345/58‑1 (MQ=255)
                                        aGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGAt                  <  1:462743/71‑1 (MQ=255)
                                        aGTTCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGAt                  <  1:1674630/71‑1 (MQ=255)
                                         gTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTc                                    <  1:2464489/52‑1 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa                                                   <  1:286851/36‑1 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCGCCGGAAACATa                                                   <  1:574340/36‑1 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:3797925/1‑43 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:1249907/1‑43 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:1789061/1‑43 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:405130/1‑43 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:1881157/1‑43 (MQ=255)
                                          ttCGATATCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTg                                            >  1:188440/1‑43 (MQ=255)
                                              atatCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCg             <  1:3982755/70‑1 (MQ=255)
                                                atCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTg          >  1:3326903/1‑71 (MQ=255)
                                                 tCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGt                                             <  1:1441922/35‑1 (MQ=255)
                                                 tCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGt                                             <  1:2511775/35‑1 (MQ=255)
                                                 tCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCg             <  1:896342/67‑1 (MQ=255)
                                                 tCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCg             <  1:1746902/67‑1 (MQ=255)
                                                 tCCAGTCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGt            <  1:1312870/68‑1 (MQ=255)
                                                      tCCGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt                                      <  1:335603/37‑1 (MQ=255)
                                                      tCCGCGCTTGCCACCGGAAAAATAAATAGTGGCAAAt                                      <  1:3291729/37‑1 (MQ=38)
                                                       ccGCGCTTGCCACCGGAAACATAAATAGTGGCAAAt                                      <  1:17040/36‑1 (MQ=255)
                                                       ccGCGCTTGCCACCGGAAACATAAATAGTGGCAAATTc                                    >  1:2394663/1‑38 (MQ=255)
                                                        cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGc        >  1:1420698/1‑65 (MQ=255)
                                                        cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGt  >  1:2693409/1‑71 (MQ=255)
                                                        cgcgCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGt  >  1:2344861/1‑71 (MQ=255)
                                                            cTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGAc    >  1:1972735/1‑65 (MQ=255)
                                                            cTTGCCACCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGAc    >  1:1454264/1‑65 (MQ=255)
                                                                  |                                                            
AAATTTGGCATCGAGGATCTTTGCCAGATCGCCTGCCGCCAGTTCGATATCCAGTCCGCGCTTGCCGCCGGAAACATAAATAGTGGCAAATTCTTGTGCGGGGGCGTCGATAATCGTTGGCAGACGT  >  minE/335096‑335222

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: