Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,294,615 T→C 23.7% R210R (CGA→CGG nac ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,294,6150TC23.7% 98.0 / 25.8 82R210R (CGA→CGGnacDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base T (16/45);  new base C (7/12);  total (23/59)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

CCGCTGGCAATCGCTGCGGTAAGCGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGA  >  minE/1294550‑1294679
                                                                 |                                                                
ccGCTGGCAATCGCTGCGGTAAGCGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGc                                                               >  1:1369065/1‑69 (MQ=255)
ccGCTGGCAATCGCTGCGGTAAGCGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGc                                                               >  1:187022/1‑69 (MQ=255)
           cGCTGCGGTAAGCGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGaaaa                                                    >  1:2971001/1‑69 (MQ=255)
            gCTGCGGTAAGCGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGaaaa                                                    >  1:1400116/1‑68 (MQ=255)
                     aGCGTGGCAATAGACTAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:1301551/69‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:205061/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:3140762/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:3338243/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:2620791/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:1061650/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:2357422/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:2013503/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:1670021/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:572622/68‑1 (MQ=255)
                       cGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCaa                                         <  1:1558851/68‑1 (MQ=255)
                                aGACTCAATTTCCCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               >  1:2912131/1‑69 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1388017/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:2443931/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:252155/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:2526403/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               >  1:1368765/1‑69 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:769209/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:3431156/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1384988/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:372162/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1907580/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:3058253/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               >  1:527784/1‑69 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1802428/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1404593/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:676449/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1656418/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:859412/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:613447/69‑1 (MQ=255)
                                aGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1429416/69‑1 (MQ=255)
                                 gACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAgg                                                         >  1:608015/1‑42 (MQ=255)
                                 gACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACg                                     >  1:1384493/1‑62 (MQ=255)
                                 gACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACg                                     >  1:2847459/1‑62 (MQ=255)
                                 gACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACg                                     >  1:3422953/1‑62 (MQ=255)
                                    tCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:3065352/65‑1 (MQ=255)
                                    tCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:1851662/65‑1 (MQ=255)
                                    tCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:119881/65‑1 (MQ=255)
                                    tCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:2238902/65‑1 (MQ=255)
                                    tCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:2135045/65‑1 (MQ=255)
                                    tCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:2803054/65‑1 (MQ=255)
                                     cAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                          <  1:3423058/69‑1 (MQ=255)
                                     cAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                          <  1:732322/69‑1 (MQ=255)
                                     cAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                          <  1:3430730/69‑1 (MQ=255)
                                     cAATTTCACCAATAACTTTTGCCGTGAGACGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAg                          <  1:3300280/69‑1 (MQ=255)
                                      aaTTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGc                         >  1:2606110/1‑69 (MQ=255)
                                      aaTTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGc                         >  1:2929786/1‑69 (MQ=255)
                                        tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                       <  1:2904942/69‑1 (MQ=255)
                                        tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                       <  1:2462302/69‑1 (MQ=255)
                                        tttCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCAc                       <  1:242047/69‑1 (MQ=255)
                                        tttCACCAATAACTTTTGCCGTGAGGCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCt                               <  1:2848716/61‑1 (MQ=255)
                                            aCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGt                    <  1:2770705/68‑1 (MQ=255)
                                            aCCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGt                    <  1:1560295/68‑1 (MQ=255)
                                               aaTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                <  1:2582921/69‑1 (MQ=255)
                                               aaTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                <  1:3008200/69‑1 (MQ=255)
                                                aTAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATc                <  1:3279490/68‑1 (MQ=255)
                                                   aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                   >  1:1292873/1‑62 (MQ=255)
                                                   aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                   >  1:2392970/1‑62 (MQ=255)
                                                   aCTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTa                   >  1:2288677/1‑62 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:508264/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:1148902/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:80109/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:1239962/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:622801/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:1837619/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:429097/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:3314151/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:315232/67‑1 (MQ=255)
                                                     ttttGCCGTGAGCAGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTCACAGCGCTGTAATCACTg            <  1:193745/67‑1 (MQ=255)
                                                     ttttGCCGAGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:1413046/67‑1 (MQ=255)
                                                     ctttGCCGTGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTg            <  1:2061566/66‑1 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:3129163/1‑67 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:509537/1‑67 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:1667411/1‑67 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:1540754/1‑67 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:1410566/1‑67 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:2658558/1‑67 (MQ=255)
                                                             tGAGCCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAACAGCGCTGTAATCACTGGGGAGaaaaa  >  1:2829861/1‑67 (MQ=255)
                                                                 |                                                                
CCGCTGGCAATCGCTGCGGTAAGCGTGGCAATAGACTCAATTTCACCAATAACTTTTGCCGTGAGTCGCCGTAGGGAAAAAGCCTCATCAACACGAAGTCTAATAGCACTGTAATCACTGGGGAGAAAGA  >  minE/1294550‑1294679

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: