Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,528,498 C→T 24.6% D83N (GAT→AAT)  dfp ← fused 4'‑phosphopantothenoylcysteine decarboxylase and phosphopantothenoylcysteine synthetase, FMN‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,528,4980CT24.6% 116.3 / 37.6 61D83N (GAT→AAT) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase and phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (26/20);  new base T (5/10);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

TTCGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCGG  >  minE/2528434‑2528563
                                                                |                                                                 
ttCGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCa                                                             >  1:1958821/1‑71 (MQ=255)
ttCGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCa                                                             >  1:1438951/1‑71 (MQ=255)
  cGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCAtt                                                           >  1:4240056/1‑71 (MQ=255)
   gCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAg                                                                 >  1:166959/1‑64 (MQ=255)
   gCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCAtt                                                           <  1:3806394/70‑1 (MQ=255)
          cGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGc                                                   >  1:1671108/1‑71 (MQ=255)
            gcagcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTc                                                 <  1:3252569/71‑1 (MQ=255)
            gcagcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTc                                                 <  1:3762045/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCaa                                               <  1:405589/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCaa                                               <  1:3350347/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCaa                                               <  1:1121420/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCaa                                               <  1:4249336/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCaa                                               <  1:1869007/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCaa                                               <  1:1716111/71‑1 (MQ=255)
              agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCaa                                               <  1:1755502/71‑1 (MQ=255)
                            tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTa                                                         <  1:1111345/47‑1 (MQ=255)
                            tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTa                                                         <  1:1490542/47‑1 (MQ=255)
                            tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCg                                 <  1:1780547/71‑1 (MQ=255)
                            tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCg                                 >  1:3436097/1‑71 (MQ=255)
                            tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCg                                 <  1:2817920/71‑1 (MQ=255)
                              aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGc                                  <  1:3919517/68‑1 (MQ=255)
                              aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGc                                  <  1:2463705/68‑1 (MQ=255)
                              aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                               >  1:667881/1‑71 (MQ=255)
                              aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                               >  1:2555667/1‑71 (MQ=255)
                                     ccGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGc                                                   <  1:2735566/44‑1 (MQ=255)
                                     ccGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGc                                                   <  1:1019683/44‑1 (MQ=255)
                                     ccGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGc                                                   <  1:2082936/44‑1 (MQ=255)
                                     ccGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGc                                                   <  1:1821563/44‑1 (MQ=255)
                                     ccGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCgg                          <  1:4234392/69‑1 (MQ=255)
                                     ccGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCt                        >  1:3857636/1‑71 (MQ=255)
                                       gTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGc                      >  1:2375494/1‑71 (MQ=255)
                                         ggCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCg                    >  1:446440/1‑71 (MQ=255)
                                         ggCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCg                    >  1:272661/1‑71 (MQ=255)
                                         ggCAGGGGAGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATgg                                          >  1:2876178/1‑49 (MQ=255)
                                             ggggCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGcc                     <  1:1184835/66‑1 (MQ=255)
                                             ggggCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGcc                     <  1:2834456/66‑1 (MQ=255)
                                              gggCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTcc               <  1:2776113/71‑1 (MQ=255)
                                                gCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcag             >  1:1606252/1‑71 (MQ=255)
                                                  gagaATCACTAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcag             >  1:3095709/1‑69 (MQ=255)
                                                  gagaATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcagca           <  1:3207455/71‑1 (MQ=255)
                                                  gagaATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcagca           <  1:1854257/71‑1 (MQ=255)
                                                   agaATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcagcag          <  1:907257/71‑1 (MQ=255)
                                                   agaATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcagcag          <  1:1870713/71‑1 (MQ=255)
                                                    gaATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGa         >  1:3630238/1‑71 (MQ=255)
                                                     aaTCACTAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGAc        >  1:2473724/1‑71 (MQ=255)
                                                       tCACTAAATCAGCCCATTTACCCAGCTCAATATgg                                          >  1:3552627/1‑35 (MQ=255)
                                                       tCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTg      >  1:2466393/1‑71 (MQ=255)
                                                       tCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTg      >  1:1381647/1‑71 (MQ=255)
                                                        cACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCAt                                     >  1:2570526/1‑39 (MQ=255)
                                                            aaaTTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCgg  >  1:806655/1‑70 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCgg                          >  1:1681866/1‑46 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGc                         >  1:3269266/1‑47 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGc                         >  1:1330258/1‑47 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGc                         >  1:2844517/1‑47 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGc                         >  1:378548/1‑47 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGc                         >  1:1860883/1‑47 (MQ=255)
                                                            aaaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCgg  >  1:3019657/1‑70 (MQ=255)
                                                             aaTCAGCCCATTTACCCAGCTCAATATGGCCCATAg                                   >  1:2930132/1‑36 (MQ=255)
                                                             aaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGcc                     <  1:1743278/50‑1 (MQ=255)
                                                             aaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGcc                     <  1:675650/50‑1 (MQ=255)
                                                             aaTCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCgg                   >  1:4058196/1‑52 (MQ=255)
                                                                |                                                                 
TTCGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCGG  >  minE/2528434‑2528563

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: