Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,486,205	C→A	T29K (ACG→AAG)	crp →	DNA‑binding transcriptional dual regulator CRP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,486,205	0	C	A	100.0%	233.6 / NA	67	T29K (ACG→AAG)	crp	DNA‑binding transcriptional dual regulator CRP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (34/33);  total (34/33)

ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485959‑3486451                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        acataaaTACATTGATGTACTGCATGTAGGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GGGACAGCGGCGTTATCTTGCTCTGGAGAAAGCTTATAACAGAGGAGAACCGCGCATGGTGCTTGGCAAACCCCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCAAGCt                                                                                                                                                                                                                                                      <  2:171193/245‑1 (MQ=255) aCAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCt                                                                                                                                                                                                                                                      <  1:117253/250‑1 (MQ=255)  cAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCtt                                                                                                                                                                                                                                                     >  1:388898/1‑250 (MQ=255)       aTACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAc                                                                                                                                                                                                                                               >  2:178225/1‑251 (MQ=255)          cATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCAcc                                                                                                                                                                                                                                              >  1:530562/1‑249 (MQ=255)          cATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGAAAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGTCAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAg                                                                                                                                                                                                                                            <  2:177295/251‑1 (MQ=255)           aTTGAGGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTTATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCTTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGTCCCGACTCTCGAATGGTTCTTGTCTCATGGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAgg                                                                                                                                                                                                                                           <  1:445479/251‑1 (MQ=255)             tGATGTACTGCTTGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAggg                                                                                                                                                                                                                                          <  2:254877/250‑1 (MQ=255)                 gTACGGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGACGCATATTCCGGCAAGCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAAAGGATAACCGCGCATGGTGCTTGGCAACCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTTCCCATCCAAGAGCAAGCTTATTCACCAGGGTGaa                                                                                                                                                                                                                                      <  1:17926/250‑1 (MQ=255)                       cATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTCGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCgg                                                                                                                                                                                                                               >  1:76343/1‑251 (MQ=255)                           tATGCAAAGGACGTCACATTACCGTGCAGTACAGATGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCACTCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGAGAAAAAGCGGAAAc                                                                                                                                                                                                                           >  1:155158/1‑251 (MQ=255)                            aTGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATACCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAAc                                                                                                                                                                                                                           <  1:623360/250‑1 (MQ=255)                            aTGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACg                                                                                                                                                                                                                          <  1:178225/251‑1 (MQ=255)                             tGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCCCGCCTGGTGCTTGGCAAACCGCAAACAGACCCGACCCTCGCATCGTTCTTGTCCCATTCCCACATTCATACGTACCCATCCAACAGCAAGCTTATTCACCacgctgccacagcggacacgc                                                                                                                                                                                                                         >  2:284089/1‑231 (MQ=255)                                          acaTTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGt                                                                                                                                                                                                            >  1:135324/1‑251 (MQ=255)                                             ttACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTa                                                                                                                                                                                                          >  1:505210/1‑250 (MQ=255)                                                ccGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAgg                                                                                                                                                                                                      >  1:52810/1‑251 (MQ=255)                                                     caggaCAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCACACCGCAAACAGACCCGACTCTAGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAACGCGGACACGCTGTACTACATCGTAAAAGGCTCtg                                                                                                                                                                                                 >  2:60069/5‑251 (MQ=255)                                                      aGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCtgt                                                                                                                                                                                                <  1:124977/251‑1 (MQ=255)                                                           aGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                           <  1:259003/251‑1 (MQ=255)                                                           aGTTGAGAGCCCCTTCCCAGGTAGCGTGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCGGGCTCGGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTGGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCTCCAGGGTGAAAAATCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAg                                                                                                                                                                                           <  1:39864/251‑1 (MQ=255)                                                            gTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGt                                                                                                                                                                                          >  1:439210/1‑251 (MQ=255)                                                                aTAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTg                                                                                                                                                                                      >  2:530123/1‑251 (MQ=255)                                                                  aGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGAt                                                                                                                                                                                    <  1:102244/251‑1 (MQ=255)                                                                  aGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGAt                                                                                                                                                                                    >  1:79141/1‑251 (MQ=255)                                                                   gCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATc                                                                                                                                                                                   <  1:488769/251‑1 (MQ=255)                                                                     cccTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCa                                                                                                                                                                                  >  1:379886/1‑250 (MQ=255)                                                                                  gCGGGAAGCATATTTCGACAATCCA‑GAGACAGCGGCTGTATCTGGCTCGGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGt                                                                                                                                                                    <  2:23965/251‑1 (MQ=255)                                                                                          aatatTTCGGCAATCCA‑GAGACAGCGGAGTGATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑t                                                                                                                                                           <  2:228865/250‑1 (MQ=255)                                                                                           atatTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑Tg                                                                                                                                                          <  1:598022/251‑1 (MQ=255)                                                                                           atatTTCGGCAA‑CCACGAGACAGCGGCGTTATCGGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGAAGGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑Tg                                                                                                                                                          <  2:112668/251‑1 (MQ=255)                                                                                             atTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGAt                                                                                                                                                        <  1:625686/251‑1 (MQ=255)                                                                                                     aaTCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCt                                                                                                                                                >  2:98517/1‑251 (MQ=255)                                                                                                     aaTCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCt                                                                                                                                                >  1:417261/1‑251 (MQ=255)                                                                                                         cA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCACACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGCGCAAGCTTATTCACCAAGGTGAAACAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGCTCACAGACGAAGCGGGTACATAAC‑TGATCCTCTCCTATCt                                                                                                                                            >  2:248772/1‑251 (MQ=255)                                                                                                         cA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCt                                                                                                                                            >  1:518901/1‑251 (MQ=255)                                                                                                          a‑gagACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTg                                                                                                                                           <  1:163262/251‑1 (MQ=255)                                                                                                                cAGCGGCGTTATCGGGCTCGGGAGAAAGTTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCa                                                                                                                                      <  2:381436/251‑1 (MQ=255)                                                                                                                  gcggcgTTATCTGGCTCTGGAGAAAGCTTATATGATAGGATAACCGCGCATGTTGCTTGACAAAGCAGAAACAGAACCGACTCTCGAATGGTTCTTGTCTCAGTGCGACATTCATAAGGACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTTTGGCAGTTCGGATCAAAGACGAAGAGGGTAATTAAA‑TTATCCTCTCCTATCGGAATCAgg                                                                                                                                    <  2:553870/251‑1 (MQ=255)                                                                                                                   gggcgTTATCTGGCTCTGGAGAAAGCATAAAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAggg                                                                                                                                   <  2:520032/250‑1 (MQ=255)                                                                                                                       gTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGAt                                                                                                                               <  1:98517/251‑1 (MQ=255)                                                                                                                                      agaAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTg                                                                                                                <  2:52810/251‑1 (MQ=255)                                                                                                                                             ttATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGACGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGt                                                                                                         >  2:38355/1‑251 (MQ=255)                                                                                                                                              tataCCAGAGGAAACCCGCGCATGGTGCTTCGCACACCGCAACCCGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTCAAGGCTCTGTTGCCGTGCTGATCAAAGACGACGAGGGTAAAGAAA‑TGATCCTCTCCTACCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTg                                                                                                          <  1:2618/249‑1 (MQ=255)                                                                                                                                              tataACAGAGGATAACCGCGTATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGAGCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCAAACTGGGACTGt                                                                                                         <  2:57575/250‑1 (MQ=255)                                                                                                                                               ataAAAGAGGATAACCGAGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGtt                                                                                                        <  2:344907/250‑1 (MQ=255)                                                                                                                                                   cAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCTCATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACAGCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAag                                                                                                   <  2:155158/251‑1 (MQ=255)                                                                                                                                                        gATAGCCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGc                                                                                              <  1:145856/251‑1 (MQ=255)                                                                                                                                                        gATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGCATGGGTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTAACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGATTCAGGGTGATTTTATTGGCGAACGGGGACTGTTTGAAGAGGGc                                                                                              >  2:130142/1‑251 (MQ=255)                                                                                                                                                          tAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGcc                                                                                             >  1:320403/1‑250 (MQ=255)                                                                                                                                                              cgcgCATGGGGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTGCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGaa                                                                                        <  2:561970/251‑1 (MQ=255)                                                                                                                                                               gcgcATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAAc                                                                                       >  1:25929/1‑251 (MQ=255)                                                                                                                                                                       tGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGAAgcgc                                                                                >  1:226126/1‑250 (MQ=255)                                                                                                                                                                                aaCCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACgg                                                                      >  2:503138/1‑250 (MQ=255)                                                                                                                                                                                          aGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGAACGTGCGAAAAc                                                              >  2:630727/1‑249 (MQ=255)                                                                                                                                                                                                gACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCtgtg                                                      >  2:511350/1‑251 (MQ=255)                                                                                                                                                                                                      cGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACACATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGGGCTGGTCAAAGAAGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCTAACTGGGCCTGTTTGTAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAACACCGCCTGTGAAGTgg                                                >  2:291143/1‑251 (MQ=255)                                                                                                                                                                                                        aaTGGTTTTTGTCTCATTGCCACATTCATAAGTAACCATCCAAGAGCAAGCTTGTTCACCAGGGTGAACAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGACGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGCGCTTTTATTGGCGAACTGGTCCTGTTTGAAGAGACCCAGGATCGTAGCGCATGGGGACGTGGGAAAACCGCCTGTGAAGTAGCt                                              >  2:250563/1‑251 (MQ=255)                                                                                                                                                                                                                ttGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTc                                      >  2:474457/1‑251 (MQ=255)                                                                                                                                                                                                                                  aaaaagtccccaccaaaGAGCAAGCTTATTCCCCAGGGGAAAAAAGCGGAAACGCTGTATTACATCGTTAAGGGCTCGGTGGCAGTGCTGATCAAAGACGAAG‑GGGTAAAAAAAATGACCCCCCCCTATCTGAATCAGGTTGATTTTATTGGCGACCTGGCCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTCCGTGCAAAACCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGcc                     <  1:316182/235‑1 (MQ=255)                                                                                                                                                                                                                                     agagccccATCCAAGAGCAAGTTTATTCCCCAGGTTGAAAAAGCGGAAGCGCGGTACTATATCGTTAAAGGCTCTGTGGCAGTTCTGATAAGAGACGAAGTGGGTAAAGAAA‑TTATCCTCCCCTATCGGAACCGGGTTGTTTTTAGTGGGAAACGGGGCCGTTTTAAAGAGGGCTGGGAACGGACCGCATGGGTACGTCCGAAATCCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAAt                  <  1:611119/245‑1 (MQ=255)                                                                                                                                                                                                                                          cccATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGGGGCTGAAATTTCGTACACAAAATTTCGCCAATTGATTc            >  1:380869/1‑251 (MQ=255)                                                                                                                                                                                                                                            cATCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTc            >  2:53030/1‑249 (MQ=255)                                                                                                                                                                                                                                              tCCAAGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTAAggg        >  2:27569/1‑250 (MQ=255)                                                                                                                                                                                                                                                 aaGAGCAAGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTACTATCTGAATCCGGGTGGTTTTATTGGCGAACTGGGCCTGTTTGAAGAGAGCCAGGAACGTAGCGCATGGGTACGTGCGCAAAACGCCTGTAACGGGGCGGAAATTTCGTCAAAAAAATTTCGCCAATTGATTCatggaa      >  2:68096/1‑245 (MQ=255)                                                                                                                                                                                                                                                  agagCAAGCTTATTAAACAGGGAGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGATGATCAAAGACGAAGAGGGTCAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAAGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAAATTCGTACAAAAAATTTCGCAAATTGATTCAGGTAAAc    >  2:514532/1‑251 (MQ=255)                                                                                                                                                                                                                                                    agCAAGCTTATTCACCAGGGGGAAAAGGCGGAAACGCTGTACTACATCGTAAAAGGCTCGGTGGCAGTGCGGATCAAAGACGAAGAGGGTAAAGAAA‑GGATCCTCTCCTATCTAAATCAGGGTGATTTTATTGGGGAACTGGGCCTTTTTGAAGAGGGCAAGGAACGGAGCGCATGGGTACGGGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAAccc  <  2:79141/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        ACAGTAATACATTGATGTACTGCATGTATGCAAAGGACGTCACATTACCGTGCAGTACAGTTGATAGCCCCTTCCCAGGTAGCGGGAAGCATATTTCGGCAATCCA‑GAGACAGCGGCGTTATCTGGCTCTGGAGAAAGCTTATAACAGAGGATAACCGCGCATGGTGCTTGGCAAACCGCAAACAGACCCGACTCTCGAATGGTTCTTGTCTCATTGCCACATTCATAAGTACCCATCCAAGAGCACGCTTATTCACCAGGGTGAAAAAGCGGAAACGCTGTACTACATCGTTAAAGGCTCTGTGGCAGTGCTGATCAAAGACGAAGAGGGTAAAGAAA‑TGATCCTCTCCTATCTGAATCAGGGTGATTTTATTGGCGAACTGGGCCTGTTTGAAGAGGGCCAGGAACGTAGCGCATGGGTACGTGCGAAAACCGCCTGTGAAGTGGCTGAAATTTCGTACAAAAAATTTCGCCAATTGATTCAGGTAAACCC  >  NC_000913/3485959‑3486451

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