Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 582,148 Δ1 bp intergenic (+51/+4) tfaD → / ← ybcY pseudogene, DLP12 prophage, tail fiber assembly protein family,Phage or Prophage Related/pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913582,1470G.100.0% 234.6 / NA 56intergenic (+50/+5)tfaD/ybcYpseudogene, DLP12 prophage, tail fiber assembly protein family,Phage or Prophage Related/pseudogene, DLP12 prophage, methyltransferase homology,Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (31/25);  total (31/25)

GTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTT  >  NC_000913/582028‑582280
                                                                                                                       |                                                                                                                                     
gtTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGaa                                                                                                                   >  2:691189/1‑139 (MQ=255)
  tGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGc                                                                                                                 >  1:48067/1‑139 (MQ=255)
        aCCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           <  2:518575/139‑1 (MQ=255)
           gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  1:658753/139‑1 (MQ=255)
           gtgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  1:302047/139‑1 (MQ=255)
            tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       >  2:169300/1‑139 (MQ=255)
            tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       >  2:268323/1‑139 (MQ=255)
            tgtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAtt                                                                                                       >  1:379483/1‑139 (MQ=255)
             gtTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                      <  2:507808/139‑1 (MQ=255)
                      tCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTAc                                                                                             >  1:543189/1‑139 (MQ=255)
                          cTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTg                                                                                         >  1:573470/1‑139 (MQ=255)
                           tGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGa                                                                                        >  1:237023/1‑139 (MQ=255)
                              aCCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAAcc                                                                                     <  2:410626/139‑1 (MQ=255)
                               ccTGGTATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCt                                                                                    >  2:114781/1‑139 (MQ=255)
                                        gAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTaa                                                                           <  1:590931/139‑1 (MQ=255)
                                         aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgttc                                                                                                                                 >  2:600035/1‑79 (MQ=255)
                                         aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTgatgttc                                                                                                                                 <  1:600035/85‑7 (MQ=255)
                                         aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAAc                                                                          >  2:238434/1‑139 (MQ=255)
                                         aGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAAc                                                                          <  1:114781/139‑1 (MQ=255)
                                             gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                      >  1:493356/1‑139 (MQ=255)
                                             gCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtc                                                                      >  2:720101/1‑139 (MQ=255)
                                              ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCtt                                                                                          <  1:251631/118‑1 (MQ=255)
                                              ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCtt                                                                                          >  2:251631/1‑118 (MQ=255)
                                              ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCtt                                                                                   >  2:312486/1‑125 (MQ=255)
                                              ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCtt                                                                                   <  1:312486/125‑1 (MQ=255)
                                              ccTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATtct                                                                     >  2:278011/1‑139 (MQ=255)
                                                 aCGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGa                                                                  <  1:268323/139‑1 (MQ=255)
                                                       ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                            <  1:720101/139‑1 (MQ=255)
                                                       ccTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGct                                                            >  2:609224/1‑139 (MQ=255)
                                                         tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                      >  1:761860/1‑95 (MQ=255)
                                                         tGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATta                                                                                                      <  2:761860/95‑1 (MQ=255)
                                                          gTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATg                                                              >  1:112095/1‑134 (MQ=255)
                                                          gTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATg                                                              <  2:112095/134‑1 (MQ=255)
                                                            cAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTgag                                                       >  2:285909/1‑139 (MQ=255)
                                                               ggAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          >  2:389910/1‑133 (MQ=255)
                                                               ggAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  1:389910/133‑1 (MQ=255)
                                                                 aGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtacta                                                                                                   <  2:110914/90‑1 (MQ=255)
                                                                 aGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATtacta                                                                                                   >  1:110914/1‑90 (MQ=255)
                                                                  gTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 <  2:635820/139‑1 (MQ=255)
                                                                  gTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGATGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGtata                                                 <  2:515216/139‑1 (MQ=255)
                                                                     aTCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATAtgt                                              >  2:755812/1‑139 (MQ=255)
                                                                          tGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTaa                                         >  2:84064/1‑139 (MQ=255)
                                                                             gATTATGCCGCAGCACGTCTTAAGCAAGAACATACTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGcc                                      >  2:108134/1‑139 (MQ=255)
                                                                                           aCGTCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCtt                        <  1:285909/139‑1 (MQ=255)
                                                                                              tCTTAAGCAAGAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTtg                     <  2:692820/139‑1 (MQ=255)
                                                                                                        gAACATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGaaaa           >  1:155933/1‑139 (MQ=255)
                                                                                                         aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           >  2:527380/1‑42 (MQ=37)
                                                                                                         aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           >  2:559512/1‑42 (MQ=37)
                                                                                                         aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           <  1:527380/42‑1 (MQ=37)
                                                                                                         aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGaaaa                                                                                                           <  1:559512/42‑1 (MQ=37)
                                                                                                         aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        <  2:584712/45‑1 (MQ=255)
                                                                                                         aaCATGCTGCGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACAt                                                                                                        >  1:584712/1‑45 (MQ=37)
                                                                                                              gctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                        >  2:163447/1‑88 (MQ=255)
                                                                                                              gctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTga                                                        <  1:163447/88‑1 (MQ=255)
                                                                                                              gctgcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGcc     >  2:305521/1‑139 (MQ=255)
                                                                                                                 gcGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCttt  <  2:573470/139‑1 (MQ=255)
                                                                                                                  cGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          <  1:247122/82‑1 (MQ=255)
                                                                                                                  cGGTTGATGTTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGctct                                                          >  2:247122/1‑82 (MQ=255)
                                                                                                                       |                                                                                                                                     
GTTGCTGAACCGTGTTGATACATCAACTGCACCTGATATTGAGTGGCCTACGAACCCTGTCAGGGAGTAATCATTGGGATTATGCCGCAGCACGTCTTAAGCAAGAACATGCTGCGGTTGGATGCTATTTTTTTCCTGAAGCGGAAAACATTACTACAGTACCTTGAACCTTGGTTTTAACATTCTCGAAATGCTCTGAGAGTATATGTGTTAAGCCTTCTTCGGAATCTTTTGTGTTTGAAAAGATGCCTTT  >  NC_000913/582028‑582280

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