Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,399,914 G→A 52.1% D239D (GAC→GAT 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,9140GA52.1% ‑3.0 / 40.8 48D239D (GAC→GATuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base A (16/9);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGAT  >  minE/2399860‑2399982
                                                      |                                                                    
gcgcCATCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:2225/1‑57 (MQ=37)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:1795127/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:2558791/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:246978/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:2413018/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:213105/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:2850536/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:3004999/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:1520159/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:3136046/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:1397561/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:137875/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:3360753/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:1312417/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:352519/1‑57 (MQ=38)
gcgcCAGCCGTAGATTGACTGGTCGTCATAGCCAACGATCATCACTTTGCCAGTATc                                                                    >  1:2603160/1‑57 (MQ=255)
       ccGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATc                                                 <  1:2337523/69‑1 (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:92951/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:521089/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:410985/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:2690590/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:2785107/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:2985619/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:944340/68‑1 (MQ=255)
              ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                           <  1:2989401/68‑1 (MQ=255)
              ttGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCAcgcg                                           <  1:1019526/68‑1 (MQ=255)
               tGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCAcgcg                                           <  1:613475/67‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTTTCCTgg                   <  1:434955/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                   <  1:2312378/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                   <  1:76924/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                   <  1:3187805/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                   <  1:1585506/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                   <  1:1986383/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                   <  1:1623109/68‑1 (MQ=255)
                                      tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTCTTGGTATCCTgg                   <  1:3304181/68‑1 (MQ=255)
                                       catcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGaa                 >  1:1106116/1‑69 (MQ=255)
                                       catcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGaa                 >  1:1199722/1‑69 (MQ=255)
                                           aCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCg             >  1:407780/1‑69 (MQ=255)
                                              ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc          >  1:86416/1‑69 (MQ=255)
                                              ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc          >  1:882514/1‑69 (MQ=255)
                                              ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc          >  1:94018/1‑69 (MQ=255)
                                              ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc          >  1:1317840/1‑69 (MQ=255)
                                              ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc          >  1:2947996/1‑69 (MQ=255)
                                                   ggTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAg     <  1:3009008/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:3380091/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:3278725/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:1445531/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:1111435/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:1032954/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:1466941/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:2081380/69‑1 (MQ=255)
                                                     tgtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGa   <  1:2081626/69‑1 (MQ=255)
                                                      gtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGat  >  1:2437741/1‑69 (MQ=255)
                                                      gtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGat  >  1:718807/1‑69 (MQ=255)
                                                      gtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGat  <  1:1994070/69‑1 (MQ=255)
                                                      gtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGat  <  1:321087/69‑1 (MQ=255)
                                                      gtCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTAGTCCACCAGGat  >  1:371367/1‑69 (MQ=255)
                                                      |                                                                    
GCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAGGAT  >  minE/2399860‑2399982

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: