Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,107 A→C 49.0% E14A (GAA→GCA)  dnaQ → DNA polymerase III epsilon subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb236,1070AC49.0% 123.4 / inf 196E14A (GAA→GCA) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (45/55);  new base C (46/50);  total (91/105)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 204 total aligned reads displayed.

AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  W3110S.gb/236038‑236176
                                                                     |                                                                     
aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                      <  1:1666722/71‑1 (MQ=255)
 caactCAAGTCTGACATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATAGTTCTCGCTACCGcaa                                                                     >  1:3698407‑M2/6‑71 (MQ=255)
 caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                     >  1:1104171‑M2/6‑71 (MQ=255)
       aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                                <  1:1023851/70‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:4225817/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:663140/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                              <  1:5332769/71‑1 (MQ=255)
        aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                              <  1:776647/71‑1 (MQ=255)
         gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:2384784/71‑1 (MQ=255)
         gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:1283759/71‑1 (MQ=255)
            tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAAGCACCGGTa                                                            >  1:1171248/1‑69 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                         <  1:3895586/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                         <  1:2701764/70‑1 (MQ=255)
              aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                         <  1:1338160/70‑1 (MQ=255)
                aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                       >  1:48590/1‑70 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                     <  1:783994/71‑1 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                     <  1:4269405/71‑1 (MQ=255)
                 tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                              <  1:874379/62‑1 (MQ=255)
                  aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                                 >  1:4257101/1‑58 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                            >  1:103763/1‑62 (MQ=255)
                   aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                                    <  1:384981/70‑1 (MQ=255)
                      gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                >  1:2382316/1‑71 (MQ=255)
                      gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                                >  1:4101269/1‑71 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                               >  1:1953321/1‑71 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                               >  1:428507/1‑71 (MQ=255)
                       aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                               >  1:4420558/1‑71 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:4107474/53‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:5006649/53‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:959100/53‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:975683/53‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:182089/53‑1 (MQ=255)
                           cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                             <  1:1990488/53‑1 (MQ=255)
                             aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:4248912/1‑71 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                              <  1:109184/61‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                              <  1:5104902/61‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                              <  1:2297810/61‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                              <  1:2606879/61‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2114352/71‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:3457083/71‑1 (MQ=255)
                                  cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:5209071/71‑1 (MQ=255)
                                    cTGCAATTACACGTCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                     <  1:4072130/52‑1 (MQ=255)
                                     tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                                   <  1:3818411/69‑1 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                        >  1:4474558/1‑47 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:3829542/1‑62 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:4387636/1‑62 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:633921/1‑62 (MQ=255)
                                      gCAATTACACGCCAGATCGTTCTCGAAACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                         >  1:484162/1‑62 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                                 >  1:4601291/1‑37 (MQ=38)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                     >  1:2262080/1‑49 (MQ=39)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                                  >  1:3365078/1‑52 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               <  1:2860232/71‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               <  1:3555857/71‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               <  1:363238/71‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               <  1:2752216/71‑1 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                               >  1:1611789/1‑55 (MQ=255)
                                       cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                               >  1:1571302/1‑55 (MQ=255)
                                        aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                               <  1:913020/70‑1 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2827075/53‑1 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1464965/53‑1 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:895581/53‑1 (MQ=255)
                                                    gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2776822/53‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:1188375/71‑1 (MQ=255)
                                                     aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:3530714/71‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAGACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2137268/52‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1422107/52‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:4903378/52‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:2515536/52‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:4414247/52‑1 (MQ=255)
                                                     aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:70027/71‑1 (MQ=255)
                                                      tcattCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga                 <  1:474702/67‑1 (MQ=255)
                                                      tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:5491852/1‑68 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4519021/1‑68 (MQ=255)
                                                      tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2435110/1‑68 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4165671/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:1552032/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:3060835/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:3732704/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:5118256/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2991266/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4898276/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4891876/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:2909040/1‑67 (MQ=255)
                                                       cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4804798/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:4022527/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:1613878/1‑67 (MQ=255)
                                                       cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt                   >  1:1780325/1‑67 (MQ=255)
                                                         ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt             >  1:5573772/1‑71 (MQ=255)
                                                           ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                       >  1:3318303/1‑59 (MQ=255)
                                                               ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1405264/41‑1 (MQ=255)
                                                               ctACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                                    <  1:1717808/41‑1 (MQ=255)
                                                                tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg      <  1:2447877/71‑1 (MQ=255)
                                                                tACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCg      <  1:4326762/71‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag                <  1:3940931/59‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag                <  1:1947358/59‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag                <  1:454179/59‑1 (MQ=255)
                                                                  ccGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGTGCCGTTg   <  1:3135261/71‑1 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    >  1:5239375/1‑71 (MQ=255)
                                                                  ccGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGtt    >  1:4926408/1‑71 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGa  <  1:1677494/71‑1 (MQ=255)
                                                                    gAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAACATCATTGAGATTGGTGCCGTTGa  <  1:2913664/71‑1 (MQ=255)
                                                                     |                                                                     
AAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCCGTTGA  >  W3110S.gb/236038‑236176

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: