Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,662 Δ1 bp coding (322/1155 nt) yagA ← CP4‑6 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,6600G.100.0% 120.4 / NA 29coding (324/1155 nt)yagACP4‑6 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (14/15);  total (14/15)

GCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGA  >  NC_000913/281531‑281789
                                                                                                                                 |                                                                                                                                 
gCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGgtgt                                                                                                                         >  2:558678/1‑139 (MQ=255)
     tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGcccc                                                                                                                    <  1:939313/139‑1 (MQ=255)
              tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                                <  2:709063/119‑3 (MQ=255)
              tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                                >  1:709063/1‑117 (MQ=255)
                  ccAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGAAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGcc                                                                                                       <  2:598455/139‑1 (MQ=255)
                   cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCa                                                                                                      <  2:947101/139‑1 (MQ=255)
                                                  ggCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                                           <  2:807060/135‑1 (MQ=255)
                                                  ggCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                                           >  1:807060/1‑135 (MQ=255)
                                                     ccGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    <  1:680110/139‑1 (MQ=255)
                                                         ggCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGcccc                                                                                                                    <  1:236125/87‑1 (MQ=255)
                                                         ggCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGcccc                                                                                                                    >  2:236125/1‑87 (MQ=255)
                                                          gCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTc                                                               <  2:786515/139‑1 (MQ=255)
                                                            gggAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCAt                                                             <  1:622527/139‑1 (MQ=255)
                                                                  gCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    <  1:293162/126‑1 (MQ=255)
                                                                  gCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    >  2:293162/1‑126 (MQ=255)
                                                                         gAAGCGCCCGGCAGCGGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACgcag                                                >  2:777694/1‑139 (MQ=255)
                                                                             cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTGTGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAgg                                            >  1:677156/1‑139 (MQ=255)
                                                                                   gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                      <  1:558678/139‑1 (MQ=255)
                                                                                         ggCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATgtcg                                >  2:775610/1‑139 (MQ=255)
                                                                                          gCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    <  2:960587/102‑1 (MQ=255)
                                                                                          gCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    >  1:960587/1‑102 (MQ=255)
                                                                                               ggCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                                           >  2:10551/1‑90 (MQ=255)
                                                                                               ggCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCg                                                                           <  1:10551/90‑1 (MQ=255)
                                                                                                    gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    <  1:931700/92‑1 (MQ=255)
                                                                                                    gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    <  1:788353/92‑1 (MQ=255)
                                                                                                    gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    >  2:788353/1‑92 (MQ=255)
                                                                                                    gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                    >  2:931700/1‑92 (MQ=255)
                                                                                                                      gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                       >  2:462252/1‑119 (MQ=255)
                                                                                                                      gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                       <  1:462252/119‑1 (MQ=255)
                                                                                                                       tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGa  >  2:210688/1‑139 (MQ=255)
                                                                                                                       tGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGa  >  1:659149/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                 
GCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGA  >  NC_000913/281531‑281789

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