Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 464,749 C→T 54.2% S202S (AGC→AGT sdhB → succinate dehydrogenase, FeS subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE464,7490CT54.2% 6.9 / 59.7 72S202S (AGC→AGTsdhBsuccinate dehydrogenase, FeS subunit
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base T (23/16);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTG  >  minE/464682‑464786
                                                                   |                                     
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGc                                       >  1:2697698/1‑68 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:3060049/1‑69 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:1869110/1‑69 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:723965/1‑69 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:2868837/1‑69 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:884478/1‑69 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:1380973/1‑69 (MQ=255)
gTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCg                                      >  1:2712403/1‑69 (MQ=255)
  ttCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTa                                    <  1:2366740/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGTTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:3080154/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTcc                                >  1:675847/1‑68 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:2042867/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:3093714/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:679753/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:2330642/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:2435890/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:664723/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:1021211/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:2789528/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:2898700/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:2905543/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:3347036/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:3065800/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:1846458/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:99595/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1084930/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1100571/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:120843/1‑69 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1337722/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1338593/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1342559/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:81857/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1884480/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:737234/69‑1 (MQ=255)
       gatcgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               >  1:1765616/1‑69 (MQ=255)
          tgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCt                             >  1:168392/1‑68 (MQ=255)
          tgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCt                             >  1:3053589/1‑68 (MQ=255)
          tgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCt                             >  1:1751092/1‑68 (MQ=255)
          tgatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCt                             >  1:635211/1‑68 (MQ=255)
           gatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGtc                                    >  1:1306879/1‑59 (MQ=255)
           gatAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTa                                    >  1:2668222/1‑60 (MQ=255)
                       aCCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATg               <  1:810780/69‑1 (MQ=255)
                       aCCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATg               <  1:2470856/69‑1 (MQ=255)
                       aCCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATg               <  1:342022/69‑1 (MQ=255)
                       aCCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATg               <  1:1551299/69‑1 (MQ=255)
                        ccGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGa              >  1:2811198/1‑69 (MQ=255)
                        ccGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGa              >  1:449132/1‑69 (MQ=255)
                        ccGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGa              >  1:562543/1‑69 (MQ=255)
                         cGAGACTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCg                               <  1:1714586/51‑1 (MQ=255)
                            gaCTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTg          <  1:2824989/69‑1 (MQ=255)
                            gaCTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTg          <  1:2926174/69‑1 (MQ=255)
                            gaCTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTg          <  1:2313843/69‑1 (MQ=255)
                            gaCTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTg          <  1:3079215/69‑1 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGCGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTg          >  1:2781731/1‑68 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:2901185/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:1078231/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:823325/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:1633795/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:1748919/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:1782189/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:1855452/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:2785528/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:468186/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:2886572/1‑69 (MQ=255)
                             aCTGACAGCCGCCTCGACGGTTTGAGCGATGCATTCAGTGTATTCCGCTGTCACAGCATCATGAACTGc         >  1:3038218/1‑69 (MQ=255)
                                    gccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:2645746/69‑1 (MQ=255)
                                    gccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:2540669/69‑1 (MQ=255)
                                    gccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:2275934/69‑1 (MQ=255)
                                    gccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:683344/69‑1 (MQ=255)
                                    gccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:445293/69‑1 (MQ=255)
                                    gccgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:1362683/69‑1 (MQ=255)
                                      cgccTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:495882/67‑1 (MQ=255)
                                                                cAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAgtg  <  1:2704365/41‑1 (MQ=255)
                                                                   |                                     
GTTTCCTGATTGATAGCCGTGATACCGAGACTGACAGCCGCCTCGACGGTTTGAGTGATGCATTCAGCGTATTCCGCTGTCACAGCATCATGAACTGCGTCAGTG  >  minE/464682‑464786

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: