Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,738,791 G→A 59.6% C340C (TGC→TGT lepA ← GTP binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,738,7910GA59.6% 25.2 / 54.6 47C340C (TGC→TGTlepAGTP binding membrane protein
Reads supporting (aligned to +/- strand):  ref base G (4/15);  new base A (5/23);  total (9/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTC  >  Exported/1738722‑1738854
                                                                     |                                                               
cagatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCg                                                                 <  1:1386624/70‑1 (MQ=255)
  gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAg                                                              <  1:987585/71‑1 (MQ=255)
  gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAg                                                              <  1:2879695/71‑1 (MQ=255)
  gatccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAg                                                              <  1:3100228/71‑1 (MQ=255)
   atccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGc                                                             >  1:751354/1‑71 (MQ=255)
   atccagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGc                                                             >  1:254531/1‑71 (MQ=255)
      cagatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGa                                                          <  1:3172500/71‑1 (MQ=255)
       agatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGaa                                                         <  1:2164931/71‑1 (MQ=255)
        gatcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGaaa                                                        <  1:698903/71‑1 (MQ=255)
         atcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAAc                                                       <  1:3167543/71‑1 (MQ=255)
          tcGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAAcc                                                      <  1:3031611/71‑1 (MQ=255)
            gTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCaa                                                    <  1:1033688/71‑1 (MQ=255)
                  aCGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGccc                                               <  1:3184090/70‑1 (MQ=255)
                    gTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAgc                                            <  1:1777260/71‑1 (MQ=255)
                     ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCa                                                     <  1:724521/61‑1 (MQ=255)
                     ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCa                                                     <  1:696135/61‑1 (MQ=255)
                     ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCa                                                     <  1:3262663/61‑1 (MQ=255)
                     ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCa                                                     <  1:1431088/61‑1 (MQ=255)
                     ttCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCa                                                     <  1:859341/61‑1 (MQ=255)
                        cAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAgcgcgc                                        <  1:1589788/71‑1 (MQ=255)
                         aGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAgcgcg                                         <  1:2080742/69‑1 (MQ=255)
                         aGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCt                                       <  1:2999730/71‑1 (MQ=255)
                           aCGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCTGCAGCGGAAACCAAAGCCCAGCGCGCTg                                      >  1:154783/1‑70 (MQ=255)
                             gTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGTAAACCAAAGcc                                                <  1:2902098/58‑1 (MQ=255)
                             gTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGcc                                                <  1:2781404/58‑1 (MQ=255)
                             gTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAg                                   >  1:3173840/1‑71 (MQ=255)
                                 cTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCttt                               <  1:3308757/71‑1 (MQ=255)
                                 cTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCtt                                <  1:3022965/70‑1 (MQ=255)
                                    gatgatCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTcc                             <  1:2827335/70‑1 (MQ=255)
                                    gatgatCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTc                              <  1:3031562/69‑1 (MQ=255)
                                    gatgatCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCg                            <  1:7762/71‑1 (MQ=255)
                                    gatgatCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCg                            >  1:2228454/1‑71 (MQ=255)
                                     atgatCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCgg                           <  1:2859492/71‑1 (MQ=255)
                                        atCTCCATGTGCAGCATGCCGAGGAAGCCACAGCGGAAAcc                                                      <  1:654564/41‑1 (MQ=37)
                                        atCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAAcc                                                      <  1:2781730/41‑1 (MQ=255)
                                        atCTCCATGTGCAGCAGGCCGAGGAAGCCACAGCGGAAAcc                                                      <  1:3513699/41‑1 (MQ=255)
                                            ccATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAg                                             >  1:1390548/1‑46 (MQ=255)
                                             cATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGc                          <  1:3561832/64‑1 (MQ=255)
                                              aTGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGaa                  <  1:1803710/71‑1 (MQ=255)
                                               tgtgCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGaa                  >  1:1448838/1‑70 (MQ=255)
                                                 tgCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                <  1:2582610/70‑1 (MQ=255)
                                                 tgCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCAAGCGCGCTGGAGCTTTCCGGCTCATAGAACa                <  1:2953186/70‑1 (MQ=255)
                                                   cagcagGCCGAGGAAGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTg             <  1:1991650/71‑1 (MQ=255)
                                                        ggCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCa        <  1:2069787/71‑1 (MQ=255)
                                                        ggCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCa        <  1:3151291/71‑1 (MQ=255)
                                                           cGAGGCAGCCACAGCGGAACCCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGc         >  1:1910420/1‑67 (MQ=255)
                                                              ggAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTc  <  1:2331631/71‑1 (MQ=255)
                                                                aaGCCACAGCGGAAACCAAAGCCCAGCGCGCTGGa                                    >  1:1378041/1‑35 (MQ=255)
                                                                     |                                                               
CAGATCCAGATCGTATTCACGTTCCAGACGTTCCTGGATGATCTCCATGTGCAGCAGGCCGAGGAAGCCGCAGCGGAAACCAAAGCCCAGCGCGCTGGAGCTTTCCGGCTCATAGAACAGTGAGGCATCGTTC  >  Exported/1738722‑1738854

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: