Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,392,402 A→G 42.3% I116I (ATT→ATC recQ ← ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,392,4020AG42.3% 12.2 / 37.4 52I116I (ATT→ATCrecQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base A (7/23);  new base G (13/9);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 1.14e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTCGCGGGTTTGCGTCGAG  >  minE/2392356‑2392462
                                              |                                                            
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:3230051/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1922180/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1530825/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:429842/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1267342/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:386287/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1165246/68‑1 (MQ=255)
gAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1072307/68‑1 (MQ=255)
 aaaGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1161409/67‑1 (MQ=255)
    gTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:3134855/69‑1 (MQ=255)
    gTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:2125211/69‑1 (MQ=255)
    gTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:829647/69‑1 (MQ=255)
    gTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:582400/69‑1 (MQ=255)
    gTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:1834088/69‑1 (MQ=255)
       aTCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:2966266/66‑1 (MQ=255)
            gCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCTGTGCGGCAGCCTGTCATCACTTCa                             <  1:2667319/68‑1 (MQ=255)
            gCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCa                             <  1:1374844/68‑1 (MQ=255)
              aTCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGt                          <  1:1898191/69‑1 (MQ=255)
                     gTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCt                                         <  1:1521219/47‑1 (MQ=255)
                        ccGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGtt                         >  1:103139/1‑60 (MQ=255)
                        ccGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTcgc                >  1:1839030/1‑69 (MQ=255)
                        ccGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTcg                 >  1:1438693/1‑68 (MQ=255)
                            ggCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGtcat                                    <  1:3183976/45‑1 (MQ=255)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:625862/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:1709328/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:850148/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:2216265/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:2297142/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:533546/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:2469293/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:125194/44‑1 (MQ=37)
                                   taAAGCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTc                              <  1:1051625/44‑1 (MQ=37)
                                   taAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTg                     <  1:3291360/53‑1 (MQ=255)
                                   taAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTg                     <  1:1027922/53‑1 (MQ=255)
                                   taAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTg                     <  1:3150245/53‑1 (MQ=255)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:932623/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:3272475/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:915762/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:446186/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:2834064/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:2637894/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:2300373/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:2073011/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:1773536/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:1663919/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:166365/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:1549638/1‑43 (MQ=37)
                                       gCAGGCGGATTTGCCCGGTGCGGCAGCCTGTCATTACTTCAAgc                          >  1:1199863/1‑43 (MQ=37)
                                        cAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTCGCGGGTTTGCGTCGa   >  1:2872380/1‑66 (MQ=255)
                                        cAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTCGCGGGTTTGCGTCGAg  >  1:3048752/1‑67 (MQ=255)
                                        cAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTCGCGGGTTTGCGTCGAg  >  1:3283738/1‑67 (MQ=255)
                                        cAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTCGCGGGTTTGCGTCGAg  >  1:2780692/1‑67 (MQ=255)
                                              |                                                            
GAAAGTTATCCAGCATCAGGCGTTCCGGGGCGATATAAAGCAGACGAATTTGCCCGGTGCGGCAGCCTGTCATCACTTCAAGTTGCTGTTCGCGGGTTTGCGTCGAG  >  minE/2392356‑2392462

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