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%	291.0 / NA	82	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 (41/41);  total (41/41)

CTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGT  >  NC_000913/4510058‑4510319                                                                                                                                   |                                                                                                                                    cTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCCGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAg                                                                                                                             >  2:634422/1‑139 (MQ=255)     ggTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGaaaa                                                                                                                         >  2:378541/1‑139 (MQ=255)                tGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCa                                                                                                              >  2:1170733/1‑139 (MQ=255)                    acacCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaaa                                                                                                          <  2:1439995/139‑1 (MQ=255)                     cacCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGaaaa                                                                                                                         <  1:1146044/123‑1 (MQ=255)                     cacCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGaaaa                                                                                                                         >  2:1146044/1‑123 (MQ=255)                        cGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      >  1:661405/1‑139 (MQ=255)                                          ttGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                    <  1:870370/107‑1 (MQ=255)                                          ttGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGa                                                                                                                    >  2:870370/1‑107 (MQ=255)                                             aTAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATgg                                                                                 >  2:416176/1‑139 (MQ=255)                                                   aTCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  1:193084/112‑1 (MQ=255)                                                   aTCGATCCATCAACTGGCCTGTATAACGCGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      >  2:193084/1‑112 (MQ=255)                                                            tCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAG‑AGCAAACGCACGCTTTTTAGAAAAATCGATGGt                                                                                                                <  2:1196219/92‑1 (MQ=255)                                                            tCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAG‑AGCAAACGCACGCTTTTTAGAAAAATCGATGGt                                                                                                                >  1:1196219/1‑92 (MQ=255)                                                               aCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               <  1:416176/139‑1 (MQ=255)                                                                  ggCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGc                                                            >  1:716129/1‑139 (MQ=255)                                                                    ccTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                          <  2:946481/139‑1 (MQ=255)                                                                      tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                          <  1:924582/137‑1 (MQ=255)                                                                      tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGt                                                          >  2:924582/1‑137 (MQ=255)                                                                      tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        <  2:716129/139‑1 (MQ=255)                                                                      tGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCt                                                        <  2:371905/139‑1 (MQ=255)                                                                           aaCGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCa                                                   <  2:661405/139‑1 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      >  2:1024779/1‑86 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  1:1024779/86‑1 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                                 <  2:879702/91‑1 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                                 >  1:879702/1‑91 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               <  2:1176011/125‑1 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGtt                                                               >  1:1176011/1‑125 (MQ=255)                                                                             cGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCAtt                                                 >  2:1389652/1‑139 (MQ=255)                                                                              gTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGCAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                            >  1:900714/1‑111 (MQ=255)                                                                              gTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAACTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                            <  2:900714/111‑1 (MQ=255)                                                                                ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                              <  1:1389652/139‑1 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGaaa                                                                                                                          >  1:1353906/1‑61 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGaaa                                                                                                                          <  2:1353906/61‑1 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCa                                                                                                              <  2:175896/73‑1 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCa                                                                                                              >  1:175896/1‑73 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            <  2:626495/91‑1 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            >  1:626495/1‑91 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                            >  1:775127/1‑139 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                            >  2:1457183/1‑139 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                            >  1:1401707/1‑139 (MQ=255)                                                                                  cTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAttt                                            >  1:195073/1‑139 (MQ=255)                                                                                     ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 >  2:989420/1‑115 (MQ=255)                                                                                     ccGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 <  1:989420/115‑1 (MQ=255)                                                                                             gaaGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAg                                 <  1:378541/138‑1 (MQ=255)                                                                                                gAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGccc                              <  1:957905/139‑1 (MQ=255)                                                                                                 aaaGCTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCa                             <  2:300335/139‑1 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      <  1:1378070/62‑1 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTgg                                                                                                      >  2:1378070/1‑62 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                            <  2:735116/88‑1 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCt                                                                            >  1:735116/1‑88 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 <  1:1030112/99‑1 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACg                                                                 >  2:1030112/1‑99 (MQ=255)                                                                                                     cTGCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCt                         >  1:577462/1‑139 (MQ=255)                                                                                                       gCACCTACCTAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCa                       <  2:597722/139‑1 (MQ=255)                                                                                                         acctacctAAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCAtt                     <  1:634422/139‑1 (MQ=255)                                                                                                               ctaAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                           <  2:154024/47‑1 (MQ=255)                                                                                                               ctaAGTAGTAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGaa                                                                                                           >  1:154024/1‑47 (MQ=255)                                                                                                                  agtagtAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGaa            <  2:1401707/139‑1 (MQ=255)                                                                                                                  agtagtAGCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGaa            >  1:1477180/1‑139 (MQ=255)                                                                                                                    tagtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAAcc          <  1:1457183/139‑1 (MQ=255)                                                                                                                      gtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                                 >  1:1471199/1‑50 (MQ=255)                                                                                                                      gtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAg                                                                                                 <  2:1471199/50‑1 (MQ=255)                                                                                                                      gtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGt                                                                <  1:560779/83‑1 (MQ=255)                                                                                                                      gtagCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGt                                                                >  2:560779/1‑83 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGTTGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                           <  1:1173765/118‑1 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGTTGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                           >  2:1173765/1‑118 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTc                                                                                      >  2:847013/1‑59 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTc                                                                                      <  1:847013/59‑1 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCtt                                                                           <  2:1212688/70‑1 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCtt                                                                           >  1:1212688/1‑70 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCt      <  2:195073/139‑1 (MQ=255)                                                                                                                        agCAAACGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCt      >  2:1230403/1‑139 (MQ=255)                                                                                                                            aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCtt                                                                           >  1:758456/1‑66 (MQ=255)                                                                                                                            aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCtt                                                                           <  2:758456/66‑1 (MQ=255)                                                                                                                            aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                              <  2:593853/95‑1 (MQ=255)                                                                                                                            aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGAt                                              >  1:593853/1‑95 (MQ=255)                                                                                                                            aaCGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGt  >  2:1382727/1‑139 (MQ=255)                                                                                                                              cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            >  1:1215879/1‑47 (MQ=255)                                                                                                                              cGCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAAtt                                                                                            <  2:1215879/47‑1 (MQ=255)                                                                                                                               gCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCa                                                   >  1:1164386/1‑87 (MQ=255)                                                                                                                               gCACGCTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCa                                                   <  2:1164386/87‑1 (MQ=255)                                                                                                                                   |                                                                                                                                    CTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGT  >  NC_000913/4510058‑4510319

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