Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,188	A→G	intergenic (+55/+502)	yjhV → / ← fecE	pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,188	0	A	G	100.0%	177.3 / NA	49	intergenic (+55/+502)	yjhV/fecE	pseudogene, KpLE2 phage‑like element/ferric citrate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/22);  total (27/22)

GCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCT  >  NC_000913/4510067‑4510315                                                                                                                          |                                                                                                                                gCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                >  2:158163/1‑139 (MQ=255)     ggTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAGGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGt                                                                                                            <  2:198664/139‑1 (MQ=255)     ggTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGt                                                                                                            >  1:402245/1‑139 (MQ=255)                     aTGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGc                                                                                            >  2:91249/1‑139 (MQ=255)                         aGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                        <  2:402245/139‑1 (MQ=255)                                  tGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                <  2:757227/106‑1 (MQ=255)                                  tGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                >  1:757227/1‑106 (MQ=255)                                        cAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTc                                                                                                           <  2:587611/105‑1 (MQ=255)                                        cAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTc                                                                                                           >  1:587611/1‑105 (MQ=255)                                         aaTCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                        <  1:158163/139‑1 (MQ=255)                                                        tGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTgg                                                         <  1:107284/139‑1 (MQ=255)                                                          gCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCg                                                       >  2:165695/1‑139 (MQ=255)                                                           ccTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                      <  2:58882/139‑1 (MQ=255)                                                             tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                             >  1:558494/1‑98 (MQ=255)                                                             tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                             <  2:558494/98‑1 (MQ=255)                                                             tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                    >  2:40053/1‑139 (MQ=255)                                                                    cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                        >  1:680148/1‑112 (MQ=255)                                                                    cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                        <  2:680148/112‑1 (MQ=255)                                                                    cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                             <  2:387462/123‑1 (MQ=255)                                                                    cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                             >  1:387462/1‑123 (MQ=255)                                                                    cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCAtt                                             >  1:232321/1‑139 (MQ=255)                                                                       ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                          <  2:232321/139‑1 (MQ=255)                                                                       ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                          <  1:91249/139‑1 (MQ=255)                                                                       ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                          >  1:334610/1‑139 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                      >  2:385146/1‑125 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                      <  1:385146/125‑1 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                        >  2:80997/1‑139 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                        >  2:788201/1‑139 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                        >  2:442960/1‑139 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                        >  1:45710/1‑139 (MQ=255)                                                                         cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                        >  2:63204/1‑139 (MQ=255)                                                                                       gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                          <  1:788201/139‑1 (MQ=255)                                                                                            cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCt                     >  2:24125/1‑139 (MQ=255)                                                                                               cACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAATTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGAtt                                                                 <  1:687054/92‑1 (MQ=255)                                                                                               cACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAATTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGAtt                                                                 >  2:687054/1‑92 (MQ=255)                                                                                                    acctaAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGtt                                                                                   <  2:618093/69‑1 (MQ=255)                                                                                                    acctaAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGtt                                                                                   >  1:618093/1‑69 (MQ=255)                                                                                                               agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                          <  1:9488/99‑1 (MQ=255)                                                                                                               agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                          >  2:9488/1‑99 (MQ=255)                                                                                                               agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCa                         <  2:69845/116‑1 (MQ=255)                                                                                                               agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCa                         >  1:69845/1‑116 (MQ=255)                                                                                                               agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCt  >  2:775615/1‑139 (MQ=255)                                                                                                               agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCt  <  2:297923/139‑1 (MQ=255)                                                                                                                   aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAAcc      <  2:707657/131‑1 (MQ=255)                                                                                                                   aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAAcc      >  1:707657/1‑131 (MQ=255)                                                                                                                     cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                      >  2:743055/1‑81 (MQ=255)                                                                                                                     cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                      <  1:743055/81‑1 (MQ=255)                                                                                                                     cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                 >  2:635600/1‑102 (MQ=255)                                                                                                                     cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                 <  1:635600/102‑1 (MQ=255)                                                                                                                          |                                                                                                                                GCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCT  >  NC_000913/4510067‑4510315

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