Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 53.5% 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,1010AC53.5% 144.9 / 157.7 146D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (33/35);  new base C (45/33);  total (78/68)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 149 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  W3110S.gb/236033‑236169
                                                                    |                                                                    
cgcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                    <  1:243449‑M2/60‑1 (MQ=255)
 tGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                   >  1:239280/1‑71 (MQ=255)
 gcgtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACAGGCCAGATCGTTCTCGCTAc                                                                   >  1:283862‑M2/11‑71 (MQ=255)
   tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                  >  1:3035584/1‑70 (MQ=255)
   tAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                  <  1:4155398/70‑1 (MQ=255)
     aaaCGACAAGTCTGAGATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                <  1:2375268/70‑1 (MQ=255)
     aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                               <  1:1705855/71‑1 (MQ=255)
     aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                <  1:167560/70‑1 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                    >  1:3269768‑M2/6‑65 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                               >  1:84422‑M2/6‑68 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                               <  1:3851497‑M2/65‑3 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                              >  1:1366635‑M2/6‑71 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                              >  1:209695‑M2/6‑71 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                              >  1:1112107‑M2/6‑71 (MQ=255)
       aCGACAAGTCTGACATAAATGACCTCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                             >  1:1847908/1‑71 (MQ=255)
       aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                             >  1:1063403/1‑71 (MQ=255)
       aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                             >  1:520847/1‑71 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:2365073/1‑71 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        >  1:2554213/1‑71 (MQ=255)
              gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                       >  1:2341911/1‑70 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                       <  1:1408901/69‑1 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                     <  1:3809013/71‑1 (MQ=255)
                 tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                      >  1:2239092/1‑52 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                  <  1:2142388/70‑1 (MQ=255)
                    cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                >  1:360545/1‑71 (MQ=255)
                     aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:1226518/1‑70 (MQ=255)
                     aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                               <  1:4053820/71‑1 (MQ=255)
                      tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:1973888/71‑1 (MQ=255)
                       aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:1286816/1‑68 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                                >  1:1793985/1‑67 (MQ=255)
                         aagACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                    <  1:1900233/60‑1 (MQ=255)
                          tGACCGCTATGAGCTCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                          <  1:3005111/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                        <  1:1896371/71‑1 (MQ=255)
                             ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                         <  1:3017172/69‑1 (MQ=255)
                              cGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTg                                      <  1:664919/71‑1 (MQ=255)
                               gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                     >  1:1825235/1‑71 (MQ=255)
                               gCTATGAGCACTGCAATAACCCGCCAGATCGTTCTCGATACCGa                                                                <  1:2442933/44‑1 (MQ=38)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                              >  1:710701/1‑45 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:3260674/1‑57 (MQ=255)
                                   tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                            <  1:678658/44‑1 (MQ=255)
                                   tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                             <  1:3982920/43‑1 (MQ=255)
                                    gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                          <  1:1644388/45‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCTCTACCGCAACCACCg                                                        <  1:2288696/46‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                        <  1:1458234/46‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                               >  1:1620798/1‑71 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                              >  1:741316/1‑40 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                              >  1:4231578/1‑40 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                       >  1:1507265/1‑47 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                               >  1:2390509/1‑71 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                          >  1:4134932/1‑58 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:2036713/71‑1 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:3263064/71‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:3746791/70‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             >  1:2883076/1‑70 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                            <  1:3586562/71‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCTTTCTCGCTACCGCAACCACCGCTATGAACCAGATTGGTGCGCACTATg                             >  1:1222860/1‑69 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                  <  1:3778021/48‑1 (MQ=39)
                                          tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                             <  1:1175565/68‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                             >  1:2574201/1‑35 (MQ=38)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGGGCACTATGa                            <  1:1792154/68‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                            <  1:1277665/68‑1 (MQ=255)
                                            cAATTACACTCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:3632208/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                >  1:85816/1‑47 (MQ=39)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                              <  1:2517281/49‑1 (MQ=39)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        >  1:650545/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGc                        >  1:2661233/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGAGCACTATGAAGGc                        >  1:4130371/1‑71 (MQ=255)
                                             aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                        <  1:3484295/54‑1 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATaaa                           >  1:1191461/1‑66 (MQ=255)
                                               ttACACGCCATATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                  <  1:394020/58‑1 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     <  1:4016643/71‑1 (MQ=255)
                                                tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                      >  1:3644913/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:3304016/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                 >  1:4235805/1‑57 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                                 >  1:3167327/1‑57 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                         >  1:1015378/1‑65 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:2135480/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     >  1:2478864/1‑69 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                     <  1:4081675/68‑1 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                  <  1:1627731/71‑1 (MQ=255)
                                                   acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa                 <  1:3580955/71‑1 (MQ=255)
                                                    cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                           <  1:476588/44‑1 (MQ=38)
                                                    cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                   >  1:72655/1‑68 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                           <  1:2181314/43‑1 (MQ=38)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:2332990/1‑70 (MQ=255)
                                                         gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                        <  1:4033696/58‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:880054/71‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:1744828/71‑1 (MQ=255)
                                                           tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          <  1:3593744/70‑1 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                  >  1:2123008/1‑45 (MQ=38)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt            >  1:4200531/1‑67 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga          >  1:1492889/1‑69 (MQ=255)
                                                             gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt       <  1:3527093/71‑1 (MQ=255)
                                                              ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                              <  1:1330540/47‑1 (MQ=255)
                                                               tctcGATACCGAAACCACCGGTATGAACCAGATTg                                         <  1:2092692/35‑1 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                   >  1:389933/1‑56 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg    >  1:940019/1‑71 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg     >  1:4025790/1‑70 (MQ=255)
                                                                  cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg  >  1:2162181/1‑71 (MQ=255)
                                                                   gATACCGAAACCACCGGTATCAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                <  1:2053177/56‑1 (MQ=255)
                                                                    |                                                                    
CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTG  >  W3110S.gb/236033‑236169

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: