Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,238	T→G	intergenic (+445/+452)	insO → / ← fecE	IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,510,238	0	T	G	100.0%	162.7 / NA	50	intergenic (+445/+452)	insO/fecE	IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (22/28);  total (22/28)

TGCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTT  >  NC_000913/4509991‑4510467                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      tGCTTTGCGCTTTCCATCTGGCGATAAAGAGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTCGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGt                                                                                                                                                                                                                                    >  2:490705/1‑251 (MQ=255)     tCGCGCTTTCCCTCTG‑CGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGc                                                                                                                                                                                                                                 <  1:214271/249‑1 (MQ=255)      cgcgCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGTCCCGGCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAATAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACGTGTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCt                                                                                                                                                                                                                                <  2:50722/250‑1 (MQ=255)       gcgcTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCttt                                                                                                                                                                                                                              <  1:138793/251‑1 (MQ=255)               ccTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCTTTCACTTTATGTCAGCTTCGGTAAGCTTGGTGAGACACCGACTCATGAAGATATCAGTGATAACAATCGATCCAGCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCGCCTACCTAAGTATTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACGGGATTAGCAATTCCGTTCCCGGGTTGCTTTTGATTTAc                                                                                                                                                                                                                      <  2:554391/251‑1 (MQ=255)                      cGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACACTCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAACCGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCg                                                                                                                                                                                                               >  1:362428/1‑251 (MQ=255)                           aaGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTg                                                                                                                                                                                                           >  2:97959/1‑250 (MQ=255)                                   aGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCAATGAt                                                                                                                                                                                                  <  1:179645/251‑1 (MQ=255)                                        gggCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATc                                                                                                                                                                                             >  1:451538/1‑251 (MQ=255)                                         ggCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATcc                                                                                                                                                                                            >  1:240438/1‑251 (MQ=255)                                         ggCCCGACTGCCGTCACTGTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATGGATAACAATCGATCCATCAACGGGCCTGTATAACGTGGCTGCCGGTGAGTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTTGCGTCTGATCATTGATTTATcc                                                                                                                                                                                            <  1:409439/251‑1 (MQ=255)                                            ccGTCTGACTTCACATTATGACAGCTTCGGCAAGCTGGGTGAGACACCGACTCATGAAGATATCCTTGATAACAATCGATCCATCAACTGGCCTGTATAACGAGGCTGCCGGTGATTAAGAAAGATGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCCTTGATTTATCCTc                                                                                                                                                                                          >  2:363929/1‑250 (MQ=255)                                             cGTCTGCCTTCACTGTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAAGCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGGTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCaa                                                                                                                                                                                        <  1:299111/251‑1 (MQ=255)                                              gTCTGACTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGACGATATCATTGCTAACCATCGATCCATCAACTGGCCGGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAAAGCACACTTTTTAGAAAAATCGCTGGTCAGAAACTGGATTAGCACTTCCGTTCCAGGGTTGCTTTTGAATTACGTTGGCGTCTGATCATTGCTTTATCCTCaa                                                                                                                                                                                        >  2:241277/1‑250 (MQ=255)                                                       cACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAAc                                                                                                                                                                               <  2:418291/250‑1 (MQ=255)                                                             aTGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGTTCATTGATTTATCCTCAAAAGCCCAACCTCCTTg                                                                                                                                                                        >  1:7361/1‑251 (MQ=255)                                                                ccAGCTCCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTa                                                                                                                                                                     <  1:335229/251‑1 (MQ=255)                                                                           aaGCTGGGTGAGACACGGACTCATGAAGATATCATTGATAAGAAAAGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAACGCTGCACCTACCTAAGAAGTAGCAAACGAACACTTTTTAGACAGATCGATGGTCAGAAACTGGATTAGCAATTCAGTTCCAGGGTTGCTTTTGATTTACGTTGGCGGCTGATCATTGATTTATCCTCAGAAGCCAAACCTCGTTGGTCATGAACCAggt                                                                                                                                                          >  1:209487/1‑249 (MQ=255)                                                                                  gTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGa                                                                                                                                                    >  2:432744/1‑250 (MQ=255)                                                                                       acacCGACTCATGAAGATATCATTGATAAAAATCAAGCCAACAACTGGCCGGTATAACGGGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCGTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCACAAAAGCCCAACCTCATTGGTAATGAACCAGCCCCGTGAATGTcc                                                                                                                                              <  2:378155/251‑1 (MQ=255)                                                                                          ccGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGct                                                                                                                                           >  1:527782/1‑251 (MQ=255)                                                                                             aCTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTg                                                                                                                                        <  2:451538/251‑1 (MQ=255)                                                                                                      gATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGc                                                                                                                               <  1:453597/251‑1 (MQ=255)                                                                                                       atatCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCg                                                                                                                              <  1:432744/251‑1 (MQ=255)                                                                                                         atCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGa                                                                                                                             <  1:225099/250‑1 (MQ=255)                                                                                                                aTAACAATCGATCCATCAACTGGACTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAAtttttt                                                                                                                     >  2:250840/1‑251 (MQ=255)                                                                                                                 tAACACTCGATACATCAACTGGACTGTATAACGTGGCTGACGGAGATTAAGACAGCTGCACCTACCTAAGTAGAAGCACACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGGTTAGCAATTCCGTTCCAGGGATGCTTTTGATTTACGTTGGCGGCTGATCATTGCTTTATCCTCAAACGACCAACCTCATTGGTAATGAACCAGCTCAGTGAATGACCGCTCTTGCACCGAGCGAAATTTTTTg                                                                                                                    >  2:517096/1‑251 (MQ=255)                                                                                                                  aaCAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGa                                                                                                                   <  1:398939/251‑1 (MQ=255)                                                                                                                    caaacGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACATAGCGAAATTTTTTGAtc                                                                                                                 <  1:208076/247‑1 (MQ=255)                                                                                                                           tccatccaCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCCAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCCCAGAGCGAAATTTTTTGATCTcccccc                                                                                                          >  2:539023/1‑251 (MQ=255)                                                                                                                                  aCTGGCCTGTATAACGTGGCTGCCGGTGATTAATAAAGCTGCACCTACCTAAGTAGTAGCAAACGCTCACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATc                                                                                                   <  1:103577/251‑1 (MQ=255)                                                                                                                                      gCCTGTATAACGGGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACAATTTTTAGAAAAATAGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGATACTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCAAACCTCATTGGCAATGAACCAGCTCCGTGAATGACCGCTCTGGCACAGAGCGAAATTTTTTGCTCTCCCCCCTGAAATCTaa                                                                                                >  2:46667/1‑250 (MQ=255)                                                                                                                                       ccTGTATAACGTGGCTGCCGGTGATTAAAAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATGGATTTATCCTCAAAAGCCCAACCGCATTGGTAATGATCCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAAc                                                                                              <  2:361777/251‑1 (MQ=255)                                                                                                                                        cTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCCAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACt                                                                                             >  1:81566/1‑251 (MQ=255)                                                                                                                                           tataACGTGGCTGCCGGAGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACtt                                                                                            >  1:173414/1‑249 (MQ=255)                                                                                                                                                          ggTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCTCAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGttttt                                                                           <  2:145118/251‑1 (MQ=255)                                                                                                                                                                 aaGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTc                                                                    <  1:175421/251‑1 (MQ=255)                                                                                                                                                                            ccctacctaAGTAGTAGCAAATGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCGAGTCCGTTCCAGGGTGGCTTTTTATTTACGTTGGCGGTGGATCAGTGAGTTAGCCTCAAAAGCCCAACCGCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTTCCGCCGGAAATGTAAACGTGGTCATGTCACGTTGTTGGGTTTCGAAAATTTTag                                                         <  2:405345/250‑2 (MQ=255)                                                                                                                                                                                   taagacgtagCAAGCGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGAGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTcgcg                                                  >  2:485777/7‑251 (MQ=255)                                                                                                                                                                                    aaGTAGTAGCAAACGCACACTTTGTAGAAAAATCGATGGTCAGAAACTGGATTATCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCAGTGATTTATCCTCAAAATCCCAACCGCAGTGGTAATGAACCATCTCCGTGAATGTCCGCTCTGGCACAGATCGAAATTTTTTGATCTCCCCCCTGAAATCTATACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATGTTAACGTCGCGt                                                 <  1:578065/251‑1 (MQ=255)                                                                                                                                                                                       tagtagCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTACACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGtttt                                              >  1:283324/1‑251 (MQ=255)                                                                                                                                                                                            gCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCCTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAAGCTAAACTTAGTCATGGCGCGTTTTTGGGTTTCTAAAATTTTAACGTCGCGTTTTTCGtt                                         <  2:241637/251‑1 (MQ=255)                                                                                                                                                                                               aaCGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGcc                                      >  1:168667/1‑251 (MQ=255)                                                                                                                                                                                                         ttttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                             >  2:458217/1‑250 (MQ=255)                                                                                                                                                                                                         ttttAGAAAAAGCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCTCGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGtt                             <  1:424786/250‑1 (MQ=255)                                                                                                                                                                                                          tttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCGCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTta                            <  1:292074/250‑1 (MQ=255)                                                                                                                                                                                                                      gATGGTCAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTcc               <  1:458217/251‑1 (MQ=255)                                                                                                                                                                                                                       atgttcataaaCTGGATTAGCAATTCCGTTCCAGGCTTGCTCTGGATTTCCGTTGGCGTCTGATCATTGCTTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCTGCTCCGTGAATGACCGCTCTGGCACAGAGCGAAAGTGTTTGATCTCCCCCCTTAAATCTAAACTTAGTTATGTCTCGTTTTTGGGTTTCTAAAATTTTCACTTCGCTTTTTTCGTTTCCGTAAGAGTTATTCAGAATTGTc                <  1:329686/241‑1 (MQ=255)                                                                                                                                                                                                                            cAGAAACTGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACATAAATGTCCGTTaa          <  1:469954/250‑1 (MQ=255)                                                                                                                                                                                                                                   tGGATTAGCAATTCCGTTCCAGGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGtt  >  2:308505/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                      TGCTTCGCGCTTTCCCTCTGGCGAAAAAGGGAGAAGGACGGGCCCGTCTGCCTTCACTTTATGCCAGCTTCGGTAAGCTGGGTGAGACACCGACTCATGAAGATATCATTGATAACAATCGATCCATCAACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAGAGTT  >  NC_000913/4509991‑4510467

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