Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,616,051 T→C 12.9% 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,0510TC12.9% 218.3 / 23.0 85intergenic (+366/‑27)prfC/osmYpeptide chain release factor RF‑3/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (6/5);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

TTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTGCTGGCT  >  W3110S.gb/4615984‑4616119
                                                                   |                                                                    
ttGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACg                                                                   <  1:3099390/71‑1 (MQ=255)
 tGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACg                                                                   >  1:3487715/1‑70 (MQ=255)
  gACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGa                                                                  <  1:2587857/70‑1 (MQ=255)
  gACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAg                                                                 >  1:6881/1‑71 (MQ=255)
       cTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             <  1:3043867/70‑1 (MQ=255)
       cTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             <  1:875691/70‑1 (MQ=255)
       cTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             <  1:1029892/70‑1 (MQ=255)
       cTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             <  1:1983294/70‑1 (MQ=255)
          tAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGACCAAAc                                                            <  1:1772048/68‑1 (MQ=255)
             agGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAg                                                      <  1:488360/71‑1 (MQ=255)
                  gAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAAt                                                  >  1:3877956/1‑70 (MQ=255)
                   aaGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACAt                                                          >  1:2761068/1‑61 (MQ=255)
                    aGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAgg                                                     >  1:3458926/1‑65 (MQ=255)
                      tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATa                                                 <  1:940250/67‑1 (MQ=255)
                      tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATa                                                 <  1:146264/67‑1 (MQ=255)
                      tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATc                                             <  1:3948955/71‑1 (MQ=255)
                      tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATc                                             <  1:3271574/71‑1 (MQ=255)
                      tctctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATc                                             <  1:1161876/71‑1 (MQ=255)
                        tctcCCTATTCAGCAATGCAACCTCGTGTTGCTAGGCTCAAATTACGAGCAAACATAc                                                        >  1:4033585/1‑58 (MQ=255)
                        tctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACg                                                                   >  1:2893110/1‑47 (MQ=255)
                        tctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAAt                                              >  1:832508/1‑68 (MQ=255)
                         ctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          >  1:3542997/1‑71 (MQ=255)
                         ctcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          >  1:3973937/1‑71 (MQ=255)
                          tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCg                                            <  1:2147908/68‑1 (MQ=255)
                          tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          >  1:4503222/1‑70 (MQ=255)
                            ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                          <  1:4519484/68‑1 (MQ=255)
                                ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTa                                     >  1:1806029/1‑69 (MQ=255)
                                ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTAt                                    <  1:1075973/70‑1 (MQ=255)
                                ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTAt                                    >  1:2434318/1‑70 (MQ=255)
                                  cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAAc                                                            >  1:41687/1‑44 (MQ=255)
                                  cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGa                                  >  1:576117/1‑70 (MQ=255)
                                  cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGa                                  >  1:3792379/1‑70 (MQ=255)
                                  cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGa                                  >  1:3338236/1‑70 (MQ=255)
                                   aGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                         <  1:1348560/46‑1 (MQ=255)
                                   aGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGAc                                 >  1:2108506/1‑70 (MQ=255)
                                   aGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACa                                <  1:4445710/71‑1 (MQ=255)
                                   aGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACa                                <  1:4364124/71‑1 (MQ=255)
                                     cAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             >  1:1099376/1‑40 (MQ=255)
                                     cAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             >  1:4033772/1‑40 (MQ=255)
                                     cAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAg                              <  1:529126/71‑1 (MQ=255)
                                     cAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAg                              <  1:173045/71‑1 (MQ=255)
                                      aaTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACAt                                                          >  1:1712350/1‑42 (MQ=255)
                                      aaTGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                             >  1:1791047/1‑71 (MQ=255)
                                        tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                         <  1:2781909/41‑1 (MQ=255)
                                        tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          <  1:1385227/56‑1 (MQ=255)
                                        tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          <  1:3705902/56‑1 (MQ=255)
                                        tGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          <  1:3700759/56‑1 (MQ=255)
                                         gCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                         >  1:3820437/1‑40 (MQ=255)
                                         gCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                           >  1:2209563/1‑54 (MQ=255)
                                         gCAACCACGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                             >  1:152789/1‑68 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATTACGAGCaaa                                                             >  1:2134890/1‑35 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGa                                           >  1:3110529/1‑53 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          <  1:2112536/54‑1 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                        >  1:854530/1‑56 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGa                                        >  1:712004/1‑56 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                          <  1:397545/70‑1 (MQ=255)
                                          cAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                          >  1:1639391/1‑70 (MQ=255)
                                           aaCCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGaa                        <  1:228810/71‑1 (MQ=255)
                                            aCCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                               >  1:1356334/1‑63 (MQ=255)
                                            aCCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                          >  1:3310742/1‑68 (MQ=255)
                                             ccTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAgg                                                     <  1:780693/40‑1 (MQ=255)
                                             ccTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATa                                                 <  1:3272644/44‑1 (MQ=255)
                                                cGTGTTGCCAGGCTCAATTTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                               <  1:4633521/59‑1 (MQ=255)
                                                 gtgtTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTc                  >  1:2476370/1‑71 (MQ=255)
                                                  tgtTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCg                 <  1:3816522/71‑1 (MQ=255)
                                                    ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGaa                        <  1:606061/62‑1 (MQ=255)
                                                    ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaa               <  1:2811000/71‑1 (MQ=255)
                                                    ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGa                >  1:4259491/1‑70 (MQ=255)
                                                      gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaaaa             >  1:4234213/1‑71 (MQ=255)
                                                      gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaaaa             >  1:1475855/1‑71 (MQ=255)
                                                           gCTCAAATTACGAGCAAACATACAGGAATAAATCGAt                                          <  1:983957/37‑1 (MQ=255)
                                                           gCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaa               <  1:1312194/64‑1 (MQ=255)
                                                            cTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                               <  1:4422926/47‑1 (MQ=255)
                                                             tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                               <  1:742175/46‑1 (MQ=255)
                                                             tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGaa                        <  1:1523613/53‑1 (MQ=255)
                                                             tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGaa                        <  1:1910000/53‑1 (MQ=255)
                                                               aaaTTACGAGCAAACATACAGGAATAAATCGATGa                                        >  1:2314920/1‑35 (MQ=255)
                                                               aaaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGAt                     <  1:429467/54‑1 (MQ=255)
                                                               aaaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaaaa             >  1:2601822/1‑62 (MQ=255)
                                                               aaaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaaaa             >  1:3722646/1‑62 (MQ=255)
                                                               aaaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTgctgg    <  1:3567076/71‑1 (MQ=255)
                                                                aaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAActc          <  1:3394589/64‑1 (MQ=255)
                                                                 aTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGa                <  1:4117956/57‑1 (MQ=255)
                                                                 aTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTgctggct  >  1:722802/1‑71 (MQ=255)
                                                                 aTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTgctggct  >  1:1087654/1‑71 (MQ=255)
                                                                   tACGAGCAAAAATACAGGAATAAATCGATGAATATGACAAGa                             >  1:1526700/1‑42 (MQ=37)
                                                                   |                                                                    
TTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTGCTGGCT  >  W3110S.gb/4615984‑4616119

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: