Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,491,595 C→T 16.3% P96P (CCG→CCA nuoJ ← NADH:ubiquinone oxidoreductase, membrane subunit J

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,491,5950CT16.3% 128.4 / 16.5 55P96P (CCG→CCAnuoJNADH:ubiquinone oxidoreductase, membrane subunit J
Reads supporting (aligned to +/- strand):  ref base C (25/21);  new base T (3/6);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTGAGC  >  minE/1491525‑1491664
                                                                      |                                                                     
cGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGcc                                                                       >  1:122383/1‑71 (MQ=255)
  aTACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCgg                                                                     <  1:980933/71‑1 (MQ=255)
    aCCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGAc                                                                   >  1:69971/1‑71 (MQ=255)
       cTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCaa                                                                <  1:3451511/71‑1 (MQ=255)
        tGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGa                                                                    >  1:3569426/1‑66 (MQ=255)
        tGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGAcc                                                                  >  1:385904/1‑68 (MQ=255)
        tGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAAt                                                               >  1:3125172/1‑71 (MQ=255)
          aTCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAAt                                                               <  1:2510808/69‑1 (MQ=255)
           tCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCca                                                            <  1:3120986/71‑1 (MQ=255)
            cGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCTGGACCAATCca                                                            >  1:3623154/1‑70 (MQ=255)
            cGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCca                                                            >  1:2288610/1‑70 (MQ=255)
            cGTTAACACCGAGGATGGCGTAAACAATCACCACAAGCATGATGGCCGACAAAATTGCCGGACCAATCca                                                            >  1:1591219/1‑70 (MQ=255)
             gTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCcac                                                           <  1:12248/70‑1 (MQ=255)
             gTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCcac                                                           <  1:716481/70‑1 (MQ=255)
                aaCACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCt                                                       >  1:1193852/1‑71 (MQ=255)
                   acCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCg                                                    >  1:2440895/1‑71 (MQ=255)
                   acCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCg                                                    <  1:306768/71‑1 (MQ=255)
                    ccGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCg                                                                      <  1:591312/52‑1 (MQ=255)
                    ccGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCg                                                                      <  1:2245418/52‑1 (MQ=255)
                    ccGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCgg                                                   >  1:1809796/1‑71 (MQ=255)
                     cGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGt                                                  <  1:1924035/71‑1 (MQ=255)
                           tGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGc                                            >  1:417647/1‑71 (MQ=255)
                            ggCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGt                                                  <  1:303482/64‑1 (MQ=255)
                             gCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGcc                                           <  1:3237340/70‑1 (MQ=255)
                              cGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCAc                                         >  1:2816328/1‑71 (MQ=255)
                                tAAACAATCACCACCAGCATGATGGCCGACAAAATTGcc                                                                       >  1:1427425/1‑39 (MQ=255)
                                tAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCa                                              <  1:144521/64‑1 (MQ=255)
                                  aaCAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGc                                            >  1:1585157/1‑64 (MQ=255)
                                  aaCAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTgg                                      <  1:2556557/70‑1 (MQ=255)
                                  aaCAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGc                                     >  1:1343029/1‑71 (MQ=255)
                                  aaCAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGc                                     <  1:1754186/71‑1 (MQ=255)
                                   aCAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGttt                                                >  1:518990/1‑59 (MQ=255)
                                    cAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGt                                   <  1:752769/71‑1 (MQ=255)
                                         accaccAGCATGATGGCCGACAAAATTGCCGGACCAATCcac                                                           >  1:179946/1‑42 (MQ=255)
                                             ccaGCATGATGGCCGACAAAATTGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTc                           >  1:862239/1‑70 (MQ=255)
                                             ccaGCATGATGGCCGACAAAATTGCCGGACCAATc                                                              >  1:1588951/1‑35 (MQ=255)
                                               aGCATGATGGCCGACAAAATTGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGAt                        <  1:3744530/71‑1 (MQ=255)
                                                gCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTc                           >  1:1938213/1‑67 (MQ=255)
                                                   tgatgGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCt                    <  1:1419664/71‑1 (MQ=255)
                                                     atgGCCGACAAAATTGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTg                   <  1:864474/70‑1 (MQ=255)
                                                     atgGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTg                   <  1:407337/70‑1 (MQ=255)
                                                      tgGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTg                                       <  1:1581825/49‑1 (MQ=255)
                                                        gCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAccc               >  1:2219231/1‑71 (MQ=255)
                                                         ccGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCg              >  1:3200140/1‑71 (MQ=255)
                                                           gACAAAATTGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGcc            <  1:1280334/71‑1 (MQ=255)
                                                           gACAAAATTGCCGGACCAATCCACACCTGCGGttt                                                <  1:1615825/35‑1 (MQ=255)
                                                            aCAAAATTGCCGGACCAATCCACACCTGCGGTTTCa                                              >  1:3868855/1‑36 (MQ=255)
                                                            aCAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGcc            >  1:1241996/1‑70 (MQ=255)
                                                             cAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGc                                            >  1:2336437/1‑37 (MQ=255)
                                                              aaaaTTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTc                           >  1:2639246/1‑53 (MQ=255)
                                                               aaaTTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGt                             <  1:2478043/50‑1 (MQ=255)
                                                                aaTTGCCGGACCAATCCACACCTGCGGTTTCAGCCAc                                         >  1:3193480/1‑37 (MQ=255)
                                                                aaTTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGt                             >  1:1387939/1‑49 (MQ=255)
                                                                 aTTGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGtt      <  1:749075/71‑1 (MQ=255)
                                                                 aTTGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGtt      <  1:753657/71‑1 (MQ=255)
                                                                  ttGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTg     <  1:255019/71‑1 (MQ=255)
                                                                   tGCTGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGtt      >  1:211306/1‑69 (MQ=255)
                                                                     ccGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTGAGc  <  1:352635/71‑1 (MQ=255)
                                                                      |                                                                     
CGATACCCTGATCGTTAACACCGAGGATGGCGTAAACAATCACCACCAGCATGATGGCCGACAAAATTGCCGGACCAATCCACACCTGCGGTTTCAGCCACTGGCGTTCCTGTTCGATTTCTGAACCGCCCAGGTTGAGC  >  minE/1491525‑1491664

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: