Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,764,608 T→C I182M (ATA→ATG yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,764,6080TC100.0% 181.5 / NA 51I182M (ATA→ATGyfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/22);  total (29/22)

CATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCATCGAT  >  NC_000913/2764478‑2764742
                                                                                                                                  |                                                                                                                                      
catcaGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                                <  1:925542/139‑1 (MQ=255)
         tgtgAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                       >  2:736616/1‑139 (MQ=255)
                  gTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                                <  1:1167367/121‑1 (MQ=255)
                  gTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCtt                                                                                                                                >  2:1167367/1‑121 (MQ=255)
                   tCGAGAATGCTGGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGGTCCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGt                                                                                                             >  1:447121/1‑139 (MQ=255)
                     gagaATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTc                                                                                                           >  2:559127/1‑139 (MQ=255)
                                                  tCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCt                                                                                                                         >  1:1101189/1‑96 (MQ=255)
                                                  tCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCt                                                                                                                         <  2:1101189/96‑1 (MQ=255)
                                                   cAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGc                                                                             >  1:554646/1‑139 (MQ=255)
                                                     aGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                       >  2:472645/1‑95 (MQ=255)
                                                     aGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                       <  1:472645/95‑1 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATc                                                                                                                            >  1:1089369/1‑83 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATc                                                                                                                            <  2:1089369/83‑1 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                       >  2:714170/1‑88 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAg                                                                                                                       <  1:714170/88‑1 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                         <  2:557293/134‑1 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                         >  1:557293/1‑134 (MQ=255)
                                                            ggCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGAt                                                                    <  2:963669/139‑1 (MQ=255)
                                                              cACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                                                                                            >  2:952513/1‑113 (MQ=255)
                                                              cACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                                                                                            <  1:952513/113‑1 (MQ=255)
                                                                      gTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                          <  2:1094399/139‑1 (MQ=255)
                                                                        cTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAAc                                                        >  2:1101828/1‑139 (MQ=255)
                                                                             gCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTta                                                   >  2:530819/1‑139 (MQ=255)
                                                                                       cGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTc                                         <  1:559127/139‑1 (MQ=255)
                                                                                           aCTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGtt                                     <  1:530819/139‑1 (MQ=255)
                                                                                              gagaATGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTaaa                                  >  1:472726/1‑139 (MQ=255)
                                                                                                 aaTGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                          >  1:785853/1‑112 (MQ=255)
                                                                                                 aaTGACCCCATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                          <  2:785853/112‑1 (MQ=255)
                                                                                                         cATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAAc                                                        >  1:125524/1‑106 (MQ=255)
                                                                                                         cATGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAAc                                                        <  2:125524/106‑1 (MQ=255)
                                                                                                           tGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                                                                                            >  1:570937/1‑68 (MQ=255)
                                                                                                           tGTAGATCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGaaaa                                                                                            <  2:570937/68‑1 (MQ=255)
                                                                                                                aTCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTc                                                                                      <  1:669761/69‑1 (MQ=255)
                                                                                                                aTCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTc                                                                                      >  2:669761/1‑69 (MQ=255)
                                                                                                                aTCTTGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTttg                >  1:459281/1‑139 (MQ=255)
                                                                                                                    tGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGa            >  1:255773/1‑139 (MQ=255)
                                                                                                                    tGGGAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGa            >  2:881110/1‑139 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                      >  2:1145483/1‑48 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCa                                                                                                      <  1:1145483/48‑1 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                         <  2:640296/77‑1 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                         >  1:640296/1‑77 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAg           >  2:363576/1‑139 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAg           >  2:96840/1‑139 (MQ=255)
                                                                                                                     gggAGGCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTCATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAg           >  2:384409/1‑139 (MQ=255)
                                                                                                                         ggCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTCATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCa       <  1:384409/139‑1 (MQ=255)
                                                                                                                          gCCAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCAt      >  1:930006/1‑139 (MQ=255)
                                                                                                                           ccAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                         <  2:1194893/71‑1 (MQ=255)
                                                                                                                           ccAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCa                                                                         >  1:1194893/1‑71 (MQ=255)
                                                                                                                            cAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                          >  2:99659/1‑85 (MQ=255)
                                                                                                                            cAGCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGaa                                                          <  1:99659/85‑1 (MQ=255)
                                                                                                                              gCCGCATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCATCGAt  <  1:1101828/139‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                      
CATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTCACAGTATCAAGGATGCGGCACCAAATGTCTTCGGCGACGGTTTCGGAACTGAGAATGACCCCATGTAGATCTTGGGAGGCCAGCCGTATGTTCTTGATCCCTAGAAGTTCCAGTTCGCGCAGGGATGAAAACACTTCCCATGTAGCCGCATCGATCCAATCTGAAACGCTTATACCTGAGTCTGTTAAAAAATTGTCTGTGCGTTTGTTGAGGTCATCGAT  >  NC_000913/2764478‑2764742

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: