Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,891 A→G R31R (CGT→CGC)  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,8910AG100.0% 193.1 / NA 54R31R (CGT→CGC) yagACP4‑6 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/33);  total (21/33)

CGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAT  >  NC_000913/281758‑282021
                                                                                                                                     |                                                                                                                                  
cgtcAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGAt                                                                                                                               >  1:156869/1‑139 (MQ=255)
 gtcAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATg                                                                                                                              <  1:934632/139‑1 (MQ=255)
       tGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcc                                                                                                                        <  1:479448/139‑1 (MQ=255)
       tGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcc                                                                                                                        <  2:1005384/139‑1 (MQ=255)
         aGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcccc                                                                                                                      >  2:74168/1‑139 (MQ=255)
         aGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcccc                                                                                                                      >  1:159290/1‑139 (MQ=255)
         aGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcc                                                                                                                        >  1:138987/1‑137 (MQ=255)
         aGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcc                                                                                                                        <  2:138987/137‑1 (MQ=255)
           cGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGt                                                                                                                             <  1:1337257/130‑1 (MQ=255)
           cGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGt                                                                                                                             >  2:1337257/1‑130 (MQ=255)
                        ggTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaa                                                                                                       >  1:978391/1‑139 (MQ=255)
                         gTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaac                                                                                                      <  2:813401/139‑1 (MQ=255)
                          tGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcc                                                                                                                        >  2:552098/1‑120 (MQ=255)
                          tGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcc                                                                                                                        <  1:552098/120‑1 (MQ=255)
                             ggAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaacaaaa                                                                                                  <  1:74168/139‑1 (MQ=255)
                              gAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaacaaaac                                                                                                 <  2:159290/139‑1 (MQ=255)
                              gAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaacaaaac                                                                                                 <  1:271273/139‑1 (MQ=255)
                                 tGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAacaaaacaaa                                                                                              >  1:1294684/1‑139 (MQ=255)
                                   cgcgGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcccc                                                                                                                      >  2:712419/1‑113 (MQ=255)
                                   cgcgGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGcccc                                                                                                                      <  1:712419/113‑1 (MQ=255)
                                     cgggcggTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCa                                                                                          <  2:1121207/139‑1 (MQ=255)
                                      gggcggTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAg                                                                                         <  2:156545/139‑1 (MQ=255)
                                      gggcggTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAg                                                                                         <  2:107108/139‑1 (MQ=255)
                                      gggcggTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAg                                                                                         <  1:409090/139‑1 (MQ=255)
                                          ggTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTGCTGCGAGGCGAACAAAACAAACTCAGTACg                                                                                     <  1:705373/139‑1 (MQ=255)
                                             ccTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTaa                                                                                  <  1:177297/139‑1 (MQ=255)
                                                 aaGACCTGCGGCAGCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGTCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGAc                                                                              <  1:511943/139‑1 (MQ=255)
                                                      cggcggCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaac                                                                                                      <  2:1154807/110‑1 (MQ=255)
                                                      cggcggCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGaac                                                                                                      >  1:1154807/1‑110 (MQ=255)
                                                      cggcggCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGt                                                                         >  1:1320922/1‑139 (MQ=255)
                                                       ggcggcACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTa                                                                        <  1:677682/139‑1 (MQ=255)
                                                              ccTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGGCCCGTCCTGCGAGGCGAACAAAACAACCTCAGTACGTAATGACAGGGTAACTCTCg                                                                 <  1:1467543/139‑1 (MQ=255)
                                                              ccTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCg                                                                 <  1:359629/139‑1 (MQ=255)
                                                                         ccAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCg                                                                 <  2:937741/128‑1 (MQ=255)
                                                                         ccAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCg                                                                 >  1:937741/1‑128 (MQ=255)
                                                                                       aCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAg                                                                                         >  2:1468679/1‑90 (MQ=255)
                                                                                       aCTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAg                                                                                         <  1:1468679/90‑1 (MQ=255)
                                                                                            tAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACg                                   <  2:1051945/139‑1 (MQ=255)
                                                                                              gCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACt                                                                                            <  1:63971/80‑1 (MQ=255)
                                                                                              gCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACt                                                                                            >  2:63971/1‑80 (MQ=255)
                                                                                                   tGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATccc                                                           <  1:727088/108‑1 (MQ=255)
                                                                                                   tGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATccc                                                           >  2:727088/1‑108 (MQ=255)
                                                                                                   tGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCtt                            >  2:1167421/1‑139 (MQ=255)
                                                                                                           gAAATGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTa                    >  1:659563/1‑139 (MQ=255)
                                                                                                             aaTGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATg                                                                                <  1:1285513/77‑1 (MQ=255)
                                                                                                             aaTGCCGAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATg                                                                                >  2:1285513/1‑77 (MQ=255)
                                                                                                                   gAAGCGATGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAg            <  1:1474578/139‑1 (MQ=255)
                                                                                                                   gAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCTAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTAt                                                                       <  1:69385/80‑1 (MQ=255)
                                                                                                                   gAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAATCAAACTCAGTACGTAATGACATGGTAt                                                                       >  2:69385/1‑80 (MQ=255)
                                                                                                                   gAAGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAg            <  2:978391/139‑1 (MQ=255)
                                                                                                                     aGCGACGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATACCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAgtg          >  1:238224/1‑139 (MQ=255)
                                                                                                                         aCGGCAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTa      >  2:775974/1‑139 (MQ=255)
                                                                                                                             cAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAt  <  2:726069/139‑1 (MQ=255)
                                                                                                                             cAGAGGGAGCGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAt  <  1:1167421/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                  
CGTCAGATGAGCGGTTCGGGGAATGGTGCGGAATGCGCGGGCGGTCCTGAAGACCGGCGGCACCTTCCTGAGCCCAGCGCTGGAGCCACTTGTAGCCGGTGGCAGGTGAAATGCCGAAGCGACGGCAGAGGGAACGGATGTTCGCCCCGTCCTGCGAGGCGAACAAAACAAACTCAGTACGTAATGACATGGTATCTCTCGCATCCCAGGGCATAAGCGACTCCATAAACGGGTTCTTATGCCTTAGTTGTAAGTGTCTACCAT  >  NC_000913/281758‑282021

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