Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 775,607 C→A 8.9% D410E (GAC→GAA icd → isocitrate dehydrogenase, specific for NADP+

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE775,6070CA 8.9% 221.6 / 5.1 67D410E (GAC→GAAicdisocitrate dehydrogenase, specific for NADP+
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base A (3/3);  total (34/33)
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.74e-01

AACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTT  >  minE/775538‑775677
                                                                     |                                                                      
aaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                                                       >  1:1380692/1‑71 (MQ=255)
aaCCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                                                       >  1:3461826/1‑71 (MQ=255)
   tgTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcgcg                                                                     <  1:3375972/69‑1 (MQ=255)
     tAACCTATGACTTCGAGAGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                                                   >  1:2277514/1‑70 (MQ=255)
       aCCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatca                                                                 <  1:942420/70‑1 (MQ=255)
        ccTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                                                                    >  1:1216677/1‑66 (MQ=255)
            tGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGa                                                                    <  1:179072/62‑1 (MQ=255)
            tGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa                                                            <  1:3517652/70‑1 (MQ=255)
            tGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGa                                                             <  1:2779013/69‑1 (MQ=255)
              aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatca                                                                 <  1:2323281/63‑1 (MQ=255)
              aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatc                                                                  >  1:2573086/1‑62 (MQ=255)
              aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaaaa                                                          >  1:1024636/1‑70 (MQ=255)
              aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAc                                                         <  1:2762455/71‑1 (MQ=255)
              aCTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAAc                                                         <  1:1865233/71‑1 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGAcg                                                                       >  1:648959/1‑53 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                                                   >  1:3083934/1‑57 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                      <  1:3236137/70‑1 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                      >  1:2410303/1‑70 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                      >  1:2405560/1‑70 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                      >  1:3347349/1‑70 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                      >  1:3995278/1‑70 (MQ=255)
                  cGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATg                                                      >  1:3953634/1‑70 (MQ=255)
                  cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATg                                                      >  1:964617/1‑70 (MQ=255)
                  cGAACGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATg                                                      >  1:2338659/1‑70 (MQ=255)
                   gaacGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTa                                                    <  1:3776972/68‑1 (MQ=255)
                   gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt                                                     <  1:1592153/70‑1 (MQ=255)
                   gAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGt                                                     <  1:2778753/70‑1 (MQ=255)
                    aacGTCTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGTa                                                    <  1:2672830/68‑1 (MQ=255)
                      cGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                                               >  1:1855742/1‑57 (MQ=255)
                         cTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGa                                                             >  1:1873766/1‑56 (MQ=255)
                         cTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCg                                              >  1:1897193/1‑71 (MQ=255)
                         cTGATGGAAGGCGCTAAGCTGCTGAAATGTTCAGAGTTTGGTGAAGCGatcat                                                                >  1:2997723/1‑53 (MQ=255)
                          tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                                                   <  1:633284/49‑1 (MQ=255)
                          tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGat                                                                   <  1:3160492/49‑1 (MQ=255)
                          tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg                                                 >  1:1065761/1‑67 (MQ=255)
                          tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATg                                                 >  1:1770847/1‑67 (MQ=255)
                          tGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCg                                              >  1:810503/1‑70 (MQ=255)
                            atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                                                   >  1:309759/1‑63 (MQ=255)
                            atggatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAg                                           >  1:341330/1‑71 (MQ=255)
                               gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACAt                                                       <  1:3061178/56‑1 (MQ=255)
                               gatggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGttt                                        >  1:136183/1‑71 (MQ=255)
                                atggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                                                   >  1:3669567/1‑59 (MQ=255)
                                 tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGt                                      <  1:1641257/71‑1 (MQ=255)
                                 tggCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGt                                      >  1:3609091/1‑71 (MQ=255)
                                   gcgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTa                                    <  1:1266873/71‑1 (MQ=255)
                                    cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGatcatc                                                               >  1:2916913/1‑43 (MQ=255)
                                    cgcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGtt                                     <  1:1202046/69‑1 (MQ=255)
                                     gcTAAACTGCTGAAATGTTCAGAGTTTGGTGCCGCGATCATCGaaa                                                           >  1:4145209/1‑46 (MQ=255)
                                     gcTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTaaa                                  <  1:648918/71‑1 (MQ=255)
                                      cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGaa                                                            <  1:1873003/44‑1 (MQ=255)
                                      cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAAt                                                  >  1:3191364/1‑54 (MQ=255)
                                      cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt                                          <  1:331815/62‑1 (MQ=255)
                                      cTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGt                                          <  1:1975351/62‑1 (MQ=255)
                                        aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                                                   >  1:1282353/1‑51 (MQ=255)
                                        aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                >  1:1124964/1‑70 (MQ=255)
                                        aaaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                >  1:3135333/1‑70 (MQ=255)
                                         agctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGaa                                                            <  1:2597352/39‑1 (MQ=255)
                                         aaCTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                >  1:2634726/1‑69 (MQ=255)
                                           ctgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATttatt                            <  1:2562336/71‑1 (MQ=255)
                                            tgctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                <  1:419521/66‑1 (MQ=255)
                                              ctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                                                   >  1:3545576/1‑45 (MQ=255)
                                              ctgAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTaa                                                   >  1:4046172/1‑45 (MQ=255)
                                                 aaaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                <  1:1884725/61‑1 (MQ=255)
                                                 aaaTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                <  1:4035568/61‑1 (MQ=255)
                                                   aTGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAg                    >  1:4001833/1‑71 (MQ=255)
                                                   aTGTTCAGAGTTTGGTGAAGCGATCATCGAAAACATGt                                                     <  1:3778614/38‑1 (MQ=255)
                                                      ttCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGc                   <  1:3102705/69‑1 (MQ=255)
                                                      ttCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGt                 <  1:2178723/71‑1 (MQ=255)
                                                          gagTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAtt                                <  1:837131/52‑1 (MQ=255)
                                                                    aCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTa                                    <  1:1063899/38‑1 (MQ=255)
                                                                    aCGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAAt                                 <  1:3997334/41‑1 (MQ=255)
                                                                     cgcgATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCgttgtt  <  1:2387288/71‑1 (MQ=255)
                                                                     |                                                                      
AACCGTAACCTATGACTTCGAGCGTCTGATGGATGGCGCTAAACTGCTGAAATGTTCAGAGTTTGGTGACGCGATCATCGAAAACATGTAATGCCGTAGTTTGTTAAATTTATTAACGGGAGCGTAACGCTCCCGTTGTT  >  minE/775538‑775677

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: