Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,399,893 G→A 48.2% V246V (GTC→GTT uvrD ← DNA‑dependent ATPase I and helicase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,399,8930GA48.2% 0.9 / 46.7 56V246V (GTC→GTTuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (18/11);  new base A (18/9);  total (36/20)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GAAACGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCG  >  minE/2399830‑2399941
                                                               |                                                
gAAACGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACcgacga                                              >  1:2585958/1‑68 (MQ=255)
gAAACGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGat                                             >  1:1031458/1‑69 (MQ=255)
gAAACGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGat                                             >  1:2673898/1‑69 (MQ=255)
gAAACGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGat                                             >  1:3237309/1‑69 (MQ=255)
 aaaCGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGatc                                            >  1:1318700/1‑69 (MQ=255)
 aaaCGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGat                                             >  1:523266/1‑68 (MQ=255)
      cTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCAc                                       <  1:2238317/69‑1 (MQ=255)
      cTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCAc                                       <  1:1500292/69‑1 (MQ=255)
      cTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCAc                                       <  1:1603919/69‑1 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:2535383/1‑69 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:3047816/1‑69 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:385326/1‑69 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:55108/1‑69 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:736563/1‑69 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:1037452/1‑69 (MQ=255)
           tatTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGc                                  >  1:908462/1‑69 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGcc                                 <  1:1637879/68‑1 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGcc                                 <  1:34626/68‑1 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGcc                                 <  1:519158/68‑1 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGcc                                 <  1:3169812/68‑1 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGcc                                 <  1:2674067/68‑1 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCg                                >  1:2831923/1‑69 (MQ=255)
             ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCg                                >  1:1019677/1‑69 (MQ=255)
                              gcgcCATCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:2225/1‑57 (MQ=37)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAgtat                            >  1:9093/1‑54 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAgtat                            >  1:3131803/1‑54 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:1397561/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:352519/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:1795127/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:3360753/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:2558791/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:1520159/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:3136046/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:213105/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:3004999/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:137875/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:2850536/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:1312417/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:2413018/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                           >  1:246978/1‑57 (MQ=38)
                              gcgcCAGCCGTAGATTGACTGGTCGTCATAGCCAACGATCATCACTTTGCCAGTATc                           >  1:2603160/1‑57 (MQ=255)
                                     ccGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATc        <  1:2337523/69‑1 (MQ=255)
                                       gTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGcc                                 >  1:1367643/1‑42 (MQ=255)
                                       gTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTcgcc                        >  1:1862475/1‑51 (MQ=255)
                                       gTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCa      >  1:3225176/1‑69 (MQ=255)
                                       gTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCa      >  1:2613707/1‑69 (MQ=255)
                                            ttGACTGGTCGTCATCGGCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:2000434/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:410985/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:2989401/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:521089/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:2985619/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:2785107/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:2690590/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:92951/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg  <  1:944340/68‑1 (MQ=255)
                                            ttGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCAcgcg  <  1:1019526/68‑1 (MQ=255)
                                             tGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCAcgcg  <  1:613475/67‑1 (MQ=255)
                                                               |                                                
GAAACGCTGAATATTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCG  >  minE/2399830‑2399941

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: