Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 23,942 G→A 14.3% G518S (GGC→AGC)  ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb23,9420GA14.3% 266.6 / 30.4 98G518S (GGC→AGC) ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (49/35);  new base A (9/5);  total (58/40)
Fisher's exact test for biased strand distribution p-value = 7.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

GTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTA  >  W3110S.gb/23876‑24008
                                                                  |                                                                  
gTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGcga                                                                >  1:1805303/1‑71 (MQ=255)
gTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGcga                                                                >  1:1563777/1‑71 (MQ=255)
gTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGcga                                                                >  1:3888577/1‑71 (MQ=255)
  ggCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGcgac                                                               >  1:1442315/1‑70 (MQ=255)
   gCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGcgacga                                                             <  1:4867030/71‑1 (MQ=255)
   gCAAAACGCGTTGAAGTCGATAGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGcgacga                                                             <  1:2151616/71‑1 (MQ=255)
      aaaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGc                                                          >  1:2292051/1‑71 (MQ=255)
      aaaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGc                                                          >  1:3317682/1‑71 (MQ=255)
      aaaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGc                                                          >  1:1846126/1‑71 (MQ=255)
      aaaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGc                                                          >  1:788757/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCt                                                         >  1:131815/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCt                                                         >  1:137861/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCt                                                         >  1:329628/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCt                                                         <  1:1840098/71‑1 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGCGACGAAGCt                                                         >  1:4149519/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGCGACGAAGCt                                                         >  1:491541/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGCGACGAAGCt                                                         >  1:2178063/1‑71 (MQ=255)
       aaCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGCGACGAAGCt                                                         >  1:1035070/1‑71 (MQ=255)
        aCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTg                                                        <  1:3458634/71‑1 (MQ=255)
        aCGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTg                                                        <  1:3549555/71‑1 (MQ=255)
         cgcgTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGa                                                       >  1:4547125/1‑71 (MQ=255)
         cgcgTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGa                                                       >  1:3157804/1‑71 (MQ=255)
           cgTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATc                                                     <  1:2418432/71‑1 (MQ=255)
           cgTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATc                                                     <  1:3503339/71‑1 (MQ=255)
            gTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGaa                                                            <  1:3414499/63‑1 (MQ=255)
            gTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGaa                                                            <  1:903906/63‑1 (MQ=255)
              tGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAg                                                   >  1:1480667/1‑70 (MQ=255)
               gAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTa                                                 >  1:2203733/1‑71 (MQ=255)
               gAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTa                                                 >  1:1490563/1‑71 (MQ=255)
               gAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTa                                                 >  1:2483623/1‑71 (MQ=255)
                aaGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTAc                                                <  1:2170648/71‑1 (MQ=255)
                aaGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTAc                                                <  1:2112571/71‑1 (MQ=255)
                aaGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTAc                                                <  1:2927556/71‑1 (MQ=255)
                    cGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTg                                             >  1:5019523/1‑70 (MQ=255)
                    cGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTg                                             >  1:4579355/1‑70 (MQ=255)
                    cGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTg                                             >  1:2684413/1‑70 (MQ=255)
                     gATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGa                                            <  1:3335234/70‑1 (MQ=255)
                     gATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGa                                            <  1:1508907/70‑1 (MQ=255)
                     gATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGa                                            <  1:3590463/70‑1 (MQ=255)
                     gATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGa                                            >  1:2718902/1‑70 (MQ=255)
                         gCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAg                                                           >  1:3947015/1‑51 (MQ=255)
                          cATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGc                                      <  1:627672/71‑1 (MQ=255)
                          cATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTAAGTGAAAGTGc                                      <  1:5086727/71‑1 (MQ=255)
                          cATCCAGGCGGGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGc                                      <  1:2619924/71‑1 (MQ=255)
                            tCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGaaa                                          >  1:429953/1‑65 (MQ=255)
                            tCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGaaa                                          >  1:3343599/1‑65 (MQ=255)
                             ccAGGCGTGGTGGGATCTCGATGGGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCg                                    >  1:2164724/1‑70 (MQ=255)
                             ccAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGCCGAAGCTGATCAGTACGTGAAAGTGcc                                     >  1:1400945/1‑69 (MQ=255)
                             ccAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGa                                                       >  1:2537273/1‑51 (MQ=255)
                              cAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGa                                  >  1:2339102/1‑71 (MQ=255)
                              cAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGa                                  >  1:2206267/1‑71 (MQ=255)
                              cAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGa                                  >  1:2745426/1‑71 (MQ=255)
                                ggCGTGGTGGGATCTCGATGCGAAAGAGATCCGCGGCGACGAAGCTGATCAGTACGTGa                                            <  1:2219665/59‑1 (MQ=255)
                                 gCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGAt                                                      <  1:4195039/48‑1 (MQ=255)
                                 gCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGaca                                >  1:2412563/1‑70 (MQ=255)
                                    tggtggGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGaca                                <  1:2869282/67‑1 (MQ=255)
                                    tggtggGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGaca                                <  1:2387392/67‑1 (MQ=255)
                                      gtggGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGAcaca                              >  1:3323204/1‑67 (MQ=255)
                                      gtggGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTg                           <  1:1327256/70‑1 (MQ=255)
                                       tggGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATc                                                     <  1:2489295/43‑1 (MQ=255)
                                       tggGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGaa                                           <  1:4574025/53‑1 (MQ=255)
                                         gtaTCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGAGACATTGGATg                       <  1:240290/69‑1 (MQ=255)
                                             ctcGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGa                                            <  1:2169366/46‑1 (MQ=255)
                                             ctcGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAg                                         >  1:4385957/1‑49 (MQ=255)
                                             ctcGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAg                                         >  1:1671754/1‑49 (MQ=255)
                                             ctcGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATg                   >  1:1599772/1‑71 (MQ=255)
                                             ctcGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATg                   >  1:121019/1‑71 (MQ=255)
                                                gATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTa                                                 >  1:5078616/1‑38 (MQ=255)
                                                gATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTa                                                 >  1:3756936/1‑38 (MQ=255)
                                                gATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGa                         >  1:4848015/1‑62 (MQ=255)
                                                gATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATg                   >  1:2380827/1‑68 (MQ=255)
                                                gATGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATg                   <  1:3109606/68‑1 (MQ=255)
                                                gATGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATg                   <  1:3473663/68‑1 (MQ=255)
                                                gATGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATg                   <  1:2774740/68‑1 (MQ=255)
                                                 aTGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGttt               <  1:278040/71‑1 (MQ=255)
                                                 aTGCGAAAGAGATCCTCAGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGttt               <  1:2877759/71‑1 (MQ=255)
                                                  tGCGAAAGAGATCCTCGGCGGCGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGttt               >  1:2794098/1‑70 (MQ=255)
                                                   gCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGaa                                           >  1:704827/1‑41 (MQ=255)
                                                   gCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGaa                                           >  1:4810984/1‑41 (MQ=255)
                                                    cGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGtt                >  1:4258840/1‑67 (MQ=255)
                                                    cGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGtt                >  1:3841058/1‑67 (MQ=255)
                                                    cGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGAc            <  1:4746937/71‑1 (MQ=255)
                                                    cGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGAc            >  1:2271225/1‑71 (MQ=255)
                                                       aaGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTg                                       >  1:3207912/1‑41 (MQ=255)
                                                       aaGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGaca                                >  1:2870924/1‑48 (MQ=255)
                                                       aaGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGaca                                >  1:4725900/1‑48 (MQ=255)
                                                        agagATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGCCTCCg        >  1:1582625/1‑71 (MQ=255)
                                                         gagaTCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCgg       <  1:2211795/71‑1 (MQ=255)
                                                          agaTCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGa      >  1:1433529/1‑71 (MQ=255)
                                                            aTCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTg                           <  1:2086089/48‑1 (MQ=255)
                                                            aTCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTg                           <  1:1049203/48‑1 (MQ=255)
                                                            aTCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATc    <  1:1149607/71‑1 (MQ=255)
                                                             tCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCt   <  1:440072/71‑1 (MQ=255)
                                                              ccTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTa  >  1:2235304/1‑71 (MQ=255)
                                                              ccTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTa  >  1:2943094/1‑71 (MQ=255)
                                                               cTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCg        <  1:1171119/64‑1 (MQ=255)
                                                                tCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATg                       <  1:4734012/48‑1 (MQ=255)
                                                                 cGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGAcaca                              >  1:685678/1‑40 (MQ=255)
                                                                  |                                                                  
GTGGCAAAACGCGTTGAAGTCGATGGCATCCAGGCGTGGTGGGATCTCGATGCGAAAGAGATCCTCGGCGACGAAGCTGATCAGTACGTGAAAGTGCCGGACACATTGGATGTATGGTTTGACTCCGGATCTA  >  W3110S.gb/23876‑24008

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: