Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 236,101 A→C 54.4% 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,1010AC54.4% 112.8 / 109.8 103D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (22/25);  new base C (22/34);  total (44/59)
Fisher's exact test for biased strand distribution p-value = 5.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236033‑236167
                                                                    |                                                                  
cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                  >  1:276579/1‑71 (MQ=255)
cTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATa                                                                  >  1:961069/1‑71 (MQ=255)
    tacaactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACcgc                                                              <  1:491490‑M2/64‑2 (MQ=255)
     aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              <  1:2499849/70‑1 (MQ=255)
     aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              <  1:2105311/70‑1 (MQ=255)
     aaaCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGa                                                              <  1:1328592/70‑1 (MQ=255)
      caactCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTAccgca                                                             <  1:486074‑M2/65‑3 (MQ=255)
       aCGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAac                                                           <  1:198854/71‑1 (MQ=255)
           cAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                                       >  1:447510/1‑71 (MQ=255)
                cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTa                                                   >  1:1817781/1‑70 (MQ=255)
                cTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  >  1:1814185/1‑71 (MQ=255)
                   aCATAAATGCCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:47374/1‑71 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  <  1:1017208/68‑1 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTAt                                                  <  1:2273727/68‑1 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                               >  1:2330018/1‑71 (MQ=255)
                   aCATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                <  1:1178316/70‑1 (MQ=255)
                    cATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:212168/1‑71 (MQ=255)
                      tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:2363355/71‑1 (MQ=255)
                      tAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCa                                            <  1:2209958/71‑1 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:1805544/1‑67 (MQ=255)
                        aaTGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAAc                                              >  1:876464/1‑67 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:2266683/71‑1 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:261642/71‑1 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:1734298/71‑1 (MQ=255)
                          tGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:1480896/71‑1 (MQ=255)
                           gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:1902729/1‑44 (MQ=255)
                           gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:1058766/1‑44 (MQ=255)
                           gACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                                  >  1:567581/1‑44 (MQ=255)
                            aCCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                       <  1:1201157/70‑1 (MQ=255)
                             ccGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGa                                                <  1:1177561/60‑1 (MQ=255)
                               gCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTgc                                   >  1:961242/1‑71 (MQ=255)
                                cTATGAGCACTGCAATTACACTCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:2318437/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:2318661/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:2007878/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:1810935/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:2039960/71‑1 (MQ=255)
                                cTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTgcg                                  <  1:377901/71‑1 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCTATACCGAAACCACCGGTATg                                                 >  1:1041411/1‑51 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                 >  1:1801684/1‑35 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATAc                                                                 >  1:253171/1‑35 (MQ=255)
                                     aGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTa                             >  1:1023266/1‑71 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa                          <  1:359908/71‑1 (MQ=255)
                                        aCTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:422676/70‑1 (MQ=255)
                                          tGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:649002/68‑1 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1209640/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1908115/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1887440/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1772938/1‑71 (MQ=255)
                                            cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:1303534/1‑71 (MQ=255)
                                            cAATTACACGCAAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGc                      >  1:724322/1‑71 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:336256/64‑1 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2130813/64‑1 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:117042/64‑1 (MQ=255)
                                              aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2324087/64‑1 (MQ=255)
                                               ttACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   <  1:1103415/71‑1 (MQ=255)
                                                tatacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1672210/4‑64 (MQ=255)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:2481494/1‑64 (MQ=255)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1356951/1‑64 (MQ=255)
                                                tACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa                         >  1:1889097/1‑64 (MQ=255)
                                                 acacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcac                   >  1:397898/1‑69 (MQ=255)
                                                 acacGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                               >  1:2555826/1‑57 (MQ=255)
                                                   acGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa               <  1:1373364/71‑1 (MQ=255)
                                                   acGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGAtt                                        <  1:1394147/46‑1 (MQ=255)
                                                        aGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          <  1:2314137/71‑1 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2542705/1‑70 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:555240/1‑70 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2178627/1‑70 (MQ=255)
                                                         gATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2455222/1‑70 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATg                           <  1:2231734/52‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCcaca                  <  1:872544/61‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2423550/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:648706/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1372308/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:337836/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:334762/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:257931/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2266368/71‑1 (MQ=255)
                                                          aTCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:2269428/71‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:2330154/41‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:1414040/41‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:1491259/41‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTgg                                      <  1:1626944/41‑1 (MQ=255)
                                                          aTCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1299425/71‑1 (MQ=255)
                                                           tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAt           >  1:743278/1‑67 (MQ=255)
                                                           tCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:648491/70‑1 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2210960/1‑67 (MQ=255)
                                                            cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:2640062/1‑67 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTg                                       >  1:2524759/1‑38 (MQ=255)
                                                            cGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCAtt          >  1:1130309/1‑67 (MQ=255)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:1860619/71‑1 (MQ=255)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:1940365/71‑1 (MQ=255)
                                                             gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     <  1:136122/71‑1 (MQ=255)
                                                             gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     >  1:1681434/1‑71 (MQ=255)
                                                             gTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGAt     >  1:1777640/1‑71 (MQ=255)
                                                              ttCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACaa                 <  1:1175271/58‑1 (MQ=255)
                                                              ttCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1975726/67‑1 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:1124258/1‑71 (MQ=255)
                                                                ctcGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTgg  >  1:977232/1‑71 (MQ=255)
                                                                ctcGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTg   <  1:1821912/70‑1 (MQ=255)
                                                                    aTACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTga        <  1:1792326/61‑1 (MQ=255)
                                                                    |                                                                  
CTGTAAAACGACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGATCATTGAGATTGG  >  W3110S.gb/236033‑236167

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: