Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,696,316 A→G 12.0% G77G (GGT→GGC tadA ← tRNA‑specific adenosine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,696,3160AG12.0% 119.9 / 8.9 50G77G (GGT→GGCtadAtRNA‑specific adenosine deaminase
Reads supporting (aligned to +/- strand):  ref base A (16/28);  new base G (4/2);  total (20/30)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

ACATTACACATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCC  >  W3110S.gb/2696246‑2696384
                                                                      |                                                                    
acaTTACACATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGacca                                                                      <  1:2370468/71‑1 (MQ=255)
     acacATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTg                                                                 >  1:1705403/1‑71 (MQ=255)
     acacATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTg                                                                 >  1:1706906/1‑71 (MQ=255)
       acaTGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTccc                                                               >  1:4056560/1‑71 (MQ=255)
         aTGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGc                                                             <  1:3622797/71‑1 (MQ=255)
          tGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCg                                                              <  1:47057/69‑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)
                          tACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGaccac                                                                     <  1:3939353/46‑1 (MQ=255)
                          tACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCAt                                                     <  1:3010541/62‑1 (MQ=255)
                           aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                         >  1:3297314/1‑57 (MQ=255)
                           aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                         >  1:4295890/1‑57 (MQ=255)
                           aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                         >  1:1402392/1‑57 (MQ=255)
                           aCAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCGCCCTGCCGTAggg                                                         >  1:3661631/1‑57 (MQ=255)
                            cAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGacca                                                                      >  1:2355238/1‑43 (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:2370685/71‑1 (MQ=255)
                                        gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGAt                              <  1:2687500/71‑1 (MQ=255)
                                        gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGAt                              <  1:1026095/71‑1 (MQ=255)
                                        gATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGACCATGATTTCTGCAt                                          >  1:1836395/1‑59 (MQ=255)
                                          tCAGACGATAATTTTGCATCACCAGACCACCCTGcc                                                               <  1:2816347/36‑1 (MQ=255)
                                           cAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGAt                                                  <  1:471081/48‑1 (MQ=255)
                                           cAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCAt                                          >  1:3321545/1‑56 (MQ=255)
                                            aGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCAt                           <  1:3162480/70‑1 (MQ=255)
                                            aGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCAt                           <  1:1795135/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:1978768/1‑71 (MQ=255)
                                                     ttttGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAAt                 <  1:1187049/71‑1 (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)
                                                                      |                                                                    
ACATTACACATGGTTCAAGCGTGACATACAACGTGGCGTCGATCAGACGATAATTTTGCATCACCAGACCACCCTGCCGCAGGGCCATGATTTCTGCATGTGCGGTGGGATCATGGCGACCAATCGGGCGGTTCCAGCC  >  W3110S.gb/2696246‑2696384

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: