Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,652,018 T→C 56.1% I84I (ATT→ATC glgC → glucose‑1‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,652,0180TC56.1% 16.7 / 49.2 66I84I (ATT→ATCglgCglucose‑1‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (16/13);  new base C (31/6);  total (47/19)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

AACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAG  >  minE/2651953‑2652086
                                                                 |                                                                    
aaCTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAg                                                                   >  1:696793/1‑69 (MQ=255)
aaCTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAg                                                                   >  1:1826941/1‑69 (MQ=255)
aaCTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAg                                                                   >  1:2309937/1‑69 (MQ=255)
            cgtcgtATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGtaa                                                       <  1:2428650/69‑3 (MQ=255)
            cgtcgtATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCa                                                       <  1:1756374/69‑1 (MQ=255)
                           aTCACCGAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaagaa                                        >  1:439251/1‑69 (MQ=255)
                           aTCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaagaa                                        >  1:1116377/1‑69 (MQ=255)
                           aTCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaagaa                                        >  1:3235705/1‑69 (MQ=255)
                           aTCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaagaa                                        >  1:3185006/1‑69 (MQ=255)
                           aTCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaagaa                                        >  1:2469215/1‑69 (MQ=255)
                           aTCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaagaa                                        >  1:1693368/1‑69 (MQ=255)
                           aTCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATgaaga                                         >  1:970537/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2918498/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2152700/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2181579/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2170270/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:110235/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2927282/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2932503/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:2936141/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:3054695/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:3135469/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:3220804/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:383208/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:5143/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:1649043/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:1590369/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:1782094/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:1956797/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:1137480/1‑68 (MQ=255)
                             cACCCAGTACCAGTCCCACACTCTGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGaaa                                       >  1:1423511/1‑68 (MQ=255)
                                   gTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:2161598/69‑1 (MQ=255)
                                   gTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:305244/69‑1 (MQ=255)
                                   gTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:2947770/69‑1 (MQ=255)
                                   gTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:347598/69‑1 (MQ=255)
                                   gTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:2079395/69‑1 (MQ=255)
                                   gTACCAGTCCCACACTCTGGTGCAGCACATTAAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:1247316/69‑1 (MQ=255)
                                    taacAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAg                               <  1:2412654/66‑1 (MQ=255)
                                           cccACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCg                        <  1:478741/69‑1 (MQ=255)
                                           cccACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCg                        >  1:811799/1‑69 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:1338156/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:527759/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:3375761/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:720599/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:1531278/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:776671/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:105928/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:1890125/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:1776291/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:1843615/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:227564/1‑35 (MQ=255)
                                              acacTCTGGTGCAGCACATCCAGCGCGGCTGGTCa                                                       >  1:2693298/1‑35 (MQ=255)
                                                 ctctGGTGCAGTACATCCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:1387822/69‑1 (MQ=255)
                                                 ctctGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  >  1:321469/1‑69 (MQ=255)
                                                 ctctGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:407206/69‑1 (MQ=255)
                                                 ctctGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:2950935/69‑1 (MQ=255)
                                                 ctctGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:1711543/69‑1 (MQ=255)
                                                 ctctGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:814019/69‑1 (MQ=255)
                                                 ctctGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:2884190/69‑1 (MQ=255)
                                                  tctGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgct                 <  1:3407668/69‑1 (MQ=255)
                                                  tctGGTGCAGCACATCCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATctgc                  <  1:606070/68‑1 (MQ=255)
                                                         cagcaCATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGa                                <  1:2460901/47‑1 (MQ=255)
                                                         cagcaCATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCcag            >  1:345245/1‑67 (MQ=255)
                                                         cagcaCATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCcag            >  1:3149976/1‑67 (MQ=255)
                                                             acaTTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCcagcac         >  1:2693155/1‑66 (MQ=255)
                                                              caTTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCcagcacag       >  1:2031106/1‑67 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:2113214/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:811135/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:2281354/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:2355651/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:670894/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:1288320/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:381485/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:2847559/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:3244031/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:3150704/69‑1 (MQ=255)
                                                                 tCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCagag  <  1:1691259/69‑1 (MQ=255)
                                                                 |                                                                    
AACTCCGGGATCCGTCGTATGGGCGTGATCACCCAGTACCAGTCCCACACTCTGGTGCAGCACATTCAGCGCGGCTGGTCATTCTTCAATGAAGAAATGAACGAGTTTGTCGATCTGCTGCCAGCACAGCAGAG  >  minE/2651953‑2652086

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