Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,696,325 C→T 11.5% L74L (CTG→CTA tadA ← tRNA‑specific adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,696,3250CT11.5% 142.6 / 8.2 52L74L (CTG→CTAtadAtRNA‑specific adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base C (17/29);  new base T (4/2);  total (21/31)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCCTTCGCC  >  W3110S.gb/2696255‑2696390
                                                                      |                                                                 
aTGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGc                                                                   <  1:3622797/71‑1 (MQ=255)
      cAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAggg                                                               >  1:3985758/1‑69 (MQ=255)
       aaGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGcc                                                             <  1:3143841/70‑1 (MQ=255)
        aGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCAt                                                           <  1:1783590/71‑1 (MQ=255)
        aGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCAt                                                           <  1:2113897/71‑1 (MQ=255)
           gTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGAt                                                        <  1:4409675/71‑1 (MQ=255)
           gTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGAt                                                        <  1:1887225/71‑1 (MQ=255)
              aCATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGAttt                                                      <  1:2824627/70‑1 (MQ=255)
                 tACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCAt                                                           <  1:3010541/62‑1 (MQ=255)
                  aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                               >  1:4295890/1‑57 (MQ=255)
                  aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                               >  1:1402392/1‑57 (MQ=255)
                  aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                               >  1:3661631/1‑57 (MQ=255)
                  aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                               >  1:3297314/1‑57 (MQ=255)
                   cAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTg                                                   >  1:965/1‑68 (MQ=255)
                   cAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCa                                                 >  1:2020975/1‑70 (MQ=255)
                   cAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCa                                                 >  1:2818913/1‑70 (MQ=255)
                   cAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCAt                                                >  1:1410841/1‑71 (MQ=255)
                   cAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCAt                                                >  1:3225622/1‑71 (MQ=255)
                        tGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATg                                                          <  1:2601543/56‑1 (MQ=255)
                        tGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCg                                           <  1:159153/71‑1 (MQ=255)
                          gCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGt                                         <  1:1637944/71‑1 (MQ=255)
                             tCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTg                                                   <  1:924854/58‑1 (MQ=255)
                               gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGAt                                    <  1:1026095/71‑1 (MQ=255)
                               gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGAt                                    <  1:2687500/71‑1 (MQ=255)
                               gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGAt                                    <  1:2370685/71‑1 (MQ=255)
                               gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGACCATGATTTCTGCAt                                                >  1:1836395/1‑59 (MQ=255)
                                  cAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGAt                                                        <  1:471081/48‑1 (MQ=255)
                                  cAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCAt                                                >  1:3321545/1‑56 (MQ=255)
                                   aGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCAt                                 <  1:1795135/70‑1 (MQ=255)
                                   aGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCAt                                 <  1:3162480/70‑1 (MQ=255)
                                     aCGATAATTTTGCATCACCAGAACACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGc                              <  1:1185858/71‑1 (MQ=255)
                                      cGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCg                             <  1:809719/71‑1 (MQ=255)
                                       gATAATTTTGCATCACGAGACCACCGTGCCGCAGGGCCATGATTTCTGCATGTGCGGt                                         <  1:3627464/58‑1 (MQ=39)
                                       gATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCg                                           <  1:477797/56‑1 (MQ=255)
                                         tAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGa                                                         >  1:347736/1‑40 (MQ=255)
                                         tAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCAtgtg                                             >  1:1376759/1‑52 (MQ=255)
                                          aaTTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGAt                                                        >  1:192250/1‑40 (MQ=255)
                                            ttttGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAAt                       <  1:1187049/71‑1 (MQ=255)
                                            ttttGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAAt                       >  1:1978768/1‑71 (MQ=255)
                                            ttttGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGATCATGGCGACCAATc                      <  1:2198455/71‑1 (MQ=255)
                                                gCATCTCCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGAt                                    <  1:3846857/54‑1 (MQ=255)
                                                  aTCACCAGACCGCCCTGCCGTAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGTCGACCAATCGGACg                 <  1:3948464/71‑1 (MQ=255)
                                                   tCACCAGACCGCCCTGCCGTAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGTCGACCAATCGGACg                 <  1:2429601/70‑1 (MQ=255)
                                                    cACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGt               >  1:402648/1‑71 (MQ=255)
                                                     aCCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGt               <  1:4042151/70‑1 (MQ=255)
                                                          accaccCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGc         >  1:1881376/1‑71 (MQ=255)
                                                           ccaccCTTCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGcc        <  1:2788670/71‑1 (MQ=255)
                                                           ccaccCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTACAGcc        <  1:1939323/71‑1 (MQ=255)
                                                            caccCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGAcc                          <  1:3731863/52‑1 (MQ=255)
                                                              cccTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCCTTc     <  1:346208/71‑1 (MQ=255)
                                                                 tGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCCTTCGcc  >  1:3564029/1‑71 (MQ=255)
                                                                    cGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCCtt      >  1:3717930/1‑64 (MQ=255)
                                                                    cGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCCtt      >  1:2293189/1‑64 (MQ=255)
                                                                      |                                                                 
ATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCCTTCGCC  >  W3110S.gb/2696255‑2696390

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: