Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,510,361	C→T	intergenic (+568/+329)	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,361	0	C	T	100.0%	157.1 / NA	51	intergenic (+568/+329)	insO/fecE	IS911B regulator fragment/ferric citrate ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/27);  total (24/27)

AACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGT  >  NC_000913/4510119‑4510600                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 aaCTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAAt                                                                                                                                                                                                                                          <  1:298525/251‑1 (MQ=255)          gTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACtt                                                                                                                                                                                                                                  >  2:322463/1‑250 (MQ=255)            ataACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGt                                                                                                                                                                                                                               >  1:635529/1‑251 (MQ=255)               aCGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCAt                                                                                                                                                                                                                            <  1:263811/251‑1 (MQ=255)                   ggCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTc                                                                                                                                                                                                                         <  2:353043/250‑1 (MQ=255)                    gCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCAc                                                                                                                                                                                                                       <  1:231937/251‑1 (MQ=255)                         cGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATGTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTTGCTCAGAGCGAAATTTTTTGATCTCCCTCCTTAAATCTAAACTGAGTCATGTCACGtttt                                                                                                                                                                                                                  <  1:372991/251‑1 (MQ=255)                                    aaaGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTaa                                                                                                                                                                                                       <  1:446622/251‑1 (MQ=255)                                          gCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGACCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAAttt                                                                                                                                                                                                  <  1:129421/250‑1 (MQ=255)                                             cctacctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCTATGGTCAGAAACTGGATTAGCAATTACGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGCTCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAAc                                                                                                                                                                                              <  2:111732/251‑1 (MQ=255)                                                acctaAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTc                                                                                                                                                                                           <  2:567240/251‑1 (MQ=255)                                                        agtagCAAACGAACACTTTTTAGAAAAATCGATGGTCAGCAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGttttt                                                                                                                                                                                   >  2:658896/1‑251 (MQ=255)                                                                           ttAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGCTGCTTTTGATTTACGTTGGCGTCTGCTCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTtat                                                                                                                                                                 >  1:65738/1‑250 (MQ=255)                                                                               aaaaaTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGa                                                                                                                                                            <  1:541456/251‑1 (MQ=255)                                                                                   aTCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCACAAGACCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATATCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATg                                                                                                                                                        >  2:578151/1‑251 (MQ=255)                                                                                            cAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAg                                                                                                                                               <  2:526838/251‑1 (MQ=255)                                                                                                         aGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTAATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTCATGAACCAGCTCCTTGACTGTCCGCTCTGGCACAGCGCGCAATTTTTTGCTCTCCCTCCTGAAATCTACACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACACACATGTCCGTTAAGCAG‑AGTTCAAAat                                                                                                                                 >  2:266191/1‑251 (MQ=255)                                                                                                              ttCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAattgatt                                                                                                                            >  2:499250/1‑251 (MQ=255)                                                                                                                  ggaccATGGTTGCTTTTGATGTACATTGGCGTCTTATCATTTATTTATCCCCAAAAGACCAATCTCATTGTTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCTCAGAGCGAAATTTTTTGATCTCCCTCCTTAAATCTAAACGTTGGCATGTCCCGTTTTTGTTTTTCTAAAATTTTATCTTCGCATTGTTCGGTGCCGTTATGGTTATACATAAATGTTCGTTAAGCATTAGT‑CAAAATTGATTGCCg                                                                                                                        <  2:51639/248‑1 (MQ=255)                                                                                                                  gTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCg                                                                                                                        <  2:309089/251‑1 (MQ=255)                                                                                                                   ttCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGACCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAAATTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑CGTTCAAAATTGATTGCCGt                                                                                                                       >  2:143315/1‑251 (MQ=255)                                                                                                                      cATGGTTGCTTTTGATTTACGGTGGCGGCTGATCATTGATTTATCCTCAAAAGCCCAACCTAATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCAGAAAGCTAAACTTAGGCATGTAAAGTTTTTGGGTTTCTAAAATTTTAACTGCGCGTTTTTCGTTACCGTAAGGGGTATACAGAAAAGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCGGTGAt                                                                                                                    <  1:266191/251‑1 (MQ=255)                                                                                                                        tGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCCCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑CGTTCAAAATTGATTGCCGTGATc                                                                                                                   >  2:305875/1‑250 (MQ=255)                                                                                                                             gCNTTTGATTTACGTTGGCGTATGATCATTGATTTAGCCTCAAAAACCCAACCTCATGGGTAATGAACCAGCAACGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTGTGATCGCCCTCCTGAAAGCTAAACTTAGTCATGTCACGTTTGTGGGTTTCTAAAATTTTAACGTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAAGGTCCATTAAGCAG‑AGTTCAAAATTGATTGCCGtcgccacgac                                                                                                              <  2:600469/250‑10 (MQ=255)                                                                                                                                tttGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTgg                                                                                                           <  1:349994/250‑1 (MQ=255)                                                                                                                                       tACGTTGGCGTCTGATAATTGATTTATCCTCAAAAGCCCAACCTCATTGGTACTGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGCGCGCAATTTTTTCATCTCCCTCCTGAAATCTAAACTTAGTCATGTCAAGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGGAAGGGTTATACAGAAATGACCGTTAAGCAG‑AGTTCAAAATTGATTGCGGTGATCACGACTGGTTTGaa                                                                                                     >  2:641623/1‑249 (MQ=255)                                                                                                                                                     aTCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCt                                                                                     >  1:68311/1‑251 (MQ=255)                                                                                                                                                            aTTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGc                                                                              <  1:40914/251‑1 (MQ=255)                                                                                                                                                             tttATCCTCAAAAGCCCATCCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGTGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAatt                                                                                                                                <  1:469929/200‑1 (MQ=255)                                                                                                                                                             tttATACTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCATCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAatt                                                                                                                                >  2:469929/1‑200 (MQ=255)                                                                                                                                                                    tCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTATGGCACAGAGCGAACTTTTTTGATCTCCCTCCTGAAATCTAAACTTAGACATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCAAAGCCTGTACAGCTCTGGttt                                                                      >  1:447905/1‑251 (MQ=255)                                                                                                                                                                     cAAAAGCCCAACCTCATTGGGAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGttt                                                                      <  1:349003/250‑1 (MQ=255)                                                                                                                                                                        aaGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTTCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGt                                                                  <  1:26869/251‑1 (MQ=255)                                                                                                                                                                        aaGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGt                                                                  >  1:608646/1‑251 (MQ=255)                                                                                                                                                                          ccccAACCTCATTGGTAATGAACCAGCTCCTTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCCCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGTTCACGACTTGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTg                                                                <  2:300174/250‑1 (MQ=255)                                                                                                                                                                                        gTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGt                                                                                                                                                                          <  1:638053/132‑1 (MQ=255)                                                                                                                                                                                        gTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGt                                                                                                                                                                          >  2:638053/1‑132 (MQ=255)                                                                                                                                                                                               aCCAGCTCCGTGAATGTCCGCTCTGGCACAGCGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGAGTTGATTATGAACCTGTCAGCCTa                                            >  1:561288/1‑250 (MQ=255)                                                                                                                                                                                                 cAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCag                                                                                                                                            >  1:611103/1‑153 (MQ=255)                                                                                                                                                                                                 cAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCag                                                                                                                                            <  2:611103/153‑1 (MQ=255)                                                                                                                                                                                                  atcTCCGTGAATGTCCGCTCGGGCACAGAGAGAAATTTGTTGATCTCCCTCCTTAAATCTAAACTTAGTCATGTCACTTTTTTGGGTGTCTAAAATTTTAACTTCGCGTTGTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTGAAGCAG‑AGTTAAAAATTGATTGCGGTGATCACGACTGGTTTGAAAGCGGCTCCAAAGCCTGTACAGCTCTGGTGTGCGTGGACTATGAGCCTGTCAGCCTaaa                                          <  2:391217/247‑1 (MQ=255)                                                                                                                                                                                                     tCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagca                                       >  1:101348/1‑249 (MQ=255)                                                                                                                                                                                                      ccGTGAAGTTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTATTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAaagcaa                                      <  1:487732/249‑1 (MQ=255)                                                                                                                                                                                                            aTGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATgg                              >  1:192471/1‑251 (MQ=255)                                                                                                                                                                                                                  gCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATg                        >  1:309085/1‑251 (MQ=255)                                                                                                                                                                                                                  gCTCTGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATg                        >  2:650441/1‑251 (MQ=255)                                                                                                                                                                                                                   ctctGGCACAGAGCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATg                        >  1:365261/1‑250 (MQ=255)                                                                                                                                                                                                                             gagCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTaa             <  1:602183/251‑1 (MQ=255)                                                                                                                                                                                                                               gCGAAATTTTTTGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATc           >  2:108315/1‑251 (MQ=255)                                                                                                                                                                                                                                       ttttGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCaga    <  2:309085/250‑1 (MQ=255)                                                                                                                                                                                                                                        tttGATCTCCCTCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGt  >  2:635996/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 AACTGGCCTGTATAACGTGGCTGCCGGTGATTAAGAAAGCTGCACCTACCTAAGTAGTAGCAAACGCACACTTTTTAGAAAAATCGATGGTCAGAAACTGGATTAGCAATTCCGTTCCATGGTTGCTTTTGATTTACGTTGGCGTCTGATCATTGATTTATCCTCAAAAGCCCAACCTCATTGGTAATGAACCAGCTCCGTGAATGTCCGCTCTGGCACAGAGCGAAATTTTTTGATCTCCCCCCTGAAATCTAAACTTAGTCATGTCACGTTTTTGGGTTTCTAAAATTTTAACTTCGCGTTTTTCGTTGCCGTAAGGGTTATACAGAAATGTCCGTTAAGCAG‑AGTTCAAAATTGATTGCCGTGATCACGACTGGTTTGAAAGCCGCGCCCAAGCCTGTACAGCTCTGGTTTGCGTTGATTATGAACCTGTCAGCCTAAAGCAAGCGGATGGACGATGAGTATTGGTAATCTTTCAGAGT  >  NC_000913/4510119‑4510600

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