Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,978,723	C→T	K56K (AAG→AAA)	unknown ←	RTX toxin transporter, ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,978,723	0	C	T	100.0%	120.3 / NA	44	K56K (AAG→AAA)	unknown	RTX toxin transporter, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/26);  total (18/26)

CAATAGATGCGCGTGATGCAATTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGAATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTTAATAAAATAAAGTGTTCACCATCATCACGCCAAGCCAATGCCGGTAAATGTAAAAAAGGCAATCTGTCAACTGTTTTGTTCGCAGTTCTTACCTTTAATCCAAGTGATTTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCCGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGGGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCGCACACACAAATATGCGGAATCACTAAAAAGCCAATAACTCATTAAGCA  >  1/1978497‑1978965                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   cAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATc                                                                                                                                                                                                                                             >  1:19746/1‑234 (MQ=255) cAATAGATGCGCGTGACGCAACTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTtgctg                                                                                                                                                                                                                              <  1:339897/249‑1 (MQ=255)                          aaTCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAgcc                                                                                                                                                                                                   <  1:183974/251‑1 (MQ=255)                           aTCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAgccc                                                                                                                                                                                                  <  1:622500/251‑1 (MQ=255)                                cTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAgcccatgcc                                                                                                                                                                                             <  1:217913/251‑1 (MQ=255)                                                 cttcAAATTCA‑TTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTAt                                                                                                                                                                           >  1:611760/1‑251 (MQ=255)                                                          cATTTTTGTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTAtt                                                                                                                                                                   >  1:875279/1‑251 (MQ=255)                                                                 gTCTAAAACAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTc                                                                                                                                                                                                       >  1:313598/1‑208 (MQ=255)                                                                      aaaCAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAgcccatgcc                                                                                                                                                                                             >  1:148087/1‑213 (MQ=255)                                                                         cAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATc                                                                                                                                                                                                                                             >  1:348643/1‑162 (MQ=255)                                                                         cAATAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAAtt                                                                                                                                                    >  1:379943/1‑251 (MQ=255)                                                                          aaTAGGATTTTTTTGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTg                                                                                                                                                   <  1:28126/251‑1 (MQ=255)                                                                                     tttGCATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGa                                                                                                                                        <  1:456801/251‑1 (MQ=255)                                                                                               aTCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGc                                                                                                                                <  1:899261/249‑1 (MQ=255)                                                                                                tCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTatat                                                                                                                                                                                    <  1:843933/196‑1 (MQ=255)                                                                                                tCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTAAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAAtt                                                                                                                                                    <  1:320068/228‑1 (MQ=255)                                                                                                   aaaaTGAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTAc                                                                                                                          <  1:146660/251‑1 (MQ=255)                                                                                                        gAGGTAACGATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAgcccatgcc                                                                                                                                                                                             <  1:300891/179‑1 (MQ=255)                                                                                                                gATCGGTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCTTTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGa                                                                                                             >  1:566830/1‑251 (MQ=255)                                                                                                                     gTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTc                                                                                                                                                                                                       <  1:240188/156‑1 (MQ=255)                                                                                                                     gTTTCTTGATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATc                                                                                                                                                                                         <  1:889672/170‑1 (MQ=255)                                                                                                                            gATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATc                                                                                                                                                                                         <  1:889977/163‑1 (MQ=255)                                                                                                                            gATCTATTTTCAATAAAATAAAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCta                                                                                                                                                                                       >  1:799754/1‑165 (MQ=255)                                                                                                                                       aataaaataaAGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCtg                                                                                      <  1:310667/251‑1 (MQ=255)                                                                                                                                           aaataaaGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAg                                                                                                           >  1:481630/1‑226 (MQ=255)                                                                                                                                             ataaaGTGGTCACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTatat                                                                                                                                                                                    <  1:502307/151‑1 (MQ=255)                                                                                                                                                       cACCATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACatat                                                                      >  1:262410/1‑251 (MQ=255)                                                                                                                                                         ccATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATcc                                                                    <  1:677991/251‑1 (MQ=255)                                                                                                                                                         ccATCATCACGCCAAGCCAATACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATc                                                                     >  1:201576/1‑250 (MQ=255)                                                                                                                                                                       agccaaTACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGatat                                                                                                                                     >  1:745969/1‑172 (MQ=255)                                                                                                                                                                       agccaaTACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTAc                                                                          <  1:470863/231‑1 (MQ=255)                                                                                                                                                                       agccaaTACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGAt                                                      >  1:480436/1‑251 (MQ=255)                                                                                                                                                                        gccaaTACCGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTAttt                                                                                                                                                                  >  1:198015/1‑142 (MQ=255)                                                                                                                                                                               ccGGTAAATGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAatat                                                                                                                                             <  1:421318/156‑1 (MQ=255)                                                                                                                                                                                      aTGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGcc                                                                                                                <  1:56080/178‑1 (MQ=255)                                                                                                                                                                                      aTGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCACACATACAAat                                       <  1:421052/251‑1 (MQ=255)                                                                                                                                                                                       tGTAAAAAAGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCACACATACAAata                                      <  1:501623/251‑1 (MQ=255)                                                                                                                                                                                          aaaaaaGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCACACATACAAATATGc                                   <  1:31942/251‑1 (MQ=255)                                                                                                                                                                                          aaaaaaGGTAATCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCACACATACAAATATGc                                   <  1:7619/251‑1 (MQ=255)                                                                                                                                                                                               aggtaaTCTATCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAAt                                                                                                                                                                               >  1:245187/5‑106 (MQ=255)                                                                                                                                                                                                     tctatCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAAc                                                            <  1:129914/211‑1 (MQ=255)                                                                                                                                                                                                       tatCAATTGTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAat                                                                                                                                               <  1:858693/128‑1 (MQ=255)                                                                                                                                                                                                               gTTTTGTTCGTAGTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCaa                                                                                                                    >  1:391153/1‑149 (MQ=255)                                                                                                                                                                                                                           gTTCTTACTTTTAATCCAAGTGATCTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCAGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGAGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCACACATACAAATATGCGGAATCACTAAAAAGCCAATAACTCATTAAGCa  >  1:184872/1‑251 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   CAATAGATGCGCGTGATGCAATTAAAATCACCTTGCTTTGATAACGTTCTTCAAATTCATTTTTGTCTAAAACAATAGGATTTTTTTGAATTAGATCAAAAATGAGGTAACGATCGGTTTCTTGATCTATTTTTAATAAAATAAAGTGTTCACCATCATCACGCCAAGCCAATGCCGGTAAATGTAAAAAAGGCAATCTGTCAACTGTTTTGTTCGCAGTTCTTACCTTTAATCCAAGTGATTTTGCTGCAAGCAACCATTTAGTTTGGTTCAGCCCATGCCCATCTATATCAAATTTATGTTTTATTTCTTCCGGGTTAATTGAAATATTATGATATTGAGCCAGTACTTCCAATGCCTGAAGGGCAAGATTAGTATTTTTCTGTGAGTCCATTACATATCCTGAAAAACAGTGATTCCGCACACACAAATATGCGGAATCACTAAAAAGCCAATAACTCATTAAGCA  >  1/1978497‑1978965

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