Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,915,239 A→G 12.3% E48G (GAG→GGG)  glyQ → glycine tRNA synthetase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,915,2390AG12.3% 167.6 / 13.4 65E48G (GAG→GGG) glyQglycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base G (5/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTAC  >  AP009048/3915173‑3915306
                                                                  |                                                                   
ccATTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCt                                                                  <  1:3067803/70‑1 (MQ=255)
ccATTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTg                                                                 <  1:539991/71‑1 (MQ=255)
  aTTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCt                                                                  <  1:2949987/68‑1 (MQ=255)
        cAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAg                                                          <  1:4039037/70‑1 (MQ=255)
        cAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAg                                                          <  1:827801/70‑1 (MQ=255)
         aaCCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTgg                                                                >  1:4566463/1‑63 (MQ=255)
         aaCCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGaa                                                        >  1:4212920/1‑71 (MQ=255)
          aCCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAAc                                                       >  1:2979235/1‑71 (MQ=255)
           ccATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAAcc                                                      >  1:2463782/1‑71 (MQ=255)
           ccATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAAcc                                                      >  1:3531490/1‑71 (MQ=255)
           ccATTGGACATGAAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAAcc                                                      >  1:4102757/1‑71 (MQ=255)
            cATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCg                                                     <  1:2840036/71‑1 (MQ=255)
                ggACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATg                                                  <  1:2016114/70‑1 (MQ=255)
                ggACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATgg                                                 >  1:21625/1‑71 (MQ=255)
                ggACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATgg                                                 >  1:3532797/1‑71 (MQ=255)
                ggACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATg                                                  >  1:3503423/1‑70 (MQ=255)
                  aCATGGAAGTCGGCGCGTGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATggcg                                               <  1:923754/71‑1 (MQ=255)
                   cATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATggcgg                                              >  1:3469975/1‑71 (MQ=255)
                      ggAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgct                                            >  1:3425323/1‑70 (MQ=255)
                        aaGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTg                                                                 <  1:1483888/47‑1 (MQ=255)
                          gTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGAt                                                   <  1:160581/59‑1 (MQ=255)
                                gcgGGAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGtt                                   >  1:4245245/1‑69 (MQ=255)
                                   ggAAGCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATg                                                  <  1:4608193/51‑1 (MQ=255)
                                   ggAACCTCTCACCCAATGACCTGTCTGCGCGGGCTGGGGCCAg                                                          <  1:2454484/43‑1 (MQ=255)
                                    gAACCTCTCACCCAATGACCTGTCTGCGCGATCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCt                             >  1:3468450/1‑71 (MQ=255)
                                    gAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGtt                                   <  1:1317128/65‑1 (MQ=255)
                                      aCCTCTCACCCAATGACCTTTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCa                                 >  1:1266343/1‑65 (MQ=255)
                                       ccTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAAc                                                       <  1:3008339/42‑1 (MQ=255)
                                       ccTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGAt                                                   >  1:1903607/1‑46 (MQ=255)
                                         tctcACCCAATGACCTGTCTGCGCGATCTGGGGCCAGa                                                         <  1:3917498/38‑1 (MQ=255)
                                          ctcACCCAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCg                                                     >  1:3745343/1‑41 (MQ=255)
                                          ctcACCCAATGACCTGTCTGCGCGATCTGGGGCCAGAACCGATGGCGGCTGCTTa                                       <  1:1056570/55‑1 (MQ=255)
                                          ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCtt                            >  1:2116824/1‑66 (MQ=255)
                                          ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcgt                       >  1:1770257/1‑71 (MQ=255)
                                          ctcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcg                        <  1:1368518/70‑1 (MQ=255)
                                           tcACCCAATGACCTGTCTGCGCGATCTGGGGCCAGAACCGATGGCGGCTGCTTa                                       <  1:4432742/54‑1 (MQ=255)
                                           tcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGAt                                                   >  1:2513126/1‑42 (MQ=255)
                                           tcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATggcgg                                              >  1:1968851/1‑47 (MQ=255)
                                           tcACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcg                        <  1:84763/69‑1 (MQ=255)
                                            cACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcgtcg                     >  1:657681/1‑71 (MQ=255)
                                            cACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcgtcg                     >  1:3986464/1‑71 (MQ=255)
                                             aCCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCtcgtc                      <  1:3300468/69‑1 (MQ=255)
                                             aCCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGc                    >  1:3682105/1‑71 (MQ=255)
                                              cccAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGgctgct                                         <  1:2818571/49‑1 (MQ=255)
                                              cccAATGACCTGGCTGCGCGAGCGGGGGCCAGAACCGATgg                                                 <  1:1688040/41‑1 (MQ=25)
                                               ccAATGACCTGTCTGCGCGGGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGca                              >  1:957578/1‑58 (MQ=255)
                                               ccAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGccc                  >  1:899472/1‑71 (MQ=255)
                                               ccAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGcc                   >  1:1215509/1‑70 (MQ=255)
                                                cAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGa                                                    >  1:1106977/1‑36 (MQ=255)
                                                cAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGa                                                    >  1:460688/1‑36 (MQ=255)
                                                cAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTa                                       >  1:2743991/1‑49 (MQ=255)
                                                  aTGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTAt                                      <  1:3477871/48‑1 (MQ=255)
                                                  aTGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCccgac               >  1:4368824/1‑71 (MQ=255)
                                                   tGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATggc                                                <  1:6238/37‑1 (MQ=255)
                                                     aCCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCccgaccg             <  1:4249214/70‑1 (MQ=255)
                                                       cTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATg          <  1:4508480/71‑1 (MQ=255)
                                                       cTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATg          <  1:906411/71‑1 (MQ=255)
                                                         gTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATg                                     >  1:1740452/1‑42 (MQ=255)
                                                          tCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATgg         <  1:3585449/69‑1 (MQ=255)
                                                            tGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATg          >  1:4278549/1‑66 (MQ=255)
                                                            tGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATg          <  1:226996/66‑1 (MQ=255)
                                                              cgcgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGc     <  1:2935541/69‑1 (MQ=255)
                                                              cgcgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCt    <  1:3188956/70‑1 (MQ=255)
                                                              cgcgAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTa   <  1:2725159/71‑1 (MQ=255)
                                                                 gAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTAc  <  1:204347/69‑1 (MQ=255)
                                                                  aGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCccga                >  1:3539610/1‑54 (MQ=255)
                                                                  |                                                                   
CCATTGTTCAACCATTGGACATGGAAGTCGGCGCGGGAACCTCTCACCCAATGACCTGTCTGCGCGAGCTGGGGCCAGAACCGATGGCGGCTGCTTATGTTCAGCCTTCTCGTCGCCCGACCGATGGTCGCTAC  >  AP009048/3915173‑3915306

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: