Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,909,126	C→T	H254H (CAC→CAT)	unknown →	Spermidine/putrescine import ABC transporter permease protein PotC (TC 3.A.1.11.1)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,909,126	0	C	T	100.0%	149.2 / NA	52	H254H (CAC→CAT)	unknown	Spermidine/putrescine import ABC transporter permease protein PotC (TC 3.A.1.11.1)
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/32);  total (20/32)

CGTAAAATCATCCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACACAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTCACCGCAAGTATGATCGCTGTTGCACTTACCGGTTGTAACGATAAAGACAACAAGCAAGCGCAAGCAACGCCGGAAGCC  >  1/1908881‑1909370                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     cGTAAAATCATCCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGa                                                                                                                                                                                                                                                 <  1:173756/251‑1 (MQ=255)   tAAAATCATCCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGAc                                                                                                                                                                                                                                                <  1:702408/250‑1 (MQ=255)   tAAAATCATCCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGAc                                                                                                                                                                                                                                                <  1:641782/250‑1 (MQ=255)          atcCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGat                                                                                                                                                                                                                                        <  1:18567/251‑1 (MQ=255)          atcCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGat                                                                                                                                                                                                                                        <  1:513728/251‑1 (MQ=255)            ccTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGatat                                                                                                                                                                                                                                      <  1:25185/251‑1 (MQ=255)                ccGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAGCCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGGCATCTAGATataat                                                                                                                                                                                                                                   <  1:41649/250‑1 (MQ=255)                 cGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATAtaataa                                                                                                                                                                                                                                 <  1:368834/251‑1 (MQ=255)                    ttGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGc                                                                                                                                                                                                                              <  1:40306/251‑1 (MQ=255)                    ttGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGc                                                                                                                                                                                                                              <  1:647485/251‑1 (MQ=255)                    ttGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGc                                                                                                                                                                                                                              <  1:881783/251‑1 (MQ=255)                    ttGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGc                                                                                                                                                                                                                              <  1:466774/251‑1 (MQ=255)                      gCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATataat                                                                                                                                                                                                                                   <  1:661632/244‑1 (MQ=255)                                  gtggtgTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAAttt                                                                                                                                                                                                                >  1:400127/1‑251 (MQ=255)                                       gtCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAg                                                                                                                                                                                                           >  1:605607/1‑251 (MQ=255)                                        tCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAg                                                                                                                                                                                                           <  1:595898/250‑1 (MQ=255)                                        tCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAg                                                                                                                                                                                                           <  1:827362/250‑1 (MQ=255)                                           ggttggttGTTAAGTTTCGCCATTTCACTGGATGACGTTGTGGCGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGat                                                                                                                                                                                                                                        <  1:60920/218‑1 (MQ=255)                                           ggttggttGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGat                                                                                                                                                                                                                                        <  1:465110/218‑1 (MQ=255)                                              tggttgTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATAc                                                                                                                                                                                                    <  1:64251/251‑1 (MQ=255)                                                gttgttAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATataat                                                                                                                                                                                                                                   >  1:778938/1‑218 (MQ=255)                                                  tgttAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATa                                                                                                                                                                                                <  1:851258/251‑1 (MQ=255)                                                    ttAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATAc                                                                                                                                                                                               <  1:829841/250‑1 (MQ=255)                                                       aGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGAc                                                                                                                                                                                                                                                <  1:307885/198‑1 (MQ=255)                                                             aCCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGc                                                                                                                                                                                     <  1:559212/251‑1 (MQ=255)                                                                   tCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGa                                                                                                                                                                                                                                                 >  1:616570/1‑185 (MQ=255)                                                                             aCGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGAc                                                                                                                                                                                                                                                >  1:309406/1‑176 (MQ=255)                                                                             aCGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGc                                                                                                                                                                                             >  1:794363/1‑227 (MQ=255)                                                                                    ggtgTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGc                                                                                                                                                                                                                              >  1:878299/1‑187 (MQ=255)                                                                                     gtgtCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATa                                                                                                                                                                                                >  1:99919/1‑216 (MQ=255)                                                                                     gtgtCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATAc                                                                                                                                                              >  1:466566/1‑250 (MQ=255)                                                                                       gtCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGat                                                                                                                                                                                                                                        <  1:29916/174‑1 (MQ=255)                                                                                            cctTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATGTAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCAttt                                                                                                                                                      >  1:437085/1‑251 (MQ=255)                                                                                            cctTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAAt                                                                                                                                                                                                 <  1:684724/208‑1 (MQ=255)                                                                                                tgtgAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATg                                                                                                                                                  <  1:588177/251‑1 (MQ=255)                                                                                                        gtgtCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATAtaataa                                                                                                                                                                                                                                 >  1:533576/1‑164 (MQ=255)                                                                                                            cAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGat                                                                                                                                                                                                                                        >  1:69799/1‑153 (MQ=255)                                                                                                                      aTCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTAATACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATAtaataa                                                                                                                                                                                                                                 >  1:656061/1‑150 (MQ=255)                                                                                                                               ttGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACgg                                                                                                                     <  1:517657/249‑1 (MQ=255)                                                                                                                                                                         gTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATataat                                                                                                                                                                                                                                   <  1:42235/97‑1 (MQ=255)                                                                                                                                                                                cccTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATc                                                                                                                                                                               >  1:246850/1‑142 (MQ=255)                                                                                                                                                                                cccTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTTTACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGAt                                                                  >  1:822744/1‑251 (MQ=255)                                                                                                                                                                                cccTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGAt                                                                  <  1:875810/251‑1 (MQ=255)                                                                                                                                                                                   ttGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGc                                                                                                                                                                                     >  1:788239/1‑133 (MQ=255)                                                                                                                                                                                          catcatGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATAc                                                                                                                                                                                                    <  1:89928/111‑1 (MQ=255)                                                                                                                                                                                          catcatGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGttttt                                                                                                                                      <  1:884144/173‑1 (MQ=255)                                                                                                                                                                                          catcatGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGATCGCTGTTGCa                                                        >  1:353962/1‑251 (MQ=255)                                                                                                                                                                                          catcatGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGATCGCTGTTGCa                                                        >  1:507269/1‑251 (MQ=255)                                                                                                                                                                                            tcatGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATa                                                                                                                                                                                                <  1:487660/113‑1 (MQ=255)                                                                                                                                                                                                              tattaCTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGATCGCTGTTGCACTTACCGGTTGTAACGATaa                                    >  1:50197/1‑251 (MQ=255)                                                                                                                                                                                                                ttaCTAGTCATTGCCAGTCAATTTATTGGCAAGAAACATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGATCGCTGTTGCACTTACCGGTTGTAACGATAAAg                                  >  1:460545/1‑251 (MQ=255)                                                                                                                                                                                                                                                 aaaCATAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTTACCGCAAGTATGATCGCTGTTGCACTTACCGGTTGTAACGATAAAGACAACAAGCAAGCGCAAGCAACGCCGGAAGcc  <  1:886257/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                     CGTAAAATCATCCTGCCGCTTGCCTTACCGGCGGTGGTGTCCGGTTGGTTGTTAAGTTTCACCATTTCACTGGATGACGTTGTGGTGTCCTCCTTTGTGAGCGGTGTCAGCTACGAAATCCTGCCGTTGCGTATCTTCTCCTTGGTGAAAACCGGAGTGACGCCGGAAGTGAATGCCCTTGCAACCATCATGATTGTGTTGTCATTATTACTAGTCATTGCCAGTCAATTTATTGGCAAGAAACACAAAGACATCTAGATATAATAATGCTTCACAAATAATTTAATAGCACATACAATACGCCTTGACGCAAAATCAAGGCGTTTTTTTATACCTGCATTTCATGCAGGTGGTTTTTTACCCCCTCTTTTACGGAGAACAAACAAAAATGAAAAAATTTGCCGGTTTATTCACCGCAAGTATGATCGCTGTTGCACTTACCGGTTGTAACGATAAAGACAACAAGCAAGCGCAAGCAACGCCGGAAGCC  >  1/1908881‑1909370

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