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.C88.1% 125.1 / 8.5 43intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (5/0);  new base C (13/24);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 1.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTA  >  NC_000913/2767303‑2767566
                                                                                                                                         |                                                                                                                                         
ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑c                                                                                                                                      >  2:775886/1‑139 (MQ=255)
 ctctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑cc                                                                                                                                     <  2:427019/139‑1 (MQ=255)
  tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑aaaatca                                                                                                                                       >  1:132297/1‑134 (MQ=255)
  tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑ccc                                                                                                                                    >  1:802569/1‑139 (MQ=255)
   ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑ccc‑c                                                                                                                                  >  2:39511/1‑139 (MQ=255)
   ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATAGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑ccc‑c                                                                                                                                  <  2:666021/139‑1 (MQ=255)
        ttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑c‑c‑ccc‑ccctca                                                                                                                             <  1:220676/139‑5 (MQ=255)
          gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑Tcaca                                                                                                                           <  2:802569/139‑1 (MQ=255)
           tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACaa                                                                                                                          >  2:17165/1‑139 (MQ=255)
               ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGACCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTc                                                                                                                      <  2:719534/139‑1 (MQ=255)
                cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACATCGTcc                                                                                                                     >  1:656076/1‑139 (MQ=255)
                 gatgatGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCt                                                                                                                    <  1:937329/139‑1 (MQ=255)
                          cGCCAAGTCCTGTCAATGTATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGa                                                                                                           <  1:538359/139‑1 (MQ=255)
                          cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGa                                                                                                           <  2:61556/139‑1 (MQ=255)
                            ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAac                                                                                                         <  1:41594/139‑1 (MQ=255)
                            ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAac                                                                                                         >  1:995577/1‑139 (MQ=255)
                                   cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  2:995577/139‑1 (MQ=255)
                                   cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  1:39511/139‑1 (MQ=255)
                                      tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               <  2:7923/139‑1 (MQ=255)
                                               ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:656076/139‑1 (MQ=255)
                                               ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:958839/139‑1 (MQ=255)
                                                tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑CCCAAC‑GGAAAACGTAATGCATAGGAAACACGTATTctctgtttcactc                                                                                     >  1:1049502/1‑127 (MQ=255)
                                                     aaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc                                                                                <  1:775886/139‑1 (MQ=255)
                                                                ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑ccc                                                                  <  1:720100/139‑1 (MQ=255)
                                                                ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGAATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑ccc                                                                  <  2:670825/139‑1 (MQ=255)
                                                                        aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CCCCC‑C‑cttta                                                        <  1:17165/139‑3 (MQ=255)
                                                                                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTc                                                                                                                      >  1:223201/1‑70 (MQ=255)
                                                                                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTc                                                                                                                      <  2:223201/70‑1 (MQ=255)
                                                                                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑CAACC‑CC‑TCCCAACGTTCTGTCTAGGGAACAAGGTCTCCCCTGTGCTCCAACGACTACAT‑TT‑C‑CCCCCAC‑TTTAACCTGCTCaata                                             >  1:1022695/1‑139 (MQ=255)
                                                                                       gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               >  2:138880/1‑90 (MQ=255)
                                                                                       gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               <  1:138880/90‑1 (MQ=255)
                                                                                                       atTAACTTGTCCAATTTGATCAAATAGCATTAAA‑C‑C‑CGC‑CC‑TCACAACGACCTGCATAGGAAACACGTTTTACACTATGCACCAACGACTAAAT‑TT‑CTCCCCC‑C‑TTTAAACTTCTCCATAATGTCACCAATTTTCTCGa                          >  2:706418/1‑139 (MQ=255)
                                                                                                               gTCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACTATTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                  >  2:423356/1‑139 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑ccccc‑ccttta                                                         <  1:975251/99‑6 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑CTCCCCC‑C‑TTTa                                                         >  2:975251/1‑99 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACACACGACTAAAT‑TT‑C‑CCCCC‑CATTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCTTatg                 >  2:209109/1‑139 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATT‑AA‑CAA‑CAC‑GC‑TCACACCCTCCTCAATAACAAACCCATTTCCCCCTCTGCACCCACCACTAACTATA‑C‑CAACA‑C‑TGTAAACTTCTCAAAAATGTCAACAACCTTATGTAAACCTTatg                 >  1:411289/1‑139 (MQ=255)
                                                                                                                accaagctgatgaAGTAGCTTT‑AACC‑C‑CCC‑CC‑TCACATCGTCCTGCATAGGGAACACGTTTTCCCCGTTGCACCCACGTCTAAGT‑TTCC‑CCCCC‑C‑TTTACACCTCTCAATAATGTCTCCAATTTTCTGGAAATCCTatg                 <  2:621399/132‑1 (MQ=255)
                                                                                                                         aTCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc        <  2:970999/139‑1 (MQ=255)
                                                                                                                          tCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               <  1:840319/55‑1 (MQ=255)
                                                                                                                          tCAAATAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               >  2:840319/1‑55 (MQ=255)
                                                                                                                          tCAAATAGCATT‑AA‑C‑C‑CCC‑CCGTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTAAAACTTCCCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc       >  2:515137/1‑139 (MQ=255)
                                                                                                                             aaTAGCATT‑AATC‑T‑CCC‑CC‑TCACAACGTCCGCCCCAGGGAACACGTTTTCCCCCGCTCAGCCACGACTAAAT‑TT‑C‑CCCCCAC‑TTCAGTCTTCCTAGTCGCGTCTCCAATTTTCTGTTATACCCATGATGCCTCTccctg    >  1:391244/1‑137 (MQ=255)
                                                                                                                              aTAGCATT‑AA‑C‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa  >  2:668324/1‑139 (MQ=255)
                                                                                                                               tAGCATT‑AACC‑C‑CCC‑CC‑TCACACCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg     >  2:392623/1‑136 (MQ=255)
                                                                                                                               tAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAAT‑TTAC‑CCCCA‑C‑TATAAACTTCTCAATAATGCCACCAATTTTCTGGAAATCCTATGTTTCCTTTACCGGTa  >  2:482860/1‑139 (MQ=255)
                                                                                                                               tAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg     <  1:392623/136‑1 (MQ=255)
                                                                                                                               tAGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa  >  2:204340/1‑139 (MQ=255)
                                                                                                                                aGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACtt                                                    <  1:323194/88‑1 (MQ=255)
                                                                                                                                aGCATT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TTCC‑CCCCC‑C‑TTTAAACtt                                                    >  2:323194/1‑88 (MQ=255)
                                                                                                                                   aTT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑c‑ccccc‑ccttta                                                         <  1:841527/80‑6 (MQ=255)
                                                                                                                                   aTT‑AACC‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CCCCC‑C‑cttta                                                        >  2:841527/1‑78 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATT‑AA‑C‑C‑CCC‑CC‑TCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAAT‑TT‑C‑CCCCC‑C‑TTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTA  >  NC_000913/2767303‑2767566

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: ‑