Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 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
*W3110S.gb236,1010AC89.2% 2511.0 / 138.6 778D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base C (378/315);  total (422/356)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 782 total aligned reads displayed.

ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  W3110S.gb/236032‑236169
                                                                     |                                                                    
   tgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAcc                                                                  <  1:1935666‑M2/62‑1 (MQ=255)
      acaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                                >  1:4047681‑M2/7‑69 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                               <  1:1031986‑M2/65‑3 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                               <  1:3608960‑M2/65‑3 (MQ=255)
       caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                              >  1:3684428‑M2/6‑71 (MQ=255)
         actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                            <  1:2467942‑M2/68‑1 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                         >  1:514885/1‑71 (MQ=255)
            cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                              >  1:3944948/1‑66 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:3422963/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:2876267/1‑71 (MQ=255)
             aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:2147471/1‑71 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                       <  1:4312600/71‑1 (MQ=255)
              aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                       <  1:2530571/71‑1 (MQ=255)
                tCTTACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                               <  1:212197/61‑3 (MQ=255)
                  tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    >  1:2713080/1‑70 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCTCTACCGCAACCACCGGTAt                                                    <  1:970246/68‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                    <  1:2046831/68‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                   <  1:88409/69‑1 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                 >  1:568677/1‑71 (MQ=255)
                    aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                    <  1:2497957/68‑1 (MQ=255)
                     cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:1079357/1‑71 (MQ=255)
                      aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                          >  1:629439/1‑60 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                  <  1:1447752/67‑1 (MQ=255)
                       tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:258365/71‑1 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:1614789/1‑68 (MQ=255)
                        aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                >  1:4097257/1‑68 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:3594180/1‑67 (MQ=255)
                         aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                             <  1:94184/70‑1 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                           <  1:3077387/71‑1 (MQ=255)
                          aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                           <  1:782479/71‑1 (MQ=255)
                           tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTAGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:427434/71‑1 (MQ=255)
                            gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                             >  1:3263829/1‑51 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:4182890/70‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:4075833/70‑1 (MQ=255)
                             aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                         <  1:1969873/70‑1 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                             >  1:714582/1‑48 (MQ=255)
                               cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:4467250/63‑1 (MQ=255)
                                gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                     <  1:2276099/71‑1 (MQ=255)
                                 cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                             <  1:4148199/46‑1 (MQ=255)
                                    tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                <  1:1140328/56‑1 (MQ=255)
                                     gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:2838718/1‑47 (MQ=255)
                                     gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                 >  1:3758819/1‑70 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                               >  1:3873054/1‑71 (MQ=255)
                                      aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                               <  1:1317942/71‑1 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                  <  1:2911674/51‑1 (MQ=255)
                                       gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              <  1:1570526/71‑1 (MQ=255)
                                        cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              >  1:1545368/1‑70 (MQ=255)
                                         aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                       <  1:3900076/44‑1 (MQ=255)
                                         aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                <  1:3642284/51‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                               <  1:1674815/67‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                           <  1:2996335/71‑1 (MQ=255)
                                          cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             >  1:4405403/1‑69 (MQ=255)
                                           tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:3290307/56‑1 (MQ=255)
                                            gCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                         <  1:2326867/39‑1 (MQ=38)
                                             cAATTAGACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:589856/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCTCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:137487/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                         >  1:1503426/1‑70 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:1855605/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:65953/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:1784586/1‑71 (MQ=255)
                                             cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:2415404/1‑71 (MQ=255)
                                               aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                      >  1:4214316/1‑71 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                <  1:4493436/44‑1 (MQ=38)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                            <  1:353559/64‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     <  1:1194505/71‑1 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                      >  1:4542424/1‑70 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                      >  1:2763586/1‑70 (MQ=255)
                                                ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                              <  1:84078/46‑1 (MQ=255)
                                                 tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              >  1:2197803/1‑61 (MQ=255)
                                                  acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         >  1:3326400/1‑49 (MQ=39)
                                                   cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                  <  1:3463367/71‑1 (MQ=255)
                                                   cacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                  <  1:750494/71‑1 (MQ=255)
                                                    acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 <  1:812841/71‑1 (MQ=255)
                                                     ctccAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                <  1:3832838/69‑1 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                <  1:2081285/71‑1 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                <  1:3138305/71‑1 (MQ=255)
                                                     cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                <  1:469397/71‑1 (MQ=255)
                                                      gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                            >  1:1765837/1‑58 (MQ=255)
                                                        cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt             >  1:1808781/1‑71 (MQ=255)
                                                         aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:2670876/54‑1 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         >  1:4444760/1‑41 (MQ=38)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:1216325/1‑70 (MQ=255)
                                                          gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:2949056/1‑70 (MQ=255)
                                                          gATCGTTCTCGATACAGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 <  1:1422213/65‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:3566827/40‑1 (MQ=38)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:4349899/52‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc                       <  1:3200892/58‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:3996346/71‑1 (MQ=255)
                                                           aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:2804750/71‑1 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:3501226/38‑1 (MQ=37)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:2016888/39‑1 (MQ=37)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:4243160/1‑68 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag         <  1:1000192/71‑1 (MQ=255)
                                                            tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:2376663/70‑1 (MQ=255)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         >  1:2003997/1‑38 (MQ=37)
                                                             cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGCCACAAGATCAtt            >  1:2631597/1‑67 (MQ=255)
                                                             cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                           >  1:3791563/1‑52 (MQ=255)
                                                              gTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          <  1:338291/36‑1 (MQ=37)
                                                               ttCTCGCTACCGCAATCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                    <  1:1249725/57‑1 (MQ=255)
                                                                 ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                          >  1:1716473/1‑49 (MQ=39)
                                                                     |                                                                    
ACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  W3110S.gb/236032‑236169

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: