Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,522,857 G→A 22.6% intergenic (‑257/+33) yicG ← / ← dinD conserved inner membrane protein/DNA‑damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,522,8570GA22.6% 88.8 / 17.7 66intergenic (‑257/+33)yicG/dinDconserved inner membrane protein/DNA‑damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (39/12);  new base A (4/11);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 1.25e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

GAAAAGAAAGGTGCAATATTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTG  >  minE/2522791‑2522917
                                                                  |                                                            
gAAAAGAAAGGTGCAATATTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGgtgt                                                             >  1:2263482/1‑68 (MQ=255)
 aaaaGAAAGGTGCAATGTTTAGCCATTGTGTTCACCGCACCTGGACTTTGGCTTACGCTTCGGGTAtc                                                            >  1:3116391/1‑68 (MQ=255)
 aaaaGAAAGGTGCAATGTTTAGCCATTGTGTTCACCGCACCTGGACTTTGGCTTACGCATCGGGTAtc                                                            >  1:649754/1‑68 (MQ=255)
 aaaaGAAAGGTGCAATGTTTAGCCATTGTGTTCACCGCACCTGGACTTTGGCTTACGCATCGGGTAtc                                                            >  1:3382842/1‑68 (MQ=255)
 aaaaGAAAGGTGCAATGTTTAGCCATTGTGTTCACCGCACCTGGACTTTGGCTTACGCATCGGGTAtc                                                            >  1:3298027/1‑68 (MQ=255)
    agaaagGTGCAATATTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCtca                                                         >  1:1277971/1‑68 (MQ=255)
    agaaagGTGCAATATTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCtca                                                         >  1:1132680/1‑68 (MQ=255)
    agaaagGTGCAATATTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCtca                                                         >  1:2158911/1‑68 (MQ=255)
                 atTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAagag                                           >  1:1245025/1‑69 (MQ=255)
                 atTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAagag                                           >  1:858198/1‑69 (MQ=255)
                 atTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAagag                                           >  1:3089638/1‑69 (MQ=255)
                 atTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAagag                                           >  1:1511296/1‑69 (MQ=255)
                 atTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAagag                                           >  1:974707/1‑69 (MQ=255)
                 atTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAagag                                           >  1:1730886/1‑69 (MQ=255)
                  tttAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATGAACAAAAAGAGt                                          >  1:708693/1‑69 (MQ=255)
                  tttAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGt                                          >  1:332061/1‑69 (MQ=255)
                  tttAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGt                                          >  1:1176281/1‑69 (MQ=255)
                  tttAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGt                                          >  1:390607/1‑69 (MQ=255)
                  tttAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGt                                          >  1:1624591/1‑69 (MQ=255)
                  tttAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGt                                          >  1:3357291/1‑69 (MQ=255)
                            gtgtTCACTGCATCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTc                                >  1:3038572/1‑69 (MQ=255)
                            gtgtTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGttt                                 >  1:3032581/1‑68 (MQ=255)
                            gtgtTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGttt                                 >  1:2745786/1‑68 (MQ=255)
                            gtgtTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGttt                                 >  1:582190/1‑68 (MQ=255)
                            gtgtTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGttt                                 >  1:2932994/1‑68 (MQ=255)
                            gtgtTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTc                                >  1:3161488/1‑69 (MQ=255)
                                     gCACCTGGACCTTGGTTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATttc                       >  1:539093/1‑69 (MQ=255)
                                     gCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACaaaa                                               >  1:2298076/1‑45 (MQ=255)
                                     gCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTAttt                        >  1:271873/1‑68 (MQ=255)
                                     gCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTAttt                        >  1:1012821/1‑68 (MQ=255)
                                     gCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTAttt                        >  1:370227/1‑68 (MQ=255)
                                     gCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATttc                       >  1:2259518/1‑69 (MQ=255)
                                     gCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATttc                       >  1:1143303/1‑69 (MQ=255)
                                        ccTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATttctt                     <  1:1930281/68‑1 (MQ=255)
                                        ccTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATttctt                     <  1:1310102/68‑1 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:382640/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:423585/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:568030/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:3081365/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:1071366/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:1107570/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:1583656/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:1257998/1‑68 (MQ=255)
                                                  ggCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtt           >  1:1275857/1‑68 (MQ=255)
                                                    cTTACGCATGGGGTATCTCAATAAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:1884264/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtttt         >  1:636251/1‑68 (MQ=255)
                                                    cTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAAtttt         >  1:3094068/1‑68 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:747482/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:1593352/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:1594550/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:1614182/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:2900375/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:228189/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:2303015/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:2627336/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:2710698/52‑1 (MQ=255)
                                                    cTTACGCATCGGGTATCTCAATCAACAAAAAGAGTTTTTATTTTCTGCTAtt                         <  1:2859030/52‑1 (MQ=255)
                                                          cATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:1799502/69‑1 (MQ=255)
                                                          cATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:828228/69‑1 (MQ=255)
                                                          cATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:678664/69‑1 (MQ=255)
                                                          cATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:1912044/69‑1 (MQ=255)
                                                          cATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:2388761/69‑1 (MQ=255)
                                                          cATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:2393439/69‑1 (MQ=255)
                                                              gggTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTTATTTTTACACTg  <  1:1903425/65‑1 (MQ=255)
                                                              gggTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:1317888/65‑1 (MQ=255)
                                                              gggTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:73684/65‑1 (MQ=255)
                                                              gggTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTg  <  1:3054090/65‑1 (MQ=255)
                                                                  |                                                            
GAAAAGAAAGGTGCAATATTTAGCCATTGTGTTCACTGCACCTGGACCTTGGCTTACGCATCGGGTGTCTCAATCAACAAAAAGAGTTTTTAGTTTCTGCTATTTCTTCTCTGTAATTTTTACACTG  >  minE/2522791‑2522917

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: