Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 74.8% 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,1010AC74.8% 414.6 / 89.1 170D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (24/19);  new base C (62/65);  total (86/84)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 171 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  W3110S.gb/236033‑236171
                                                                    |                                                                      
   gtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCtacag                                                                   >  1:490947‑M2/9‑69 (MQ=255)
   gtacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                                   <  1:1289611‑M2/63‑1 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                                 <  1:1082445‑M2/65‑3 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                                                >  1:1534398‑M2/6‑71 (MQ=255)
         ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAacca                                                             >  1:2081888‑M2/3‑71 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          <  1:138897/71‑1 (MQ=255)
            aaGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAGCGAAACCACCg                                                          <  1:359241/71‑1 (MQ=255)
              gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        <  1:1877129/71‑1 (MQ=255)
              gTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         >  1:1229164/1‑70 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGt                                                        >  1:1579273/1‑70 (MQ=255)
               tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                       <  1:714669/71‑1 (MQ=255)
                 tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                       <  1:93999/69‑1 (MQ=255)
                 tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:1638290/1‑70 (MQ=255)
                 tGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                      >  1:46055/1‑70 (MQ=255)
                 tGACATAAATGACCGCTATGAGCACAGCAATTACACGCCAGATCGTTCTCGATACCGa                                                                  <  1:623428/58‑1 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           <  1:1451202/63‑1 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                     <  1:67567/69‑1 (MQ=255)
                    cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:424707/1‑71 (MQ=255)
                    cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  <  1:1122217/71‑1 (MQ=255)
                     aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATg                                                     <  1:880138/67‑1 (MQ=255)
                      tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 <  1:174595/54‑3 (MQ=255)
                       aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:1043842/1‑68 (MQ=255)
                       aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:1636306/1‑68 (MQ=255)
                       aaaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCgg                                                         >  1:1111476/1‑61 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGa                                                    >  1:1129544/1‑65 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:1021700/1‑67 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGa                                              <  1:1256620/71‑1 (MQ=255)
                         aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                             >  1:228549/1‑71 (MQ=255)
                         aTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                             <  1:689144/71‑1 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAtt                                            <  1:1567866/71‑1 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTgg                                          <  1:1186257/71‑1 (MQ=255)
                             ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                                 >  1:1286890/1‑45 (MQ=255)
                             ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATAcc                                                                    >  1:551222/1‑44 (MQ=255)
                             ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGaaa                                                                >  1:268727/1‑48 (MQ=255)
                               gCTATTAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                       >  1:349353/1‑71 (MQ=255)
                               gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                       >  1:750570/1‑71 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:429098/1‑59 (MQ=255)
                                 tATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               <  1:2024350/61‑1 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCg                                                          >  1:1943535/1‑49 (MQ=255)
                                  aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                       >  1:1378416/1‑52 (MQ=255)
                                   tGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:519784/71‑1 (MQ=255)
                                    gAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:584581/70‑1 (MQ=255)
                                    gAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAcc                                                 <  1:1900510/56‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgc                                       >  1:250042/1‑65 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                                  <  1:1054267/70‑1 (MQ=255)
                                      gCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                <  1:1061847/71‑1 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:1299183/71‑1 (MQ=255)
                                       cACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTAt                                >  1:286174/1‑70 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:129356/70‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGAACTATg                               >  1:376056/1‑70 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCa                                    <  1:1522516/65‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:823890/69‑1 (MQ=255)
                                         cTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:2060984/69‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATg                                                     <  1:1736713/46‑1 (MQ=39)
                                          tGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAg                            >  1:900744/1‑71 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                    <  1:14603/47‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                               <  1:1088053/68‑1 (MQ=255)
                                           gCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAgg                           >  1:789887/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                                   >  1:1970680/1‑46 (MQ=39)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAAc                                                  >  1:1846998/1‑47 (MQ=39)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                          <  1:1191915/71‑1 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGaa                                                   >  1:194415/1‑46 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                                 <  1:1937689/64‑1 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                           <  1:1152686/36‑1 (MQ=38)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                                             <  1:209011/50‑1 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCca                        >  1:826558/1‑71 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:316924/71‑1 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                            >  1:1276762/1‑48 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:1504924/1‑69 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    <  1:1837429/55‑1 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       >  1:1835964/1‑68 (MQ=255)
                                                  cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAg                    <  1:612708/71‑1 (MQ=255)
                                                   acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                                               >  1:1808753/1‑43 (MQ=38)
                                                    cGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                  <  1:1263929/71‑1 (MQ=255)
                                                     gCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                                   <  1:1392823/53‑1 (MQ=255)
                                                       cAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCa                <  1:1032249/70‑1 (MQ=255)
                                                        aGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              <  1:1309064/71‑1 (MQ=255)
                                                        aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1069660/1‑71 (MQ=255)
                                                        aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              <  1:458408/71‑1 (MQ=255)
                                                         gATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:1170472/1‑70 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCTGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                       <  1:1598243/60‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAATATCATTga            <  1:221830/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:2038327/1‑69 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGAt                  <  1:1667104/65‑1 (MQ=255)
                                                           tCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCa                                    >  1:1310904/1‑46 (MQ=39)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt              >  1:416813/1‑67 (MQ=255)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         <  1:1831978/71‑1 (MQ=255)
                                                              ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAtt        >  1:1510811/1‑71 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt         <  1:142055/69‑1 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       >  1:1270405/1‑71 (MQ=255)
                                                               tctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg       >  1:597706/1‑71 (MQ=255)
                                                               tctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                             >  1:1550265/1‑49 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTgg                                          >  1:662507/1‑35 (MQ=37)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTgg                                          >  1:1757800/1‑35 (MQ=37)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      >  1:54771/1‑71 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg      >  1:659007/1‑71 (MQ=255)
                                                                ctcGATACCGAAACCACCGGTATGAACCAGATTgg                                          >  1:924695/1‑35 (MQ=255)
                                                                 tcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTgag           <  1:1328989/65‑1 (MQ=255)
                                                                  cGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTg    >  1:814885/1‑71 (MQ=255)
                                                                    aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGcc  <  1:899496/71‑1 (MQ=255)
                                                                    |                                                                      
CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGGTGCC  >  W3110S.gb/236033‑236171

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: