Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,419,159 G→A *163* (TAG→TAA) 
S51S (AGC→AGT) 		 sieB →
ydaF ←		 phage superinfection exclusion protein, Rac prophage
uncharacterized protein, Rac prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,419,1590GA100.0% 210.1 / NA 65*163* (TAG→TAA) 
S51S (AGC→AGT) 		sieB
ydaF		phage superinfection exclusion protein, Rac prophage
uncharacterized protein, Rac prophage
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/35);  total (30/35)

GAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACA  >  NC_000913/1419039‑1419274
                                                                                                                        |                                                                                                                   
gAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTcc                                                                                                   <  2:1161848/139‑1 (MQ=255)
           ttCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        <  1:1256322/139‑1 (MQ=255)
            tCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcgc                                                                                       <  1:1476384/139‑1 (MQ=255)
                cATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        >  2:920586/1‑134 (MQ=255)
                cATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        <  1:920586/134‑1 (MQ=255)
                      agTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                             >  1:481175/1‑139 (MQ=255)
                          aTAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGt                                                                         >  1:1393878/1‑139 (MQ=255)
                           tAGCGCAAACTGTTCTTATTCTCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACt                                                                                                               >  1:948410/1‑100 (MQ=255)
                           tAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACt                                                                                                               <  2:948410/100‑1 (MQ=255)
                             gcgcAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGtttt                                                                      >  1:1199786/1‑139 (MQ=255)
                                  aaCTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                          <  1:902516/98‑1 (MQ=255)
                                  aaCTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                          >  2:902516/1‑98 (MQ=255)
                                  aaCTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATg                                                                 >  2:746221/1‑139 (MQ=255)
                                         cTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                          >  2:661367/1‑139 (MQ=255)
                                         cTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                          <  1:1201412/139‑1 (MQ=255)
                                           tattatCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                        <  1:746221/139‑1 (MQ=255)
                                               atCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGt                                                    >  1:72512/1‑139 (MQ=255)
                                                tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTATACGTTTTAATCGTTAACTTATTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   <  2:1322124/139‑1 (MQ=255)
                                                tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   <  1:661367/139‑1 (MQ=255)
                                                tctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   <  2:1199786/139‑1 (MQ=255)
                                                 ctcGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGtat                                                  >  2:624466/1‑139 (MQ=255)
                                                   cGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   >  2:1013462/1‑136 (MQ=255)
                                                   cGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCAGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   <  1:1013462/136‑1 (MQ=255)
                                                    gTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        >  1:10529/1‑98 (MQ=255)
                                                    gTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        <  2:10529/98‑1 (MQ=255)
                                                    gTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACg                                               >  1:999007/1‑139 (MQ=255)
                                                      cACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTc                                             >  2:1299412/1‑139 (MQ=255)
                                                      cACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTT‑ATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        >  1:1613955/1‑95 (MQ=255)
                                                      cACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTT‑ATCGTTAACTTACTGTGTGCTTCTCCAACCATCGgcg                                                                                        <  2:1613955/95‑1 (MQ=255)
                                                          gAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATgg                                         >  1:750074/1‑139 (MQ=255)
                                                                  cTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTc                                                                            >  1:1587326/1‑96 (MQ=255)
                                                                  cTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTc                                                                            <  2:1587326/96‑1 (MQ=255)
                                                                   tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCATTAACTTACTGTGTGCTTCTCCAACCATCg                                                                                           >  2:711101/1‑80 (MQ=255)
                                                                   tCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCATTAACTTACTGTGTGCTTCTCCAACCATCg                                                                                           <  1:711101/80‑1 (MQ=255)
                                                                          ttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                        <  1:1582576/108‑1 (MQ=255)
                                                                          ttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                        >  2:1582576/1‑108 (MQ=255)
                                                                          ttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGt                                                    <  1:912520/112‑1 (MQ=255)
                                                                          ttATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGt                                                    >  2:912520/1‑112 (MQ=255)
                                                                             tGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                        <  1:159726/105‑1 (MQ=255)
                                                                             tGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCttt                                                        >  2:159726/1‑105 (MQ=255)
                                                                               aaGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCAc                                                                                   >  2:967621/1‑76 (MQ=255)
                                                                               aaGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCAc                                                                                   <  1:967621/76‑1 (MQ=255)
                                                                                 gTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   <  2:1249979/106‑1 (MQ=255)
                                                                                 gTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGta                                                   >  1:1249979/1‑106 (MQ=255)
                                                                                 gTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCa                                       <  1:801859/118‑1 (MQ=255)
                                                                                 gTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCa                                       >  2:801859/1‑118 (MQ=255)
                                                                                 gTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTg                  >  1:978716/1‑139 (MQ=255)
                                                                                    aTTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                          >  2:157860/1‑96 (MQ=255)
                                                                                    aTTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCt                                                          <  1:157860/96‑1 (MQ=255)
                                                                                       cTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                             <  2:1026562/74‑1 (MQ=255)
                                                                                       cTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGttt                                                                             >  1:1026562/1‑74 (MQ=255)
                                                                                       cTGGAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAac            <  1:624466/139‑1 (MQ=255)
                                                                                          gAACAGCAGGAGTAGACGTTTTAATCGTTAACTTACTGGGTGCTTTTCCAACCTTCGGGGCGCACGAGTTTTGTTTGTAAATGTTTAGCTTTTGGTTTACGTCATGGCAGTTAACGTTCCATCCTGGTTGGGGAACAcg         <  1:582287/139‑1 (MQ=255)
                                                                                             cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg      <  2:72512/139‑1 (MQ=255)
                                                                                             cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg      <  2:750074/139‑1 (MQ=255)
                                                                                             cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg      <  2:481175/139‑1 (MQ=255)
                                                                                             cagcagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACAcgccg      <  1:778615/139‑1 (MQ=255)
                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                          >  2:661112/1‑37 (MQ=255)
                                                                                               gcagGAGTAGACGTTTTAATCGTTAACTTACtgtgtg                                                                                                          <  1:661112/37‑1 (MQ=255)
                                                                                                cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCAt                                                                                             <  1:309335/49‑1 (MQ=255)
                                                                                                cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCAt                                                                                             >  2:309335/1‑49 (MQ=255)
                                                                                                cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGtt                              <  1:629756/112‑1 (MQ=255)
                                                                                                cagGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGtt                              >  2:629756/1‑112 (MQ=255)
                                                                                                 agGAGTAGACGTTTTAATCGTTAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGcaca  <  1:1231974/139‑1 (MQ=255)
                                                                                                                     ttAACTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGtt                                                                              <  1:1508214/43‑1 (MQ=255)
                                                                                                                        |                                                                                                                   
GAAGGGCATCATTCAACATCAGAGTGATAGCGCAAACTGTTCTTATTATCTCGTCACCGAAAAATACTCCCATTTTATGAAGTTATTCTGGAACAGCAGGAGTAGACGTTTTAATCGTTAGCTTACTGTGTGCTTCTCCAACCATCGGCGCGCACCAGTTTCGGTTTTAAATGTTTTGCTTTTGGTATACGTCATGGCAGTGAACGTTCCATCCTGGTTGGGGAACACGCCGCACA  >  NC_000913/1419039‑1419274

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