Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,429,684 A→G S212S (TCA→TCG)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,6840AG100.0% 160.3 / NA 47S212S (TCA→TCG) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/26);  total (21/26)

GTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGC  >  NC_000913/1429537‑1429827
                                                                                                                                                   |                                                                                                                                               
gtcgtcGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCg                                                                                                                                              <  2:90600/151‑1 (MQ=255)
                  gTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGcc                                                                                                                            >  1:319894/1‑151 (MQ=255)
                      tCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCca                                                                                                                           <  2:479619/148‑1 (MQ=255)
                       ccAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACAt                                                                                                                        <  2:511219/150‑1 (MQ=255)
                          gcgcAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGc                                                                                                                    >  2:520299/1‑151 (MQ=255)
                           cgcAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAGTCGCGTCGGGCGTCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCTTCACTACAATCAGCTGTCACATCTGCa                                                                                                                   <  2:716078/151‑1 (MQ=255)
                                  aCGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACcg                                                                                                            <  1:626502/151‑1 (MQ=255)
                                   cgccATCAGCAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGCGTCCTCAAAAAGCGCGGCGCCCACCAGTGCCGGTGCGGCGAAAACGTCAGAACCGAATGCTCCGGCGTCACTACAATCAGCAGCCCCATCTGCATCCAcc                                                                                                             >  2:420396/4‑149 (MQ=255)
                                       aTCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGAcc                                                                                                       >  2:477964/1‑151 (MQ=255)
                                        tCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGAcc                                                                                                       >  2:471187/1‑150 (MQ=255)
                                            caaAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGa                                                                                                   >  2:735633/1‑150 (MQ=255)
                                            caaAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGa                                                                                                   >  2:297644/1‑150 (MQ=255)
                                             aaaGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGa                                                                                                   >  1:444674/1‑149 (MQ=255)
                                               aGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGc                                                                                               <  2:100548/151‑1 (MQ=255)
                                                ggCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGc                                                                                               >  2:627789/1‑150 (MQ=255)
                                                 gCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCAt                                                                                             >  2:469215/1‑151 (MQ=255)
                                                  ccACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATc                                                                                            >  2:637574/1‑151 (MQ=255)
                                                    aCTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAg                                                                                          <  2:419600/151‑1 (MQ=255)
                                                               aaaaaGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGAcc                                                                               <  1:373937/151‑1 (MQ=255)
                                                                 aaaGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTc                                                                             >  1:200488/1‑151 (MQ=255)
                                                                  aaGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCg                                                                            <  1:71831/151‑1 (MQ=255)
                                                                  aaGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCg                                                                            <  2:469924/151‑1 (MQ=255)
                                                                          ccGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGa                                                                    <  2:301851/151‑1 (MQ=255)
                                                                                      ccTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTc                                                        <  1:164705/151‑1 (MQ=255)
                                                                                        tCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCa                                                       >  2:220149/1‑150 (MQ=255)
                                                                                         cAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCaaa                                                     <  2:761661/151‑1 (MQ=255)
                                                                                          aaaaaGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCaaa                                                     >  2:264424/1‑150 (MQ=255)
                                                                                           aaaaGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCaaa                                                     >  1:4113/1‑149 (MQ=255)
                                                                                                cgcgGGGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCCGCCTCAAaag                                                   >  2:239827/1‑146 (MQ=255)
                                                                                                      ggcCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCaaaa                                        <  2:683906/151‑1 (MQ=255)
                                                                                                             aGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAg                                 <  2:719052/151‑1 (MQ=255)
                                                                                                                gCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGaaa                              <  1:57593/151‑1 (MQ=255)
                                                                                                                       cggcgAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGc                       >  2:513012/1‑151 (MQ=255)
                                                                                                                        ggcgAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGca                      <  1:742998/151‑1 (MQ=255)
                                                                                                                           gAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGcatc                    >  2:30105/1‑150 (MQ=255)
                                                                                                                           gAAAACGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGcatc                    <  1:700719/150‑1 (MQ=255)
                                                                                                                            aaaaCGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGcatcat                  >  1:747977/1‑151 (MQ=255)
                                                                                                                               aCGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCaa               <  1:269208/151‑1 (MQ=255)
                                                                                                                                cGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAg              <  2:200480/151‑1 (MQ=255)
                                                                                                                                cGTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAg              <  2:368103/151‑1 (MQ=255)
                                                                                                                                 gTCAGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAgc             <  1:532486/151‑1 (MQ=255)
                                                                                                                                    aGAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGcc          >  1:65516/1‑151 (MQ=255)
                                                                                                                                     gAAACGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCa         <  1:71294/151‑1 (MQ=255)
                                                                                                                                       aaCGAATGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCagt       <  1:613460/151‑1 (MQ=255)
                                                                                                                                           aaTGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTg   <  2:194134/151‑1 (MQ=255)
                                                                                                                                            aTGCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGc  <  2:644512/151‑1 (MQ=255)
                                                                                                                                              gCTTCGGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGc  >  1:718427/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTCGTCGGCTCAGTCAGCGTCTTCCAGCGCAGGAACGGCATCAACAAAGGCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGCAAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGC  >  NC_000913/1429537‑1429827

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