Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,281,039 G→A 31.1% R94R (CGC→CGT fabR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,281,0390GA31.1% 56.6 / 29.3 64R94R (CGC→CGTfabRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (18/26);  new base A (2/18);  total (20/44)
Fisher's exact test for biased strand distribution p-value = 1.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

AATGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTA  >  minE/2280980‑2281100
                                                           |                                                             
aaTGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGt                                                      >  1:2189451/1‑69 (MQ=255)
       agaCCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGACGCATGAGTTGGCGt                                                >  1:988930/1‑68 (MQ=255)
       agaCCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGACGCATGAGTTGGCGt                                                >  1:974727/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGc                                             >  1:1134227/1‑67 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1326402/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1330122/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1379187/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:2917290/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:2537994/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:2299473/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:679544/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1993618/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1667177/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1924645/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCa                                            >  1:1723164/1‑68 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCAt                                           <  1:2592151/69‑1 (MQ=255)
           cGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTCGCGCATGAGTTGGCGTAGCAt                                           <  1:979247/69‑1 (MQ=255)
                     aTCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGctc                                 <  1:3060361/69‑1 (MQ=255)
                     aTCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGctc                                 <  1:2605892/69‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTTGCGCATGAGTTGGCGTAGCATTAAACCGCTCtcatca                           >  1:1840820/1‑67 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:2324652/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:2922494/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:3237898/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:3368783/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:1441516/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:1512040/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:2393960/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:1920954/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:1604799/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTa                                         <  1:1642443/53‑1 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCtcatca                           >  1:2249411/1‑67 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCtcatca                           >  1:1417513/1‑67 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCtcatca                           >  1:142064/1‑67 (MQ=255)
                             cccGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAAc                         >  1:455373/1‑69 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCGTCAACCATGGTCa                 <  1:2011705/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:689111/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:1129323/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:894985/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:3297680/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:658578/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:389695/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:72984/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:1193542/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:581776/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:1365740/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:260443/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:1690806/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:1827193/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:3135493/67‑1 (MQ=255)
                                       gCGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGGCa                 <  1:2100287/67‑1 (MQ=255)
                                        cGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:271370/66‑1 (MQ=255)
                                        cGATACGCTGACGCGCCTGACGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCa                 <  1:2028599/66‑1 (MQ=255)
                                            aCGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGt                   <  1:3406484/60‑1 (MQ=255)
                                            aCGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGt                   <  1:850265/60‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:1297728/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:1313901/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:1586278/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:2048979/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:2382801/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:4057/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:2505973/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:358629/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:2843015/69‑1 (MQ=255)
                                                    gcgcCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTa  <  1:3150889/69‑1 (MQ=255)
                                                           |                                                             
AATGTGGAGACCGAGGTGCGGATCACACTCCCGCCTTTGGCGATACGCTGACGCGCCTGGCGCATGAGTTGGCGTAGCATTAAACCGCTCTCATCAACCATGGTCAGACCCAGTTCGTCTA  >  minE/2280980‑2281100

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