Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,679 T→C D102G (GAC→GGC)  yagA ← CP4‑6 prophage; integrase core domain‑containing protein YagA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,6790TC100.0% 165.5 / NA 47D102G (GAC→GGC) yagACP4‑6 prophage; integrase core domain‑containing protein YagA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/28);  total (19/28)

CGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACC  >  NC_000913/281569‑281812
                                                                                                              |                                                                                                                                     
cgtGTTCGTACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCCGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCt                                                                                                          >  2:419967/1‑140 (MQ=255)
cgtGTTCGAACCGGCCCGTGGCGNNAATGCCCGGTGGAGCGCCCGGCAGGAGGCCATGGCGGGCCATCAGGGTATGGACGGGGTTGGAGGCGGGCCTGGTGTGCCCCTGGGCCCTGAGGCCACGCTTTATTTTTTGGGCt                                                                                                          >  2:1152941/1‑140 (MQ=255)
              cccGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATg                                                                                            >  1:273481/1‑140 (MQ=255)
               ccGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGa                                                                                           <  2:87841/140‑1 (MQ=255)
                  tGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGGGCTGGAGGCGGGGATGGGGTGCCCCCGGCCCCCGAGCCCGCGCTTTATCTTGCGGGCTCCCCAGCGGTCCTGACgg                                                                                        >  2:442646/1‑140 (MQ=255)
                   ggCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCCTCAGGTTATGGAAGGTGCTGAAGGCGGGCATGGGGTGCCCCCGGCCCCCGAGCCAGCGCCTAATCTTGCGGGCTCCCCAGCGTTCATGACggg                                                                                       >  2:260686/1‑139 (MQ=255)
                      gggAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCAt                                                                                    <  1:681136/140‑1 (MQ=255)
                              ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGcc                                                                               <  2:372616/137‑1 (MQ=255)
                              ccGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGcc                                                                               >  1:372616/1‑137 (MQ=255)
                                  tGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACgca                                                                        >  1:328928/1‑140 (MQ=255)
                                   gAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACgcag                                                                       >  1:147607/1‑140 (MQ=255)
                                        gcCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAggg                                                                  <  1:377017/140‑1 (MQ=255)
                                            ggCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGt                                                              >  1:179396/1‑140 (MQ=255)
                                             gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                                             <  1:154604/140‑1 (MQ=255)
                                             gcagcaGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                                             <  1:911869/140‑1 (MQ=255)
                                             gcagcaGGCCATGGCGGGCCATCAGGTTATGGAAGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTg                                                             <  1:260686/140‑1 (MQ=255)
                                               agcagGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGAt                                                           <  1:890525/140‑1 (MQ=255)
                                                        tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGAt                                                           <  1:624457/131‑1 (MQ=255)
                                                        tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGAt                                                           >  2:624457/1‑131 (MQ=255)
                                                        tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGa                                                  >  2:60322/1‑140 (MQ=255)
                                                        tGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGa                                                  >  2:138202/1‑140 (MQ=255)
                                                            gggCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                                              <  2:328928/140‑1 (MQ=255)
                                                              gCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                                              >  1:1160875/1‑138 (MQ=255)
                                                              gCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAg                                              <  2:1160875/138‑1 (MQ=255)
                                                              gCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc                                             >  2:74526/1‑139 (MQ=255)
                                                              gCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGc                                             <  1:74526/139‑1 (MQ=255)
                                                               ccATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCgg                                           >  1:79667/1‑140 (MQ=255)
                                                                cATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGt                                          <  2:179396/140‑1 (MQ=255)
                                                                cATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGt                                          <  2:988799/140‑1 (MQ=255)
                                                                   cAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCg                                       <  1:88172/140‑1 (MQ=255)
                                                                     ggTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCggg                                     <  1:208336/140‑1 (MQ=255)
                                                                       ttATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa                                   <  1:175948/140‑1 (MQ=255)
                                                                       ttATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa                                   <  1:419967/140‑1 (MQ=255)
                                                                       ttATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGa                                   <  1:60322/140‑1 (MQ=255)
                                                                        tATGGACGGTGCTGAAGGCGGGCATGGTGTACCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGaa                                  <  1:969187/140‑1 (MQ=255)
                                                                         aTGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAAt                                 <  1:1152941/140‑1 (MQ=255)
                                                                                gTGCTGAAGGCGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCGGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCgg                          >  2:803429/1‑140 (MQ=255)
                                                                                          cGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  2:389443/140‑1 (MQ=255)
                                                                                          cGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  2:273481/140‑1 (MQ=255)
                                                                                          cGGGCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  2:147607/140‑1 (MQ=255)
                                                                                           gggCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                >  1:607287/1‑139 (MQ=255)
                                                                                           gggCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                >  1:8535/1‑139 (MQ=255)
                                                                                           gggCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  2:607287/139‑1 (MQ=255)
                                                                                           gggCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  2:1098988/139‑1 (MQ=255)
                                                                                           gggCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  2:8535/139‑1 (MQ=255)
                                                                                           gggCATGGTGTGCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                >  1:1098988/1‑139 (MQ=255)
                                                                                                   tgtgCCCCTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGa       >  2:555820/1‑140 (MQ=255)
                                                                                                       ccccTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                >  2:312392/1‑127 (MQ=255)
                                                                                                       ccccTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCgcggg                <  1:312392/127‑1 (MQ=255)
                                                                                                        cccTGGCCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGAcc  >  1:762642/1‑140 (MQ=255)
                                                                                                              |                                                                                                                                     
CGTGTTCGAACCGGCCCGTGGCGGGAATGCCCGGTGAAGCGCCCGGCAGCAGGCCATGGCGGGCCATCAGGTTATGGACGGTGCTGAAGGCGGGCATGGTGTGCCCCTGGTCCTCGAGCCAGCGCTTAATCTTGCGGGCTCCCCAGCGTTCATGACGGTCATGGGCCATACGCAGCAGGGCCGTGATGTCGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACC  >  NC_000913/281569‑281812

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