Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 214,925 A→C 89.2% D12A (GAT→GCT)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported214,9250AC89.2% 2965.5 / 183.9 895D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (38/59);  new base C (422/376);  total (460/435)
Fisher's exact test for biased strand distribution p-value = 1.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 907 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  Exported/214857‑214991
                                                                    |                                                                  
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:3688484/1‑59 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:432205/1‑59 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:3439344/1‑59 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                        >  1:2805912/1‑69 (MQ=255)
             aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACATGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:5034342/71‑1 (MQ=255)
             aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:1000195/1‑64 (MQ=255)
             aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:3774484/1‑64 (MQ=255)
             aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     <  1:1465359/71‑1 (MQ=255)
              gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                     >  1:5083846/1‑70 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                        <  1:3379756/66‑1 (MQ=255)
                cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAc                                                                 >  1:4518616/1‑56 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:1294832/1‑52 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:760053/69‑1 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2410332/1‑71 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2899224/1‑71 (MQ=255)
                     aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:224452/1‑70 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACCCGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                              >  1:1151800/1‑67 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                             <  1:1922721/52‑3 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                            >  1:861099/1‑53 (MQ=255)
                         aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:1649540/68‑1 (MQ=255)
                         aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:4772421/68‑1 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:3991811/1‑52 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:2861363/1‑52 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:2443672/1‑52 (MQ=255)
                               gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                   <  1:552096/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:1102015/56‑1 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGCTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcgc                                 >  1:645405/1‑70 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:1390041/1‑44 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:4471251/1‑44 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:2502209/1‑63 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                >  1:5159803/1‑71 (MQ=255)
                                     aGCACTGCAATTACATGCCAGATCGTTCTCGCTACCGCAac                                                           >  1:1792054/1‑41 (MQ=37)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                         <  1:2766493/59‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1619523/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1723705/1‑71 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:4968466/71‑1 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:907735/71‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:457841/71‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:3089495/71‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                 <  1:2531947/47‑1 (MQ=39)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:978788/1‑56 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3340090/69‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:4820613/47‑1 (MQ=39)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:1083858/47‑1 (MQ=39)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                <  1:4057963/63‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:4358515/68‑1 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:828822/1‑46 (MQ=39)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:3277713/1‑46 (MQ=39)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                          >  1:1415613/1‑51 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                         >  1:1029955/1‑52 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                         >  1:2928019/1‑52 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:1974550/66‑1 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                       <  1:4853969/70‑1 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGAAaccac                                                        >  1:1897495/1‑37 (MQ=255)
                                             aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                          >  1:599103/1‑50 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                    >  1:1346151/1‑71 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                      <  1:3193473/52‑1 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                          >  1:42705/1‑48 (MQ=255)
                                                 acacGCCAGATTGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:3145174/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                     >  1:961309/1‑51 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                     >  1:4903141/1‑51 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                   >  1:1499719/1‑53 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                    >  1:2611551/1‑52 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:793581/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:2575515/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:4642220/1‑69 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        >  1:4718019/1‑47 (MQ=39)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:3854449/1‑62 (MQ=255)
                                                   tcGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:4201508/45‑1 (MQ=38)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:708798/36‑1 (MQ=37)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                <  1:1065472/36‑1 (MQ=37)
                                                      ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   <  1:2040154/64‑1 (MQ=255)
                                                      ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            <  1:1687738/71‑1 (MQ=255)
                                                      ccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            <  1:5127230/71‑1 (MQ=255)
                                                       cAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2654123/1‑35 (MQ=37)
                                                       cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa            <  1:2184701/70‑1 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3635314/54‑1 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          >  1:230053/1‑55 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          >  1:1998967/1‑55 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          <  1:3034141/71‑1 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3563774/53‑1 (MQ=255)
                                                         gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:2990261/40‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCTGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2475160/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:1525132/39‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:4401673/39‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                        <  1:5189391/39‑1 (MQ=37)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:2136527/53‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:32511/1‑69 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1565910/1‑69 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2731007/1‑69 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2360532/71‑1 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGGATGAACCAGATTGGTgc                                   >  1:4673504/1‑42 (MQ=37)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:1238289/71‑1 (MQ=255)
                                                             gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 >  1:3508368/1‑59 (MQ=255)
                                                             gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     >  1:2053935/1‑71 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:3913042/48‑1 (MQ=39)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt              <  1:4948093/61‑1 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt    <  1:3824671/71‑1 (MQ=255)
                                                               tctcGCTACCGCACCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg   >  1:4138390/1‑71 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 >  1:1851273/1‑56 (MQ=255)
                                                                    |                                                                  
CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  Exported/214857‑214991

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: