Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	92,612	+A	coding (176/222 nt)	unknown →	Sialic acid‑specific 9‑O‑acetylesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	92,612	1	.	A	100.0%	119.8 / NA	45	G59D (GGT→GAT)	unknown	Sialic acid‑specific 9‑O‑acetylesterase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base A (31/14);  total (31/14)

GTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGGCATTGTATATGTCCTTTTTTATTGTCAACTTTTGCTTCTTAAAGTGCGGTGG‑TTTTTACGGCGTTTTGAAAATCACTTATGCCCGACCCGAATTGTAATATCCCTGAATGCCAGGGCATTGCGGCTAACGCCGGGCTGTACCTCAATACGCACACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGACCACCTCGGCGC  >  1/92491‑92826                                                                                                                           |                                                                                                                                                                                                                       gtttaaacatagcacgctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCt                                                                                                                             >  1:863311/22‑213 (MQ=19)  tttaaacatagcacgctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTGTATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAgg                                                                                       >  1:23162/21‑250 (MQ=255)  tttaaacatagcacgctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAgg                                                                                       <  1:88134/230‑1 (MQ=255)        catagcacgctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAg                                                                                >  1:564121/15‑251 (MQ=255)         atagcacgctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAg                                                                                >  1:92038/14‑250 (MQ=255)         atagcacgctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGc                                                                               <  1:858111/238‑1 (MQ=255)                 ctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATgtcc                                                                                                                                                                                        <  1:508031/133‑4 (MQ=21)                 ctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGcc                                                                                                                                        <  1:119711/181‑1 (MQ=25)                 ctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCTAGGCTGTCAGCCACttgtt                                                                       >  1:856684/6‑251 (MQ=255)                 ctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACttgtt                                                                       >  1:804044/6‑251 (MQ=255)                 ctccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACttgtt                                                                       >  1:651259/6‑251 (MQ=255)                   ccttttgtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACttgttg                                                                      >  1:557908/4‑250 (MQ=255)                         gtttgtttgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTTCGGCTAATGCCCGGTTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATg                                                                >  1:657370/1‑250 (MQ=255)                                tgtttgtttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGt                                                        >  1:626285/1‑251 (MQ=255)                                      tttgtttgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTg                                                  <  1:134266/251‑1 (MQ=255)                                            tgtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGt                                            <  1:603975/251‑1 (MQ=255)                                             gtttgtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTAt                                                                                                        >  1:569947/1‑190 (MQ=255)                                                 gtttgtttgtttgtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGcc                                        >  1:762435/1‑250 (MQ=255)                                                             gtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCAc                                                                            >  1:551133/1‑202 (MQ=255)                                                             gtttgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGt                                            >  1:400327/1‑234 (MQ=255)                                                                tgttt‑gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAAtg                                                                                                                        >  1:606460/1‑154 (MQ=255)                                                                  gtttgCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATc                                                                                                                                                                      <  1:882847/104‑1 (MQ=255)                                                                    ttgCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTAt                     <  1:383425/249‑1 (MQ=255)                                                                      gCATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATcc                                                                                                                                                                     <  1:853778/105‑1 (MQ=255)                                                                       cATTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGAt                  >  1:170735/1‑251 (MQ=255)                                                                        aTTGTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGAtt                 <  1:462446/251‑1 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATgtcc                                                                                                                                                                                        >  1:596978/1‑78 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGt                                                                                                                                 >  1:263632/1‑136 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGGATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGa              >  1:435741/1‑251 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATcc                                                                                                      >  1:252855/1‑163 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCgg                                                                   <  1:639106/198‑1 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAg                                             <  1:854441/220‑1 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCt                                >  1:464081/1‑233 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGAtt                 >  1:54683/1‑248 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGAtt                 >  1:343979/1‑248 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTg               >  1:227627/1‑250 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGa              >  1:4501/1‑251 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGa              >  1:851310/1‑251 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGa              >  1:563625/1‑251 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGa              >  1:356458/1‑251 (MQ=255)                                                                           gTATATGTCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGa              >  1:894519/1‑251 (MQ=255)                                                                                  tCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGACCACCt        >  1:727645/1‑250 (MQ=255)                                                                                  tCCTTTTTTATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGACCACCTc       <  1:547041/251‑1 (MQ=255)                                                                                      tttttATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGACCACCTCGgcg   >  1:776912/1‑251 (MQ=255)                                                                                       ttttATTGTCAAATTTTGATTATAAAAATGCAGTGGATTTTTACGGCGTTTTGAAAATCACTTATGTCCCACCCGAATTGTAATATCCCTAAATGCCAGCGCGTTGCGGCTAATGCCCGGCTGTACCTCAATGCGCAGACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGACCACCTCGgcgc  <  1:552294/251‑1 (MQ=255)                                                                                                                           |                                                                                                                                                                                                                       GTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGTTTGGCATTGTATATGTCCTTTTTTATTGTCAACTTTTGCTTCTTAAAGTGCGGTGG‑TTTTTACGGCGTTTTGAAAATCACTTATGCCCGACCCGAATTGTAATATCCCTGAATGCCAGGGCATTGCGGCTAACGCCGGGCTGTACCTCAATACGCACACGAATTTTATCCACATTTTGCGCGAGGCTGTCAGCCACTTGTTGCGGATGCCAATCGTCAGTTGTCCAGTCGCCCAAAAGCTCGGATGGGTATGATTTGACCACCTCGGCGC  >  1/92491‑92826

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