Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,291,126 G→A 66.3% 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,1260GA66.3% 96.6 / 84.6 86intergenic (+241/+212)yeeN/asnWconserved hypothetical protein/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base G (19/10);  new base A (40/17);  total (59/27)
Fisher's exact test for biased strand distribution p-value = 8.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAT  >  minE/1291056‑1291191
                                                                      |                                                                 
aaTTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTg                                                                   <  1:1265411/71‑1 (MQ=255)
aaTTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTg                                                                   <  1:1475287/71‑1 (MQ=255)
  ttACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                                 >  1:1547676/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATGCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:408330/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:667276/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:2694441/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:2135621/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:1965672/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:1823299/1‑71 (MQ=255)
    aCGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACg                                                               >  1:495848/1‑71 (MQ=255)
        aTAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                                 >  1:1317389/1‑65 (MQ=255)
        aTAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                                 >  1:2382735/1‑65 (MQ=255)
        aTAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTAc                                                           >  1:3643790/1‑71 (MQ=255)
        aTAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTAc                                                           >  1:748699/1‑71 (MQ=255)
        aTAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTAc                                                           >  1:1974617/1‑71 (MQ=255)
              tCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATg                                                     >  1:3244514/1‑71 (MQ=255)
              tCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATg                                                     >  1:975913/1‑71 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTg                                                                   <  1:792748/50‑1 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTg                                                                   <  1:3634673/50‑1 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAAc                                                 <  1:1493524/68‑1 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACa                                                >  1:1305607/1‑69 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACa                                                >  1:1300628/1‑69 (MQ=255)
                     aaaGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACa                                                >  1:3846695/1‑69 (MQ=255)
                        gACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGtt                                           <  1:1773113/71‑1 (MQ=255)
                            cgATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAAc                                                 <  1:1123017/61‑1 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTtata                                                                 >  1:910537/1‑41 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTtata                                                                 >  1:62312/1‑41 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCt                                   >  1:2315080/1‑71 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCt                                   >  1:2417961/1‑71 (MQ=255)
                                aGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCt                                   >  1:964480/1‑71 (MQ=255)
                                 gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGaa                                                  <  1:4127337/55‑1 (MQ=255)
                                 gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGaa                                                  <  1:4049130/55‑1 (MQ=255)
                                 gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGaa                                                  <  1:3956810/55‑1 (MQ=255)
                                 gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGaa                                                  <  1:118097/55‑1 (MQ=255)
                                 gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACa                                                >  1:4028246/1‑57 (MQ=255)
                                 gCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACag                                               <  1:570976/58‑1 (MQ=255)
                                   gACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACa                                                >  1:2468263/1‑55 (MQ=255)
                                      aTCCCGCATAAGGCATTTTTCTCTTTATCTTtata                                                                 <  1:1338151/35‑1 (MQ=255)
                                      aTCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                                 >  1:2251041/1‑35 (MQ=255)
                                      aTCCCGCATAAGGCATTTTTCTCTTTATCTTTGTa                                                                 >  1:3149422/1‑35 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACaga                                              >  1:693382/1‑51 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:1651620/1‑59 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:4162963/1‑59 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:1549680/1‑59 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATc                           <  1:3036371/70‑1 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATc                           <  1:3586956/70‑1 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATc                           <  1:2425107/70‑1 (MQ=255)
                                         ccGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCg                          >  1:465579/1‑71 (MQ=255)
                                          cGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACag                                               <  1:2997257/49‑1 (MQ=255)
                                          cGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                         <  1:782466/71‑1 (MQ=255)
                                          cGCATAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                         <  1:306034/71‑1 (MQ=255)
                                           gCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                         >  1:2315708/1‑70 (MQ=255)
                                            cATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACa                                                <  1:4071059/46‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                         <  1:3576574/67‑1 (MQ=255)
                                              tAAGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGt                         <  1:1464529/67‑1 (MQ=255)
                                               aaGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATg                     <  1:888699/70‑1 (MQ=255)
                                                aGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa                       >  1:1002186/1‑67 (MQ=255)
                                                aGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGa                       >  1:50942/1‑67 (MQ=255)
                                                aGGCATTTTTCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTATGACCTTGCGAATc                           >  1:567466/1‑63 (MQ=255)
                                                 ggCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                              >  1:2562156/1‑43 (MQ=255)
                                                 ggCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                              >  1:814322/1‑43 (MQ=255)
                                                 ggCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                              >  1:1852100/1‑43 (MQ=255)
                                                 ggCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACaga                                              >  1:1314578/1‑43 (MQ=255)
                                                     tttttCTCTTTATCTTTATACTGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg              <  1:813682/71‑1 (MQ=255)
                                                     tttttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg              <  1:1107682/71‑1 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:3289075/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:3284273/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:2320189/1‑45 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:3786699/1‑45 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:3782755/1‑45 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGa                                      >  1:335906/1‑45 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:2831977/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:3783608/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:450471/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:1116773/1‑71 (MQ=255)
                                                       tttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTgg             >  1:3270754/1‑70 (MQ=255)
                                                        ttCTCTTTATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            >  1:273033/1‑70 (MQ=255)
                                                            ctTTATCTTTGTTCGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGt                <  1:3293675/62‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            <  1:428406/64‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            <  1:3148782/64‑1 (MQ=255)
                                                              ttATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTggg            <  1:344838/64‑1 (MQ=255)
                                                               tATCTTTGTACGGTACTTCATGGAACAGAGttttt                                        >  1:1064843/1‑35 (MQ=255)
                                                               tATCTTTATACGGTACTTCATGGAACAGAGTTTTTGAc                                     >  1:638229/1‑38 (MQ=255)
                                                               tATCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGAc     >  1:3212911/1‑70 (MQ=255)
                                                                aTCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACaa   >  1:1280867/1‑71 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTg              <  1:3957600/59‑1 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACaa   <  1:2180286/70‑1 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACaa   <  1:188226/70‑1 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAt  >  1:3834142/1‑71 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAt  >  1:3432385/1‑71 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAt  >  1:3342984/1‑71 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAt  >  1:2192742/1‑71 (MQ=255)
                                                                 tCTTTATACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAt  >  1:1343734/1‑71 (MQ=255)
                                                                      |                                                                 
AATTACGCATAGTTTCCCGATAAAGACGCGATAGCGACATCCCGCATAAGGCATTTTTCTCTTTATCTTTGTACGGTACTTCATGGAACAGAGTTTTTGACCTTGCGAATCGTGATGTCTGTTGGGGAGGGACAAT  >  minE/1291056‑1291191

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: