Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,986,078	C→T	R48K (AGG→AAG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,986,078	0	C	T	100.0%	200.9 / NA	71	R48K (AGG→AAG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/41);  total (30/41)

AAGTGCGGTCGTCTTACAAAACGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGAAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACTTTATTCATCGTGCCGCAACCGAAAATAATCTGATATTGTCCTGCTACGGCAAACTGCACTTTTACACCCTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTTGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTTCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGGCGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCCGTGCACCGAAT  >  1/1985834‑1986323                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      aaGTGCGGTCGTCTTACAAAACGTTTTTCATTATCAGAAAGGTGCCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAAtg                                                                                                                                                                                                                                                 <  1:346536/251‑1 (MQ=255) aaGTGCGGTCGTCTTACAAAACGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAAtg                                                                                                                                                                                                                                                 >  1:796887/1‑251 (MQ=255) aaGTGCGGTCGTCTTACAAAACGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAAtg                                                                                                                                                                                                                                                 <  1:656002/251‑1 (MQ=255) aaGTGCGGTCGTCTTACAAAACGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAAtg                                                                                                                                                                                                                                                 <  1:182422/251‑1 (MQ=255)        gtcgtcTTACAAAACGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCa                                                                                                                                                                                                                                           <  1:802281/250‑1 (MQ=255)                    aaCGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAAc                                                                                                                                                                                                                              <  1:778064/251‑1 (MQ=255)                       gTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACttt                                                                                                                                                                                                                           >  1:451681/1‑251 (MQ=255)                                         aGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGaa                                                                                                                                                                                                          >  1:297441/1‑250 (MQ=255)                                           tgccAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGaacaa                                                                                                                                                                                                       <  1:320466/249‑1 (MQ=255)                                                 ccTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACgcgc                                                                                                                                                                                                  >  1:239768/1‑250 (MQ=255)                                                         cAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCAt                                                                                                                                                                                          <  1:130361/250‑1 (MQ=255)                                                           gATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATgg                                                                                                                                                                                        >  1:132615/1‑250 (MQ=255)                                                            aTTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGgc                                                                                                                                                                                       <  1:399774/250‑1 (MQ=255)                                                                   cTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTggg                                                                                                                                                                               >  1:629246/1‑251 (MQ=255)                                                                   cTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTggg                                                                                                                                                                               >  1:740367/1‑251 (MQ=255)                                                                   cTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGCGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTggg                                                                                                                                                                               >  1:441028/1‑251 (MQ=255)                                                                     ttAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTa                                                                                                                                                                             >  1:802919/1‑251 (MQ=255)                                                                      tAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGt                                                                                                                                                                              <  1:676829/249‑1 (MQ=255)                                                                           gAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCtt                                                                                                                                                                                                                       >  1:651220/1‑203 (MQ=255)                                                                           gAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGaac                                                                                                                                                                                                         >  1:215402/1‑217 (MQ=255)                                                                                  tCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAatat                                                                                                                                                                <  1:50502/251‑1 (MQ=255)                                                                                                      tttGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACttt                                                                                                                                                                                                                           <  1:32886/172‑1 (MQ=255)                                                                                                      tttGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGgcgc                                                                                                                                                                                     <  1:834626/210‑1 (MQ=255)                                                                                                      tttGAAGGGAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGtttt                                                                                                                                            >  1:702419/1‑251 (MQ=255)                                                                                                            gggAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTg                                                                                                                                                                         >  1:191264/1‑216 (MQ=255)                                                                                                             ggAGTGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATaa                                                                                                                                     >  1:677390/1‑251 (MQ=255)                                                                                                                 tGGAAAGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGgc                                                                                                                                                                                       <  1:95511/197‑1 (MQ=255)                                                                                                                  ggAAAGACAGAAGGAAATAAAGAAAGCAGCATACCGGGCTATTTAATCAATTCGGCGTGGACCTTATTCATCGTTACGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACCCCTTCAAtg                                                                                                                                                                                                                                                 <  1:768614/138‑1 (MQ=255)                                                                                                                      aGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATAtttcct                                                                                                                                                           >  1:26748/1‑217 (MQ=255)                                                                                                                      aGACTGAAGGAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATAtttc                                                                                                                                                             >  1:627086/1‑217 (MQ=255)                                                                                                                            aaggaaaTAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGAttt                                                                                                                      <  1:151209/247‑1 (MQ=255)                                                                                                                                 aaTAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTaaga                                                                                                                 <  1:383584/251‑1 (MQ=255)                                                                                                                                 aaTAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTaaga                                                                                                                 <  1:142420/251‑1 (MQ=255)                                                                                                                                    aaagaaagCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATAtttcct                                                                                                                                                           >  1:735056/1‑203 (MQ=255)                                                                                                                                           gcagcaTACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACGCCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGa                                                                                                                         <  1:791399/233‑1 (MQ=255)                                                                                                                                               caTACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCAtatt                                                                                                      <  1:695531/248‑1 (MQ=255)                                                                                                                                                 tACCGGGCCATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCg                                                                                                 >  1:149741/1‑251 (MQ=255)                                                                                                                                                        ccATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAg                                                                                          <  1:869287/251‑1 (MQ=255)                                                                                                                                                        ccATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAg                                                                                          <  1:303686/251‑1 (MQ=255)                                                                                                                                                         cATTTAATCAATTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAgt                                                                                         >  1:67364/1‑251 (MQ=255)                                                                                                                                                                  aaTTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAaat                                                                                <  1:237041/251‑1 (MQ=255)                                                                                                                                                                  aaTTCGGCGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAaat                                                                                <  1:795451/251‑1 (MQ=255)                                                                                                                                                                         cGTGGACCTTATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGt                                                                         <  1:710246/251‑1 (MQ=255)                                                                                                                                                                             gaccttATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAat                                                                                                                                                                  >  1:87765/5‑158 (MQ=255)                                                                                                                                                                              accttATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGc                                                                                                                                                                                                                                       >  1:400115/4‑88 (MQ=255)                                                                                                                                                                              accttATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTcac                                                                                                                                                                     >  1:17084/4‑154 (MQ=255)                                                                                                                                                                              accttATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTaaga                                                                                                                 <  1:304805/203‑1 (MQ=255)                                                                                                                                                                                cttATTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTc                                                                  <  1:111413/250‑1 (MQ=255)                                                                                                                                                                                   aTTCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTctaa                                                               >  1:203718/1‑251 (MQ=255)                                                                                                                                                                                     tCATCGTTCCGCAACCGAAAATAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCaa                                                                                                                                        <  1:860712/176‑1 (MQ=255)                                                                                                                                                                                                      aaaaTAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAg                                             >  1:835646/1‑250 (MQ=255)                                                                                                                                                                                                       aaaTAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAAt                                                                    >  1:788255/1‑226 (MQ=255)                                                                                                                                                                                                       aaaTAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGcc                                           <  1:449686/251‑1 (MQ=255)                                                                                                                                                                                                       aaaTAATCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGcc                                           <  1:547760/251‑1 (MQ=255)                                                                                                                                                                                                            atCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAaat                                                                                <  1:530129/209‑1 (MQ=255)                                                                                                                                                                                                            atCTGATACTGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCgg                                      >  1:816334/1‑251 (MQ=255)                                                                                                                                                                                                                 atactGTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAaat                                                                                <  1:773625/200‑1 (MQ=255)                                                                                                                                                                                                                      gTCCTGCTACGGCAAATTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTctaact                                                             <  1:368371/218‑1 (MQ=255)                                                                                                                                                                                                                      gTCCTGCTACGGCAAACTGCACTTTTACACCTTCAAtgttgt                                                                                                                                                                                                                                             <  1:590114/42‑1 (MQ=255)                                                                                                                                                                                                                      gTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTctaact                                                             <  1:750162/218‑1 (MQ=255)                                                                                                                                                                                                                      gTCCTGCTACGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTg                             >  1:342508/1‑250 (MQ=255)                                                                                                                                                                                                                              aCGGCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTaaga                                                                                                                 <  1:606588/158‑1 (MQ=255)                                                                                                                                                                                                                                 gCAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACt                  <  1:854517/250‑1 (MQ=255)                                                                                                                                                                                                                                  cAAACTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTc                 <  1:66957/250‑1 (MQ=255)                                                                                                                                                                                                                                   aaaCTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCAt               >  1:798346/1‑251 (MQ=255)                                                                                                                                                                                                                                      cTGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGa                                                                                                                         <  1:339294/142‑1 (MQ=255)                                                                                                                                                                                                                                       tGCACTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATctgc           <  1:892912/251‑4 (MQ=255)                                                                                                                                                                                                                                           cTTTTACACCTTCAATGTTGTTAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAaataa                                                                              >  1:92332/1‑180 (MQ=255)                                                                                                                                                                                                                                           cTTTTACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATctgcgca        >  1:326309/1‑244 (MQ=255)                                                                                                                                                                                                                                             tttACACCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACaa                                                                                                                                                                   <  1:118668/93‑1 (MQ=255)                                                                                                                                                                                                                                                acacCTTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTCGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTCCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGACGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCTGCGCACCGAAt  <  1:575442/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      AAGTGCGGTCGTCTTACAAAACGTTTTTCATTATCAGAAAAGTACCAACCTGTTTGCAGATTCGCGCTTTAATGGAACAAATCTAAATCACATTTTTAGCATTTGAAGGGAGTGGAAAGACTGAAGAAAATAAAGAAAGCAGCATACCGGGCCATTTAATCAATTCGGCGTGGACTTTATTCATCGTGCCGCAACCGAAAATAATCTGATATTGTCCTGCTACGGCAAACTGCACTTTTACACCCTCAATGTTGTCAATGCCTTGTCAACTTTGCTTTTGTTATTCAGAACAACGCGCCAACGCATGGCGCAGTGGGTAGCTGTCACAATATTTTCTTTGTCGCCAAGTTTTTCAATAACCTGTTGGGCGATTTTAAGAAGATTCATATTATTCGCCGGAAGTGTCTTAAATAATAAGTAGAATTCTAACTAAAAGGCTTCGCTAAGCCAGCGGAATTTCGTGCACACTTAACTCATCCGTGCACCGAAT  >  1/1985834‑1986323

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