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% 196.7 / NA 55V945A (GTG→GCG) stfRRac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/28);  total (27/28)

CAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGAC  >  NC_000913/1431750‑1432011
                                                                                                                                    |                                                                                                                                 
caGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAACCTTGCAGCCGCTTATCCGTCAGGCGCGATccc                                                                                                                             <  1:916327/139‑1 (MQ=255)
                 ggCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                            <  1:1538770/139‑1 (MQ=255)
                         aaTCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAg                                                                                                    >  2:1598662/1‑139 (MQ=255)
                                     aTCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGcc                                                                                        >  1:385572/1‑139 (MQ=255)
                                               gTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                              <  2:943106/139‑1 (MQ=255)
                                                ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGc                                                                             >  1:1483796/1‑139 (MQ=255)
                                                ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGc                                                                             <  1:917134/139‑1 (MQ=255)
                                                    gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                         >  2:880456/1‑139 (MQ=255)
                                                    gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                         >  2:498824/1‑139 (MQ=255)
                                                     tCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGat                                                                                                                         >  1:1463384/1‑90 (MQ=255)
                                                     tCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGat                                                                                                                         <  2:1463384/90‑1 (MQ=255)
                                                     tCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTAt                                                                        >  2:1533138/1‑139 (MQ=255)
                                                        ggTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGt                                                                     >  1:1386568/1‑139 (MQ=255)
                                                        ggTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGt                                                                     >  2:485680/1‑139 (MQ=255)
                                                                     tGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                        <  2:806392/139‑1 (MQ=255)
                                                                         ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGAtt                                                                                                       <  2:926553/88‑1 (MQ=255)
                                                                         ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGAtt                                                                                                       >  1:926553/1‑88 (MQ=255)
                                                                            gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTgg                                                                                                                 <  1:1162715/75‑1 (MQ=255)
                                                                            gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTgg                                                                                                                 >  2:1162715/1‑75 (MQ=255)
                                                                            gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCAtt                                                 >  1:1136204/1‑139 (MQ=255)
                                                                              aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATc                                                                                                                               <  2:27977/59‑1 (MQ=255)
                                                                              aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATc                                                                                                                               >  1:27977/1‑59 (MQ=255)
                                                                              aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                      >  2:277490/1‑68 (MQ=255)
                                                                              aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                      <  1:277490/68‑1 (MQ=255)
                                                                              aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                        >  1:1087834/1‑130 (MQ=255)
                                                                              aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                        <  2:1087834/130‑1 (MQ=255)
                                                                               ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                         <  2:524650/112‑1 (MQ=255)
                                                                               ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                         >  1:524650/1‑112 (MQ=255)
                                                                               ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTaaa                                              >  1:398697/1‑139 (MQ=255)
                                                                                gCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAAt                                             >  1:1551514/1‑139 (MQ=255)
                                                                                    ttGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTc                                                                                 >  2:809405/1‑99 (MQ=255)
                                                                                    ttGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTc                                                                                 <  1:809405/99‑1 (MQ=255)
                                                                                      gACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATa                                       <  2:398697/139‑1 (MQ=255)
                                                                                            tCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGAt                                                                                                                                >  1:658696/1‑44 (MQ=255)
                                                                                            tCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGAt                                                                                                                                <  2:658696/44‑1 (MQ=255)
                                                                                               gCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcg                              <  1:495428/139‑1 (MQ=255)
                                                                                                 ttACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgcc                            <  2:385572/139‑1 (MQ=255)
                                                                                                    cccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagc                         <  1:498824/139‑1 (MQ=255)
                                                                                                     ccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                  <  1:1291058/97‑1 (MQ=255)
                                                                                                     ccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                  >  2:1291058/1‑97 (MQ=255)
                                                                                                      cGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagcgc                       <  2:1136204/139‑1 (MQ=255)
                                                                                                              gCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTa               <  2:1551514/139‑1 (MQ=255)
                                                                                                               cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                      <  2:339308/35‑1 (MQ=17)
                                                                                                               cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATg                                                                                                                      >  1:339308/1‑35 (MQ=17)
                                                                                                               cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgcca                           >  2:1273362/1‑126 (MQ=255)
                                                                                                               cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgcca                           <  1:1273362/126‑1 (MQ=255)
                                                                                                                  ccGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGa           >  2:1256209/1‑139 (MQ=255)
                                                                                                                   cGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGt                                                                          <  1:259054/75‑1 (MQ=255)
                                                                                                                   cGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGt                                                                          >  2:259054/1‑75 (MQ=255)
                                                                                                                     cTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAg            <  1:774808/135‑1 (MQ=255)
                                                                                                                     cTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAg            >  2:774808/1‑135 (MQ=255)
                                                                                                                     cTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGAttt        >  2:254313/1‑139 (MQ=255)
                                                                                                                         tCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTgggg    <  1:880456/139‑1 (MQ=255)
                                                                                                                         tCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTgggg    <  1:1256209/139‑1 (MQ=255)
                                                                                                                           cGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGAc  <  2:108479/139‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
CAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGAC  >  NC_000913/1431750‑1432011

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: ‑