Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,291,126 G→A 30.1% intergenic (+241/+212) yeeN → / ← asnW conserved hypothetical protein/tRNA‑Asn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,291,1260GA30.1% 122.9 / 67.5 83intergenic (+241/+212)yeeN/asnWconserved hypothetical protein/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base A (14/11);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

AATTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGG  >  minE/1291056‑1291181
                                                                      |                                                       
aaTTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTg                                                         <  1:1578804/71‑1 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACg                                                     >  1:1839434/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACg                                                     >  1:1351895/1‑71 (MQ=255)
         tAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACt                                                >  1:3950507/1‑71 (MQ=255)
              tCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCAt                                            <  1:2168246/70‑1 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACa                                      >  1:2174315/1‑69 (MQ=255)
                      aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                       >  1:65463/1‑51 (MQ=255)
                      aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                       >  1:2038542/1‑51 (MQ=255)
                      aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                       >  1:1944756/1‑51 (MQ=255)
                      aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACagag                                   >  1:2290692/1‑71 (MQ=255)
                      aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACagag                                   >  1:3525265/1‑71 (MQ=255)
                      aaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACagag                                   >  1:1813838/1‑71 (MQ=255)
                          cgcgATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGtttt                               <  1:1813433/71‑1 (MQ=255)
                              aTAGCGACATCCCGCATAAGTCATTTTTCTCTTTATGTTTGTACgg                                                    <  1:2095329/46‑1 (MQ=38)
                              aTAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACgg                                                    <  1:3382369/46‑1 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCt                         >  1:2561528/1‑71 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCt                         >  1:2537364/1‑71 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCt                         >  1:2153876/1‑71 (MQ=255)
                                    aCATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCg                     <  1:2978425/71‑1 (MQ=255)
                                       tCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTc                                              <  1:3131207/43‑1 (MQ=255)
                                       tCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTc                                              <  1:3050810/43‑1 (MQ=255)
                                       tCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTc                                              <  1:373787/43‑1 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                    >  1:334757/1‑51 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                    >  1:3958086/1‑51 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                    >  1:2388721/1‑51 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATc                 <  1:1032458/70‑1 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCg                >  1:2919797/1‑71 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCg                <  1:2130092/71‑1 (MQ=255)
                                          cGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:844462/71‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTc                                              <  1:526011/36‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATg           <  1:1537061/71‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATg           <  1:1139450/71‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACa                                      <  1:3969702/44‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACa                                      <  1:743652/44‑1 (MQ=255)
                                               aaGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACa                                      >  1:2583324/1‑43 (MQ=255)
                                               aaGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa             >  1:1781609/1‑68 (MQ=255)
                                               aaGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa             >  1:594269/1‑68 (MQ=255)
                                                  gCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTg       >  1:795521/1‑71 (MQ=255)
                                                  gCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTg       >  1:781663/1‑71 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCt                         <  1:1463682/50‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCt                         <  1:1257617/50‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCt                         <  1:816074/50‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCTAATCGTGATGTCTg       <  1:1009241/68‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTg       <  1:3648756/68‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTg       <  1:1948443/68‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTCCGAATCGTGATGTCTg       <  1:4129088/68‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGttttt                              <  1:2616622/45‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGttttt                              <  1:1231775/45‑1 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGttttt                              >  1:295317/1‑43 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGaa                   >  1:2039172/1‑54 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGaa                   >  1:1294825/1‑54 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGaa                   >  1:2144294/1‑54 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGaa                   >  1:148240/1‑54 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTg       >  1:3320159/1‑66 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg  >  1:1919531/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg  >  1:542329/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg  >  1:1814360/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   >  1:923750/1‑70 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   >  1:9553/1‑70 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   >  1:2312524/1‑70 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg    >  1:3984047/1‑69 (MQ=255)
                                                       tttCTCGTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   >  1:530072/1‑70 (MQ=255)
                                                        ttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa             >  1:877855/1‑59 (MQ=255)
                                                        ttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa             >  1:1594480/1‑59 (MQ=255)
                                                          ctctTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGaa                   >  1:2517301/1‑51 (MQ=255)
                                                             tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGa                            <  1:3881688/39‑1 (MQ=255)
                                                             tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGa                            <  1:1489079/39‑1 (MQ=255)
                                                             tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTTCGAATCGTGATGTCTGTTggg  <  1:2644757/65‑1 (MQ=255)
                                                             tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg  <  1:3885850/65‑1 (MQ=255)
                                                             tttATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg  <  1:1719792/65‑1 (MQ=255)
                                                             tttATCTTTATACGGTACTTCATGGAACAGAGttttt                              <  1:1253489/37‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   <  1:3755763/63‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   <  1:2840139/63‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   <  1:1625122/63‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg   <  1:3573052/63‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGTAACAGAGTTTTTGACCTTGCGAATCGt               <  1:2111514/46‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:2927902/46‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:580745/46‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:725430/46‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:2941569/46‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:493553/46‑1 (MQ=255)
                                                                   tttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:1048301/46‑1 (MQ=255)
                                                                    ttGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt               <  1:1446008/45‑1 (MQ=255)
                                                                     tGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa             >  1:2384039/1‑46 (MQ=255)
                                                                      |                                                       
AATTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGG  >  minE/1291056‑1291181

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