Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,541,815	T→C	L218L (TTA→TTG)	unknown ←	Putative DNA‑binding protein in cluster with Type I restriction‑modification system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,541,815	0	T	C	100.0%	200.1 / NA	70	L218L (TTA→TTG)	unknown	Putative DNA‑binding protein in cluster with Type I restriction‑modification system
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (27/43);  total (27/43)

AAAGTATTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAATAACGCCGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTAT  >  1/1541570‑1542057                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   aaaGTATTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGc                                                                                                                                                                                                                                               <  1:136495/251‑1 (MQ=255)      aTTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCg                                                                                                                                                                                                                                          <  1:879739/251‑1 (MQ=255)             ttAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTcccc                                                                                                                                                                                                                                   <  1:901150/251‑1 (MQ=255)                  cAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCttttt                                                                                                                                                                                                                              >  1:596048/1‑251 (MQ=255)                       gTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCa                                                                                                                                                                                                                           >  1:587726/1‑249 (MQ=255)                           tACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACg                                                                                                                                                                                                                     <  1:201599/251‑1 (MQ=255)                                  cTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGACGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGttt                                                                                                                                                                                                              <  1:456341/251‑1 (MQ=255)                                     gCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTaa                                                                                                                                                                                                           >  1:810608/1‑251 (MQ=255)                                          tAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGTCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTc                                                                                                                                                                                                                             <  1:451803/228‑1 (MQ=255)                                             gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGTCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAAt                                                                                                                                                                                                          <  1:344599/244‑1 (MQ=255)                                             gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                   >  1:335900/1‑251 (MQ=255)                                             gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                   >  1:78679/1‑251 (MQ=255)                                             gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                   >  1:623770/1‑251 (MQ=255)                                              tgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAAt                                                                                                                                                                                                  <  1:249776/251‑1 (MQ=255)                                                         ccGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATa                                                                                                                                                                                       >  1:698017/1‑251 (MQ=255)                                                            cAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGt                                                                                                                                                                                    <  1:3103/251‑1 (MQ=255)                                                             aaaTACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTc                                                                                                                                                                                   <  1:832347/251‑1 (MQ=255)                                                             aaaTACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGCCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTc                                                                                                                                                                                   <  1:672193/251‑1 (MQ=255)                                                                 aCCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCt                                                                                                                                                                               <  1:612487/251‑1 (MQ=255)                                                                  ccATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTa                                                                                                                                                                              >  1:456425/1‑251 (MQ=255)                                                                      aaaaaCAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAacac                                                                                                                                                                          <  1:552502/251‑1 (MQ=255)                                                                      aaaaaCAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAacac                                                                                                                                                                          <  1:417930/251‑1 (MQ=255)                                                                      aaaaaCAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAacac                                                                                                                                                                          <  1:213860/251‑1 (MQ=255)                                                                           cAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACa                                                                                                                                                                     >  1:489334/1‑251 (MQ=255)                                                                           cAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACa                                                                                                                                                                     >  1:652429/1‑251 (MQ=255)                                                                              gAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAg                                                                                                                                                                    <  1:463568/249‑1 (MQ=255)                                                                               aaaGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAAc                                                                                                                                                                 <  1:767028/251‑1 (MQ=255)                                                                                       cGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCAtattat                                                                                                                                                         <  1:654778/251‑1 (MQ=255)                                                                                           cGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                       >  1:23182/1‑249 (MQ=255)                                                                                           cGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTg                                                                                                                                                     >  1:891714/1‑251 (MQ=255)                                                                                            gTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTg                                                                                                                                                     >  1:313659/1‑250 (MQ=255)                                                                                                  tGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGACGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGc                                                                                                                                              <  1:497704/251‑1 (MQ=255)                                                                                                    ccATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGc                                                                                                                                                                                                                    <  1:406308/179‑1 (MQ=255)                                                                                                    ccATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAAc                                                                                                                                                                 >  1:466187/1‑230 (MQ=255)                                                                                                              aCGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTGATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGtt                                                                                                                                  >  1:569589/1‑251 (MQ=255)                                                                                                                               ataaAATACAATAAATTCGCAGCTCGGCTAGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTGCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCGCTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                       <  1:185487/213‑1 (MQ=255)                                                                                                                                   aaTACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGTAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGcc                                                                                                             <  1:453119/251‑1 (MQ=255)                                                                                                                                   aaTACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGcc                                                                                                             <  1:634262/251‑1 (MQ=255)                                                                                                                                       cAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAAc                                                                                                                                                                      <  1:215737/190‑1 (MQ=255)                                                                                                                                        aaTAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                   >  1:706422/1‑160 (MQ=255)                                                                                                                                             aTTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAaac                                                                                                   <  1:877711/251‑1 (MQ=255)                                                                                                                                                  ccgcTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCg                                                                                                                                    <  1:467632/213‑1 (MQ=255)                                                                                                                                                     cTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAaactccaactc                                                                                           <  1:796490/251‑1 (MQ=255)                                                                                                                                                               aaCATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAAt                                                                                 <  1:24877/251‑1 (MQ=255)                                                                                                                                                                 cATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATc                                                                                >  1:768504/1‑250 (MQ=255)                                                                                                                                                                       gATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGcac                                                                          <  1:848744/250‑1 (MQ=255)                                                                                                                                                                          aaaaaaGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCacaaca                                                                      >  1:274131/1‑251 (MQ=255)                                                                                                                                                                              aaGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAAc                                                                                                                                                                      <  1:703540/151‑1 (MQ=255)                                                                                                                                                                                  gtgtTACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGAtt                                                                                                                                         <  1:21824/176‑1 (MQ=255)                                                                                                                                                                                     tCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGc                                                           <  1:269498/251‑1 (MQ=255)                                                                                                                                                                                            aTCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATcc                                                                               >  1:543255/1‑224 (MQ=255)                                                                                                                                                                                             tCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAGCTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCtt                                                   >  1:527664/1‑251 (MQ=255)                                                                                                                                                                                             tCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCacaac                                                                       <  1:890863/231‑1 (MQ=255)                                                                                                                                                                                                ttGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTg                                                                                                                                                                                                        >  1:898018/1‑99 (MQ=255)                                                                                                                                                                                                ttGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCg                                                                                                                                    >  1:840134/1‑167 (MQ=255)                                                                                                                                                                                                ttGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGt                                                >  1:792670/1‑251 (MQ=255)                                                                                                                                                                                                 tGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTc                                               <  1:272966/251‑1 (MQ=255)                                                                                                                                                                                                       aTCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAAc                                                                                                       <  1:569253/189‑1 (MQ=255)                                                                                                                                                                                                            ggCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGc                                                           <  1:210457/228‑1 (MQ=255)                                                                                                                                                                                                             gCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAg                                    <  1:710472/250‑1 (MQ=255)                                                                                                                                                                                                                aTGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAgattg                                <  1:763339/251‑1 (MQ=255)                                                                                                                                                                                                                  ggCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTccc                                                                                                                                                                                                                                    <  1:547567/53‑1 (MQ=255)                                                                                                                                                                                                                        ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCTCTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGAc                                                                                                                              >  1:119044/1‑149 (MQ=255)                                                                                                                                                                                                                        ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                       <  1:289670/124‑1 (MQ=255)                                                                                                                                                                                                                        ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCa                        >  1:153483/1‑251 (MQ=255)                                                                                                                                                                                                                        ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCa                        >  1:318235/1‑251 (MQ=255)                                                                                                                                                                                                                             aaCCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTc                   <  1:277146/251‑1 (MQ=255)                                                                                                                                                                                                                                      ttCGTTGTTGGCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTTATGTAACACCAACAGAACCATAGTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACACCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGc          <  1:431641/251‑1 (MQ=255)                                                                                                                                                                                                                                               gtCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGt                                                                                                                                                                                    <  1:363127/72‑1 (MQ=255)                                                                                                                                                                                                                                               gtCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTAt  >  1:74566/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                   AAAGTATTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAATAACGCCGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTAT  >  1/1541570‑1542057

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