Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,616,051 T→C 11.3% 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,0510TC11.3% 158.0 / 11.8 62intergenic (+366/‑27)prfC/osmYpeptide chain release factor RF‑3/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (33/22);  new base C (3/4);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTGCTGGCTG  >  W3110S.gb/4615982‑4616120
                                                                     |                                                                     
cATTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTa                                                                      >  1:915776/1‑71 (MQ=255)
cATTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTa                                                                      >  1:666464/1‑71 (MQ=255)
 aTTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTAc                                                                     <  1:1989786/71‑1 (MQ=255)
    gACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAg                                                                  >  1:1931167/1‑71 (MQ=255)
    gACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAg                                                                  >  1:2456289/1‑71 (MQ=255)
    gACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAg                                                                  >  1:145124/1‑71 (MQ=255)
             agagGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          <  1:41934/70‑1 (MQ=255)
              gagGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:448975/1‑69 (MQ=255)
                gggTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAgg                                                      <  1:79971/71‑1 (MQ=255)
                gggTGAAGTCTATCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAgg                                                      <  1:1171314/71‑1 (MQ=255)
                  gTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAg                                                                  <  1:1758000/57‑1 (MQ=255)
                  gTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAg                                                                  <  1:1145581/57‑1 (MQ=255)
                    gAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACg                                                                    >  1:1257643/1‑53 (MQ=255)
                    gAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAAAATACAGGAAt                                                   >  1:704879/1‑70 (MQ=255)
                            tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:378206/1‑55 (MQ=255)
                            tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:766786/1‑55 (MQ=255)
                            tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:685808/1‑55 (MQ=255)
                            tcCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACa                                                        >  1:2425285/1‑57 (MQ=255)
                              ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGAc                                        <  1:353126/71‑1 (MQ=255)
                              ccTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGAt                                           <  1:495015/68‑1 (MQ=255)
                               cTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGAc                                        >  1:1166524/1‑70 (MQ=255)
                                  ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          <  1:2476180/49‑1 (MQ=255)
                                  ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATg                                    <  1:2166764/71‑1 (MQ=255)
                                  ttCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTAt                                     <  1:2606981/70‑1 (MQ=255)
                                  ttCAGCAATCCAACCTCGTTTTGCCAGGCTCAAATTACGAGCAAACATa                                                          <  1:251615/49‑1 (MQ=255)
                                    cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTAc                                                                     >  1:1479217/1‑36 (MQ=255)
                                    cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGa                                   <  1:2734630/70‑1 (MQ=255)
                                    cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGa                                            >  1:108291/1‑61 (MQ=255)
                                    cAGCAATGCAACCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGa                                         >  1:109154/1‑64 (MQ=255)
                                     aGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGAc                                  <  1:2218755/70‑1 (MQ=255)
                                     aGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGAc                                  <  1:122161/70‑1 (MQ=255)
                                        aaTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATc                                              <  1:937362/55‑1 (MQ=255)
                                        aaTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                              >  1:454258/1‑71 (MQ=255)
                                        aaTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                              >  1:1868757/1‑71 (MQ=255)
                                        aaTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                              >  1:1510558/1‑71 (MQ=255)
                                        aaTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                              >  1:2610975/1‑71 (MQ=255)
                                         aTGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATc                                              <  1:2242762/54‑1 (MQ=255)
                                           gCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACt                            >  1:1432282/1‑70 (MQ=255)
                                           gCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACt                            >  1:1357232/1‑70 (MQ=255)
                                            cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGa                                   >  1:2108693/1‑62 (MQ=255)
                                            cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                           >  1:1832538/1‑70 (MQ=255)
                                            cAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                           >  1:1189906/1‑70 (MQ=255)
                                            cAACCTCGTGTTGCCAGGCTAAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTg                           >  1:1741922/1‑70 (MQ=255)
                                              aCCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAg                        <  1:1245842/71‑1 (MQ=255)
                                              aCCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAg                        <  1:1223861/71‑1 (MQ=255)
                                              aCCTCGTGTTGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGAc                                  <  1:294601/61‑1 (MQ=255)
                                               ccTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:2515518/1‑36 (MQ=255)
                                               ccTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:1570349/1‑36 (MQ=255)
                                               ccTCGTGTTGCCAGGCTCAAATTACGAGCAAACATa                                                          >  1:886774/1‑36 (MQ=255)
                                                 tCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGAtt                     >  1:946166/1‑71 (MQ=255)
                                                      ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                                <  1:295580/55‑1 (MQ=255)
                                                      ttGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCTaa                <  1:2792246/71‑1 (MQ=255)
                                                      ttGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGAc                                        <  1:386950/47‑1 (MQ=255)
                                                      ttGCCAGGCTCAAATCACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGa                 >  1:913991/1‑70 (MQ=255)
                                                        gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGa                              >  1:2126697/1‑55 (MQ=255)
                                                        gCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGaaaa              >  1:1615552/1‑71 (MQ=255)
                                                         ccAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAAc             <  1:1391285/71‑1 (MQ=255)
                                                               tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGa                                   >  1:2721279/1‑43 (MQ=255)
                                                               tCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgct       >  1:330661/1‑71 (MQ=255)
                                                                cAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgctg      >  1:2393398/1‑71 (MQ=255)
                                                                 aaaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTctgctg      <  1:2046175/70‑1 (MQ=255)
                                                                  aaTTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAActct          <  1:2164125/65‑1 (MQ=255)
                                                                    ttACGAGCAAACATACAGGAATAAATCGATGACTATGACaa                                <  1:93254/41‑1 (MQ=255)
                                                                    ttACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTgctggctg  >  1:2363070/1‑71 (MQ=255)
                                                                     |                                                                     
CATTGACTTCTGTAGAGGGTGAAGTCTCTCCCTATTCAGCAATGCAACCTCGTGTTGCCAGGCTCAAATTACGAGCAAACATACAGGAATAAATCGATGACTATGACAAGACTGAAGATTTCGAAAACTCTGCTGGCTG  >  W3110S.gb/4615982‑4616120

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: