Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,616,051 T→C 9.7% intergenic (+366/‑27) prfC → / → osmY peptide chain release factor RF‑3/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,616,0510TC 9.7% 289.2 / 19.1 108intergenic (+366/‑27)prfC/osmYpeptide chain release factor RF‑3/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (55/42);  new base C (2/8);  total (58/50)
Fisher's exact test for biased strand distribution p-value = 4.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

Only 100 of 108 total aligned reads displayed.

TGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTGCTG  >  W3110S.gb/4615985‑4616116
                                                                  |                                                                 
tGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTAc                                                                 >  1:2549272/1‑69 (MQ=255)
   cTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCa                                                            <  1:1700759/71‑1 (MQ=255)
   cTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCa                                                            <  1:3380493/71‑1 (MQ=255)
      cTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCaaa                                                          <  1:869022/70‑1 (MQ=255)
      cTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCaaa                                                          <  1:4065766/70‑1 (MQ=255)
                tGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                               <  1:605096/71‑1 (MQ=255)
                tGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                               <  1:2074054/71‑1 (MQ=255)
                 gAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                               <  1:3214203/70‑1 (MQ=255)
                  aaGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGc                                                             >  1:3082977/1‑55 (MQ=255)
                  aaGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                               >  1:1643913/1‑69 (MQ=255)
                  aaGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                               >  1:1648427/1‑69 (MQ=255)
                    gTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGATCAAACATAc                                                     <  1:677284/61‑1 (MQ=255)
                     tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                          <  1:3130752/55‑1 (MQ=255)
                     tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                          <  1:2062398/55‑1 (MQ=255)
                         tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                       >  1:1199283/1‑70 (MQ=255)
                          cccTATTCAGCAATGCAACCTCGTATTGCCAGGCTCAAATAACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:394533/1‑71 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:1123990/52‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:1246450/52‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:1805545/52‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:2958349/52‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       <  1:3933771/68‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                       <  1:3064874/68‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                       <  1:355515/68‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                       <  1:1482279/68‑1 (MQ=255)
                           ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                       <  1:1940260/68‑1 (MQ=255)
                              aTTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:3346369/1‑64 (MQ=255)
                              aTTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:1474704/1‑64 (MQ=255)
                              aTTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:2021416/1‑64 (MQ=255)
                              aTTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:2026988/1‑64 (MQ=255)
                              aTTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:2969652/1‑64 (MQ=255)
                               ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATaaa                                            >  1:1866103/1‑59 (MQ=255)
                               ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:3767559/1‑63 (MQ=255)
                                   gCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTa                                  <  1:1839984/65‑1 (MQ=255)
                                   gCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTa                                  <  1:2997753/65‑1 (MQ=255)
                                   gCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                            <  1:1627482/71‑1 (MQ=255)
                                   gCAATGCAACCTCGTGTTGCCAGGCTCAAATCACTAGCAAATATACAGGAATAAATCGATGACTATGACaa                            <  1:23900/71‑1 (MQ=255)
                                      aTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                          <  1:3674339/38‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGCGCAAACATAc                                                     >  1:761093/1‑42 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACa                                                        <  1:1395502/39‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:1491804/40‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:278022/40‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                       <  1:3443148/40‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       <  1:634993/56‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       <  1:2000552/56‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       <  1:2701442/56‑1 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGa                               >  1:1946434/1‑64 (MQ=255)
                                       tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGa                               >  1:2106980/1‑64 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                      >  1:2660061/1‑39 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                      >  1:460987/1‑39 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATaaa                                            >  1:1567994/1‑49 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATaaa                                            >  1:2952378/1‑49 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       >  1:2507553/1‑54 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       >  1:1849360/1‑54 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       >  1:54559/1‑54 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                       >  1:3532113/1‑54 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:2939321/1‑56 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:2944007/1‑56 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:3097639/1‑56 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                       >  1:3747428/1‑70 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                       >  1:2080188/1‑70 (MQ=255)
                                         cAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                       >  1:3460271/1‑70 (MQ=255)
                                            ccTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                            <  1:1966332/62‑1 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATaaa                                            >  1:471206/1‑44 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATa                                              <  1:3729243/42‑1 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATa                                              <  1:1988774/42‑1 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATc                                          <  1:1922543/46‑1 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:3169803/1‑51 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:1732046/1‑51 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                          >  1:1881201/1‑62 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                          >  1:3343315/1‑62 (MQ=255)
                                              tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGa                      <  1:3571777/66‑1 (MQ=255)
                                                gtgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:3294814/1‑49 (MQ=255)
                                                gtgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:3952509/1‑49 (MQ=255)
                                                 tgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTAt                                 <  1:17107/52‑1 (MQ=255)
                                                 tgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTAt                                 <  1:3713470/52‑1 (MQ=255)
                                                 tgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTAt                                 <  1:2674979/52‑1 (MQ=255)
                                                 tgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTAt                                 <  1:2726437/52‑1 (MQ=255)
                                                   ttGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                               <  1:652422/36‑1 (MQ=255)
                                                   ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGAc                                    <  1:856015/47‑1 (MQ=255)
                                                   ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGAc                                    <  1:311862/47‑1 (MQ=255)
                                                   ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGa             >  1:2603135/1‑70 (MQ=255)
                                                     gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                        >  1:2372777/1‑41 (MQ=255)
                                                     gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTa                                  >  1:3855377/1‑47 (MQ=255)
                                                     gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTa                                  >  1:3921236/1‑47 (MQ=255)
                                                     gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTa                                  >  1:2548521/1‑47 (MQ=255)
                                                        aGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGAc                              <  1:1897936/48‑1 (MQ=255)
                                                         ggCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAActct      >  1:1371429/1‑71 (MQ=255)
                                                         ggCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAActct      >  1:179703/1‑71 (MQ=255)
                                                          gCTCAAATTACGAGCAAACATACAGGAATAAATCGATGAc                                    <  1:3897066/40‑1 (MQ=255)
                                                          gCTCAAATTACGAGCAAACATACAGGAATAAATCGATGAc                                    <  1:3460301/40‑1 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATg                                      >  1:1526043/1‑36 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:2614802/1‑37 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:2067586/1‑37 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgct   >  1:3762837/1‑71 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgct   >  1:3735534/1‑71 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgct   >  1:304686/1‑71 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgct   >  1:3453576/1‑71 (MQ=255)
                                                            tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgct   >  1:1792059/1‑71 (MQ=255)
                                                             cAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgctg  >  1:2264095/1‑71 (MQ=255)
                                                              aaaTTACGAGCAAACATACAGGAATAAATCGATGa                                     >  1:4036630/1‑35 (MQ=255)
                                                                  |                                                                 
TGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTGCTG  >  W3110S.gb/4615985‑4616116

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: