Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,514,287	G→C	L8V (CTT→GTT)	fecB ←	ferric citrate ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,514,287	0	G	C	100.0%	177.1 / NA	56	L8V (CTT→GTT)	fecB	ferric citrate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (27/29);  total (27/29)

GCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGT  >  NC_000913/4514041‑4514521                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            gCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACaaaa                                                                                                                                                                                                                                         <  1:25671/251‑1 (MQ=255)    gcgcgcgcGTTCCGACGGACTACCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGGTTTGCATCGGTATCGGCGGCAATCCCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACgg                                                                                                                                                                                                                                      <  2:404468/251‑1 (MQ=255)       ggcgcgTTCCGCCGGACTGCCACGGTTTCAGGTGCGCACGCCCTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATa                                                                                                                                                                                                                                   <  2:107089/250‑1 (MQ=255)            gTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCACAGGCGTGGCTGATCACCAGCAGCACGCCTGCAAAAACAAAACGGATAAATGc                                                                                                                                                                                                                              >  2:326556/1‑251 (MQ=255)                cGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGTATAAATGCCAAc                                                                                                                                                                                                                          <  1:111608/251‑1 (MQ=255)                     gACTGCCACGGTTTCAGGTGCGCACGCACTTCTGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCCGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGGGTGTAAACGTGCCGTATTCGTCCTGAAACGTGCCGGCCAAGGCGTGGCTGTTCACCAGCAGCAGGCCTGCAAAAACACAACGGGTAAAAGCCAACA‑TAAt                                                                                                                                                                                                                    >  2:223017/1‑251 (MQ=255)                         gCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATcac                                                                                                                                                                                                                 <  1:326556/250‑1 (MQ=255)                             cgggttCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGGGTTTTTTCGAGTGTAAACGTGCCGGGTTCGTCCTGAACCGGGGCGGCAAAGGCGAGGCTTATCACCAGCAGCACGCCTGCACAAACAAACAGGATAAATGCCaaaaaaaacccctacc                                                                                                                                                                                                             >  2:360286/5‑237 (MQ=255)                                tttCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGc                                                                                                                                                                                                         >  1:14090/1‑251 (MQ=255)                                      gTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTaaa                                                                                                                                                                                                     >  1:373083/1‑249 (MQ=255)                                          gcacgcacTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCg                                                                                                                                                                                               <  2:421829/251‑1 (MQ=255)                                                acTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGc                                                                                                                                                                                         >  1:385302/1‑251 (MQ=255)                                                     gggCAGGATGCGTTTAGCATCGGTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGATCCAGCACCACAATCCGTTGTGGAGTTTTTTCGAGTGTAAACGGGCCGTGTTCGTCCTGAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCAAGCTAAAAGCCCGGCAAGACGGGc                                                                                                                                                                                    >  1:470310/1‑251 (MQ=255)                                                            aTGCGGTTTGCATCGTTCTCGTCGGCAATACCGATCGCGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑GAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAaca                                                                                                                                                                             <  2:4252/251‑1 (MQ=255)                                                               cGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCAAAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGACGGCCTGCAAAAACCAAACGGAAAACTGCCA‑CACTAATCACATTCCAGCTAAAAACCCGGCAAGCCGGGCGTTAACACAt                                                                                                                                                                          >  2:392923/1‑251 (MQ=255)                                                                gTTTGGCATCGTGATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTGCAGCACCACAATCCGTTGTTGCGGTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGTCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGGCAACA‑TATTCTCATTCCAGCTAAAAGCCCGGCAAGCCGTGCGTTAACACATc                                                                                                                                                                         <  1:277983/251‑1 (MQ=255)                                                                  tttGCATCGTTATCGTCGGCAATACCGATCGGGCTGCCGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGGTCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑AAATACCCTTCCAGCTAAAAGCCCGACAAGCCGGGCGTTAAAACATc                                                                                                                                                                         >  2:170758/1‑249 (MQ=255)                                                                       aTCGTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGC‑AACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAACGGCGTGGCTGATCACCAGCATCAGGCCTGCAAAAACAAAACGGATAAATGCCAAAA‑TAATCACATTCCAACTAAAAGCCCGGCAAGCCGGGCGTTAACACATCCAAACTTc                                                                                                                                                                 >  2:73961/1‑251 (MQ=255)                                                                          gTTATCGTCGGCAATACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCa                                                                                                                                                                >  2:39436/1‑250 (MQ=255)                                                                                      aaTACCGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGcct                                                                                                                                                   >  2:255349/1‑249 (MQ=255)                                                                                          ccGATCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATAc                                                                                                                                               <  1:255349/251‑1 (MQ=255)                                                                                             aTCGGGCTGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGc                                                                                                                                            <  2:328347/251‑1 (MQ=255)                                                                                                          ccACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCCCCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGc                                                                                                                               <  1:321396/251‑1 (MQ=255)                                                                                                            ccggcggCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCt                                                                                                                             <  1:364948/250‑1 (MQ=255)                                                                                                                   ccAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACACCGTGCGCGGCTGGCCTGCATaa                                                                                                                       >  2:17865/1‑249 (MQ=255)                                                                                                                        gcATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAAAACATCAGAAATTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCt                                                                                                                 >  1:35501/1‑251 (MQ=255)                                                                                                                                gAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCCTATACAGCGTGCGCGGCTGGCCTGCATAGATTCCTttgttgtt                                                                                                         >  2:441045/1‑251 (MQ=255)                                                                                                                                gAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTttgttgtt                                                                                                         <  1:329202/251‑1 (MQ=255)                                                                                                                                       agTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAAAACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTgtcgtc                                                                                                   >  2:62756/1‑250 (MQ=255)                                                                                                                                             aGCACCCCAATCCGTTGTGGCTTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAagag                                                                                            <  1:280704/251‑1 (MQ=255)                                                                                                                                                     aaTCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATg                                                                                      >  1:155720/1‑249 (MQ=255)                                                                                                                                                             gtgGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCt                                                                            <  1:349463/251‑1 (MQ=255)                                                                                                                                                             gtgGCGTTTTTGCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCt                                                                            <  2:222468/251‑1 (MQ=255)                                                                                                                                                                  gTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCg                                                                       >  2:177503/1‑251 (MQ=255)                                                                                                                                                                         cGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATg                                                                 <  1:135350/250‑1 (MQ=255)                                                                                                                                                                            gtgtAAACGTGCCGTGGTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACGTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGtttt                                                             <  2:407285/251‑1 (MQ=255)                                                                                                                                                                                caaCGTGCCGTGTTCGGCCTGAACCGCGGCGGCAACAGCGGGCCTGATCACCAGCAGCAAGCCTGCAAAAACAAAACGGATAAATGCAAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAGCGACCCCTGCATATACAGCGTGCGCGGCTAGCCTGCATCGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGCAGTCATGGTCGAAGATGTTTTTcac                                                         >  1:179721/2‑251 (MQ=255)                                                                                                                                                                                aaaCGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTcac                                                         >  1:179727/1‑251 (MQ=255)                                                                                                                                                                                     tGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCg                                                     >  1:143031/1‑250 (MQ=255)                                                                                                                                                                                        cGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATg                                                 <  2:155720/251‑1 (MQ=255)                                                                                                                                                                                           gtTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCa                                              <  2:143031/251‑1 (MQ=255)                                                                                                                                                                                               gTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGTCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGtt                                          <  1:436280/251‑1 (MQ=255)                                                                                                                                                                                               gTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGtt                                          <  1:155788/251‑1 (MQ=255)                                                                                                                                                                                                 cccgAACCGGGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCa                                        <  2:19416/248‑1 (MQ=255)                                                                                                                                                                                                 ccTGAACCGGGGCGGCCAAGGTGCTGTTGATCACCACCAGCAGGTCTGCAAAAACACAACGGATAAATGCCAACA‑TAATCACATACGAGCTAAAAGCCCGGCAAGTGCGGCGTTAACACATCAGAACTTCTTCGGTCCCTGCATATACAGCGTACGCGGCTGGCCTGCATATATGCCTTTGTTGTTTGCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGGTGTTTTTCTCACCGAATGCTAGGTTc                                         <  2:429626/250‑1 (MQ=255)                                                                                                                                                                                                    gAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTc                                         >  1:414925/1‑247 (MQ=255)                                                                                                                                                                                                      aCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCt                                   >  1:413913/1‑251 (MQ=255)                                                                                                                                                                                                      aCCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCt                                   <  1:39436/251‑1 (MQ=255)                                                                                                                                                                                                          gggcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCa                               <  2:277189/250‑1 (MQ=255)                                                                                                                                                                                                            gcggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATc                             >  1:406469/1‑251 (MQ=255)                                                                                                                                                                                                              ggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGCGGAAGTAGACCTGGTCGAAGATGTTTTTCACACCGAATGCCGGGTTCAGATCTGCCCTCTg                           >  2:19763/1‑251 (MQ=255)                                                                                                                                                                                                              ggcAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTg                           <  1:73961/251‑1 (MQ=255)                                                                                                                                                                                                                               tCACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACg          >  1:425993/1‑251 (MQ=255)                                                                                                                                                                                                                                cACCAGCAGCAGGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACg          <  1:62756/250‑1 (MQ=255)                                                                                                                                                                                                                                       agcagGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGt  <  2:373083/251‑1 (MQ=255)                                                                                                                                                                                                                                       agcagGCCTGCAAAAACAAAACGGATAAATGCCAACA‑TAATCACATTACAGCTAAAAGCCCGACAAGCCGGGCGATAACACATCAGAACTTCAACGACCCCAGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGCAGTAGTCCTGGTCGAAGATGTTTTTCACACAGCATGCCAAGATAAAATCTGCCATCTGCGGGCCTAAGTCCTACGcccagcct  >  2:250293/1‑245 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            GCTGCGCGCGCGTTCCGACGGACTGCCACGGTTTCAGGTGCGCACGCACTTCGGGCAGGATGCGTTTTGCATCGTTATCGTCGGCAATACCGATCGGGATGACGTCCACGGCGGCCAGCGCATCGGCGAACGAGAGTTCCAGCACCACAATCCGTTGTGGCGTTTTTTCGAGTGTAAACGTGCCGTGTTCGTCCTGAACCGTGGCGGCAAAGGCGTGGCTGATCACCAGCAGCAGGCCTGCAAAAAGAAAACGGATAAATGCCAACA‑TAATCACATTCCAGCTAAAAGCCCGGCAAGCCGGGCGTTAACACATCAGAACTTCAACGACCCCTGCATATACAGCGTGCGCGGCTGGCCTGCATAGATGCCTTTGTTGTTGTCGTCATAAGAGCGGATGAAGTAGTCCTGGTCGAAGATGTTTTTCACACCGAATGCCAGGTTCAGATCTGCCATCTGCGGGCCAAAGTCATACGCCACGCGT  >  NC_000913/4514041‑4514521

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