Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 977,473 C→T 14.5% D558N (GAT→AAT)  dcp ← dipeptidyl carboxypeptidase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE977,4730CT14.5% 186.4 / 21.7 69D558N (GAT→AAT) dcpdipeptidyl carboxypeptidase II
Reads supporting (aligned to +/- strand):  ref base C (23/36);  new base T (6/4);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 3.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGCGCCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTT  >  minE/977403‑977541
                                                                      |                                                                    
aGCGCCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATc                                                                      >  1:371887/1‑71 (MQ=255)
aGCGCCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATc                                                                      >  1:3263749/1‑71 (MQ=255)
   gcCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATc                                                                      >  1:1441273/1‑68 (MQ=255)
   gcCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCgag                                                                   <  1:2825844/71‑1 (MQ=255)
   gcCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCgag                                                                   <  1:2726547/71‑1 (MQ=255)
   gcCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCgag                                                                   <  1:1841752/71‑1 (MQ=255)
      cGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCgag                                                                   <  1:1047425/68‑1 (MQ=255)
      cGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCgag                                                                   <  1:2777390/68‑1 (MQ=255)
      cGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCgag                                                                   <  1:2745954/68‑1 (MQ=255)
      cGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAg                                                                <  1:762404/71‑1 (MQ=255)
       gCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGt                                                               >  1:993939/1‑71 (MQ=255)
           ttCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTgcgg                                                           >  1:2384066/1‑71 (MQ=255)
              aaaaTCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTgcggcg                                                         <  1:874905/70‑1 (MQ=255)
               aaaTCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTGCGGCGCt                                                       <  1:1479959/71‑1 (MQ=255)
               aaaTCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCt                                                       <  1:688794/71‑1 (MQ=255)
                aaTCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCt                                                       >  1:986605/1‑70 (MQ=255)
                aaTCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCt                                                       >  1:1012080/1‑70 (MQ=255)
                    atcGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAg                                                    >  1:3484672/1‑69 (MQ=255)
                    atcGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGc                                                   <  1:1864764/70‑1 (MQ=255)
                    atcGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGc                                                   <  1:2772507/70‑1 (MQ=255)
                     tcGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTgc                                                             <  1:1740731/59‑1 (MQ=255)
                          aTCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGc                                            >  1:3875417/1‑71 (MQ=255)
                          aTCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGc                                            <  1:3268568/71‑1 (MQ=255)
                          aTCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGc                                            <  1:980303/71‑1 (MQ=255)
                                 aTTGCTTCGTTTTCTTCCAGGCAATGCCAGCTCATATCGAGAAGTGCGGCGCTAAGCAGTTCGctcatctc                                     <  1:912091/71‑1 (MQ=255)
                                 aTTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTGCGGCGCTAAGCAGTTCGct                                           <  1:437310/65‑1 (MQ=255)
                                 aTTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGctcatctc                                     <  1:35668/71‑1 (MQ=255)
                                  ttGCTTCGTTTTCTTCCAGGCAATGCCAGCGCGTATTGAGAAGTGCGGCGCt                                                       <  1:1908428/52‑1 (MQ=255)
                                  ttGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTaa                                                     <  1:2145324/54‑1 (MQ=255)
                                  ttGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGct                                           <  1:770627/64‑1 (MQ=255)
                                  ttGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGAGCTaa                                                     <  1:243895/54‑1 (MQ=255)
                                   tGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGctcatctca                                    >  1:1898933/1‑70 (MQ=255)
                                    gCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATc                                                                      >  1:2863402/1‑35 (MQ=255)
                                    gCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATc                                                                      >  1:1884854/1‑35 (MQ=255)
                                    gCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTgcgg                                                           >  1:1442985/1‑46 (MQ=255)
                                    gCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTgcgg                                                           >  1:440493/1‑46 (MQ=255)
                                    gCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTgcgg                                                           >  1:701760/1‑46 (MQ=255)
                                    gCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATa                                  >  1:150078/1‑71 (MQ=255)
                                     cTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATa                                  <  1:2479146/70‑1 (MQ=255)
                                      ttCGTTTTCTTCCAGGCAATGCCAGCGCATATTGAGAAGTgccg                                                           <  1:2376546/44‑3 (MQ=39)
                                        cGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCa                                                  >  1:964789/1‑51 (MQ=255)
                                        cGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCa                                                  >  1:211271/1‑51 (MQ=255)
                                        cGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCa                                                  >  1:218216/1‑51 (MQ=255)
                                         gTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCt                              <  1:538068/70‑1 (MQ=255)
                                          ttttCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCtt                             <  1:2131178/70‑1 (MQ=255)
                                            ttcttcCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCa                                                  >  1:1334202/1‑47 (MQ=255)
                                            ttcttcCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATa                                  >  1:1673144/1‑63 (MQ=255)
                                            ttcttcCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgt                          >  1:1500598/1‑71 (MQ=255)
                                            ttcttcCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgt                          >  1:1750486/1‑71 (MQ=255)
                                                 ccAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgttg                        <  1:1413324/68‑1 (MQ=255)
                                                 ccAGGCAATGCCACCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgttg                        <  1:831529/68‑1 (MQ=255)
                                                  cAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGaa                      >  1:3939846/1‑69 (MQ=255)
                                                  cAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGaa                      >  1:1547304/1‑69 (MQ=255)
                                                   aGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGctcatct                                      <  1:372635/52‑1 (MQ=255)
                                                   aGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGctcatct                                      <  1:2901565/52‑1 (MQ=255)
                                                   aGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAg                   <  1:2242138/71‑1 (MQ=255)
                                                     gCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGc                 <  1:3706760/71‑1 (MQ=255)
                                                       aaTGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgt                          >  1:301512/1‑60 (MQ=255)
                                                        aTGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTgg              <  1:3412419/71‑1 (MQ=255)
                                                           ccAGCGCATATCGTGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt           <  1:936037/71‑1 (MQ=255)
                                                           ccAGCGCATATCGTGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt           <  1:3515422/71‑1 (MQ=255)
                                                            cAGCGCATATTGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCAt           >  1:3427294/1‑70 (MQ=255)
                                                              gcgcATATTGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTAc        >  1:3051178/1‑71 (MQ=255)
                                                                gcATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgtt                         >  1:2754856/1‑52 (MQ=255)
                                                                gcATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTttgtt                         >  1:1689727/1‑52 (MQ=255)
                                                                gcATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGc      <  1:1335307/71‑1 (MQ=255)
                                                                  atatCGAGAAGTGCGGCGCTAAGCAGTTCGctcatc                                       >  1:3759039/1‑36 (MQ=255)
                                                                  atatCGAGAAGTGCGGCGCTAAGCAGTTCGctcat                                        >  1:3142823/1‑35 (MQ=255)
                                                                  atatCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAt    <  1:2609445/71‑1 (MQ=255)
                                                                   tatCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGc                 <  1:1421104/57‑1 (MQ=255)
                                                                    atCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGc                 <  1:1923113/56‑1 (MQ=255)
                                                                    atCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGc                 <  1:3088053/56‑1 (MQ=255)
                                                                     tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAttt  >  1:1805046/1‑70 (MQ=255)
                                                                     tCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCAttt  >  1:886995/1‑70 (MQ=255)
                                                                      |                                                                    
AGCGCCCGCAATTCAAAATCATCGACATCCTGCATTGCTTCGTTTTCTTCCAGGCAATGCCAGCGCATATCGAGAAGTGCGGCGCTAAGCAGTTCGCTCATCTCATACCCTTTGTTGAACAGGCTGGCATTACGCATTT  >  minE/977403‑977541

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: