Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,683,289	G→A	I53I (ATC→ATT)	dctA ←	C4 dicarboxylate/orotate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,683,289	0	G	A	100.0%	147.6 / NA	41	I53I (ATC→ATT)	dctA	C4 dicarboxylate/orotate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/26);  total (15/26)

ACAATGCCCTGGTCTTTCGCCTGATCGGCGTAAACCGCTACCGCTTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAGATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACA  >  NC_000913/3683059‑3683529                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 aCAATGCCCTGGTCTTTCGCCTGATCGGCGTAAACCGCTACCGCTTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGatcat                                                                                                                                                                                                                              <  1:456246/251‑1 (MQ=255)   aaTGCCCTGGTCTTTCGCCTGATCGGCGTAAACCGCTACCGCTTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATAATCt                                                                                                                                                                                                                            >  1:310098/1‑251 (MQ=255)            gTCTTTCGCCTGATCGGCGTAAACCGCTACCGCTTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCt                                                                                                                                                                                                                   >  1:442360/1‑251 (MQ=255)                      ggATCGGCGTAAACCGCTACCGCTTGCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGcc                                                                                                                                                                                                         <  2:91676/250‑1 (MQ=255)                            gCGTAAACCGCTACCGCTTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGcc                                                                                                                                                                                                   >  2:28472/1‑251 (MQ=255)                                      ctaccgctTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGCTGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGCTTc                                                                                                                                                                                         >  2:339283/1‑251 (MQ=255)                                         ccgctTTCGCATCTAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCATCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGCCCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAATCCGTCGCCAAGCGGTTTCAtt                                                                                                                                                                                      <  2:86482/251‑1 (MQ=255)                                         ccgctTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCAtt                                                                                                                                                                                      <  2:374431/251‑1 (MQ=255)                                          cgctTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCAttt                                                                                                                                                                                     >  1:343017/1‑251 (MQ=255)                                                         cGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAg                                                                                                                                                                      <  2:19730/251‑1 (MQ=255)                                                             gCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTc                                                                                                                                                                                 >  1:408213/1‑236 (MQ=255)                                                             gCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTc                                                                                                                                                                                 <  2:408213/236‑1 (MQ=255)                                                              ccGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAACGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAg                                                                                                                                                                 <  1:35794/251‑1 (MQ=255)                                                                               ccGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGaa                                                                                                                                                 <  2:343017/250‑1 (MQ=255)                                                                                  gCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCATGATAGAAATGGCCAAGGAGAATAcc                                                                                                                                             >  2:36431/1‑251 (MQ=255)                                                                                         gCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCTACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCt                                                                                                                                      <  1:44482/251‑1 (MQ=255)                                                                                                gacgTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTg                                                                                                                               >  2:447546/1‑251 (MQ=255)                                                                                                         gatgatAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGa                                                                                                                         <  1:441779/248‑1 (MQ=255)                                                                                                         gatgatAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGa                                                                                                                         >  2:441779/1‑248 (MQ=255)                                                                                                         gatgatAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCt                                                                                                                      <  2:255144/251‑1 (MQ=255)                                                                                                          atgatAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTg                                                                                                                     <  2:275964/251‑1 (MQ=255)                                                                                                                           cAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTNTa                                                                                                     >  2:421348/1‑250 (MQ=255)                                                                                                                                  aTGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACagaga                                                                                             <  1:332788/251‑1 (MQ=255)                                                                                                                                  aTGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAAGTCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACagaga                                                                                             <  2:332284/251‑1 (MQ=255)                                                                                                                                       aCTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCa                                                                                                                            <  2:155537/215‑1 (MQ=255)                                                                                                                                       aCTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCa                                                                                                                            >  1:155537/1‑215 (MQ=255)                                                                                                                                          gACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCa                                                                                     <  2:293330/251‑1 (MQ=255)                                                                                                                                                                gCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAAcc                                                               <  1:351485/251‑1 (MQ=255)                                                                                                                                                                           ggTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGc                                                                                                                                                                                               <  2:190771/112‑1 (MQ=255)                                                                                                                                                                           ggTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGc                                                                                                                                                                                               >  1:190771/1‑112 (MQ=255)                                                                                                                                                                            gTACGACCGACCGCCTTCATGCTTTCCATGCCCTCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATGTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGaaa                                                                                                                  <  2:432614/188‑1 (MQ=255)                                                                                                                                                                            gTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGaaa                                                                                                                  >  1:432614/1‑188 (MQ=255)                                                                                                                                                                            gTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTa                                                   <  1:28472/251‑1 (MQ=255)                                                                                                                                                                              aCGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGGCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTAGCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTaa                                                  <  2:190392/250‑1 (MQ=255)                                                                                                                                                                                   cgaccgCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGAAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTAt                                            <  2:144579/251‑1 (MQ=255)                                                                                                                                                                                     accgCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATgg                                          <  1:436831/251‑1 (MQ=255)                                                                                                                                                                                     accgCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATgg                                          <  1:239765/251‑1 (MQ=255)                                                                                                                                                                                        gCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGATCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGACAGGACCTGAAAGCAAAGGCTTTTAAACAGAAAGGTTTTCATAGGATGGCCTTTAGTAAAACCACAGGTCTTGTACGTTTATGGGGt                                       >  2:297493/1‑251 (MQ=255)                                                                                                                                                                                             cATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAAGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCt                                   <  1:173423/250‑1 (MQ=255)                                                                                                                                                                                                                         acgacgGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACAt      >  1:410277/1‑251 (MQ=255)                                                                                                                                                                                                                             cgGTACAAAAAATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca  <  1:421348/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 ACAATGCCCTGGTCTTTCGCCTGATCGGCGTAAACCGCTACCGCTTTCGCATCAAGCGTTGCCGGATCGACGTTCATTCCGGCACCAGGCTGCACGACGTTAACGATGATAAGACCAATAATCAGCGCGATGGTACTGACAATTTCAAAGTAAAGCAGTGCGACTGCGCCGGTACGACCGACCGCCTTCATGCTTTCCATGCCCGCAATGCCCGTTACGACGGTACAAAAGATGACAGGAGCGATGATCATCTTAATGAGCTTAACGAAGCCGTCGCCAAGCGGTTTCATTTGCTCGCCTATTTCAGGATAGAAATGGCCAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACA  >  NC_000913/3683059‑3683529

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