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%	251.9 / NA	69	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 (35/34);  total (35/34)

GCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGT  >  NC_000913/4510057‑4510319                                                                                                                                    |                                                                                                                                    gCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTa                                                                                                                              >  2:371201/1‑139 (MQ=255)     cGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGaaa                                                                                                                          >  2:535816/1‑139 (MQ=255)                      cacCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGTTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAAc                                                                                                         >  2:955950/1‑139 (MQ=255)                         cGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  1:1153323/139‑1 (MQ=255)                         cGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  2:1314553/139‑1 (MQ=255)                           aCTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                           <  1:1116518/132‑1 (MQ=255)                           aCTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                           >  2:1116518/1‑132 (MQ=255)                           aCTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGAt                                                                                                    >  2:1466250/1‑139 (MQ=255)                                   aGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            <  1:535816/139‑1 (MQ=255)                                                     tCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      >  2:398681/1‑111 (MQ=255)                                                     tCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  1:398681/111‑1 (MQ=255)                                                     tCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCaa                                                                                              <  1:614772/119‑1 (MQ=255)                                                     tCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCaa                                                                                              >  2:614772/1‑119 (MQ=255)                                                        aTCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGt                                                                                        <  2:1635198/122‑1 (MQ=255)                                                        aTCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGt                                                                                        >  1:1635198/1‑122 (MQ=255)                                                         tCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGAt                                                                      >  1:1636716/1‑139 (MQ=255)                                                          ccATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGAtt                                                                     <  2:1347866/139‑1 (MQ=255)                                                              cAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 <  1:371201/139‑1 (MQ=255)                                                                aCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               <  1:1466250/139‑1 (MQ=255)                                                                    gCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                    >  2:1447024/1‑82 (MQ=255)                                                                    gCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                    <  1:1447024/82‑1 (MQ=255)                                                                     ccTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                          <  2:336276/139‑1 (MQ=255)                                                                           taACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTc                                                                                           <  1:1243092/100‑1 (MQ=255)                                                                           taACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTc                                                                                           >  2:1243092/1‑100 (MQ=255)                                                                           taACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGC‑‑TTGATTTACGTTGGCGTCTGATCa                                                   >  1:1366582/1‑138 (MQ=255)                                                                           taACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGC‑‑TTGATTTACGTTGGCGTCTGATCa                                                   <  2:1366582/138‑1 (MQ=255)                                                                              cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCg                                                                                                                     >  2:798750/1‑71 (MQ=255)                                                                              cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCg                                                                                                                     <  1:798750/71‑1 (MQ=255)                                                                              cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        >  1:611475/1‑132 (MQ=255)                                                                              cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        <  2:611475/132‑1 (MQ=255)                                                                                 ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaaa                                                                                                          <  2:774989/79‑1 (MQ=255)                                                                                 ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaaa                                                                                                          >  1:774989/1‑79 (MQ=255)                                                                                 ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTACAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                              <  2:1051428/139‑1 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCg                                                                                         >  2:518525/1‑94 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCg                                                                                         <  1:518525/94‑1 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                              >  1:661707/1‑137 (MQ=255)                                                                                   cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                              <  2:661707/137‑1 (MQ=255)                                                                                     gCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTAt                                          <  2:1636716/139‑1 (MQ=255)                                                                                      ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCa                                                   <  1:654477/129‑1 (MQ=255)                                                                                      ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCa                                                   >  2:654477/1‑129 (MQ=255)                                                                                                      cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            >  1:192526/1‑72 (MQ=255)                                                                                                      cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            <  2:192526/72‑1 (MQ=255)                                                                                                      cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGc                                >  2:168497/1‑132 (MQ=255)                                                                                                      cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGc                                <  1:168497/132‑1 (MQ=255)                                                                                                         cACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCttt                                                                          >  2:1228737/1‑87 (MQ=255)                                                                                                         cACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCttt                                                                          <  1:1228737/87‑1 (MQ=255)                                                                                                         cACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAt                      >  1:1483161/1‑139 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            >  1:1690459/1‑54 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            <  2:1690459/54‑1 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               <  1:1130892/83‑1 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               >  2:1130892/1‑83 (MQ=255)                                                                                                                        tagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGc       >  2:401466/1‑139 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      >  1:579501/1‑43 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  2:579501/43‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                          >  2:844223/1‑87 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                          <  1:844223/87‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        <  2:695636/89‑1 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        >  1:695636/1‑89 (MQ=255)                                                                                                                         agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCt      >  2:1115916/1‑139 (MQ=255)                                                                                                                          gCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGt                  <  2:1413476/126‑1 (MQ=255)                                                                                                                          gCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGt                  >  1:1413476/1‑126 (MQ=255)                                                                                                                          gCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTc     >  1:1019960/1‑139 (MQ=255)                                                                                                                             aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 >  1:661778/1‑76 (MQ=255)                                                                                                                             aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 <  2:661778/76‑1 (MQ=255)                                                                                                                             aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGt  >  2:978046/1‑139 (MQ=255)                                                                                                                              acgcacgcTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               <  1:707215/77‑1 (MQ=255)                                                                                                                              acgcacgcTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               >  2:707215/1‑77 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        >  1:346959/1‑83 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        <  2:346959/83‑1 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                     >  1:1602909/1‑102 (MQ=255)                                                                                                                               cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCa                                     <  2:1602909/102‑1 (MQ=255)                                                                                                                                    |                                                                                                                                    GCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGT  >  NC_000913/4510057‑4510319

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