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%	204.8 / NA	59	E17K (GAA→AAA)	frdB	fumarate reductase iron‑sulfur protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/31);  total (28/31)

TTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCATTGCGGAGTGCGCCAT  >  NC_000913/4380181‑4380764                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                        tttCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATTAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAgc                                                                                                                                                                                                                                                                                              <  1:9679/300‑1 (MQ=255) tttCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAgc                                                                                                                                                                                                                                                                                              <  1:42041/300‑1 (MQ=255) tttCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAg                                                                                                                                                                                                                                                                                               <  1:260883/299‑1 (MQ=255)        gAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGcaccac                                                                                                                                                                                                                                                                                       >  1:487240/1‑300 (MQ=255)             gggTCGTATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCATGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCACAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCCCTTCATAGAATGCGCTCTGCGGTACGGTACCGACTTTCGGGTTATAGTGCACCACCTCaa                                                                                                                                                                                                                                                                                  >  2:22824/1‑300 (MQ=255)             gggTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCaa                                                                                                                                                                                                                                                                                  >  2:485943/1‑300 (MQ=255)             gggTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTATTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCaa                                                                                                                                                                                                                                                                                  >  2:22788/1‑300 (MQ=255)               gTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTACCCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAAtt                                                                                                                                                                                                                                                                                >  1:697122/1‑300 (MQ=255)                                                aaGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGccttct                                                                                                                                                                                                                                                >  1:347624/1‑299 (MQ=255)                                                 aGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGccttctcc                                                                                                                                                                                                                                              <  1:188788/300‑1 (MQ=255)                                                          aCGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTg                                                                                                                                                                                                                                     <  1:640417/300‑1 (MQ=255)                                                                        aTACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCt                                                                                                                                                                                                                        <  1:668213/299‑1 (MQ=255)                                                                                 gtgtAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGc                                                                                                                                                                                                              >  1:443810/1‑300 (MQ=255)                                                                                    tAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATc                                                                                                                                                                                                           <  1:588569/300‑1 (MQ=255)                                                                                             aGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGcc                                                                                                                                                                                                  <  2:393277/300‑1 (MQ=255)                                                                                               gAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCAc                                                                                                                                                                                                <  1:542355/300‑1 (MQ=255)                                                                                                               agtgttGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGATCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGc                                                                                                                                                                                <  2:443900/296‑1 (MQ=255)                                                                                                               aGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACTGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGc                                                                                                                                                                                <  1:86283/300‑1 (MQ=255)                                                                                                                gTTTTTGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGc                                                                                                                                                                                <  2:443858/299‑1 (MQ=255)                                                                                                                           ttgttAACCATCATGCCGCAGGACCCACACATCGCCATACGGCAGGACCAGCGGTAGTTTAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTAATAGAATGCGCTATGCGTTGCGGTATCGACTTTCGGGTTATAGCGCACCAACTAAAGTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGTCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCgt                                                                                                                                                                    >  2:477461/1‑300 (MQ=255)                                                                                                                           ttgttAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGTGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGc                                                                                                                                                                      >  2:498363/1‑298 (MQ=255)                                                                                                                           ttgttAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCgt                                                                                                                                                                    >  1:8282/1‑300 (MQ=255)                                                                                                                             gttAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCgta                                                                                                                                                                   >  1:592063/1‑299 (MQ=255)                                                                                                                                                     acaAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTANTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTc                                                                                                                                            <  2:572702/298‑1 (MQ=255)                                                                                                                                                     acaAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTAcgc                                                                                                                                            <  2:16980/298‑3 (MQ=255)                                                                                                                                                     acaAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTc                                                                                                                                            <  1:78922/298‑1 (MQ=255)                                                                                                                                                     acaAACCGCCATACGGCAGGACCAGCGGTAACTCAGGTCCGGTGCCAGGTTGTCTTTGATTTAGCCCAGCGCATCCAGTAATTATGTANTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTc                                                                                                                                            <  2:572715/298‑1 (MQ=255)                                                                                                                                                        aaTCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGc                                                                                                                                       <  1:181492/300‑1 (MQ=255)                                                                                                                                                                  cGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGcc                                                                                                                             >  1:579729/1‑300 (MQ=255)                                                                                                                                                                              cGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCAACTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGACTTATCGGCTGCCTCCGCTTCGCAACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCAAGGCGAGTCGTGCAATCAGCATc                                                                                                                    >  2:282772/1‑297 (MQ=255)                                                                                                                                                                                        ggTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCATTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGa                                                                                                         >  2:660702/1‑298 (MQ=255)                                                                                                                                                                                              gtggaggggtGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGTACCACCTCAATTTTCGGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTAAGCCTTATCGGTTGCATCCGCTTCGCCACCGTAAACGCGTGTAGCTGGCGGCAGCGTAGTAATCtt                                                                                                                                                           >  2:245467/10‑242 (MQ=255)                                                                                                                                                                                              gTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTATACTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCAAGTCGCTGTACTCCAGGCGAGTCTTGCAATCAGCATCGAGGCAGGCGAGGgtgtgt                                                                                                 >  2:245874/1‑300 (MQ=255)                                                                                                                                                                                              gTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTAAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGgtgtgt                                                                                                 >  1:269268/1‑300 (MQ=255)                                                                                                                                                                                               tGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGtt                                                                                                <  1:501728/300‑1 (MQ=255)                                                                                                                                                                                               tGCCAGGTTGTCTTTGANGTAGCCCAGCGCATCCAGNAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTCTGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGNAAGGCGAGGGTGTGtt                                                                                                <  1:613959/300‑1 (MQ=255)                                                                                                                                                                                                   aGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAg                                                                                            <  2:144166/300‑1 (MQ=255)                                                                                                                                                                                                      ttGTCTTTGATGTAGCCCAGCGCATCCAGTACTGAGTTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGTTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATCGGCTGCTTCCGCCTTATCGGCTACATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGa                                                                                          <  2:353482/299‑1 (MQ=255)                                                                                                                                                                                                        gTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCAATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGt                                                                                       >  1:265780/1‑300 (MQ=255)                                                                                                                                                                                                         tCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGtt                                                                                      <  1:1567/300‑1 (MQ=255)                                                                                                                                                                                                               aTGTAGCCCAGCGCATCCATTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTAAGCCATTCGACTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAGTCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCAATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGa                                                                                    >  1:31575/1‑296 (MQ=255)                                                                                                                                                                                                                    gCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTAATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTCGAGGAAGTTGACGTCGTCAc                                                                           >  2:61034/1‑300 (MQ=255)                                                                                                                                                                                                                        aGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTAAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCt                                                                        <  1:353347/299‑1 (MQ=255)                                                                                                                                                                                                                                cAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAAc                                                               >  2:665614/1‑300 (MQ=255)                                                                                                                                                                                                                                 aGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAAcc                                                              <  1:649455/300‑1 (MQ=255)                                                                                                                                                                                                                                  ggaaTGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGGTTGAGGGAGTTGatgtaggtatggtaggggagagtt                                                             <  1:22788/298‑24 (MQ=255)                                                                                                                                                                                                                                  gTAATGAGGTATTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTTGCTGCTTCCGCCTTATCGGCTGCATCCGCTTTGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCt                                                             <  2:639046/300‑1 (MQ=255)                                                                                                                                                                                                                                  gTAATGAGGTAGTTGCGTCATAATGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGGTTGAGGAAGTTGatgtaggtatggtaggggagagtt                                                             <  1:22824/300‑24 (MQ=255)                                                                                                                                                                                                                                             gTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACg                                                  <  1:439931/300‑1 (MQ=255)                                                                                                                                                                                                                                                             cTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCANCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTCCGTCAAGACGCTGGTGCGcgccgc                                    >  2:233454/1‑298 (MQ=255)                                                                                                                                                                                                                                                                   aGAATGCGCTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTCCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGActct                            >  2:632368/1‑300 (MQ=255)                                                                                                                                                                                                                                                                        gcgcTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAATGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTAc                        >  1:90760/1‑299 (MQ=255)                                                                                                                                                                                                                                                                        gcgcTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAANCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCNATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACg                       <  2:535500/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          gcTATGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACATAAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGCCGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTaagt                      >  2:90401/1‑296 (MQ=255)                                                                                                                                                                                                                                                                             aTGCGGTGCGGTATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCa                  >  1:479619/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                      ggtATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCATTGCGGAg          <  1:62967/299‑1 (MQ=255)                                                                                                                                                                                                                                                                                       gtATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGCCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTNGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCATTGCGGAGt         >  1:516490/1‑299 (MQ=255)                                                                                                                                                                                                                                                                                        tATCGACTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCATTGCGGAGTgc       <  1:296683/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             aCTTTCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCATTGCGGAGTGCGCCAt  >  1:710871/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                                        TTTCGATGAAGTGGGTCATATCGACCACCAGATCGCGTTCAATCGGGAAGTTAGCTAACGCTTCAACCTTCATACCGTCGGTGTAATCACGCAGGAAGGTTTTACATGCCAGTTTTGGCACGTTGTTAACCATCATGCCGCAGGAACCACAAATCGCCATACGGCAGGACCAGCGGTAGCTCAGGTCCGGTGCCAGGTTGTCTTTGATGTAGCCCAGCGCATCCAGTAATGAGGTAGTTGCGTCATAAGGCACTTCATAGAATGCGCTATGCGGTGCGGTATCGACTTCCGGGTTATAGCGCACCACCTCAATTTTCAGGTTTTTCATCTCAGCCATTCGCCTTCTCCTTCTTATTGGCTGCTTCCGCCTTATCGGCTGCATCCGCTTCGCCACCGTAAACGCGTTTAGCTGGCGGCAGCGTAGTAATCTTCACGTCGCTGTACTCCAGGCGAGTCGTGCCATCAGCATCGCGGAAGGCGAGGGTGTGTTTGAGGAAGTTGACGTCGTCACGCTCGGTGCAACCTTCGTCCAGACGCTGGTGCGCGCCGCGGGACTCTTTACGTGCCATTGCGGAGTGCGCCAT  >  NC_000913/4380181‑4380764

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