Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,159 G→A *163* (TAG→TAA
S51S (AGC→AGT
sieB →
ydaF ←
phage superinfection exclusion protein, Rac prophage
uncharacterized protein, Rac prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,1590GA100.0% 79.3 / NA 27*163* (TAG→TAA
S51S (AGC→AGT
sieB
ydaF
phage superinfection exclusion protein, Rac prophage
uncharacterized protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/13);  total (14/13)

TAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATC  >  NC_000913/1419032‑1419251
                                                                                                                               |                                                                                            
tAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                   <  1:461944/139‑1 (MQ=255)
                            gagTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGtt                                                       >  1:100241/1‑139 (MQ=255)
                             agTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                      >  1:186986/1‑139 (MQ=255)
                                 aTAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGt                                                  >  2:472165/1‑139 (MQ=255)
                                     cgcAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTa                                              >  1:727936/1‑139 (MQ=255)
                                         aaCTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGATTCGGTTTTAAATg                                          >  1:128806/1‑139 (MQ=255)
                                              ttCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTg                                     >  2:489263/1‑139 (MQ=255)
                                                       tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                            <  2:727936/139‑1 (MQ=255)
                                                           gTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                            >  1:455509/1‑135 (MQ=255)
                                                           gTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                            <  2:455509/135‑1 (MQ=255)
                                                                       ttcTCCTATGTTTTGAAGTTAGTCTGTAACAGCTGGTGTAGTGGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAGCCATGGGCGCGGACCTGTTTCGGTTTGAAATGTTTTGCTTTTTGTATACGTCATGGCAGTGa            <  2:783741/137‑1 (MQ=255)
                                                                       tACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGa            <  2:100241/139‑1 (MQ=255)
                                                                          tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                   <  2:426254/65‑1 (MQ=255)
                                                                          tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                   >  1:426254/1‑65 (MQ=255)
                                                                          tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                 <  1:689714/115‑1 (MQ=255)
                                                                          tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                 >  2:689714/1‑115 (MQ=255)
                                                                                tttATGAAGTTATTCTGGAACAGCATGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCAt   <  1:71476/139‑1 (MQ=255)
                                                                                 ttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATc  <  2:321023/139‑1 (MQ=255)
                                                                                 ttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATc  >  2:613108/1‑139 (MQ=255)
                                                                                         ttATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                 >  2:296983/1‑100 (MQ=255)
                                                                                         ttATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                 <  1:296983/100‑1 (MQ=255)
                                                                                                            gTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                 <  2:4734/49‑1 (MQ=255)
                                                                                                            gTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                 >  1:4734/1‑49 (MQ=255)
                                                                                                               gACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATc  >  1:684348/1‑109 (MQ=255)
                                                                                                               gACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATc  <  2:684348/109‑1 (MQ=255)
                                                                                                                aCGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                      <  1:600027/56‑1 (MQ=255)
                                                                                                                aCGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                      >  2:600027/1‑56 (MQ=255)
                                                                                                                               |                                                                                            
TAGAACGGAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATC  >  NC_000913/1419032‑1419251

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