Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	93,358	A→C	S172A (TCG→GCA)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	93,358	0	A	C	100.0%	159.9 / NA	54	S172A (TCG→GCA)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (26/28);  total (26/28)

CTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAATGTACGGCGCGGCGACAGGCGGCGGGCAATGCGCTGCTCTGCGGCAGTGACGGATTGATGTACAGAGGTCAGCCACTGCAGGTTTTCTGTGACATCCTCGCTCCAATCCGGCATAAACGACCAATCAGTTGAGCGCGAGCCTATAATGCGCAGGGTAACTGCACTTTTACCCAAAAAATTAAAGGTAACCACACAATCAATGGTGGGCGCGCCTTGCATGCCGACGTTAACCGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCCACAACAGACACCGATAAAAGTTTGACCGCGCTGCGGTTGGCGTTCCAAATTT  >  1/93135‑93594                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             cTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATaa                                                                                                                                                                                                                    <  1:894178/250‑1 (MQ=255) cTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATaa                                                                                                                                                                                                                    <  1:418612/250‑1 (MQ=255) cTGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATaa                                                                                                                                                                                                                    <  1:36114/250‑1 (MQ=255)  tGTCGGTAAAAATCGGCATCGACCATCCGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATg                                                                                                                                                                                                                  <  1:883236/251‑1 (MQ=255)  tGTCGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATg                                                                                                                                                                                                                  <  1:742735/251‑1 (MQ=255)     cGGTAAAAATCGGCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATgcgc                                                                                                                                                                                                               <  1:524431/251‑1 (MQ=255)       gTAAAAATCGGCACCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATgcgc                                                                                                                                                                                                               <  1:392872/249‑1 (MQ=255)                ggCATCGACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAAccg                                                                                                                                                                                                     >  1:299155/1‑248 (MQ=255)                      gACCATACGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAAccg                                                                                                                                                                                                     >  1:130013/1‑242 (MQ=255)                         catacGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTAcc                                                                                                                                                                                           <  1:550401/248‑1 (MQ=255)                          atacGCGCGAGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTAccc                                                                                                                                                                                          <  1:599167/249‑1 (MQ=255)                                   aGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATaa                                                                                                                                                                                                                    >  1:248100/1‑216 (MQ=255)                                   aGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTa                                                                                                                                                                                 >  1:785298/1‑251 (MQ=255)                                   aGCCATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTa                                                                                                                                                                                 >  1:576117/1‑251 (MQ=255)                                     ccATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGGTGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAggg                                                                                                                                                                                                           <  1:607539/223‑1 (MQ=255)                                     ccATAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAAccg                                                                                                                                                                                                     >  1:40064/1‑227 (MQ=255)                                       aTAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAgg                                                                                                                                                                             <  1:152848/251‑1 (MQ=255)                                        tAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGGTTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAgg                                                                                                                                                                             >  1:614305/1‑250 (MQ=255)                                        tAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATaa                                                                                                                                                                                                                    <  1:207314/211‑1 (MQ=255)                                        tAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGt                                                                                                                                                                                                          >  1:852650/1‑221 (MQ=255)                                        tAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAgg                                                                                                                                                                             >  1:627490/1‑250 (MQ=255)                                        tAACCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGt                                                                                                                                                                            >  1:155892/1‑251 (MQ=255)                                          aCCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTaa                                                                                                                                                                          >  1:621698/1‑251 (MQ=255)                                          aCCATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTaa                                                                                                                                                                          >  1:879522/1‑251 (MQ=255)                                           ccATAGAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTATTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTaa                                                                                                                                                                          >  1:397328/1‑250 (MQ=255)                                                 aGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACaacac                                                                                                                                                                    <  1:545443/250‑1 (MQ=255)                                                   cgccgcTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACAcaat                                                                                                                                                                 <  1:190192/251‑1 (MQ=255)                                                     ccgcTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACAcaatca                                                                                                                                                               <  1:25156/251‑1 (MQ=255)                                                      cgcTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACAcaatcaa                                                                                                                                                              <  1:507922/251‑1 (MQ=255)                                                       gcTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAAcaac                                                                                                                                                                      >  1:325259/1‑239 (MQ=255)                                                         ttCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAAtgg                                                                                                                                                           <  1:396453/251‑1 (MQ=255)                                                            aaaGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAAtggtgg                                                                                                                                                        <  1:684616/251‑1 (MQ=255)                                                                         tCAACATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCAt                                                                                                                                           <  1:83562/251‑1 (MQ=255)                                                                            aCATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGcc                                                                                                                                        <  1:671204/251‑1 (MQ=255)                                                                            aCATCGCCGAAACTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGcc                                                                                                                                        <  1:479769/251‑1 (MQ=255)                                                                                     aaaCTTACTTTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAGTGGTGGGCGAGCCCTGCATGCCGACGTTAAc                                                                                                                               <  1:18224/251‑1 (MQ=255)                                                                                           acttTAAACTCAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTggg                                                                                                                                                       <  1:547628/221‑1 (MQ=255)                                                                                                     cAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAg                                                                                                                >  1:743282/1‑250 (MQ=255)                                                                                                     cAAACGTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAg                                                                                                                >  1:553731/1‑250 (MQ=255)                                                                                                         cgTACGGCGCGGCGACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAAtgg                                                                                                                                                           >  1:173771/2‑203 (MQ=255)                                                                                                               gcgcggcgACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAAtggtg                                                                                                                                                         >  1:350008/1‑199 (MQ=255)                                                                                                                      gACAAGCGGCGGGCAATGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAatat                                                                                              >  1:748436/1‑251 (MQ=255)                                                                                                                                    aaTGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGcc                                                                                <  1:689940/251‑1 (MQ=255)                                                                                                                                    aaTGCGCTGCTCTGCTGCCGTCACAGATTGATGCACGGAAGTCAGCCACTCCAGATTTTCGGTCACCTCTTGTGACCAGTCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGcc                                                                                <  1:56647/251‑1 (MQ=255)                                                                                                                                                                                                   cacctcttgtgaccagtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGt                 >  1:255976/17‑251 (MQ=255)                                                                                                                                                                                                     cctcttgtgaccagtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACtt                                                                      <  1:158018/182‑1 (MQ=255)                                                                                                                                                                                                     cctcttgtgaccagtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTg               >  1:292718/15‑251 (MQ=255)                                                                                                                                                                                                        cttgtgaccagtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTcc                                                                                                                          >  1:596024/12‑141 (MQ=255)                                                                                                                                                                                                        cttgtgaccagtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTGGCg            >  1:702135/12‑251 (MQ=255)                                                                                                                                                                                                               ccagtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTGGCggtgctta     >  1:823810/5‑244 (MQ=255)                                                                                                                                                                                                                  gtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTGGCGTTCCacactt  >  1:203569/2‑246 (MQ=255)                                                                                                                                                                                                                  gtCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTGGCGTTCCacact   >  1:412288/2‑246 (MQ=255)                                                                                                                                                                                                                   tCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTGGCGTTCCacact   <  1:218516/249‑5 (MQ=255)                                                                                                                                                                                                                   tCCGGCAGAAATGCCCAATCGGTTGAACGCGAGCCTATAATGCGCAGGGTAACCGCACTTTTACCCAAAAAATTAAAGGTAACAACACAATCAATGGTGGGCGAGCCCTGCATGCCGACGTTAACTGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCGACAACAGACACCGATAAAAGCTTGACCGCGCTGCGGTTGGCGTTCCacac    <  1:344951/248‑4 (MQ=255)                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             CTGTCGGTAAAAATCGGCATCGACCACACGCGCGAGCCATAACCATACAGCGCCGCTTCAAAGCGCTGACGCTCAACATCGCCGAAACTTACTTTAAACTCAAATGTACGGCGCGGCGACAGGCGGCGGGCAATGCGCTGCTCTGCGGCAGTGACGGATTGATGTACAGAGGTCAGCCACTGCAGGTTTTCTGTGACATCCTCGCTCCAATCCGGCATAAACGACCAATCAGTTGAGCGCGAGCCTATAATGCGCAGGGTAACTGCACTTTTACCCAAAAAATTAAAGGTAACCACACAATCAATGGTGGGCGCGCCTTGCATGCCGACGTTAACCGTCCATTTTTTAAGCGCCAGCGGGTTAAATATGCCTGATGTCGGCCCCTCTAACTTAATCCCCTCACCACCCACAACAGACACCGATAAAAGTTTGACCGCGCTGCGGTTGGCGTTCCAAATTT  >  1/93135‑93594

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