Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,445 (C)7→8 intergenic (+90/‑265) rnlB → / → yfjP CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,4381.C83.0% 134.9 / 19.9 47intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (7/1);  new base C (16/23);  total (23/24)
Fisher's exact test for biased strand distribution p-value = 2.26e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGG  >  NC_000913/2767301‑2767571
                                                                                                                                           |                                                                                                                                         
aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑aaa                                                                                                                                           >  1:1208286/1‑138 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      >  2:378434/1‑139 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      <  1:704111/139‑1 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      >  2:517013/1‑139 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      >  2:190376/1‑139 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c‑t‑ccc                                                                                                                                >  2:747730/1‑139 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      <  2:672190/135‑1 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      >  1:672190/1‑135 (MQ=255)
          ttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c‑c‑cccctca                                                                                                                            <  2:891142/139‑5 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCTGTCAATCGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTcaca                                                                                                                           <  2:297854/139‑1 (MQ=255)
              ctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAAc                                                                                                                         <  1:247002/139‑1 (MQ=255)
                           aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTc                                                                                                                      >  1:1205806/1‑129 (MQ=255)
                           aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTc                                                                                                                      <  2:1205806/129‑1 (MQ=255)
                            cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGTATT‑AACC‑C‑C‑C‑CCCTcaca                                                                                                                           >  1:795199/1‑123 (MQ=255)
                            cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTcaca                                                                                                                           <  2:795199/123‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  1:378434/139‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  2:1230534/139‑1 (MQ=255)
                                     cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACTAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑CTC‑CCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  1:747730/139‑1 (MQ=255)
                                       gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACATCATGCATAGAGAACACGTTTTcc                                                                                                >  1:565254/1‑139 (MQ=255)
                                             ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCACACAACGACCTGCATAGGGAACACGTTTTCCCCtgag                                                                                          >  2:675944/1‑137 (MQ=255)
                                                 ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  1:460471/139‑1 (MQ=255)
                                                 ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:646889/139‑1 (MQ=255)
                                                           cttcttTCTAAACGAACCTGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c                                                                                                                                          <  2:240336/80‑1 (MQ=255)
                                                           cttcttTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c                                                                                                                                          >  1:240336/1‑80 (MQ=255)
                                                                 tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                      <  2:949371/139‑1 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTc                                                                                                                      >  1:442302/1‑88 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTc                                                                                                                      <  2:442302/88‑1 (MQ=255)
                                                                            cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAAc                                                           <  1:190376/139‑1 (MQ=255)
                                                                               tataGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTcaca                                                                                                                           >  1:379302/1‑72 (MQ=255)
                                                                               tataGAATTCGCCCTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTcaca                                                                                                                           <  2:379302/72‑1 (MQ=255)
                                                                                      ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  2:1214720/83‑1 (MQ=255)
                                                                                      ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAac                                                                                                         >  1:1214720/1‑83 (MQ=255)
                                                                                        cGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGa                                                                                                           <  1:70638/79‑1 (MQ=255)
                                                                                        cGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGa                                                                                                           >  2:70638/1‑79 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCTAATAGAATT‑AA‑CAA‑C‑C‑GATTCTCAACCGTCTGAATAGGGAACACCTTTTCCCCTGTCCACGTACGACCAATTCTCCCCCCCTATATACtaccccacaatac                                              >  1:828923/1‑127 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  >  1:366268/1‑87 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  2:366268/87‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt                                              >  2:1029891/1‑139 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTAACAACGCCCTTCATAGGGAAAACTTTGTCTCCTGAGCACCCCCAAGTAAATTTCTCCCCCTTTAACCGTATCAATAATGt                                              >  2:93582/1‑139 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑CAC‑C‑CCCTCACAACGTCATGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt                                              >  2:589476/1‑139 (MQ=255)
                                                                                                tGGCATGATATCAACTTGACCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCGGCATAGAGAACACGTTTTCCCCTGTGCACCCATGACTAAAATTCCCCCCCTTTACACTTCTCACTAATGTCActcat                                       <  2:455503/139‑5 (MQ=255)
                                                                                                  gCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      >  1:118642/1‑90 (MQ=255)
                                                                                                  gCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:118642/90‑1 (MQ=255)
                                                                                                            aaCTTGTCCAATATTATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGCCCTCCATAGCGAACACGTTTTCCCCTGTGTACCCACGACTAAATTTCCCCCTCTTCAAAATTCTCAATAATGTCACCAATTTTCTGGAAATc                           <  2:714722/139‑1 (MQ=255)
                                                                                                                tGTCCAACATGATCAAATAGCATT‑AA‑CTC‑C‑C‑CCATTACCAAGTCCTGCGTAGGGAACAGGTGTTCCCCTGTACACCCAAGACTAAATTTCCCTCCCTTTCACTTTCTCAATAAAGTCACCAATTTTGTATAAATGCTag                       >  1:755788/1‑138 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTc                                                                                                                      <  2:311174/42‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTc                                                                                                                      >  1:311174/1‑42 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AA‑C‑C‑CAC‑CCCTCACAACGTCCTGCATAGGGAAAACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCACCCTTTAAACATCTCAACAATGTCACCAAGTTTCTGGAAATCCTatga                     >  2:431788/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATT‑AA‑C‑C‑C‑CACCCTGACTACTTCCTTCATAGTGAACACCATTACACCTGTCAACCCACGAGTAACTTTCACTCCCCTTAAAATACACACTAATCTCTCCAATTTTCTGGTAATCCTaaca                     >  2:449428/1‑136 (MQ=255)
                                                                                                                          gATCAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             <  1:325958/139‑1 (MQ=255)
                                                                                                                          gATCAAATACCCTTGAC‑C‑C‑C‑C‑CCCTAGCAACGTCCTGCATAGGCAACACGTTTTCCCCTGTGCACCCACGACTAAATAACCCCCCCTTTATACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             <  2:237933/139‑1 (MQ=255)
                                                                                                                             cAAATAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg          <  1:431788/139‑1 (MQ=255)
                                                                                                                               aaTAGCATT‑AA‑C‑C‑CAC‑CCCCCACAACGTCCTGCATAGGGAACACGTTTTCCCCTCTGGACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATTCCACCAATTTTATGGAAAACCTATTATACCTTTCCCgat        >  1:60518/1‑137 (MQ=255)
                                                                                                                                 tAGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaa      <  1:93582/139‑1 (MQ=255)
                                                                                                                                  aGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             >  2:231999/1‑131 (MQ=255)
                                                                                                                                  aGCATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             <  1:231999/131‑1 (MQ=255)
                                                                                                                                    cATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGCAAATCCTATGATGCCTTTcc            >  2:113603/1‑130 (MQ=255)
                                                                                                                                    cATT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc            <  1:113603/130‑1 (MQ=255)
                                                                                                                                     aTT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGCCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg  >  1:274625/1‑139 (MQ=255)
                                                                                                                                     aTT‑AACC‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCACCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg  >  2:811985/1‑139 (MQ=255)
                                                                                                                                     aTT‑AA‑C‑C‑CAC‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACt                                                          >  1:749557/1‑83 (MQ=255)
                                                                                                                                           |                                                                                                                                         
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑C‑C‑CCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGG  >  NC_000913/2767301‑2767571

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑