Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,399,911 G→A 55.1% T240T (ACC→ACT 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,9110GA55.1% 5.6 / 38.0 49T240T (ACC→ACTuvrDDNA‑dependent ATPase I and helicase II
Reads supporting (aligned to +/- strand):  ref base G (11/11);  new base A (18/9);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAG  >  minE/2399843‑2399979
                                                                    |                                                                    
ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCg                                                                      >  1:1019677/1‑69 (MQ=255)
ttCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCg                                                                      >  1:2831923/1‑69 (MQ=255)
                 gcgcCATCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:2225/1‑57 (MQ=37)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAgtat                                                                  >  1:3131803/1‑54 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAgtat                                                                  >  1:9093/1‑54 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:352519/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:3360753/1‑57 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:213105/1‑57 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:3136046/1‑57 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:1795127/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:1397561/1‑57 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:137875/1‑57 (MQ=38)
                 gcgcCAGCCGTAGATTGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATc                                                                 >  1:1312417/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:2613707/1‑69 (MQ=255)
                          gTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCa                                            >  1:3225176/1‑69 (MQ=255)
                               ttGACTGGTCGTCATCGGCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                        <  1:2000434/68‑1 (MQ=255)
                               ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                        <  1:2989401/68‑1 (MQ=255)
                               ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                        <  1:944340/68‑1 (MQ=255)
                               ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                        <  1:410985/68‑1 (MQ=255)
                               ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                        <  1:521089/68‑1 (MQ=255)
                               ttGACTGGTCGTCATCGCCAACGATCATCACTTTGCCAGTATCGCCTGCCAGCAGTCGGATCCAcgcg                                        <  1:92951/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)
                               ttGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCAcgcg                                        <  1:1019526/68‑1 (MQ=255)
                                tGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCAcgcg                                        <  1:613475/67‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTTTCCTgg                <  1:434955/68‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                <  1:3187805/68‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                <  1:2312378/68‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                <  1:1986383/68‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                <  1:76924/68‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                <  1:1623109/68‑1 (MQ=255)
                                                       tcatcaCTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTgg                <  1:1585506/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:1317840/1‑69 (MQ=255)
                                                               ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc       >  1:94018/1‑69 (MQ=255)
                                                               ttGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTcc       >  1:2947996/1‑69 (MQ=255)
                                                                    ggTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAg  <  1:3009008/69‑1 (MQ=255)
                                                                    |                                                                    
TTCTCCACCTGCGCCCCGCGCCAGCCGTAGATTGACTGGTCGTCATCACCGACGATCATCACTTTGCCGGTGTCGCCCGCCAGCAGGCGGATCCACGCGTACTGAATGTTGTTGGTATCCTGGAATTCGTCCACCAG  >  minE/2399843‑2399979

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: