Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 21,090 A→G 55.2% D341G (GAC→GGC)  ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported21,0900AG55.2% 6.3 / 34.5 29D341G (GAC→GGC) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (2/11);  new base G (9/7);  total (11/18)
Fisher's exact test for biased strand distribution p-value = 5.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGC  >  Exported/21023‑21151
                                                                   |                                                             
cTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgacga                                                            <  1:2321071/71‑1 (MQ=255)
  cGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgacgac                                                           >  1:471828/1‑70 (MQ=255)
   ggCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTa                                                         <  1:3604775/71‑1 (MQ=255)
    gCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTa                                                         <  1:617845/70‑1 (MQ=255)
         cACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACt                                                          >  1:1641356/1‑64 (MQ=255)
            gTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCgg                                                <  1:1788474/71‑1 (MQ=255)
             ttACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGggcgac                                                           >  1:1018290/1‑59 (MQ=255)
             ttACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGggcgac                                                           >  1:3354694/1‑59 (MQ=255)
             ttACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGggcgac                                                           >  1:392841/1‑59 (MQ=255)
             ttACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGggcgac                                                           >  1:2243589/1‑59 (MQ=255)
                  cTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGa                                           >  1:592569/1‑70 (MQ=255)
                     gATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGa                                                    <  1:3620259/58‑1 (MQ=255)
                     gATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGa                                                    <  1:3615921/58‑1 (MQ=255)
                       tGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGt                                               >  1:619778/1‑61 (MQ=255)
                       tGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGGCAGAAATACg                                     <  1:3364932/71‑1 (MQ=255)
                          cGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGcc                                  <  1:679003/71‑1 (MQ=255)
                          cGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGcc                                  <  1:564034/71‑1 (MQ=255)
                             tACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTg                                >  1:1640239/1‑70 (MQ=255)
                              aCCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTgg                               >  1:3401199/1‑70 (MQ=255)
                              aCCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTgg                               >  1:3427512/1‑70 (MQ=255)
                               ccGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGa                              <  1:3145380/70‑1 (MQ=255)
                                      cGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCt                       <  1:1306266/70‑1 (MQ=255)
                                          cacaCCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTg                                >  1:2221960/1‑57 (MQ=255)
                                          cacaCCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAAcc                   <  1:1553955/70‑1 (MQ=255)
                                              ccGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGtt              <  1:1029548/71‑1 (MQ=255)
                                               cgcgCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTg             >  1:2668103/1‑71 (MQ=255)
                                               cgcgCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTg             >  1:904601/1‑71 (MQ=255)
                                                 cgcCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGc           >  1:902725/1‑71 (MQ=255)
                                                   ccTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGcc          <  1:1206698/70‑1 (MQ=255)
                                                    cTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCg                 <  1:1224830/62‑1 (MQ=255)
                                                    cTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCg                 <  1:680962/62‑1 (MQ=255)
                                                      ggCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCt                       <  1:1328300/54‑1 (MQ=255)
                                                        ccACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTg                                <  1:2125285/43‑1 (MQ=255)
                                                          aCGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACgg   >  1:842849/1‑70 (MQ=255)
                                                          aCGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGc  <  1:2621759/71‑1 (MQ=255)
                                                                   acgacTATGTGATCGGTCAGAAATACGGACTGGAAACCGCTAACCCGGTTGGCCCGGAc     <  1:190913/59‑1 (MQ=255)
                                                                   |                                                             
CTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGC  >  Exported/21023‑21151

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: