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

AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGG  >  NC_000913/2767301‑2767571
                                                                                                                                           |                                                                                                                                         
aGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTaaa                                                                                                                                           >  1:1208286/1‑138 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      >  2:378434/1‑139 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      <  1:704111/139‑1 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      >  2:517013/1‑139 (MQ=255)
  ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      >  2:190376/1‑139 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccctccc                                                                                                                                >  2:747730/1‑139 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      <  2:672190/135‑1 (MQ=255)
      tGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      >  1:672190/1‑135 (MQ=255)
          ttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccccctca                                                                                                                            <  2:891142/139‑5 (MQ=255)
            gTCTCCCGATGATGTACGCCAAGTCCTGTCAATCGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                           <  2:297854/139‑1 (MQ=255)
              ctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAAc                                                                                                                         <  1:247002/139‑1 (MQ=255)
                           aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      >  1:1205806/1‑129 (MQ=255)
                           aCGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      <  2:1205806/129‑1 (MQ=255)
                            cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGTATTAACCCCCCCCTcaca                                                                                                                           >  1:795199/1‑123 (MQ=255)
                            cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                           <  2:795199/123‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  1:378434/139‑1 (MQ=255)
                              ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  2:1230534/139‑1 (MQ=255)
                                     cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACTAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCTCCCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  1:747730/139‑1 (MQ=255)
                                       gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACATCATGCATAGAGAACACGTTTTcc                                                                                                >  1:565254/1‑139 (MQ=255)
                                             ggATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCACACAACGACCTGCATAGGGAACACGTTTTCCCCtgag                                                                                          >  2:675944/1‑137 (MQ=255)
                                                 ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  1:460471/139‑1 (MQ=255)
                                                 ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:646889/139‑1 (MQ=255)
                                                           cttcttTCTAAACGAACCTGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAc                                                                                                                                          <  2:240336/80‑1 (MQ=255)
                                                           cttcttTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAc                                                                                                                                          >  1:240336/1‑80 (MQ=255)
                                                                 tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                                      <  2:949371/139‑1 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      >  1:442302/1‑88 (MQ=255)
                                                                    aaaCGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      <  2:442302/88‑1 (MQ=255)
                                                                            cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAAc                                                           <  1:190376/139‑1 (MQ=255)
                                                                               tataGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                           >  1:379302/1‑72 (MQ=255)
                                                                               tataGAATTCGCCCTTTTGGCATGATATTAACCTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                           <  2:379302/72‑1 (MQ=255)
                                                                                      ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  2:1214720/83‑1 (MQ=255)
                                                                                      ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                         >  1:1214720/1‑83 (MQ=255)
                                                                                        cGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGa                                                                                                           <  1:70638/79‑1 (MQ=255)
                                                                                        cGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGa                                                                                                           >  2:70638/1‑79 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCTAATAGAATTAACAACCGATTCTCAACCGTCTGAATAGGGAACACCTTTTCCCCTGTCCACGTACGACCAATTCTCCCCCCCTATATACtaccccacaatac                                              >  1:828923/1‑127 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  >  1:366268/1‑87 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  2:366268/87‑1 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt                                              >  2:1029891/1‑139 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTAACAACGCCCTTCATAGGGAAAACTTTGTCTCCTGAGCACCCCCAAGTAAATTTCTCCCCCTTTAACCGTATCAATAATGt                                              >  2:93582/1‑139 (MQ=255)
                                                                                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCACCCCCTCACAACGTCATGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGt                                              >  2:589476/1‑139 (MQ=255)
                                                                                                tGGCATGATATCAACTTGACCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCGGCATAGAGAACACGTTTTCCCCTGTGCACCCATGACTAAAATTCCCCCCCTTTACACTTCTCACTAATGTCActcat                                       <  2:455503/139‑5 (MQ=255)
                                                                                                  gCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      >  1:118642/1‑90 (MQ=255)
                                                                                                  gCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:118642/90‑1 (MQ=255)
                                                                                                            aaCTTGTCCAATATTATCAAATAGCATTAACCCCCCCCTCACAACGCCCTCCATAGCGAACACGTTTTCCCCTGTGTACCCACGACTAAATTTCCCCCTCTTCAAAATTCTCAATAATGTCACCAATTTTCTGGAAATc                           <  2:714722/139‑1 (MQ=255)
                                                                                                                tGTCCAACATGATCAAATAGCATTAACTCCCCCATTACCAAGTCCTGCGTAGGGAACAGGTGTTCCCCTGTACACCCAAGACTAAATTTCCCTCCCTTTCACTTTCTCAATAAAGTCACCAATTTTGTATAAATGCTag                       >  1:755788/1‑138 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      <  2:311174/42‑1 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      >  1:311174/1‑42 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCACCCCTCACAACGTCCTGCATAGGGAAAACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCACCCTTTAAACATCTCAACAATGTCACCAAGTTTCTGGAAATCCTatga                     >  2:431788/1‑139 (MQ=255)
                                                                                                                  tCCAATATGATCAAATAGCATTAACCCCACCCTGACTACTTCCTTCATAGTGAACACCATTACACCTGTCAACCCACGAGTAACTTTCACTCCCCTTAAAATACACACTAATCTCTCCAATTTTCTGGTAATCCTaaca                     >  2:449428/1‑136 (MQ=255)
                                                                                                                          gATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             <  1:325958/139‑1 (MQ=255)
                                                                                                                          gATCAAATACCCTTGACCCCCCCCTAGCAACGTCCTGCATAGGCAACACGTTTTCCCCTGTGCACCCACGACTAAATAACCCCCCCTTTATACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             <  2:237933/139‑1 (MQ=255)
                                                                                                                             cAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg          <  1:431788/139‑1 (MQ=255)
                                                                                                                               aaTAGCATTAACCCACCCCCCACAACGTCCTGCATAGGGAACACGTTTTCCCCTCTGGACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATTCCACCAATTTTATGGAAAACCTATTATACCTTTCCCgat        >  1:60518/1‑137 (MQ=255)
                                                                                                                                 tAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTaa      <  1:93582/139‑1 (MQ=255)
                                                                                                                                  aGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             >  2:231999/1‑131 (MQ=255)
                                                                                                                                  aGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc             <  1:231999/131‑1 (MQ=255)
                                                                                                                                    cATTAACCCCCCCCTCACAACGTCCTGCATATGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGCAAATCCTATGATGCCTTTcc            >  2:113603/1‑130 (MQ=255)
                                                                                                                                    cATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc            <  1:113603/130‑1 (MQ=255)
                                                                                                                                     aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGCCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg  >  1:274625/1‑139 (MQ=255)
                                                                                                                                     aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCACCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg  >  2:811985/1‑139 (MQ=255)
                                                                                                                                     aTTAACCCACCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACt                                                          >  1:749557/1‑83 (MQ=255)
                                                                                                                                           |                                                                                                                                         
AGCCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGG  >  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: