Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,269,969 G→A 5.5% A74T (GCA→ACA)  yedP → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,269,9690GA 5.5% 242.5 / 3.4 72A74T (GCA→ACA) yedPconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/24);  new base A (2/2);  total (46/26)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

ATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTATTAGCC  >  minE/1269899‑1270038
                                                                      |                                                                     
aTCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGgcg                                                                       <  1:3578257/71‑1 (MQ=255)
aTCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGgcg                                                                       <  1:683164/71‑1 (MQ=255)
      ggAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTg                                                                  >  1:2734132/1‑70 (MQ=255)
      ggAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGa                                                                 >  1:1219193/1‑71 (MQ=255)
      ggAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGa                                                                 <  1:68525/71‑1 (MQ=255)
      ggAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGa                                                                 <  1:2292135/71‑1 (MQ=255)
           tGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATcc                                                              >  1:630288/1‑69 (MQ=255)
            gCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGc                                                           >  1:871460/1‑71 (MQ=255)
            gCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGc                                                           >  1:3074721/1‑71 (MQ=255)
             cTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCt                                                          >  1:7456/1‑71 (MQ=255)
                tACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCtt                                                         <  1:3860380/69‑1 (MQ=255)
                  cTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCt                                                      >  1:3280054/1‑70 (MQ=255)
                  cTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCt                                                      >  1:1366953/1‑70 (MQ=255)
                  cTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTg                                                     <  1:1454780/71‑1 (MQ=255)
                     gCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAg                                                   >  1:3579579/1‑70 (MQ=255)
                     gCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAg                                                   >  1:266537/1‑70 (MQ=255)
                     gCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAg                                                   >  1:1457133/1‑70 (MQ=255)
                     gCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAg                                                   >  1:1980878/1‑70 (MQ=255)
                      cAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCa                                                                     >  1:3618209/1‑51 (MQ=255)
                      cAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCa                                                                     >  1:2484377/1‑51 (MQ=255)
                      cAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCa                                                                     >  1:2372094/1‑51 (MQ=255)
                        aaaaaCGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGa                                                                 <  1:2167531/53‑1 (MQ=255)
                         aaaaCGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATg                                              >  1:539934/1‑71 (MQ=255)
                             cGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAg                                          >  1:222229/1‑71 (MQ=255)
                             cGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAg                                          <  1:1584245/71‑1 (MQ=255)
                              gTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAgg                                         <  1:2731803/71‑1 (MQ=255)
                                  gggCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGgaga                                      >  1:3673649/1‑70 (MQ=255)
                                     cTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGa                                  >  1:3301009/1‑71 (MQ=255)
                                      tACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCtt                                                         <  1:2373376/47‑1 (MQ=255)
                                      tACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCtt                                                         <  1:2703312/47‑1 (MQ=255)
                                      tACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCtt                                                         <  1:1486573/47‑1 (MQ=255)
                                       aCAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCATGAGATAGACg                                >  1:85694/1‑71 (MQ=255)
                                       aCAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATa                                    <  1:637218/67‑1 (MQ=255)
                                       aCAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATa                                    <  1:1197097/67‑1 (MQ=255)
                                        cAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACgg                               <  1:2292731/71‑1 (MQ=255)
                                        cAAGGTTTACCGCTGATTGCCGAAAATGGCACAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACgg                               >  1:1323304/1‑71 (MQ=255)
                                        cAAGGTTTACCGCTGATTGCCGAAAATGGCACAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACgg                               >  1:689360/1‑71 (MQ=255)
                                         aaGGTTTACCGCTGATTGCCGATAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGa                                  >  1:2495758/1‑67 (MQ=255)
                                         aaGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGa                                  >  1:2626635/1‑67 (MQ=255)
                                         aaGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGa                                  >  1:2346876/1‑67 (MQ=255)
                                         aaGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGt                              >  1:3505969/1‑71 (MQ=255)
                                         aaGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGt                              >  1:2235761/1‑71 (MQ=255)
                                         aaGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGt                              >  1:2722383/1‑71 (MQ=255)
                                           ggTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGttt                            >  1:288714/1‑71 (MQ=255)
                                           ggTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGttt                            <  1:2627160/71‑1 (MQ=255)
                                            gTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGc                                                  <  1:1347766/48‑1 (MQ=255)
                                            gTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGc                                                  <  1:1253327/48‑1 (MQ=255)
                                            gTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGc                                                  <  1:2362431/48‑1 (MQ=255)
                                            gTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGtttt                           >  1:857117/1‑71 (MQ=255)
                                             tttATCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTc                          <  1:721262/71‑1 (MQ=255)
                                             tttACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTc                          <  1:2605183/71‑1 (MQ=255)
                                              ttACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTc                          >  1:875849/1‑70 (MQ=255)
                                               tACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCa                        >  1:305076/1‑71 (MQ=255)
                                                  cGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCAc                       >  1:2304905/1‑69 (MQ=255)
                                                    cTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcat                   >  1:2208738/1‑71 (MQ=255)
                                                    cTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcat                   >  1:482351/1‑71 (MQ=255)
                                                    cTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcat                   >  1:2725640/1‑71 (MQ=255)
                                                    cTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcat                   >  1:2314136/1‑71 (MQ=255)
                                                     tGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGc                                            >  1:1074100/1‑45 (MQ=255)
                                                      gATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcatca                 >  1:3632225/1‑71 (MQ=255)
                                                       aTTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcat                   <  1:2268728/68‑1 (MQ=255)
                                                        ttGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCatcatc               <  1:643888/71‑1 (MQ=255)
                                                         tGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACgg                               >  1:325937/1‑54 (MQ=255)
                                                         tGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACgg                               >  1:1070822/1‑54 (MQ=255)
                                                          gCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGTCAGGAGATAGACGGTTTTCCACGCATCAtct              <  1:501713/70‑1 (MQ=255)
                                                          gCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCAtct              <  1:56516/70‑1 (MQ=255)
                                                          gCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCAtct              >  1:2834133/1‑70 (MQ=255)
                                                             gAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAgg          >  1:2460879/1‑71 (MQ=255)
                                                             gAAAATGGCACAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAgg          <  1:827369/71‑1 (MQ=255)
                                                             gAAAATGGCACAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAgg          <  1:3602825/71‑1 (MQ=255)
                                                               aaaTGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGcat                   >  1:2443269/1‑60 (MQ=255)
                                                               aaaTGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTa        >  1:2534917/1‑71 (MQ=255)
                                                               aaaTGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTa        >  1:3789481/1‑71 (MQ=255)
                                                                aaTGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTAt       <  1:690861/71‑1 (MQ=255)
                                                                    gcgcAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTATTAGc   <  1:1970771/71‑1 (MQ=255)
                                                                     cgcAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACTGTTTTCCACGCTTCATCTCAGGTATTAGcc  <  1:640738/71‑1 (MQ=255)
                                                                     cgcAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTATTAGcc  <  1:2385726/71‑1 (MQ=255)
                                                                      |                                                                     
ATCAGCGGAAATGCTGTACTTGCAAAAAACGTTGGGGCTACAAGGTTTACCGCTGATTGCCGAAAATGGCGCAGTGATCCAGCTTGCTGAGCAATGGCAGGAGATAGACGGTTTTCCACGCATCATCTCAGGTATTAGCC  >  minE/1269899‑1270038

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: