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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported21,0900AG14.6% 121.4 / 13.2 48D341G (GAC→GGC) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (15/26);  new base G (2/5);  total (17/31)
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.72e-01

TCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCG  >  Exported/21021‑21145
                                                                     |                                                       
tCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgac                                                        <  1:3492836/71‑1 (MQ=255)
tCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgac                                                        <  1:1719252/71‑1 (MQ=255)
  cTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgacga                                                      <  1:4120248/71‑1 (MQ=255)
  cTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgacga                                                      <  1:4534671/71‑1 (MQ=255)
     ggCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACt                                                    <  1:2321174/70‑1 (MQ=255)
      gCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAt                                                  <  1:4701468/71‑1 (MQ=255)
       cGATCACGTTACCCTGGATGCCGGTAGCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAtg                                                 <  1:1336369/71‑1 (MQ=255)
       cGATCACGGTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAtg                                                 <  1:1571065/71‑1 (MQ=255)
         aTCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAtgtg                                               <  1:4527248/71‑1 (MQ=255)
           cACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGAt                                             >  1:600790/1‑71 (MQ=255)
                tACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAtg                                                 >  1:3590110/1‑62 (MQ=255)
                tACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAtg                                                 >  1:4141810/1‑62 (MQ=255)
                   ccTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGa                                     <  1:1079386/71‑1 (MQ=255)
                     tGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgacga                                                      >  1:644152/1‑52 (MQ=255)
                     tGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGgacga                                                      >  1:3348854/1‑52 (MQ=255)
                     tGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGaaa                                   >  1:3253573/1‑71 (MQ=255)
                        aTGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATAc                                <  1:2018067/71‑1 (MQ=255)
                         tGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGGCAGAAATAc                                <  1:3749990/70‑1 (MQ=255)
                                  cGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGaaa                      <  1:2513942/71‑1 (MQ=255)
                                  cGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGaaa                      <  1:3082428/71‑1 (MQ=255)
                                   ggTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAAc                     >  1:481543/1‑71 (MQ=255)
                                     tGCCGTTCACACCGCCCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCg                   <  1:4735041/71‑1 (MQ=255)
                                       ccGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGc                  >  1:212399/1‑70 (MQ=255)
                                       ccGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCt                 <  1:14860/71‑1 (MQ=255)
                                         gTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTaa               >  1:3128641/1‑71 (MQ=255)
                                         gTTCACACCGCGCCTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTaa               >  1:2863041/1‑71 (MQ=255)
                                         gTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGcc                            <  1:1669099/58‑1 (MQ=255)
                                             acacCGCGCCTGGCCACGGCCCGGACGACTATGTGAt                                             <  1:1536775/37‑1 (MQ=255)
                                             acacCGCGCCTGGCCACGGCCCGGACGACTATGTGAt                                             <  1:2562760/37‑1 (MQ=255)
                                             acacCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAAcc             <  1:2363249/69‑1 (MQ=255)
                                              cacCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCgg          <  1:1694309/71‑1 (MQ=255)
                                                ccGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTaa               >  1:123316/1‑64 (MQ=255)
                                                 cgcgCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTg       >  1:604077/1‑71 (MQ=255)
                                                 cgcgCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTg       >  1:1056838/1‑71 (MQ=255)
                                                 cgcgCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTg       >  1:4585780/1‑71 (MQ=255)
                                                   cgcCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCt                 <  1:3462547/59‑1 (MQ=255)
                                                   cgcCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCt                 <  1:3869450/59‑1 (MQ=255)
                                                   cgcCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCt                 <  1:210677/59‑1 (MQ=255)
                                                   cgcCTGGCAACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGc     >  1:1672451/1‑71 (MQ=255)
                                                     ccTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGccc   >  1:2686724/1‑71 (MQ=255)
                                                      cTGGCCACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCg  <  1:1239595/71‑1 (MQ=255)
                                                      cTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCg                   <  1:2071205/54‑1 (MQ=255)
                                                      cTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCg                   <  1:374248/54‑1 (MQ=255)
                                                      cTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCg  <  1:900429/71‑1 (MQ=255)
                                                       tGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAAc                     <  1:1757505/51‑1 (MQ=255)
                                                       tGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAAc                     <  1:782503/51‑1 (MQ=255)
                                                       tGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTaa               >  1:430475/1‑57 (MQ=255)
                                                       tGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCATGGAAAc                     <  1:2498692/51‑1 (MQ=255)
                                                          ccACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGc                             <  1:4453288/40‑1 (MQ=255)
                                                          ccACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGc                             <  1:3620182/40‑1 (MQ=255)
                                                          ccACGGCCCGGGCGACTATGTGATCGGTCAGAAATACGGc                             <  1:3281073/40‑1 (MQ=255)
                                                          ccACAGCCCGGGCGACTATGTGATCGGTCAGAAATACGGc                             <  1:2114790/40‑1 (MQ=255)
                                                                   ggACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCg                   >  1:2950984/1‑41 (MQ=255)
                                                                   ggACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCg                   >  1:829262/1‑41 (MQ=255)
                                                                     |                                                       
TCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCG  >  Exported/21021‑21145

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: