Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 83,694 G→A 31.1% V184V (GTG→GTA ftsW → integral membrane protein involved in stabilising FstZ ring during cell division

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE83,6940GA31.1% 87.1 / 44.2 93V184V (GTG→GTAftsWintegral membrane protein involved in stabilising FstZ ring during cell division
Reads supporting (aligned to +/- strand):  ref base G (13/51);  new base A (16/13);  total (29/64)
Fisher's exact test for biased strand distribution p-value = 1.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

CTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGTTGTTTGT  >  minE/83628‑83759
                                                                  |                                                                 
cTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:3426497/69‑1 (MQ=255)
cTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:3212628/69‑1 (MQ=255)
cTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:2417961/69‑1 (MQ=255)
cTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:1587080/69‑1 (MQ=255)
cTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:982686/69‑1 (MQ=255)
 tATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAtt                                                                >  1:2425292/1‑69 (MQ=255)
 tATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 >  1:2617732/1‑68 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:102394/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:3383124/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:3279966/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:3248239/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:419091/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:2842773/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:2777749/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:531276/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:2623797/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:2483811/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:110258/68‑1 (MQ=255)
     tGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:1889598/68‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGTGCGTGATTCTGgt                                                           <  1:192588/66‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTgg                                                            <  1:2556807/65‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTgg                                                            <  1:2713196/65‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTgg                                                            <  1:705306/65‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGgt                                                           <  1:2792252/66‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGgt                                                           <  1:1795893/66‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGgt                                                           <  1:1602821/66‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGgt                                                           <  1:1943121/66‑1 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:103583/1‑69 (MQ=255)
         gCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:588326/1‑69 (MQ=255)
              aaGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAg                                                   <  1:1660789/69‑1 (MQ=255)
              aaGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAg                                                   <  1:2376137/69‑1 (MQ=255)
                 gCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:2333450/69‑1 (MQ=255)
                 gCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:1219594/69‑1 (MQ=255)
                 gCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:2991144/69‑1 (MQ=255)
                 gCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:2913704/69‑1 (MQ=255)
                 gCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:1310816/69‑1 (MQ=255)
                 gCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:1770570/69‑1 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTAct                                            <  1:466998/69‑1 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGtt                                                        >  1:2922251/1‑57 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:3328577/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:3201900/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:3123836/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:3078373/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:3019221/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:1660104/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:2897014/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:2809999/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:1520368/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:352486/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:3354003/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:1170190/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:685668/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:538110/1‑69 (MQ=255)
                     cgaAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTGGTGTTGGCAGTGTTAct                                            >  1:1071303/1‑69 (MQ=255)
                        aGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTAATTCTg                                                             <  1:535107/49‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGTGCGTGAt                                                                 <  1:2601185/44‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:3317678/44‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:1378666/44‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:3426335/44‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:2325853/44‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGAt                                                                 <  1:2173515/44‑1 (MQ=255)
                         gTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGc                                                     >  1:2228352/1‑56 (MQ=255)
                          tACGTAATAACCTGCGCGGCTTCCTGAAACCTATGGGCGTGAt                                                                 <  1:1856804/43‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:3198012/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:2915630/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:489911/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:1243884/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:1239246/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               >  1:1704087/1‑69 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:1495625/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:1531007/69‑1 (MQ=255)
                                  aaCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               <  1:892781/69‑1 (MQ=255)
                                   aCCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCa                               >  1:1530881/1‑68 (MQ=255)
                                        cgcgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:425008/1‑38 (MQ=255)
                                        cgcgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:1035816/1‑38 (MQ=255)
                                        cgcgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:3330823/1‑38 (MQ=255)
                                        cgcgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:1212535/1‑38 (MQ=255)
                                        cgcgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:2448682/1‑38 (MQ=255)
                                        cgcgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGtt                                                        >  1:1493712/1‑38 (MQ=255)
                                          cgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTgg                       <  1:435421/69‑1 (MQ=255)
                                          cgGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTgg                       <  1:842210/69‑1 (MQ=255)
                                           ggCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTgg                       <  1:3163986/68‑1 (MQ=255)
                                              ttCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:3247294/40‑1 (MQ=255)
                                              ttCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:2832931/40‑1 (MQ=255)
                                              ttCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:2716956/40‑1 (MQ=255)
                                              ttCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:58749/40‑1 (MQ=255)
                                              ttCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:1639684/40‑1 (MQ=255)
                                              ttCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAgtgt                                                <  1:838939/40‑1 (MQ=255)
                                                               cGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGttgtttgt  <  1:1113959/69‑1 (MQ=255)
                                                               cGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGttgtttgt  <  1:1549517/69‑1 (MQ=255)
                                                               cGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGttgtttgt  <  1:2616324/69‑1 (MQ=255)
                                                               cGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGttgtttgt  <  1:2854508/69‑1 (MQ=255)
                                                               cGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGttgtttgt  <  1:360385/69‑1 (MQ=255)
                                                                  |                                                                 
CTATCTGGTGCGTAAAGGCGACGAAGTACGTAATAACCTGCGCGGCTTCCTGAAACCGATGGGCGTGATTCTGGTGTTGGCAGTGTTACTGCTGGCACAGCCAGACCTTGGTACGGTGGTGGTGTTGTTTGT  >  minE/83628‑83759

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: