Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,507 (C)7→8 intergenic (+152/‑203) 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,5001.C81.1% 136.1 / 27.1 54intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (6/4);  new base C (22/21);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATA  >  NC_000913/2767367‑2767631
                                                                                                                                            |                                                                                                                                       
ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑cc‑c                                                                                                                                   <  1:720100/139‑1 (MQ=255)
ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGAATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑cc‑c                                                                                                                                   <  2:670825/139‑1 (MQ=255)
        aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CC‑CCC‑C‑cttta                                                                                                                         <  1:17165/139‑3 (MQ=255)
                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑CAACCCC‑TCCCAACGTTCTGTCTAGGGAACAAGGTCTCCCCTGTGCTCCAACGACTACAT‑TT‑C‑CC‑CCCAC‑TTTAACCTGCTCaata                                                                                                              >  1:1022695/1‑139 (MQ=255)
                                       atTAACTTGTCCAATTTGATCAAATAGCATTAAA‑C‑C‑CGCCC‑TCACAACGACCTGCATAGGAAACACGTTTTACACTATGCACCAACGACTAAAT‑TT‑CTCC‑CCC‑C‑TTTAAACTTCTCCATAATGTCACCAATTTTCTCGa                                                                                           >  2:706418/1‑139 (MQ=255)
                                               gTCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACTATTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                                                                                   >  2:423356/1‑139 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑cc‑ccc‑ccttta                                                                                                                          <  1:975251/99‑6 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑CTCC‑CCC‑C‑TTTa                                                                                                                          >  2:975251/1‑99 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACACACGACTAAAT‑TT‑C‑CC‑CCC‑CATTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCTTatg                                                                                  >  2:209109/1‑139 (MQ=255)
                                                tCCAATATGATCAAATAGCATT‑AA‑CAA‑CACGC‑TCACACCCTCCTCAATAACAAACCCATTTCCCCCTCTGCACCCACCACTAACTATA‑C‑CA‑ACA‑C‑TGTAAACTTCTCAAAAATGTCAACAACCTTATGTAAACCTTatg                                                                                  >  1:411289/1‑139 (MQ=255)
                                                accaagctgatgaAGTAGCTTT‑AACC‑C‑CCCCC‑TCACATCGTCCTGCATAGGGAACACGTTTTCCCCGTTGCACCCACGTCTAAGT‑TTCC‑CC‑CCC‑C‑TTTACACCTCTCAATAATGTCTCCAATTTTCTGGAAATCCTatg                                                                                  <  2:621399/132‑1 (MQ=255)
                                                         aTCAAATAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                         <  2:970999/139‑1 (MQ=255)
                                                          tCAAATAGCATT‑AA‑C‑C‑CCCCCGTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTAAAACTTCCCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                                                        >  2:515137/1‑139 (MQ=255)
                                                             aaTAGCATT‑AATC‑T‑CCCCC‑TCACAACGTCCGCCCCAGGGAACACGTTTTCCCCCGCTCAGCCACGACTAAAT‑TT‑C‑CC‑CCCAC‑TTCAGTCTTCCTAGTCGCGTCTCCAATTTTCTGTTATACCCATGATGCCTCTccctg                                                                     >  1:391244/1‑137 (MQ=255)
                                                              aTAGCATT‑AA‑C‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                                                                   >  2:668324/1‑139 (MQ=255)
                                                               tAGCATT‑AACC‑C‑CCCCC‑TCACACCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                                                      >  2:392623/1‑136 (MQ=255)
                                                               tAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAAT‑TTAC‑CC‑CCA‑C‑TATAAACTTCTCAATAATGCCACCAATTTTCTGGAAATCCTATGTTTCCTTTACCGGTa                                                                   >  2:482860/1‑139 (MQ=255)
                                                               tAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg                                                                      <  1:392623/136‑1 (MQ=255)
                                                               tAGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa                                                                   >  2:204340/1‑139 (MQ=255)
                                                                aGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACtt                                                                                                                     <  1:323194/88‑1 (MQ=255)
                                                                aGCATT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACtt                                                                                                                     >  2:323194/1‑88 (MQ=255)
                                                                   aTT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑cc‑ccc‑ccttta                                                                                                                          <  1:841527/80‑6 (MQ=255)
                                                                   aTT‑AACC‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CC‑CCC‑C‑cttta                                                                                                                         >  2:841527/1‑78 (MQ=255)
                                                                    attaccc‑c‑ccccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATg                                                               <  2:132297/133‑1 (MQ=255)
                                                                     ttaacc‑c‑ccccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTAACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                                                              >  1:540801/6‑139 (MQ=255)
                                                                            c‑ccccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCTCCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                        <  2:1022695/139‑1 (MQ=255)
                                                                              ccccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTTTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                                                       <  1:558881/139‑1 (MQ=255)
                                                                               cccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGCTGCCTTTCCCGGTAAATGGt                                                             >  2:53354/1‑132 (MQ=255)
                                                                               cccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGt                                                             <  1:53354/132‑1 (MQ=255)
                                                                               cccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGa                                                      >  2:305794/1‑139 (MQ=255)
                                                                               cccc‑TCACAAAGTCCTGCCTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCCCAATAATGTCAACAATTTTCTGGAAATCCTATGATCCCTTTCCCGGTAAATGGTGCTGTGa                                                      >  1:549153/1‑139 (MQ=255)
                                                                                ccc‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                     <  1:482860/139‑1 (MQ=255)
                                                                                         aCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGaaa                                             <  2:540801/139‑1 (MQ=255)
                                                                                            tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAt                                          >  1:841582/1‑139 (MQ=255)
                                                                                               tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCt                                       <  1:265534/139‑1 (MQ=255)
                                                                                                gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                     <  1:220455/124‑1 (MQ=255)
                                                                                                gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                     >  2:220455/1‑124 (MQ=255)
                                                                                                gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                      <  2:440827/139‑1 (MQ=255)
                                                                                                 cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                     >  1:902192/1‑139 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                                                                >  1:110429/1‑107 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                                                                <  2:110429/107‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                                                       <  2:182549/116‑1 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgtg                                                       >  1:182549/1‑116 (MQ=255)
                                                                                                      ggAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCt                                >  1:306223/1‑139 (MQ=255)
                                                                                                        aaCACGTTTTCCCCTGTGCACCCACGACTAAATCTC‑C‑CC‑CCC‑C‑TGTAAACTTCTCTATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTAt                              <  2:570018/139‑1 (MQ=255)
                                                                                                           acGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTAtgtg                           >  2:180028/1‑139 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGt                                                             <  1:663862/100‑1 (MQ=255)
                                                                                                                ttCCCCTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGt                                                             >  2:663862/1‑100 (MQ=255)
                                                                                                                  ccccTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                         <  2:971386/86‑1 (MQ=255)
                                                                                                                  ccccTGTGCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc                                                                         >  1:971386/1‑86 (MQ=255)
                                                                                                                       gtgCACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAAt               >  1:970616/1‑139 (MQ=255)
                                                                                                                          cACCCACGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCtt            <  1:180028/139‑1 (MQ=255)
                                                                                                                               aCGACTAAATGTC‑C‑CC‑CCC‑C‑TTTATACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat       <  2:156011/139‑1 (MQ=255)
                                                                                                                               aCGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat       <  1:305794/139‑1 (MQ=255)
                                                                                                                               aCGACTAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat       <  2:902192/139‑1 (MQ=255)
                                                                                                                               aCGACTAAAT‑TT‑C‑CCGCCC‑C‑TTTAAACTTCTCAATAATGTAACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTTAGAAACATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat       <  1:706418/139‑1 (MQ=255)
                                                                                                                                    tAAAT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAgatagata  >  1:824966/1‑139 (MQ=255)
                                                                                                                                      aaT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCag         >  2:741273/1‑130 (MQ=255)
                                                                                                                                      aaT‑TTCC‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCag         <  1:741273/130‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                       
CTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑CCCCC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CC‑CCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATA  >  NC_000913/2767367‑2767631

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