Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 236,101 A→C 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,1010AC93.3% 1390.5 / 44.5 420D12A (GAT→GCT) dnaQDNA polymerase III epsilon subunit
Reads supporting (aligned to +/- strand):  ref base A (28/0);  new base C (199/193);  total (227/193)
Fisher's exact test for biased strand distribution p-value = 1.90e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 420 total aligned reads displayed.

GACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGA  >  W3110S.gb/236042‑236154
                                                           |                                                     
ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                     <  1:2175950‑M2/60‑1 (MQ=255)
ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                     <  1:799367‑M2/60‑1 (MQ=255)
ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                     <  1:154591‑M2/60‑1 (MQ=255)
ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                     <  1:1162502‑M2/60‑1 (MQ=255)
ctCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                     <  1:1253314‑M2/60‑1 (MQ=255)
 tCAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTa                                                     <  1:1679384‑M2/60‑1 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:1525197/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:650836/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2918090/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2875421/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:1711114/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2825161/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2706839/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2674096/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2571200/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2546138/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:2503604/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:812394/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:205377/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:906058/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:111594/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:337179/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:1186885/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:580198/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:1374192/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:1359352/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                  >  1:1434398/1‑61 (MQ=255)
    aGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCg                                                  >  1:1077311/1‑61 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                               <  1:1103459/62‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                               <  1:620337/62‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                               <  1:836699/62‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                               <  1:1866842/62‑1 (MQ=255)
      tCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGcaa                                               <  1:710061/62‑1 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:1812044/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:1515354/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:274786/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:2378509/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:1413814/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:1673318/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:777931/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:292928/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgca                                                >  1:2036999/1‑56 (MQ=255)
         gACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCgc                                                 >  1:1626134/1‑56 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                          <  1:50839/61‑1 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                          <  1:1960621/61‑1 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                          <  1:2585282/61‑1 (MQ=255)
            aTAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAaccacc                                          <  1:1599878/61‑1 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                            <  1:694933/62‑1 (MQ=255)
                         aTGAGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGAt                            <  1:1989003/62‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:2011490/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:1780094/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:1908313/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:194971/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:255465/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:252312/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:1278425/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:2808939/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:2433753/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:1727813/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:2084522/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:542649/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:1656854/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:767147/53‑1 (MQ=255)
                            aGCACTGCAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGaa                                  <  1:38356/53‑1 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:2607892/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:88179/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:2772101/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGa                             >  1:1393300/1‑51 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:1833513/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:2452150/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:2402585/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:1235749/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:2843631/1‑62 (MQ=255)
                                   cAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCAc                  >  1:1023823/1‑62 (MQ=255)
                                     aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAg                              >  1:215893/1‑48 (MQ=39)
                                     aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                 >  1:225983/1‑61 (MQ=255)
                                     aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                 >  1:1575137/1‑61 (MQ=255)
                                     aTTACACGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACt                 >  1:2339369/1‑61 (MQ=255)
                                     aTTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACt                 >  1:2318861/1‑61 (MQ=255)
                                         cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa            >  1:155781/1‑62 (MQ=255)
                                         cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa            >  1:1479242/1‑62 (MQ=255)
                                         cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa            >  1:448959/1‑62 (MQ=255)
                                         cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa            >  1:247026/1‑62 (MQ=255)
                                         cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGaa            >  1:670903/1‑62 (MQ=255)
                                         cacGCCAGATCGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGa             >  1:2838615/1‑61 (MQ=255)
                                               aGATCGTTCTCGCTACCGCAACCCCCGGTATGAACCAGATTGGTGCGCACTATGAAGGcc        >  1:730491/1‑60 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGTCCACAAGa  <  1:1000337/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:946716/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:351480/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:2895032/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:2285592/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:212776/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:204857/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:1877952/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:184458/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:1496288/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:1329660/62‑1 (MQ=255)
                                                   cGTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGACACAAGa  <  1:2027945/62‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:789194/61‑1 (MQ=255)
                                                    gTTCTCGCTACCGCAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGa  <  1:939253/61‑1 (MQ=255)
                                                           |                                                     
GACAAGTCTGACATAAATGACCGCTATGAGCACTGCAATTACACGCCAGATCGTTCTCGATACCGAAACCACCGGTATGAACCAGATTGGTGCGCACTATGAAGGCCACAAGA  >  W3110S.gb/236042‑236154

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: