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% 70.7 / NA 18coding (324/1155 nt)yagACP4‑6 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (9/9);  total (9/9)

AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGG  >  NC_000913/281525‑281788
                                                                                                                                       |                                                                                                                                
aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                               <  1:404671/139‑3 (MQ=255)
     ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCCCGTCGTGTTCGAACCGGCCCGTGGCGGTAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGgtg                                                                                                                         <  1:792596/139‑1 (MQ=255)
                   aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                               >  1:284645/1‑118 (MQ=255)
                   aTCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                               <  2:284645/120‑3 (MQ=255)
                                                     aCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGTggca                                                                                                                               >  2:696765/1‑84 (MQ=255)
                                                     aCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggc                                                                                                                                >  2:125658/1‑84 (MQ=255)
                                                     aCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggc                                                                                                                                <  1:125658/85‑2 (MQ=255)
                                                                  gggAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCAt                                                            <  2:1057375/139‑1 (MQ=255)
                                                                   ggAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATg                                                           >  2:943748/1‑139 (MQ=255)
                                                                              tGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGccc                                                                                                                    <  1:391001/71‑1 (MQ=255)
                                                                              tGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGccc                                                                                                                    >  2:391001/1‑71 (MQ=255)
                                                                                         gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                     <  1:252621/139‑1 (MQ=255)
                                                                                                  cATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCa                            >  1:315748/1‑139 (MQ=255)
                                                                                                    tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGa                          >  1:1091136/1‑139 (MQ=255)
                                                                                                     ggCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGAt                         >  2:101236/1‑139 (MQ=255)
                                                                                                     ggCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGAt                         >  2:1212230/1‑139 (MQ=255)
                                                                                                       cGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                     >  1:116026/1‑125 (MQ=255)
                                                                                                       cGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                     <  2:116026/125‑1 (MQ=255)
                                                                                                        gggCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                      <  1:1212230/139‑1 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                     >  1:900199/1‑122 (MQ=255)
                                                                                                          gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                     <  2:900199/122‑1 (MQ=255)
                                                                                                                   ttATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa           <  2:315748/139‑1 (MQ=255)
                                                                                                                   ttATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa           <  1:1250439/139‑1 (MQ=255)
                                                                                                                            gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCgg  >  1:652238/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                
AAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGG  >  NC_000913/281525‑281788

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: