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.C95.2% 136.3 / ‑0.2 42intergenic (+83/‑272)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base C (19/21);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCG  >  NC_000913/2767305‑2767563
                                                                                                                                        |                                                                                                                             
tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCCAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AA‑ccccc                                                                                                                          <  1:442069/139‑1 (MQ=255)
        gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTcaca                                                                                                                   <  1:419322/139‑1 (MQ=255)
             ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              <  1:726143/139‑1 (MQ=255)
              cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTcc                                                                                                             <  2:687740/139‑1 (MQ=255)
                           cAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCGGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTcaca                                                                                                                   <  2:916503/120‑1 (MQ=255)
                           cAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTcaca                                                                                                                   >  1:916503/1‑120 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              >  2:949315/1‑117 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              <  1:1270211/117‑1 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              >  2:1270211/1‑117 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              <  1:949315/117‑1 (MQ=255)
                                        tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTcaca                                                                                                                   <  1:1095714/107‑1 (MQ=255)
                                        tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTcaca                                                                                                                   >  2:1095714/1‑107 (MQ=255)
                                                             tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  1:1193997/114‑1 (MQ=255)
                                                             tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  2:1193997/1‑114 (MQ=255)
                                                             gtataagagacaggTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                              <  1:1009376/126‑1 (MQ=255)
                                                                taagagacagctataGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCACTC‑CC‑CCCCCCCTCACAACGTCCTGCCTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCcccccc                                                           <  2:144029/128‑1 (MQ=255)
                                                                     aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              <  1:449674/83‑1 (MQ=255)
                                                                     aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              >  2:449674/1‑83 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AA‑ccccccc                                                                                                                        <  2:895994/61‑1 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AA‑ccccccc                                                                                                                        >  1:895994/1‑61 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TTAAA‑CCCCAACTCGCAACGTCCTGAAAAGGGAATACGATCTTCACTGTGGGACCACAACTAAATTTCTCCCCCTTTGAACtttctggccctta                                      >  1:144029/1‑128 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCTTTCATAGGGAACACGATTTCCCCTGTGCACCAACGACTAAATTTCTCCCCCTTTAAACTACTCaataattt                                      >  1:391833/1‑137 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:316999/90‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:316999/1‑90 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:1278071/90‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:1278071/1‑90 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCTCATGTGAATGCACGACTAAATTTCCCCCTCACTAAAcctgtctgttatac                                      >  2:1009376/1‑126 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACt                                                                       >  2:500582/1‑106 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACt                                                                       <  1:500582/106‑1 (MQ=255)
                                                                                                             gTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg              >  2:241281/1‑139 (MQ=255)
                                                                                                             gTCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTACCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg              >  2:911981/1‑139 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTcaca                                                                                                                   >  1:809231/1‑37 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              >  2:741148/1‑42 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTc                                                                                                              <  1:741148/42‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  2:120742/1‑65 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:1000327/65‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:411589/65‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:411589/1‑65 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  1:120742/65‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:1000327/1‑65 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCt                                                         <  1:487840/95‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACAAGTTTTCCCCTCTGCACCCACGACTAAATTTCCCCCCCt                                                         >  2:487840/1‑95 (MQ=255)
                                                                                                                cAATATGATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg           <  1:454525/139‑1 (MQ=255)
                                                                                                                      gATCAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc     <  1:1343333/139‑1 (MQ=255)
                                                                                                                         cAAATAGCA‑TT‑AACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg  <  1:911981/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                             
TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCA‑TT‑AA‑CCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCG  >  NC_000913/2767305‑2767563

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