Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,431,882 T→C V945A (GTG→GCG)  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,431,8820TC100.0% 106.2 / NA 31V945A (GTG→GCG) stfRRac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/17);  total (14/17)

CCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAA  >  NC_000913/1431746‑1432014
                                                                                                                                         |                                                                                                                                    
ccaccaGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGa                                                                                                                                    <  2:271155/139‑1 (MQ=255)
    caGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATcc                                                                                                                                 <  2:136594/139‑1 (MQ=255)
                                 ccGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCa                                                                                                   <  2:842088/139‑1 (MQ=255)
                                    tGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCaaa                                                                                                 <  2:855873/138‑1 (MQ=255)
                                    tGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCaaa                                                                                                 >  1:855873/1‑138 (MQ=255)
                                                aCCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTc                                                                                    >  1:910479/1‑139 (MQ=255)
                                                     tCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                         >  1:603632/1‑97 (MQ=255)
                                                     tCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                         <  2:603632/97‑1 (MQ=255)
                                                        gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                            >  1:781999/1‑139 (MQ=255)
                                                                          tGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                           <  2:705117/139‑1 (MQ=255)
                                                                              ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCAGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGc                                                       >  2:660676/1‑139 (MQ=255)
                                                                               gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGt                                                                             >  1:350529/1‑116 (MQ=255)
                                                                               gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGt                                                                             <  2:350529/116‑1 (MQ=255)
                                                                               gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCa                                                      >  1:472721/1‑139 (MQ=255)
                                                                                 gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCAtt                                                    <  2:910479/139‑1 (MQ=255)
                                                                                     gCTTTTGACAAATCTGCTTCCCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGATGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAAt                                                <  1:660676/139‑1 (MQ=255)
                                                                                     gCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAAt                                                >  2:77072/1‑139 (MQ=255)
                                                                                         ttGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                               >  2:425259/1‑72 (MQ=255)
                                                                                         ttGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                               <  1:425259/72‑1 (MQ=255)
                                                                                                aTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCcac                                     >  1:916146/1‑139 (MQ=255)
                                                                                                    gCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcg                                 <  2:20539/139‑1 (MQ=255)
                                                                                                      ttACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                     >  1:805923/1‑101 (MQ=255)
                                                                                                      ttACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                     <  2:805923/101‑1 (MQ=255)
                                                                                                            gAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCa                         <  1:730449/139‑1 (MQ=255)
                                                                                                              aaCTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATc                       <  2:472721/139‑1 (MQ=255)
                                                                                                                    cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                           >  2:629322/1‑97 (MQ=255)
                                                                                                                    cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                           <  1:629322/97‑1 (MQ=255)
                                                                                                                                 gTCAGGCGCGATCCCTGATATGCGTAGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTTAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACg    >  2:602233/1‑139 (MQ=255)
                                                                                                                                   cAGGCGCGATCCTTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGaa  <  1:602233/139‑1 (MQ=255)
                                                                                                                                    aGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                            >  2:306138/1‑64 (MQ=255)
                                                                                                                                    aGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                            <  1:306138/64‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CCACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAA  >  NC_000913/1431746‑1432014

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