Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,815,599 C→T 5.4% A102T (GCT→ACT)  trmH ← tRNA (Guanosine‑2'‑O‑)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,815,5990CT 5.4% 306.7 / 5.1 93A102T (GCT→ACT) trmHtRNA (Guanosine‑2'‑O‑)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (41/47);  new base T (2/3);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

ATGCCCGTTTTCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGCCAA  >  W3110S.gb/3815529‑3815667
                                                                      |                                                                    
aTGCCCGTTTTCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGc                                                                      <  1:1930717/71‑1 (MQ=255)
 tGCCCGTTTTCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCg                                                                     <  1:2624752/71‑1 (MQ=255)
   cccGTTTTCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGtt                                                                   <  1:409561/71‑1 (MQ=255)
         ttCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAg                                                              >  1:1034369/1‑70 (MQ=255)
         ttCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAg                                                              >  1:1897213/1‑70 (MQ=255)
           ctcCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGa                                                             >  1:241273/1‑69 (MQ=255)
           ctcCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGa                                                             >  1:3307050/1‑69 (MQ=255)
              cTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTAt                                                                 >  1:3351157/1‑62 (MQ=255)
              cTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGaaa                                                           <  1:240297/68‑1 (MQ=255)
              cTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGaaa                                                           <  1:1646649/68‑1 (MQ=255)
              cTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGa                                                         <  1:1791963/70‑1 (MQ=255)
                 tCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATggg                                                     <  1:1396301/71‑1 (MQ=255)
                 tCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATggg                                                     <  1:1666039/71‑1 (MQ=255)
                   ccATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGt                                                                    >  1:179488/1‑54 (MQ=255)
                   ccATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGt                                                                    >  1:1813038/1‑54 (MQ=255)
                    cATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAATATGGGTTg                                                  <  1:2946003/71‑1 (MQ=255)
                     aTCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGaaa                                                           <  1:4610906/61‑1 (MQ=255)
                      tCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGtt                                                   <  1:3220812/68‑1 (MQ=255)
                      tCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGcc                                                <  1:2657149/71‑1 (MQ=255)
                      tCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGc                                                 >  1:2657506/1‑70 (MQ=255)
                      tCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGc                                                 >  1:2401113/1‑70 (MQ=255)
                      tCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGc                                                 >  1:2621403/1‑70 (MQ=255)
                      tCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGc                                                 >  1:1798212/1‑70 (MQ=255)
                        aaaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGc                                                 >  1:1225200/1‑68 (MQ=255)
                         aaaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGa                                             <  1:4067132/71‑1 (MQ=255)
                          aaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGaa                                            <  1:3048786/71‑1 (MQ=255)
                          aaTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGaa                                            >  1:50980/1‑71 (MQ=255)
                           aTGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAAt                                           >  1:309271/1‑71 (MQ=255)
                            tGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAAt                                           >  1:2273392/1‑70 (MQ=255)
                            tGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAAt                                           >  1:935102/1‑70 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCATTTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         <  1:3893644/70‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                        <  1:3900000/71‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                        <  1:165601/71‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         <  1:2852185/70‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         <  1:2846458/70‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTg                                        <  1:690970/71‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTg                                        <  1:104663/71‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTGGCCAGAATCTg                                        <  1:4050629/71‑1 (MQ=255)
                              cAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGCTGCCAGAATCTg                                        <  1:4949764/71‑1 (MQ=255)
                               aGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATGGGTTGCCAGAATCTg                                        <  1:795620/70‑1 (MQ=255)
                               aGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGc                                       <  1:3181079/71‑1 (MQ=255)
                                 gTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATc                                          <  1:3241449/66‑1 (MQ=255)
                                   cGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATggg                                                     <  1:836061/53‑1 (MQ=255)
                                   cGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGt                                                    >  1:2169631/1‑54 (MQ=255)
                                       cGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTg                                        <  1:1405264/62‑1 (MQ=255)
                                         aGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGaaa                                                           <  1:2279485/41‑1 (MQ=255)
                                           tgtAATCAATTTCGCGGAAATCGACAGCGTTATCATaaa                                                           <  1:2220691/39‑1 (MQ=255)
                                           tgtAATCAATTTCGCGGAAATCGACAGCGTTATCAGa                                                             <  1:3878933/37‑1 (MQ=255)
                                             tAATCAATTTCGCGGAAATCGACATCGTTATCAGAAAGATGGGTTGCCAGa                                             <  1:2633481/51‑1 (MQ=255)
                                             tAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGa                                             <  1:425609/51‑1 (MQ=255)
                                             tAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         >  1:3677458/1‑55 (MQ=255)
                                             tAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         >  1:1383863/1‑55 (MQ=255)
                                             tAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         >  1:3952570/1‑55 (MQ=255)
                                             tAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCt                                         >  1:2335760/1‑55 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGAt                                                        <  1:2057105/38‑1 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGAt                                                        <  1:2653805/38‑1 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGAt                                                        <  1:1400106/38‑1 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGAt                                                        <  1:44661/38‑1 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGAt                                                        <  1:4405116/38‑1 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATTGGTTGCCAGAATCTGCATGCCCTGGCCTTTga                       <  1:350550/71‑1 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCAt                                     >  1:1558642/1‑57 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCtat                         >  1:847809/1‑67 (MQ=255)
                                               aTCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTga                       <  1:3421119/71‑1 (MQ=255)
                                                tCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGAt                                                        <  1:2020700/37‑1 (MQ=255)
                                                 cAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATgg                                                      >  1:4704817/1‑38 (MQ=255)
                                                 cAATTTCGCGGAAATCGACAGTGTTATCAGAAAGATgg                                                      >  1:4543519/1‑38 (MQ=255)
                                                  aaTTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGcc                                                >  1:4110449/1‑43 (MQ=255)
                                                  aaTTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGcc                                                >  1:4553234/1‑43 (MQ=255)
                                                  aaTTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGAt                    <  1:3032613/71‑1 (MQ=255)
                                                     ttCGCGGAAATCGACAGCGTTATCAGAAAGATggg                                                     >  1:287817/1‑35 (MQ=255)
                                                     ttCGCGGAAATCGACAGCGTTATCAGAAAGATggg                                                     >  1:2229424/1‑35 (MQ=255)
                                                     ttCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGa                                             <  1:1349956/43‑1 (MQ=255)
                                                     ttCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGa                                             <  1:3182711/43‑1 (MQ=255)
                                                     ttCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGa                                             <  1:4483216/43‑1 (MQ=255)
                                                      tCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGc                >  1:4349817/1‑71 (MQ=255)
                                                          ggAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACg            >  1:3264428/1‑71 (MQ=255)
                                                           gAAATCGACAGCGTTATCTGAAAGATGGGTTGCCAGAATc                                          <  1:639859/40‑1 (MQ=255)
                                                           gAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCttt                         <  1:119948/57‑1 (MQ=255)
                                                            aaaTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACgg           >  1:3338038/1‑70 (MQ=255)
                                                            aaaTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACgg           >  1:3331875/1‑70 (MQ=255)
                                                            aaaTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACgg           >  1:2246982/1‑70 (MQ=255)
                                                            aaaTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc          <  1:986419/71‑1 (MQ=255)
                                                            aaaTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGc          <  1:2437552/71‑1 (MQ=255)
                                                              aTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTgaga                     >  1:1817869/1‑58 (MQ=255)
                                                              aTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTgaga                     >  1:1895604/1‑58 (MQ=255)
                                                              aTCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCAt        >  1:185646/1‑71 (MQ=255)
                                                                cGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCg      >  1:3857291/1‑71 (MQ=255)
                                                                cGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCg      >  1:3260265/1‑71 (MQ=255)
                                                                cGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCg      >  1:3239672/1‑71 (MQ=255)
                                                                  aCAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGc     >  1:4289162/1‑70 (MQ=255)
                                                                  aCAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGc     >  1:1978934/1‑70 (MQ=255)
                                                                  aCAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGc     >  1:94566/1‑70 (MQ=255)
                                                                    aGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGCCaa  <  1:4858168/71‑1 (MQ=255)
                                                                    aGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGCCaa  <  1:2896827/71‑1 (MQ=255)
                                                                      cGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGcc    <  1:3919083/67‑1 (MQ=255)
                                                                      cGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGcc    <  1:3913335/67‑1 (MQ=255)
                                                                      cGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGcc    <  1:3268605/67‑1 (MQ=255)
                                                                      |                                                                    
ATGCCCGTTTTCTCCTGTCCCATCAAAATGCAGGTCGGGCGAGTGTAATCAATTTCGCGGAAATCGACAGCGTTATCAGAAAGATGGGTTGCCAGAATCTGCATGCCCTGGCCTTTGAGATGAGCGACGGCATCGCCAA  >  W3110S.gb/3815529‑3815667

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: