Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	509,517	C→T	E621E (GAG→GAA)	copA ←	Cu(+) exporting P‑type ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	509,517	0	C	T	100.0%	221.9 / NA	63	E621E (GAG→GAA)	copA	Cu(+) exporting P‑type ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/31);  total (32/31)

CCATTGCCACCTGACGTCCTTCACTTTGCAGATGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCTTCAGTC  >  NC_000913/509228‑509796                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                        ccaTTGCCACCTGACGTCCTTCACTTTGCAGATGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGa                                                                                                                                                                                                                                                                               >  1:151707/1‑300 (MQ=255) ccaTTGCCACCTGACGTCCTTCACTTTGCAGATGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGa                                                                                                                                                                                                                                                                               >  2:376658/1‑300 (MQ=255)      gCCACCTGACGTCCCTCACTTTGCAGATGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCATGAGCAGTAATTTCCGCTTCGATAGCt                                                                                                                                                                                                                                                                          >  1:225285/1‑300 (MQ=255)                gTCCTTCACTTTGCAGATGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCaa                                                                                                                                                                                                                                                                >  1:594624/1‑300 (MQ=255)                                aTGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTaa                                                                                                                                                                                                                                                >  1:160732/1‑300 (MQ=255)                                 tGTTTGATCGCTTCTGCTTTACCGTCCGGCAGCACCCCTGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAAc                                                                                                                                                                                                                                               <  2:283020/300‑1 (MQ=255)                                 tGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGTTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAAc                                                                                                                                                                                                                                               >  1:391047/1‑300 (MQ=255)                                   tttGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTTGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCNCCACGCTATCACTACGCAACGGATCGNGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAg                                                                                                                                                                                                                                             <  2:462325/300‑1 (MQ=255)                                   tttGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCNCCACGCTATCACTACGCAACGGATCGNGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAg                                                                                                                                                                                                                                             <  2:462895/300‑1 (MQ=255)                                   tttGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAg                                                                                                                                                                                                                                             <  1:241291/300‑1 (MQ=255)                                   tttGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAg                                                                                                                                                                                                                                             >  2:660285/1‑300 (MQ=255)                                           cTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTNTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGAt                                                                                                                                                                                                                                     <  2:515581/300‑1 (MQ=255)                                                           cGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCAGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCATGATTGCCCAGCAATAACGc                                                                                                                                                                                                                     >  1:210655/1‑300 (MQ=255)                                                                            tCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCAc                                                                                                                                                                                                    <  1:412920/300‑1 (MQ=255)                                                                             cACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCAcc                                                                                                                                                                                                   >  2:75863/1‑300 (MQ=255)                                                                               ccTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGc                                                                                                                                                                                                 >  1:150261/1‑300 (MQ=255)                                                                                 tcatcaATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTc                                                                                                                                                                                               >  1:214876/1‑300 (MQ=255)                                                                                           cTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGcc                                                                                                                                                                                     <  1:174424/300‑1 (MQ=255)                                                                                                  tttGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCACTCACGCCCAGCCCGCGCaa                                                                                                                                                                              <  2:38690/300‑1 (MQ=255)                                                                                                          tCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCAACAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGATACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAACAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCACGTACAATGTGCGGa                                                                                                                                                                      >  2:93655/1‑300 (MQ=255)                                                                                                             cATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAAc                                                                                                                                                                   <  1:75304/300‑1 (MQ=255)                                                                                                                     gtggtTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGAcc                                                                                                                                                           >  1:652179/1‑300 (MQ=255)                                                                                                                     gtggtTGGGTTATCCCAGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCATGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGCTGCTCATTTAACAGCGCCTGATTGCCCAGTAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGAGGAAACGGTTGAcc                                                                                                                                                           >  2:57875/1‑300 (MQ=255)                                                                                                                           gggTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCTCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGt                                                                                                                                                     >  1:618037/1‑300 (MQ=255)                                                                                                                            ggTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATTGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCCTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTa                                                                                                                                                    >  2:238745/1‑300 (MQ=255)                                                                                                                                        tCAATATCACCAGACGATATCCCGCTTTATTCAGATGTTGTATCGCCTCCACGCTATCACTACGCAACGGATCGCTTACTGCCATCAGTGCTACCGCTTTCCCTTCAACCGCCAGCGGCACAGGCGTTGCCCCTTTCGATGTCTGAGCAGTGAGTGCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCGTTTAACAGCGCCTGATTGCCTAGCAATAACGCACGACCTTCAGCTTCACCGCTCACACCCAACCCGCGCAATCTGCGGAAAccttcgccctgcattctctccacctccc                                                                                                                                        <  2:247522/300‑27 (MQ=255)                                                                                                                                           aCATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCAACTg                                                                                                                                     >  2:678588/1‑300 (MQ=255)                                                                                                                                                 ccAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTAt                                                                                                                               >  1:580595/1‑300 (MQ=255)                                                                                                                                                        atatCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGTTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTACTTTATCGAGGAt                                                                                                                        <  2:369653/300‑1 (MQ=255)                                                                                                                                                        atatCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTACCATCAGTGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTTCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGAt                                                                                                                        <  2:178646/300‑1 (MQ=255)                                                                                                                                                                       cAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCgg                                                                                                         <  2:430804/300‑1 (MQ=255)                                                                                                                                                                          gCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATg                                                                                                      >  1:320065/1‑300 (MQ=255)                                                                                                                                                                                  gcgcCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAAc                                                                                              <  1:34705/300‑1 (MQ=255)                                                                                                                                                                                   cgccgcCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAAcc                                                                                             <  1:632972/300‑1 (MQ=255)                                                                                                                                                                                   cgccgcCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCACTTTCCCGTCAACCACCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTNCCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAAcc                                                                                             <  2:154375/300‑1 (MQ=255)                                                                                                                                                                                    gccgccACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTNTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTAGCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCt                                                                                            <  2:503572/300‑1 (MQ=255)                                                                                                                                                                                          aCGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTcc                                                                                      <  1:442660/300‑1 (MQ=255)                                                                                                                                                                                           cGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCACCCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGACCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATCTCACCTGCTTTATCAAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCa                                                                                     >  2:145401/1‑300 (MQ=255)                                                                                                                                                                                               aTCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCGCGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTgc                                                                                 <  2:665473/300‑1 (MQ=255)                                                                                                                                                                                                   cTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTgccgcc                                                                             >  1:160596/1‑300 (MQ=255)                                                                                                                                                                                                      cGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCNGCACAGTCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGAGAGCGGAGGGcgaagatcgtgtgccggtgtcggcgtg                                                                          <  2:518981/300‑27 (MQ=255)                                                                                                                                                                                                          aCGGATCGCGGACTGTCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGAc                                                                      <  2:225860/300‑1 (MQ=255)                                                                                                                                                                                                                       tGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTTCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTgcgc                                                         <  2:662503/300‑1 (MQ=255)                                                                                                                                                                                                                          cagcagGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTtc                                                      >  1:545248/1‑300 (MQ=255)                                                                                                                                                                                                                            gcagGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATATGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTtcat                                                    <  2:4707/300‑1 (MQ=255)                                                                                                                                                                                                                               gggCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCNAGCAATAACGCATGACCTTCAGCTTCAACGCTCACGCCCAGCCCGCGCAATGTGCGGAAACNGTTGACCTGCGGTAGCTGCATATCAACTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCaa                                                 >  2:551001/1‑300 (MQ=255)                                                                                                                                                                                                                               gggCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCaa                                                 >  1:326591/1‑300 (MQ=255)                                                                                                                                                                                                                                ggCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGNCCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAAc                                                >  1:466076/1‑300 (MQ=255)                                                                                                                                                                                                                                                tCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTc                                <  1:259056/300‑1 (MQ=255)                                                                                                                                                                                                                                                   aCCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGNGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTNTTCACt                             >  2:559013/1‑300 (MQ=255)                                                                                                                                                                                                                                                           cagcaCAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAAc                     >  1:407109/1‑300 (MQ=255)                                                                                                                                                                                                                                                            agcaCAGGCGTTGCCCCTTTCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCAGTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGGGCGGAAACCGTGGACCTGCGGGAGCTGCATATCACCTGCTTTAGCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAAcc                    <  2:209537/300‑1 (MQ=255)                                                                                                                                                                                                                                                             gcaCAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCt                   >  2:219968/1‑300 (MQ=255)                                                                                                                                                                                                                                                                 atgcggtGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCAGCAACATCAGCAAATGTTTTCACTGCGACAACCTGCgg               <  2:206335/294‑1 (MQ=255)                                                                                                                                                                                                                                                                  ggCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCANGCCCAGCCCGCGCAATGTGCGGAAACNGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCCCTGCGCCAACCTGCGGc              >  2:467545/1‑300 (MQ=255)                                                                                                                                                                                                                                                                  ggCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGc              <  2:379423/300‑1 (MQ=255)                                                                                                                                                                                                                                                                   gCGTTGCCCCTTGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCt             >  1:687007/1‑300 (MQ=255)                                                                                                                                                                                                                                                                   gCGTTGCCCATTGCGATGCATGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCACAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCAACTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCCCTGCGACCACCTGCGGCt             >  2:182159/1‑300 (MQ=255)                                                                                                                                                                                                                                                                         ccccTGGCGATGCCGGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAAGGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCtt       <  2:181338/300‑1 (MQ=255)                                                                                                                                                                                                                                                                             ttGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCTtcagt   >  1:45256/1‑300 (MQ=255)                                                                                                                                                                                                                                                                              tGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCTtcagtc  <  1:43252/300‑1 (MQ=255)                                                                                                                                                                                                                                                                              tGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCTtcagtc  <  2:302437/300‑1 (MQ=255)                                                                                                                                                                                                                                                                              tGCGATGCCTGAGCAGTAATTTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCAACTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCTtcagtc  <  1:674162/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                        CCATTGCCACCTGACGTCCTTCACTTTGCAGATGTTTGATCGCTTCGGCTTTACCGTCCGGCAGCACCCCGGCGATCACCTCATCAATCCCTGCTTCTTTGGCGATCGCATTGGCGGTGGTTGGGTTATCCCCGGTCAACATCACCAGACGATATCCCGCTTTATGCAGGCGTTGCAGCGCCGCCACGCTATCACTACGCAACGGATCGCGTACTGCCAGCAGGGCTACCGCTTTCCCGTCAACCGCCAGCAGCACAGGCGTTGCCCCTTGCGATGCCTGAGCAGTAATCTCCGCTTCGATAGCTTTGGTACCAACCTGTTGCTCATTTAACAGCGCCTGATTGCCCAGCAATAACGCATGACCTTCAGCTTCACCGCTCACGCCCAGCCCGCGCAATGTGCGGAAACCGTTGACCTGCGGTAGCTGCATATCACCTGCTTTATCGAGGATCGCTCGTGCCAGCGGATGGCTGGAACCTTGCTCCAGTGCCGCCGCCAGACGCAATGCCTGCGCTTCATCAACATCAGCAAATGTTTTCACTGCGACAACCTGCGGCTTCCCTTCAGTC  >  NC_000913/509228‑509796

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