Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,471,105	G→A	pseudogene (291/291 nt) P138L (CCT→CTT)	tfaS → yfdK ←	pseudogene, CPS‑53 (KpLE1) prophage, tail fiber assembly protein fragment,Phage or Prophage Related CPS‑53 (KpLE1) prophage, conserved protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,471,105	0	G	A	100.0%	124.8 / NA	41	pseudogene (291/291 nt) P138L (CCT→CTT)	tfaS yfdK	pseudogene, CPS‑53 (KpLE1) prophage, tail fiber assembly protein fragment,Phage or Prophage Related CPS‑53 (KpLE1) prophage, conserved protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (20/21);  total (20/21)

GGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACT  >  NC_000913/2470981‑2471238                                                                                                                             |                                                                                                                                      ggCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCt                                                                                                                         >  1:469173/1‑139 (MQ=255)       aGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCa                                                                                                                   >  2:215903/1‑139 (MQ=255)           gttGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               <  2:385902/137‑1 (MQ=255)                 ggcagAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACCGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTAt                                                                                                         >  2:500572/1‑139 (MQ=255)                  gcagAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATc                                                                                                        >  1:59973/1‑139 (MQ=255)                            cccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCa                                                                                                                               <  1:35486/106‑1 (MQ=255)                            cccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCa                                                                                                                               >  2:35486/1‑106 (MQ=255)                             ccGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCa                                                                                             <  2:396412/139‑1 (MQ=255)                              cGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAg                                                                                            <  2:127756/139‑1 (MQ=255)                               gACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGGCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGt                                                                                           <  2:744525/139‑1 (MQ=255)                                          gtgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTgg                                                                                                                       <  2:262292/100‑1 (MQ=18)                                          gtgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTgg                                                                                                                       >  1:262292/1‑100 (MQ=21)                                           tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcg                                                                                                                     <  1:42992/101‑1 (MQ=18)                                           tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGgcg                                                                                                                     >  2:42992/1‑101 (MQ=17)                                           tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                               >  2:383829/1‑139 (MQ=255)                                           tgCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATaa                                                                               >  2:428498/1‑139 (MQ=255)                                                  ttACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  2:699459/139‑1 (MQ=255)                                                     cATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        <  1:442625/136‑1 (MQ=255)                                                     cATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGcc                                                                        >  2:442625/1‑136 (MQ=255)                                                     cATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCcaca                                                                     >  1:58580/1‑139 (MQ=255)                                                          aCGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATta                                                                >  2:629106/1‑139 (MQ=255)                                                                tGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTgc                                                          >  1:539696/1‑139 (MQ=255)                                                                          aGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCt                                                <  2:539696/139‑1 (MQ=255)                                                                            aTACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               <  2:68068/74‑1 (MQ=25)                                                                            aTACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTgg                                                                                                               >  1:68068/1‑74 (MQ=255)                                                                                     ccGCGCCGGATGTCATCTTGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAga                                     <  2:58580/139‑1 (MQ=255)                                                                                      cgcgCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGcc                                                       >  1:86926/1‑120 (MQ=255)                                                                                      cgcgCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGcc                                                       <  2:86926/120‑1 (MQ=255)                                                                                        cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     >  2:413287/1‑72 (MQ=255)                                                                                        cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGAc                                                                                                     <  1:413287/72‑1 (MQ=255)                                                                                        cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                  >  2:384152/1‑139 (MQ=255)                                                                                        cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                  >  1:150134/1‑139 (MQ=255)                                                                                        cgcCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAgacga                                  >  1:330097/1‑139 (MQ=255)                                                                                                            gAACTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGc              <  2:330097/139‑1 (MQ=255)                                                                                                               cTGCCGGAGGAGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCa           <  1:215903/139‑1 (MQ=255)                                                                                                                    ggaggaGTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGtt      <  2:59973/139‑1 (MQ=255)                                                                                                                      aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa    <  2:359670/139‑1 (MQ=255)                                                                                                                      aggagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTa    <  1:384152/139‑1 (MQ=255)                                                                                                                        gagTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACt  <  2:150134/139‑1 (MQ=255)                                                                                                                         agTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACt  >  2:622325/1‑138 (MQ=255)                                                                                                                         agTAAGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACt  <  1:622325/138‑1 (MQ=255)                                                                                                                             |                                                                                                                                      GGCGGAAGCTGATGCAGGCAGAGACCTCCCGACTTAACACTGTGCTGGATTACATTGACGCGGTGACGGCAACAGATACCAGCACCGCGCCGGATGTCATCTGGCCTGAACTGCCGGAGGAGTAGGCCATTCAATATCTGGCGCACTGGAAGTATCGACCAGCTCCAGTGCGTCCAGATAATCCAGCCACAAATTATATTGCGCCAGTTCCTCACCTTTCAGACGACCAATAGCCGCTTTACCAGCCCATTGTTTACT  >  NC_000913/2470981‑2471238

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