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

CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTA  >  NC_000913/2767303‑2767566
                                                                                                                                         |                                                                                                                                         
ccTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccc                                                                                                                                      >  2:775886/1‑139 (MQ=255)
 ctctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccc                                                                                                                                     <  2:427019/139‑1 (MQ=255)
  tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTaaaatca                                                                                                                                       >  1:132297/1‑134 (MQ=255)
  tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                                    >  1:802569/1‑139 (MQ=255)
   ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccc                                                                                                                                  >  2:39511/1‑139 (MQ=255)
   ctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATAGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccc                                                                                                                                  <  2:666021/139‑1 (MQ=255)
        ttGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAcccccccctca                                                                                                                             <  1:220676/139‑5 (MQ=255)
          gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                           <  2:802569/139‑1 (MQ=255)
           tctcCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACaa                                                                                                                          >  2:17165/1‑139 (MQ=255)
               ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGACCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      <  2:719534/139‑1 (MQ=255)
                cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACATCGTcc                                                                                                                     >  1:656076/1‑139 (MQ=255)
                 gatgatGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCt                                                                                                                    <  1:937329/139‑1 (MQ=255)
                          cGCCAAGTCCTGTCAATGTATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGa                                                                                                           <  1:538359/139‑1 (MQ=255)
                          cGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGa                                                                                                           <  2:61556/139‑1 (MQ=255)
                            ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                         <  1:41594/139‑1 (MQ=255)
                            ccAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAac                                                                                                         >  1:995577/1‑139 (MQ=255)
                                   cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  2:995577/139‑1 (MQ=255)
                                   cTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGtttt                                                                                                  <  1:39511/139‑1 (MQ=255)
                                      tCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               <  2:7923/139‑1 (MQ=255)
                                               ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:656076/139‑1 (MQ=255)
                                               ttttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCAcc                                                                                      <  2:958839/139‑1 (MQ=255)
                                                tttGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCAACGGAAAACGTAATGCATAGGAAACACGTATTctctgtttcactc                                                                                     >  1:1049502/1‑127 (MQ=255)
                                                     aaCTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGAc                                                                                <  1:775886/139‑1 (MQ=255)
                                                                ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccc                                                                  <  1:720100/139‑1 (MQ=255)
                                                                ctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGAATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTcccc                                                                  <  2:670825/139‑1 (MQ=255)
                                                                        aCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCcttta                                                        <  1:17165/139‑3 (MQ=255)
                                                                                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      >  1:223201/1‑70 (MQ=255)
                                                                                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                                      <  2:223201/70‑1 (MQ=255)
                                                                                    ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCAACCCCTCCCAACGTTCTGTCTAGGGAACAAGGTCTCCCCTGTGCTCCAACGACTACATTTCCCCCCACTTTAACCTGCTCaata                                             >  1:1022695/1‑139 (MQ=255)
                                                                                       gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               >  2:138880/1‑90 (MQ=255)
                                                                                       gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               <  1:138880/90‑1 (MQ=255)
                                                                                                       atTAACTTGTCCAATTTGATCAAATAGCATTAAACCCGCCCTCACAACGACCTGCATAGGAAACACGTTTTACACTATGCACCAACGACTAAATTTCTCCCCCCTTTAAACTTCTCCATAATGTCACCAATTTTCTCGa                          >  2:706418/1‑139 (MQ=255)
                                                                                                               gTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACTATTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTat                  >  2:423356/1‑139 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccccccccttta                                                         <  1:975251/99‑6 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCTCCCCCCTTTa                                                         >  2:975251/1‑99 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACACACGACTAAATTTCCCCCCCATTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCTTatg                 >  2:209109/1‑139 (MQ=255)
                                                                                                                tCCAATATGATCAAATAGCATTAACAACACGCTCACACCCTCCTCAATAACAAACCCATTTCCCCCTCTGCACCCACCACTAACTATACCAACACTGTAAACTTCTCAAAAATGTCAACAACCTTATGTAAACCTTatg                 >  1:411289/1‑139 (MQ=255)
                                                                                                                accaagctgatgaAGTAGCTTTAACCCCCCCCTCACATCGTCCTGCATAGGGAACACGTTTTCCCCGTTGCACCCACGTCTAAGTTTCCCCCCCCTTTACACCTCTCAATAATGTCTCCAATTTTCTGGAAATCCTatg                 <  2:621399/132‑1 (MQ=255)
                                                                                                                         aTCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc        <  2:970999/139‑1 (MQ=255)
                                                                                                                          tCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               <  1:840319/55‑1 (MQ=255)
                                                                                                                          tCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                               >  2:840319/1‑55 (MQ=255)
                                                                                                                          tCAAATAGCATTAACCCCCCCGTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTAAAACTTCCCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc       >  2:515137/1‑139 (MQ=255)
                                                                                                                             aaTAGCATTAATCTCCCCCTCACAACGTCCGCCCCAGGGAACACGTTTTCCCCCGCTCAGCCACGACTAAATTTCCCCCCACTTCAGTCTTCCTAGTCGCGTCTCCAATTTTCTGTTATACCCATGATGCCTCTccctg    >  1:391244/1‑137 (MQ=255)
                                                                                                                              aTAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa  >  2:668324/1‑139 (MQ=255)
                                                                                                                               tAGCATTAACCCCCCCCTCACACCGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg     >  2:392623/1‑136 (MQ=255)
                                                                                                                               tAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACCCGTTTTCCCCTGTGCACCCACGACTAAATTTACCCCCACTATAAACTTCTCAATAATGCCACCAATTTTCTGGAAATCCTATGTTTCCTTTACCGGTa  >  2:482860/1‑139 (MQ=255)
                                                                                                                               tAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg     <  1:392623/136‑1 (MQ=255)
                                                                                                                               tAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTa  >  2:204340/1‑139 (MQ=255)
                                                                                                                                aGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACtt                                                    <  1:323194/88‑1 (MQ=255)
                                                                                                                                aGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCCTTTAAACtt                                                    >  2:323194/1‑88 (MQ=255)
                                                                                                                                   aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccccccccttta                                                         <  1:841527/80‑6 (MQ=255)
                                                                                                                                   aTTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCcttta                                                        >  2:841527/1‑78 (MQ=255)
                                                                                                                                         |                                                                                                                                         
CCTCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTA  >  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: