Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 61.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,1010AC61.5% 249.4 / 136.5 158D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base C (49/48);  total (78/80)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 158 total aligned reads displayed.

AACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTG  >  W3110S.gb/236039‑236166
                                                              |                                                                 
caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                            >  1:304750‑M2/6‑68 (MQ=255)
caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                            <  1:1072238‑M2/65‑3 (MQ=255)
 aactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                          >  1:1644982‑M2/5‑71 (MQ=255)
    tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                       >  1:778568‑M2/2‑71 (MQ=255)
    tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                       >  1:131006‑M2/2‑71 (MQ=255)
    aCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAaccac                                                       >  1:639234/1‑71 (MQ=255)
       aGTCTGACATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    >  1:2687040/1‑71 (MQ=255)
       aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacc                                                         <  1:1569766/66‑1 (MQ=255)
       aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:443998/71‑1 (MQ=255)
       aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:1310181/71‑1 (MQ=255)
           tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                 >  1:1179559/1‑70 (MQ=255)
             aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                <  1:389655/69‑1 (MQ=255)
              cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                             >  1:1445659/1‑71 (MQ=255)
                tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                       <  1:2039662/59‑1 (MQ=255)
                tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                            >  1:399619/1‑70 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                       >  1:1291837/1‑58 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACAGCAaccac                                                       >  1:1781532/1‑58 (MQ=255)
                 aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                          >  1:1345628/1‑55 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                          <  1:2975761/70‑1 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                          <  1:2637324/70‑1 (MQ=255)
                  aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                           >  1:570570/1‑69 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                <  1:1763204/62‑1 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTa                                                  >  1:2519216/1‑60 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                       <  1:3615504/71‑1 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                       <  1:355809/71‑1 (MQ=255)
                    tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                       <  1:308638/71‑1 (MQ=255)
                     gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                      >  1:1122769/1‑71 (MQ=255)
                      aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                      <  1:2810635/70‑1 (MQ=255)
                      aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                      <  1:2175756/70‑1 (MQ=255)
                      aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                              >  1:2576117/1‑62 (MQ=255)
                          cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                 >  1:1837101/1‑71 (MQ=255)
                           tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                               <  1:1986626/56‑1 (MQ=255)
                             tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:842645/49‑1 (MQ=255)
                             tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                        <  1:298307/61‑1 (MQ=255)
                             tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                              >  1:1507451/1‑71 (MQ=255)
                             tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                              <  1:601354/71‑1 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                             <  1:3652424/70‑1 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                            >  1:1430097/1‑71 (MQ=255)
                               aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                           >  1:3250861/1‑56 (MQ=255)
                                gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                         <  1:1918796/57‑1 (MQ=255)
                                gCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                           <  1:1929646/71‑1 (MQ=255)
                                 cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:3459523/71‑1 (MQ=255)
                                 cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:619702/71‑1 (MQ=255)
                                   cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          >  1:2252280/1‑69 (MQ=255)
                                    tGCAATTACACGCCAGATCGTTCTCGCTACCGCAac                                                          <  1:3221797/36‑1 (MQ=38)
                                    tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                     >  1:3315030/1‑41 (MQ=255)
                                     gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                      >  1:2928611/1‑71 (MQ=255)
                                     gCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACAAGATTGGTGCGCACTATGAAgg                      >  1:886102/1‑71 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                    >  1:3447728/1‑56 (MQ=255)
                                      cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                               >  1:3372638/1‑61 (MQ=255)
                                       aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:1719089/39‑1 (MQ=38)
                                       aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:2673917/39‑1 (MQ=38)
                                       aaTTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:887832/39‑1 (MQ=38)
                                        aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                         <  1:3589853/65‑1 (MQ=255)
                                        aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                   >  1:1299917/1‑71 (MQ=255)
                                         ttACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                    <  1:3507358/37‑1 (MQ=38)
                                         ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                           <  1:1014464/46‑1 (MQ=39)
                                         ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgcg                                 <  1:2328292/56‑1 (MQ=255)
                                          tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                              >  1:2933083/1‑42 (MQ=38)
                                          tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTa                            >  1:3013918/1‑60 (MQ=255)
                                          tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                                         >  1:1102695/1‑47 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCACCCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                   >  1:1933698/1‑68 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                        >  1:1747825/1‑63 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  >  1:797105/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  >  1:1396514/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  >  1:363724/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  >  1:2307619/1‑69 (MQ=255)
                                           acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCa                                           <  1:1704940/44‑1 (MQ=255)
                                           acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                  >  1:1258611/1‑69 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg               <  1:1649251/71‑1 (MQ=255)
                                            cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg               <  1:3177293/71‑1 (MQ=255)
                                              cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt             <  1:205411/71‑1 (MQ=255)
                                               gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc            <  1:1004515/71‑1 (MQ=255)
                                                ccAGATCGTTCTCGCTACCGCAACCCCCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc            <  1:83578/70‑1 (MQ=255)
                                                 cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                         >  1:1024491/1‑40 (MQ=38)
                                                  agctCGTTCTCGCTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         <  1:728153/68‑1 (MQ=255)
                                                  aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                        >  1:877412/1‑40 (MQ=38)
                                                  aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                         <  1:2743702/55‑1 (MQ=255)
                                                  aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         <  1:2266333/71‑1 (MQ=255)
                                                  aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa           <  1:1862823/69‑1 (MQ=255)
                                                   gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTg                                   <  1:3292664/44‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                           <  1:2901253/51‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                       >  1:3134863/1‑55 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                   <  1:2393204/59‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATc            <  1:1386426/66‑1 (MQ=255)
                                                    aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga       <  1:2396305/71‑1 (MQ=255)
                                                     tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:2363089/51‑1 (MQ=255)
                                                     tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                          <  1:2092654/51‑1 (MQ=255)
                                                     tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag      <  1:664927/71‑1 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         >  1:1912814/1‑67 (MQ=255)
                                                      cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         >  1:2150163/1‑67 (MQ=255)
                                                      cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         >  1:3671523/1‑67 (MQ=255)
                                                      cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt         >  1:2390760/1‑67 (MQ=255)
                                                       gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                        >  1:1528907/1‑51 (MQ=255)
                                                       gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgaga     <  1:3176692/70‑1 (MQ=255)
                                                        ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt   >  1:3117729/1‑71 (MQ=255)
                                                        ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt   <  1:784854/71‑1 (MQ=255)
                                                          ctctATACCGAAACCACCGGTATGAACCAGATTGATGCGCAc                              >  1:1283432/1‑42 (MQ=38)
                                                          ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                        >  1:850571/1‑48 (MQ=39)
                                                          ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg  >  1:2598438/1‑70 (MQ=255)
                                                              |                                                                 
AACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTG  >  W3110S.gb/236039‑236166

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: