Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,528,498 C→T 21.8% 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,4980CT21.8% 117.1 / 27.0 55D83N (GAT→AAT) dfpfused 4'‑phosphopantothenoylcysteine decarboxylase and phosphopantothenoylcysteine synthetase, FMN‑binding
Reads supporting (aligned to +/- strand):  ref base C (26/17);  new base T (7/5);  total (33/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCG  >  minE/2528436‑2528562
                                                              |                                                                
cGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCAtt                                                          >  1:399179/1‑71 (MQ=255)
 gCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATc                                                                  <  1:534146/62‑1 (MQ=255)
 gCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCAttt                                                         <  1:1843881/71‑1 (MQ=255)
   cATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTAc                                                       >  1:1654053/1‑71 (MQ=255)
       ccGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTa                                                        >  1:297880/1‑66 (MQ=255)
       ccGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCa                                                    >  1:1469731/1‑70 (MQ=255)
       ccGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAg                                                   >  1:995663/1‑71 (MQ=255)
       ccGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAg                                                   >  1:314257/1‑71 (MQ=255)
            agcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCaa                                              <  1:1667221/71‑1 (MQ=255)
             gcaACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTc                                                <  1:1978690/68‑1 (MQ=255)
              caACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAata                                            >  1:331886/1‑71 (MQ=255)
              caACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAata                                            >  1:920627/1‑71 (MQ=255)
              caACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAata                                            >  1:1213465/1‑71 (MQ=255)
              caACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAat                                             >  1:197468/1‑70 (MQ=255)
                acacGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATg                                          >  1:187874/1‑71 (MQ=255)
                 cacGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCAttt                                                         <  1:1184573/55‑1 (MQ=255)
                  acGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCACATTTACCCAGCTCAata                                            >  1:2194052/1‑67 (MQ=255)
                   cGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCAATATGGcc                                       <  1:585091/71‑1 (MQ=255)
                         aTCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTa                                                        <  1:558972/48‑1 (MQ=255)
                          tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCaa                                              <  1:70315/57‑1 (MQ=255)
                          tCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCa                                                            <  1:61884/43‑1 (MQ=255)
                           cAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCa                                                    >  1:957722/1‑50 (MQ=255)
                           cAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCa                                     >  1:1903493/1‑65 (MQ=255)
                           cAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCg                                >  1:909106/1‑70 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCAata                                            >  1:509682/1‑57 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATTAGCCCATTTACCCAGCTCAata                                            >  1:2214075/1‑57 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATTGcc                                       >  1:1899626/1‑62 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                              >  1:616752/1‑71 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                              >  1:1104866/1‑71 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                              >  1:2199288/1‑71 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                              >  1:1056121/1‑71 (MQ=255)
                            aaaTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGc                              >  1:1597189/1‑71 (MQ=255)
                              aTCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCtt                            <  1:2085266/71‑1 (MQ=255)
                               tCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGc                                                               >  1:1411687/1‑35 (MQ=255)
                               tCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGc                                                               >  1:847287/1‑35 (MQ=255)
                                 tGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCa                                               <  1:1626797/49‑1 (MQ=255)
                                 tGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCa                                               <  1:1593899/49‑1 (MQ=255)
                                   ccGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTAccc                                                     >  1:689082/1‑41 (MQ=255)
                                   ccGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCt                       <  1:428948/71‑1 (MQ=255)
                                     gTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCaa                                              <  1:732650/46‑1 (MQ=255)
                                     gTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGc                     <  1:1776820/71‑1 (MQ=255)
                                     gTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGc                     <  1:1187164/71‑1 (MQ=255)
                                      tGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCt                       <  1:1875935/68‑1 (MQ=255)
                                      tGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCt                       <  1:1873274/68‑1 (MQ=255)
                                      tGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGcc                    >  1:114650/1‑71 (MQ=255)
                                         cAGGGGCGAGAATCACTAAATCAGCCCATTTACACAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCggg                 <  1:671515/71‑1 (MQ=255)
                                           ggggCGAGAATCACTAAATTAGCCCATTTACCCAGCTCaa                                              <  1:1652953/40‑1 (MQ=255)
                                            gggCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTc               >  1:1393538/1‑70 (MQ=255)
                                             ggCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCca             <  1:2125876/71‑1 (MQ=255)
                                              gCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAata                                            >  1:162344/1‑39 (MQ=255)
                                              gCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAata                                            >  1:507750/1‑39 (MQ=255)
                                              gCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcag            <  1:1601090/71‑1 (MQ=255)
                                                gagaATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCcagca          <  1:555005/71‑1 (MQ=255)
                                                       aCTAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTc   >  1:260737/1‑71 (MQ=255)
                                                         tAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCg  >  1:277525/1‑70 (MQ=255)
                                                         tAAATTAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCg  >  1:2103672/1‑70 (MQ=255)
                                                              |                                                                
CGCCATTCCGGCAGCAACACGGGCAATCAAATCTGCCGTGGCAGGGGCGAGAATCACTAAATCAGCCCATTTACCCAGCTCAATATGGCCCATAGCGGCTTCGGCTGCCGGGTCCAGCAGACTGTCG  >  minE/2528436‑2528562

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: