Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,380,469	C→T	E17K (GAA→AAA)	frdB ←	fumarate reductase iron‑sulfur protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,380,469	0	C	T	100.0%	148.0 / NA	42	E17K (GAA→AAA)	frdB	fumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/17);  total (25/17)

GGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCC  >  NC_000913/4380226‑4380711                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   ggAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGtt                                                                                                                                                                                                                                              >  1:267272/1‑250 (MQ=255)         gCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGcac                                                                                                                                                                                                                                     >  2:153233/1‑251 (MQ=255)           tAACGCTTCAACCTTCATACCGTCGGTGTATTCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGcacc                                                                                                                                                                                                                                    >  1:258747/1‑250 (MQ=255)           tAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGcaccc                                                                                                                                                                                                                                   >  2:189157/1‑250 (MQ=255)                  tCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCa                                                                                                                                                                                                                              <  2:15550/249‑1 (MQ=255)                  tCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAAt                                                                                                                                                                                                                            >  1:260681/1‑251 (MQ=255)                                gTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTc                                                                                                                                                                                                               >  2:111643/1‑250 (MQ=255)                                gTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGGGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGGTGCGTCATCAGGCACTTCATAGAATGCGCTATGCGGTGCGGGATCGACTTTCGGGTTATAGCGCACCACCTCCATTTTCAGGTTTTTc                                                                                                                                                                                                               >  2:237572/1‑250 (MQ=255)                                           cACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCAtt                                                                                                                                                                                                   <  1:109525/251‑1 (MQ=255)                                           cACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCAtt                                                                                                                                                                                                   >  1:68286/1‑251 (MQ=255)                                                          tACATGCCAGTTTTGGCACGTTGTTTACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGccttctccttct                                                                                                                                                                                     <  2:176619/250‑1 (MQ=255)                                                                ccAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGgc                                                                                                                                                                              >  1:104950/1‑251 (MQ=255)                                                                                         catGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATcc                                                                                                                                                     <  1:113721/251‑1 (MQ=255)                                                                                         catGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATcc                                                                                                                                                     <  1:189157/251‑1 (MQ=255)                                                                                          atGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCg                                                                                                                                                    >  1:95968/1‑251 (MQ=255)                                                                                           tGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGGAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATTCGc                                                                                                                                                   >  2:4890/1‑251 (MQ=255)                                                                                              cGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCCCTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTc                                                                                                                                                <  1:45744/251‑1 (MQ=255)                                                                                                    aaCCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCAcc                                                                                                                                          >  1:124297/1‑251 (MQ=255)                                                                                                      ccACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGcctt                                                                                                                                                                                             >  2:272740/1‑198 (MQ=255)                                                                                                      ccACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGcctt                                                                                                                                                                                             <  1:272740/198‑1 (MQ=255)                                                                                                          aaaTCGCCATACGGCAGGACCAGCGGTAGCTAAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCCTAAGGCACTTCCTAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCCTCCGCTTCGCCACCGTAAAc                                                                                                                                    >  2:187678/1‑251 (MQ=255)                                                                                                           aaTCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAAcg                                                                                                                                   <  1:294522/251‑1 (MQ=255)                                                                                                                  aTACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGttt                                                                                                                              >  1:139508/1‑249 (MQ=255)                                                                                                                               aGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTgg                                                                                                                        >  1:145848/1‑242 (MQ=255)                                                                                                                               aGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCAGTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTgg                                                                                                                        <  2:145848/242‑1 (MQ=255)                                                                                                                                          aGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTc                                                                                                     >  1:263315/1‑250 (MQ=255)                                                                                                                                                  ggCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCGTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTg                                                                                            <  2:204950/250‑1 (MQ=255)                                                                                                                                                         ttGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGAAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGATGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTcc                                                                                      >  2:157426/1‑250 (MQ=255)                                                                                                                                                           gTCTTTGATGTAGCGCAGCGCAGCCAGTAATGAGGTAGTTGCGTCTTAAGGCACTTCATAGAATTCTCTATGCGGTGCGGTATCGGCTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTGCTCCTTCTTATTGGCTGCTTCCGCCTTATCAGCTGCATCCGCTTCGCCACCGTACACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCTCGTCGCTGTACTCCAg                                                                                    <  2:227540/250‑1 (MQ=255)                                                                                                                                                              tttGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCg                                                                                 >  2:180961/1‑250 (MQ=255)                                                                                                                                                                        cccAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCAt                                                                      >  2:181247/1‑251 (MQ=255)                                                                                                                                                                        cccAGCGCATCCAGTAATGAGGCAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCAt                                                                      >  1:51699/1‑251 (MQ=255)                                                                                                                                                                                     gTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCATTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGa                                                          >  1:205626/1‑250 (MQ=255)                                                                                                                                                                                       aaTGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCGCCACCTCAATTTTCATGTTTTTCATCTCAGCCAGTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGTCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAAGCTTCACGTCGCTGTGGGCCAGGCGAGTCGTGCCATCAGCATCGCTGAAgg                                                       <  2:126561/251‑1 (MQ=255)                                                                                                                                                                                           aGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCg                                                     >  1:133409/1‑249 (MQ=255)                                                                                                                                                                                             gTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCATGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGAGAGTCGTGCCATCAGCATCTCGGAAGGCGAggg                                                 <  1:237572/251‑1 (MQ=255)                                                                                                                                                                                                  tGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCa                                                                                                                                                                                                              >  2:218768/1‑89 (MQ=255)                                                                                                                                                                                                  tGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCa                                                                                                                                                                                                              <  1:218768/89‑1 (MQ=255)                                                                                                                                                                                                       cATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAg                                       <  1:23821/251‑1 (MQ=255)                                                                                                                                                                                                               ccTTCATAGAATGCTCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGa                               <  2:177807/250‑1 (MQ=255)                                                                                                                                                                                                                           gcgcTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGc                    <  1:228032/250‑1 (MQ=255)                                                                                                                                                                                                                                              cGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTcc  >  2:64461/1‑249 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   GGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCC  >  NC_000913/4380226‑4380711

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