Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,475 C→T 10.9% H366H (CAC→CAT icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,4750CT10.9% 150.2 / 8.0 55H366H (CAC→CATicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (4/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

AATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAA  >  minE/775408‑775545
                                                                   |                                                                      
aaTATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATg                                                                     <  1:3244265/71‑1 (MQ=255)
 atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                    >  1:2080564/1‑71 (MQ=255)
 atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                    >  1:3271018/1‑71 (MQ=255)
 atatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATgg                                                                    >  1:3924770/1‑71 (MQ=255)
  tatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                   <  1:2697726/71‑1 (MQ=255)
  tatGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                   <  1:796021/71‑1 (MQ=255)
   atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATggg                                                                   <  1:32292/70‑1 (MQ=255)
   atGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGt                                                                  >  1:2085816/1‑71 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGtt                                                                 >  1:2893131/1‑71 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGtt                                                                 >  1:88874/1‑71 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGt                                                                  <  1:734120/70‑1 (MQ=255)
    tGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGt                                                                  <  1:1839639/70‑1 (MQ=255)
       cGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGa                                                              >  1:218588/1‑71 (MQ=255)
         gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGAcc                                                            >  1:3059513/1‑71 (MQ=255)
         gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGAcc                                                            <  1:2267106/71‑1 (MQ=255)
         gTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGAcc                                                            >  1:3802079/1‑71 (MQ=255)
          tCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTCCGCCACATGGGTTGGACCg                                                           <  1:2033961/71‑1 (MQ=255)
             ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACAt                                                                      <  1:1392972/57‑1 (MQ=255)
             ggACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAg                                                        <  1:2203892/71‑1 (MQ=255)
               aCAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                      <  1:1138752/71‑1 (MQ=255)
                cAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCgg                                                     >  1:714894/1‑71 (MQ=255)
                cAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCg                                                      >  1:2753179/1‑70 (MQ=255)
                 aaaGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGc                                                    <  1:1649686/71‑1 (MQ=255)
                   aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcaca                                                                       >  1:1084858/1‑50 (MQ=255)
                   aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCcaca                                                                       >  1:3587809/1‑50 (MQ=255)
                   aGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGt                                                                  >  1:1573541/1‑55 (MQ=255)
                     tAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGa                                                              >  1:3629737/1‑57 (MQ=255)
                      aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                               >  1:1124496/1‑71 (MQ=255)
                      aaaTCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACt                                               >  1:3263982/1‑71 (MQ=255)
                           cTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTg                                                  <  1:777707/63‑1 (MQ=255)
                               ctctATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTa                                      <  1:3803664/71‑1 (MQ=255)
                                  tattatTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaa                                     <  1:3346182/69‑1 (MQ=255)
                                    ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAg                                                        >  1:1041134/1‑48 (MQ=255)
                                    ttattCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAg                                                        >  1:1626294/1‑48 (MQ=255)
                                         ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                             <  1:1377093/70‑1 (MQ=255)
                                         ctctcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                             <  1:3127539/70‑1 (MQ=255)
                                            tcCGCTGAGATGATGCTGCGCCACATGTGTTGGACCGAAg                                                        >  1:607895/1‑40 (MQ=255)
                                            tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGa                            <  1:1239885/68‑1 (MQ=255)
                                            tcCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGa                            <  1:1499814/68‑1 (MQ=255)
                                              cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAg                          <  1:2283565/68‑1 (MQ=255)
                                              cGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGgcg                       >  1:3301351/1‑71 (MQ=255)
                                                cTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATgg                             >  1:1714068/1‑63 (MQ=255)
                                                 tGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCaa                    >  1:3162991/1‑71 (MQ=255)
                                                    gatgatGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCa                 <  1:2997594/71‑1 (MQ=255)
                                                    gacgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCa                 <  1:2011454/68‑1 (MQ=255)
                                                     ctaatgCTGCGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCaa                <  1:3770846/68‑1 (MQ=255)
                                                     atgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCaa                <  1:1334831/71‑1 (MQ=255)
                                                      tgatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCaa                <  1:164725/70‑1 (MQ=255)
                                                       gatgCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTaa                                     <  1:3803803/48‑1 (MQ=255)
                                                         tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTa                                      >  1:2101086/1‑45 (MQ=255)
                                                         tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgcc            <  1:2660882/71‑2 (MQ=255)
                                                         tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgc             >  1:3095905/1‑70 (MQ=255)
                                                         tgctgcGCCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATGGAAGGCGCAATCAAtgc             >  1:1458380/1‑70 (MQ=255)
                                                           ctgcGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGaa                           <  1:236511/54‑1 (MQ=255)
                                                              cgcCATATGGGTTGGACTGAAGCGGCTGACCTGATTGTTAAAGGTATgg                             >  1:339422/1‑49 (MQ=37)
                                                              cgcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAg                                   >  1:352914/1‑43 (MQ=255)
                                                               gcCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAAc               >  1:3032286/1‑62 (MQ=255)
                                                                  acaTGGGTTGGACCGAAGCGGCTGACTTAATTGTTaaa                                    <  1:3699681/38‑1 (MQ=255)
                                                                   caTGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTaa  >  1:890819/1‑71 (MQ=255)
                                                                   |                                                                      
AATATGCCGGTCAGGACAAAGTAAATCCTGGCTCTATTATTCTCTCCGCTGAGATGATGCTGCGCCACATGGGTTGGACCGAAGCGGCTGACTTAATTGTTAAAGGTATGGAAGGCGCAATCAACGCGAAAACCGTAA  >  minE/775408‑775545

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