Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,089,079	C→T	P149L (CCG→CTG)	unknown →	Inner membrane transport permease YhhJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,089,079	0	C	T	100.0%	187.6 / NA	69	P149L (CCG→CTG)	unknown	Inner membrane transport permease YhhJ
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/40);  total (29/40)

CAACAAAGCCTGATTCCGCCGAACTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCGCCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCCACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCCGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACGGGAGCGCGGCACGCTGGAACATATTTTGGTGATGCCCGTACGCTCCA  >  1/1088835‑1089285                                                                                                                                                                                                                                                     |                                                                                                                                                                                                               caacaaAGCCTGATTCCGCCGAACTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCAc                                                                                                                                                                                                            <  1:670326/249‑1 (MQ=255)        gCCTGATTCCGCCGAACTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAg                                                                                                                                                                                                    >  1:673450/1‑250 (MQ=255)         ccTGATTCCGCCGAACTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGc                                                                                                                                                                                                   <  1:168582/250‑1 (MQ=255)              ttCCGCCGAACTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAAc                                                                                                                                                                                             <  1:131833/251‑1 (MQ=255)                      aaCTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACt                                                                                                                                                                                     <  1:844904/251‑1 (MQ=255)                          tCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAg                                                                                                                                                                                                    <  1:771660/232‑1 (MQ=255)                          tCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCa                                                                                                                                                                                         >  1:239688/1‑243 (MQ=255)                               aaaGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAcc                                                                                                                                                                            <  1:317951/251‑1 (MQ=255)                                 aGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACgg                                                                                                                                                                                           <  1:811215/234‑1 (MQ=255)                                  gTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTc                                                                                                                                                                                <  1:687576/244‑1 (MQ=255)                                  gTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAc                                                                                                                                                                             >  1:456730/1‑247 (MQ=255)                                  gTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACcca                                                                                                                                                                          >  1:7486/1‑250 (MQ=255)                                  gTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACcca                                                                                                                                                                          >  1:38761/1‑250 (MQ=255)                                        aCAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCGCCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTaa                                                                                                                                                                   <  1:157977/251‑1 (MQ=255)                                         cAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTaa                                                                                                                                                                   <  1:737075/250‑1 (MQ=255)                                          aaTTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCGATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAAtg                                                                                                                                                                 <  1:145630/251‑1 (MQ=255)                                            ttCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAc                                                                                                                                                                             <  1:133334/237‑1 (MQ=255)                                             tCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGc                                                                                                                                                              <  1:84154/251‑1 (MQ=255)                                               aTGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAcc                                                                                                                                                                            >  1:150973/1‑235 (MQ=255)                                                       caGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAAcc                                                                                                                                                     >  1:808992/1‑250 (MQ=255)                                                                  ccTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAcc                                                                                                                                                                            >  1:379171/1‑216 (MQ=255)                                                                     gATGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTa                                                                                                                                      >  1:819625/1‑251 (MQ=255)                                                                       tGGACACAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAAc                                                                                                                                    <  1:104594/251‑1 (MQ=255)                                                                          acacAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTc                                                                                                                                                                                    <  1:326392/200‑1 (MQ=255)                                                                          acacAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTgg                                                                                                                                 >  1:786421/1‑251 (MQ=255)                                                                          acacAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTgg                                                                                                                                 >  1:41641/1‑251 (MQ=255)                                                                          acacAGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTgg                                                                                                                                 >  1:591876/1‑251 (MQ=255)                                                                             caGGAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGttt                                                                                                                              >  1:726217/1‑251 (MQ=255)                                                                               ggAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAcc                                                                                                                                                                            <  1:585140/203‑1 (MQ=255)                                                                               ggAGAATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTAt                                                                                                                            >  1:32574/1‑251 (MQ=255)                                                                                  gaATATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATggg                                                                                                                         >  1:650771/1‑251 (MQ=255)                                                                                   aaTATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAg                                                                                                                                                                                                    >  1:525723/1‑175 (MQ=255)                                                                                   aaTATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATggg                                                                                                                         <  1:498288/250‑1 (MQ=255)                                                                                   aaTATACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATggg                                                                                                                         <  1:315666/250‑1 (MQ=255)                                                                                    atatACGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTtat                                                                                                                                                         <  1:666591/217‑1 (MQ=255)                                                                                      ataCGTTCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTc                                                                                                                                           >  1:770648/1‑229 (MQ=255)                                                                                            tCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAAc                                                                                                               <  1:187706/251‑1 (MQ=255)                                                                                            tCATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAAc                                                                                                               <  1:629203/251‑1 (MQ=255)                                                                                             cATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAg                                                                                                                                                                                                    >  1:565764/1‑165 (MQ=255)                                                                                             cATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACa                                                                                                              <  1:112702/251‑1 (MQ=255)                                                                                             cATTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACa                                                                                                              <  1:800254/251‑1 (MQ=255)                                                                                              aTTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAccc                                                                                                                                                                           >  1:289409/1‑189 (MQ=255)                                                                                              aTTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACAt                                                                                                             >  1:472676/1‑251 (MQ=255)                                                                                              aTTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACAt                                                                                                             >  1:268483/1‑251 (MQ=255)                                                                                              aTTTTAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACAt                                                                                                             <  1:764495/251‑1 (MQ=255)                                                                                               ttttAGATATTCCACCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACgg                                                                                                                                                                                           <  1:374735/172‑1 (MQ=255)                                                                                                                  ttATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACa                                                                                                              >  1:129695/1‑230 (MQ=255)                                                                                                                  ttATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATttat                                                                                         >  1:560/1‑251 (MQ=255)                                                                                                                       aaCGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAAcc                                                                                    <  1:837937/251‑1 (MQ=255)                                                                                                                       aaCGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAAcc                                                                                    <  1:468899/251‑1 (MQ=255)                                                                                                                       aaCGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAAc                                                                                     <  1:855194/250‑1 (MQ=255)                                                                                                                       aaCGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAAc                                                                                     <  1:280815/250‑1 (MQ=255)                                                                                                                         cGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTgg                                                                                                        >  1:383334/1‑229 (MQ=255)                                                                                                                         cGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTg                                                                                                <  1:651870/237‑1 (MQ=255)                                                                                                                          gTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATAc                                                                                                                                                                                      <  1:462951/150‑1 (MQ=255)                                                                                                                              tGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGc                                                                             <  1:128838/251‑1 (MQ=255)                                                                                                                                 aCTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCtgg                                                                          <  1:760814/251‑1 (MQ=255)                                                                                                                                 aCTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCtgg                                                                          <  1:604460/251‑1 (MQ=255)                                                                                                                                          gCGTAACCCGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATggg                                                                                                                         <  1:792106/195‑1 (MQ=255)                                                                                                                                                  cGCTGATTCAGTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATATCCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCAtt                                                         <  1:688652/251‑1 (MQ=255)                                                                                                                                                            gTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATAc                                                                                                                                                                                      >  1:760517/1‑116 (MQ=255)                                                                                                                                                            gTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACggg                                               >  1:891899/1‑251 (MQ=255)                                                                                                                                                            gTTATTATCCGACGCGACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACggg                                               >  1:412280/1‑251 (MQ=255)                                                                                                                                                                                cAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACGGGAGCGCGGCACGCTGGAACat                           <  1:368723/251‑1 (MQ=255)                                                                                                                                                                                                      ttCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACGGGAGCGCGGCACGCTGGAACATATTTTGGTGATGCCCGTACGc      <  1:686870/250‑1 (MQ=255)                                                                                                                                                                                                        ccGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACGGGAGCGCGGCACGCTGGAACATATTTTGGTGATGCCCGTACGCTcc   <  1:866219/251‑1 (MQ=255)                                                                                                                                                                                                         cGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACGGGAGCGCGGCACGCTGGAACATATTTTGGTGATGCCCGTACGCTCCa  >  1:496330/1‑251 (MQ=255)                                                                                                                                                                                                          gCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAAt                                                                                             >  1:780674/1‑159 (MQ=255)                                                                                                                                                                                                                                        aTTAACGATTTTCTGCACCTTCAAAGCAATAACGGCATACTCATTCAAc                                                                                                                                                                             <  1:564400/49‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                               CAACAAAGCCTGATTCCGCCGAACTTCCGTAAAGTGCATACAATTCATGCACAACAGGCGGATCGCCTGATGGACACAGGAGAATATACGTTCATTTTAGATATTCCGCCGAATTATGAACGTGATGTACTTGCGGGGCGTAACCCGCTGATTCAGTTATTATCCGACGCCACGGCAATGACGCAAGCCGCCATCGGTTCCGCTTATATTTCACAAATTTTTCACGGTGAAATTAACGATTTTCCGCACCTTCAAAGCAATAACGGCATACTCATTCAACCCACCATTAATGTGCTTTATAACCCGAATTTTTCCGGTAACTGGTTTATGGGCGGCATGAACATTGTGGGCTATTTGAATTTATTAACCATGTTGCTGGTGGGCGCCGCTGTCATTCGGGAACGGGAGCGCGGCACGCTGGAACATATTTTGGTGATGCCCGTACGCTCCA  >  1/1088835‑1089285

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