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

GAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAT  >  NC_000913/281524‑281781
                                                                                                                                        |                                                                                                                         
gAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggc                                                                                                                         <  2:1232950/139‑2 (MQ=255)
gAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggc                                                                                                                         <  2:113229/139‑2 (MQ=255)
 aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                        <  2:244609/139‑3 (MQ=255)
 aaaGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggca                                                                                                                        <  1:139800/139‑3 (MQ=255)
      ggCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGgtg                                                                                                                  <  2:881654/139‑1 (MQ=255)
            tAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGcccc                                                                                                            <  2:1076159/139‑1 (MQ=255)
                     tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggcatg                                                                                                                      >  2:394362/1‑117 (MQ=255)
                     tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCggcatg                                                                                                                      <  1:394362/121‑5 (MQ=255)
                     tCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCg                                                                                                   >  1:493226/1‑139 (MQ=255)
                          cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCa                                                                                              <  1:700487/139‑1 (MQ=255)
                          cAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCa                                                                                              >  1:17230/1‑139 (MQ=255)
                                   ttCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGc                                                                                                <  1:172727/128‑1 (MQ=255)
                                   ttCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGc                                                                                                >  2:172727/1‑128 (MQ=255)
                                         gcgTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCg                                                                               <  1:466199/139‑1 (MQ=255)
                                          cgtcgtGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCgg                                                                              <  2:17230/139‑1 (MQ=255)
                                                      aCCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGt                                                                  >  1:112121/1‑139 (MQ=255)
                                                         ggCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCa                                                               >  1:1137100/1‑139 (MQ=255)
                                                                      aaTGCCCGGTGAAGCCCCCGGCAGCAGGCCATGGGGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATggg                                                  <  2:1113052/139‑1 (MQ=255)
                                                                          cccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc                                                                       >  2:929604/1‑114 (MQ=255)
                                                                          cccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTcccc                                                                       <  1:929604/114‑1 (MQ=255)
                                                                           ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCAt                                                     <  1:1284248/131‑1 (MQ=255)
                                                                           ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCAt                                                     >  2:1284248/1‑131 (MQ=255)
                                                                                gAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTccc                                                                        <  1:1330805/107‑1 (MQ=255)
                                                                                gAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTccc                                                                        >  2:1330805/1‑107 (MQ=255)
                                                                                  aGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACgcagca                                      <  2:493226/139‑1 (MQ=255)
                                                                                    cgcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAgg                                    >  1:1296735/1‑139 (MQ=255)
                                                                                         ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTg                                                                                                          >  2:637887/1‑64 (MQ=255)
                                                                                         ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTg                                                                                                          <  1:637887/64‑1 (MQ=255)
                                                                                                    aTGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGc                                                                                                <  2:174891/63‑1 (MQ=255)
                                                                                                    aTGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGc                                                                                                >  1:174891/1‑63 (MQ=255)
                                                                                                     tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGa                   >  1:991423/1‑139 (MQ=255)
                                                                                                      ggCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAgcgc                                                                                          <  2:255673/67‑1 (MQ=255)
                                                                                                      ggCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAgcgc                                                                                          >  1:255673/1‑67 (MQ=255)
                                                                                                        cGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                            <  2:222154/95‑1 (MQ=255)
                                                                                                        cGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                            >  1:222154/1‑95 (MQ=255)
                                                                                                         gggCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg               <  2:1137100/139‑1 (MQ=255)
                                                                                                          ggCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc              <  2:112121/139‑1 (MQ=255)
                                                                                                           gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCg             >  2:222115/1‑139 (MQ=255)
                                                                                                           gCCATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCg             >  2:738672/1‑139 (MQ=255)
                                                                                                            ccATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg               <  2:41964/136‑1 (MQ=255)
                                                                                                            ccATCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg               >  1:41964/1‑136 (MQ=255)
                                                                                                               tCAGGTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCAGTGATGTCGTCAGATGAGCGGTTc         >  1:753735/1‑139 (MQ=255)
                                                                                                                  ggTTATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCggg      <  2:753735/139‑1 (MQ=255)
                                                                                                                    ttATGGACGGTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa    <  1:222115/139‑1 (MQ=255)
                                                                                                                             gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAt  >  1:217451/1‑132 (MQ=255)
                                                                                                                             gTGCTGAAGGCGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAt  <  2:217451/132‑1 (MQ=255)
                                                                                                                                        |                                                                                                                         
GAAAGTGGCCCTTAAAATCCATCTGCCAGAGGCGGTTCGGCGCGTCGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAT  >  NC_000913/281524‑281781

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