Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 56.9% 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,1010AC56.9% 126.7 / 103.7 102D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base C (30/28);  total (50/52)
Fisher's exact test for biased strand distribution p-value = 5.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 108 total aligned reads displayed.

AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  W3110S.gb/236031‑236171
                                                                      |                                                                      
aaCTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        <  1:437255/71‑1 (MQ=255)
       aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaa                                                                 <  1:602316/71‑1 (MQ=255)
         tcGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                               <  1:476949/70‑1 (MQ=255)
          cGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:1894217/1‑61 (MQ=255)
          cGACAAGTCTGACATAAATGAACGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:1604680/1‑61 (MQ=255)
          actCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:1016834‑M2/4‑69 (MQ=255)
            tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                            >  1:617256‑M2/2‑71 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           >  1:1826717/1‑71 (MQ=255)
             cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                            >  1:846564/1‑70 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           <  1:224803/70‑1 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccac                                                            <  1:2036180/69‑1 (MQ=255)
              aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:2152824/1‑71 (MQ=255)
               aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCgg                                                         >  1:1161628/1‑71 (MQ=255)
                gTCTGACATAAATGACCGCTATTAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGt                                                        <  1:81344/71‑1 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCACCCACCgg                                                         >  1:1174870/1‑70 (MQ=255)
                gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         >  1:14250/1‑70 (MQ=255)
                   tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:648936/1‑52 (MQ=255)
                     aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:1145759/1‑71 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         <  1:2066782/62‑1 (MQ=255)
                        tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         <  1:261036/62‑1 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGa                                                                        >  1:2090865/1‑45 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:546103/1‑67 (MQ=255)
                          aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTAt                                                      >  1:1843586/1‑63 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                             <  1:409702/71‑1 (MQ=255)
                           aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                             >  1:911496/1‑71 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           <  1:1154799/56‑1 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                            <  1:582882/71‑1 (MQ=255)
                            tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                            <  1:1573669/71‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTg                                           <  1:1519861/70‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:253225/70‑1 (MQ=255)
                              aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           <  1:1652217/70‑1 (MQ=255)
                                   tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        >  1:520069/1‑52 (MQ=255)
                                     tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:1362825/71‑1 (MQ=255)
                                      gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   >  1:531338/1‑70 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                               >  1:848738/1‑41 (MQ=255)
                                       aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                               >  1:1904894/1‑41 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:2123410/1‑65 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAt                                             >  1:583451/1‑56 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               >  1:839777/1‑70 (MQ=255)
                                          aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:1041112/70‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:83368/68‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             <  1:147901/71‑1 (MQ=255)
                                           cTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             <  1:1119111/71‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:664409/52‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                              <  1:898039/53‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                    <  1:468761/47‑1 (MQ=255)
                                            tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                    <  1:1015894/47‑1 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGt                                         >  1:2116980/1‑57 (MQ=255)
                                             gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                           >  1:2001991/1‑71 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                          >  1:2168163/1‑55 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:109742/1‑68 (MQ=255)
                                              cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          >  1:335191/1‑71 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:1893006/64‑1 (MQ=255)
                                                aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:738517/64‑1 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                              <  1:1462544/48‑1 (MQ=39)
                                                 ttACACGCCAGATCGTTCTCGATACCGAAACCACCg                                                          <  1:1929165/36‑1 (MQ=255)
                                                 ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:1018950/71‑1 (MQ=255)
                                                  tACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGacac                       >  1:2135451/1‑70 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                          >  1:555329/1‑50 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:1204659/1‑56 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:187022/1‑56 (MQ=255)
                                                   acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:1443272/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAg                                               >  1:1752967/1‑45 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:841367/1‑69 (MQ=255)
                                                   acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:816643/1‑69 (MQ=255)
                                                    cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:1679021/1‑68 (MQ=255)
                                                      cGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:2068902/1‑60 (MQ=255)
                                                       gCCAGATGGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:135107/56‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:1856432/56‑1 (MQ=255)
                                                       gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           >  1:2107829/1‑45 (MQ=255)
                                                       gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                           >  1:1341278/1‑45 (MQ=255)
                                                       gCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:1628256/57‑1 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:1744627/1‑71 (MQ=255)
                                                         cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt               >  1:1935024/1‑71 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              <  1:49692/71‑1 (MQ=255)
                                                          aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:123948/1‑71 (MQ=255)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg             <  1:655887/71‑1 (MQ=255)
                                                           gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTg             <  1:639031/71‑1 (MQ=255)
                                                           gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                  <  1:807437/66‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:858347/35‑1 (MQ=37)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                                <  1:1903107/35‑1 (MQ=37)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:1148050/71‑1 (MQ=255)
                                                            aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga            <  1:2182048/71‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACGAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:1272973/63‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCa                                                <  1:254611/35‑1 (MQ=255)
                                                            aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             <  1:1217106/54‑1 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATAGGTGCGCACTATGAAGGCCACAAGa                   >  1:1790035/1‑63 (MQ=255)
                                                             tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATAGGTGCGCACTATGAAGGCCACAAGa                   >  1:942146/1‑63 (MQ=255)
                                                             tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag           <  1:1812280/71‑1 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                     >  1:660053/1‑60 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                     >  1:1137453/1‑60 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1796270/1‑67 (MQ=255)
                                                              cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:296312/1‑67 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                      <  1:1896194/59‑1 (MQ=255)
                                                              cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1105466/1‑67 (MQ=255)
                                                               gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         <  1:1455900/71‑1 (MQ=255)
                                                                ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt        <  1:1287075/71‑1 (MQ=255)
                                                                  ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       >  1:1880411/1‑70 (MQ=255)
                                                                     gATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      <  1:754962/68‑1 (MQ=255)
                                                                      aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  <  1:2080218/71‑1 (MQ=255)
                                                                      |                                                                      
AACTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  W3110S.gb/236031‑236171

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: