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

TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCG  >  NC_000913/2767305‑2767563
                                                                                                                                        |                                                                                                                             
tctGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCCAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccc                                                                                                                          <  1:442069/139‑1 (MQ=255)
        gTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                   <  1:419322/139‑1 (MQ=255)
             ccGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              <  1:726143/139‑1 (MQ=255)
              cGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTcc                                                                                                             <  2:687740/139‑1 (MQ=255)
                           cAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCGGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                   <  2:916503/120‑1 (MQ=255)
                           cAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                   >  1:916503/1‑120 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              >  2:949315/1‑117 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              <  1:1270211/117‑1 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              >  2:1270211/1‑117 (MQ=255)
                                   gTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              <  1:949315/117‑1 (MQ=255)
                                        tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                   <  1:1095714/107‑1 (MQ=255)
                                        tGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                   >  2:1095714/1‑107 (MQ=255)
                                                             tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  1:1193997/114‑1 (MQ=255)
                                                             tctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  2:1193997/1‑114 (MQ=255)
                                                             gtataagagacaggTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTccc                                                              <  1:1009376/126‑1 (MQ=255)
                                                                taagagacagctataGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCACTCCCCCCCCCCTCACAACGTCCTGCCTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCcccccc                                                           <  2:144029/128‑1 (MQ=255)
                                                                     aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              <  1:449674/83‑1 (MQ=255)
                                                                     aaCCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              >  2:449674/1‑83 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccccc                                                                                                                        <  2:895994/61‑1 (MQ=255)
                                                                                aaTTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAccccccc                                                                                                                        >  1:895994/1‑61 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAACCCCAACTCGCAACGTCCTGAAAAGGGAATACGATCTTCACTGTGGGACCACAACTAAATTTCTCCCCCTTTGAACtttctggccctta                                      >  1:144029/1‑128 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCTTTCATAGGGAACACGATTTCCCCTGTGCACCAACGACTAAATTTCTCCCCCTTTAAACTACTCaataattt                                      >  1:391833/1‑137 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:316999/90‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:316999/1‑90 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:1278071/90‑1 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:1278071/1‑90 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCTCATGTGAATGCACGACTAAATTTCCCCCTCACTAAAcctgtctgttatac                                      >  2:1009376/1‑126 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACt                                                                       >  2:500582/1‑106 (MQ=255)
                                                                                     gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACt                                                                       <  1:500582/106‑1 (MQ=255)
                                                                                                             gTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg              >  2:241281/1‑139 (MQ=255)
                                                                                                             gTCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTACCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatg              >  2:911981/1‑139 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTcaca                                                                                                                   >  1:809231/1‑37 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              >  2:741148/1‑42 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTc                                                                                                              <  1:741148/42‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  2:120742/1‑65 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:1000327/65‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  2:411589/65‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:411589/1‑65 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       <  1:120742/65‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTccc                                                                                       >  1:1000327/1‑65 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCt                                                         <  1:487840/95‑1 (MQ=255)
                                                                                                              tCCAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACAAGTTTTCCCCTCTGCACCCACGACTAAATTTCCCCCCCt                                                         >  2:487840/1‑95 (MQ=255)
                                                                                                                cAATATGATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgatg           <  1:454525/139‑1 (MQ=255)
                                                                                                                      gATCAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTc     <  1:1343333/139‑1 (MQ=255)
                                                                                                                         cAAATAGCATTAACCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCg  <  1:911981/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                             
TCTGTTTTGTCTCCCGATGATGTACGCCAAGTCCTGTCAATGGATTTTTGAAACTCTTCTTTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCG  >  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: