Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	516,290	T→C	D186G (GAC→GGC)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	516,290	0	T	C	100.0%	211.5 / NA	74	D186G (GAC→GGC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/43);  total (31/43)

GCGGCGGATTGTCCAATTTCGATTATTTATACGGCGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGTCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCAAATAAATCGCACCTTTTGGCGCATCATAATAA  >  1/516048‑516529                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                gcggcggATTGTCCAATTTCAATAATTTATACGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTc                                                                                                                                                                                                                                          <  1:832750/250‑1 (MQ=255)     cggATTGTCCAATTTCGATTATTTATACGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATc                                                                                                                                                                                                                                     <  1:619144/251‑1 (MQ=255)                        tatttatACGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAAt                                                                                                                                                                                                                  >  1:752291/1‑251 (MQ=255)                         atttatACGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATg                                                                                                                                                                                                                 >  1:393209/1‑251 (MQ=255)                         atttatACGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATg                                                                                                                                                                                                                 >  1:883919/1‑251 (MQ=255)                          tttatACGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATg                                                                                                                                                                                                                 >  1:351456/1‑250 (MQ=255)                             ataCGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATc                                                                                                                                                                                                                                     <  1:58399/227‑1 (MQ=255)                             ataCGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCtt                                                                                                                                                                                                             >  1:211432/1‑251 (MQ=255)                             ataCGGTGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCtt                                                                                                                                                                                                             >  1:515682/1‑251 (MQ=255)                                  gtgcCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCg                                                                                                                                                                                                        >  1:47241/3‑251 (MQ=255)                                  gtgcCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCg                                                                                                                                                                                                        >  1:679210/3‑251 (MQ=255)                                   tgcCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCgg                                                                                                                                                                                                       <  1:90914/250‑1 (MQ=255)                                          ccGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGttt                                                                                                                                                                                                <  1:233213/251‑1 (MQ=255)                                              aaGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAatc                                                                                                                                                                                            <  1:727821/251‑1 (MQ=255)                                                       cgccgcATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATcc                                                                                                                                                                                   <  1:668233/251‑1 (MQ=255)                                                            cATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATcatcat                                                                                                                                                                                       <  1:893459/242‑1 (MQ=255)                                                            cATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTgg                                                                                                                                                                              >  1:753780/1‑251 (MQ=255)                                                            cATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTgg                                                                                                                                                                              >  1:473144/1‑251 (MQ=255)                                                            cATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTgg                                                                                                                                                                              >  1:233793/1‑251 (MQ=255)                                                                         gcgcCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTAccc                                                                                                                                                                                                                      <  1:364913/198‑1 (MQ=255)                                                                         gcgcCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGt                                                                                                                                                                   >  1:234980/1‑249 (MQ=255)                                                                                 gTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATcc                                                                                                                                                                                   >  1:512148/1‑225 (MQ=255)                                                                                         gATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCg                                                                                                                                                                                 >  1:144951/1‑219 (MQ=255)                                                                                                   gCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATcc                                                                                                                                                                                   <  1:533176/207‑1 (MQ=255)                                                                                                       aaTATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGGTTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAg                                                                                                                                   >  1:513720/1‑251 (MQ=255)                                                                                                       aaTATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGAtttttt                                                                                                                                                                       <  1:285082/215‑1 (MQ=255)                                                                                                              cGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGA‑TTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGat                                                                                                                           <  1:306604/251‑1 (MQ=255)                                                                                                                 cTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATcc                                                                                                                                                                                   <  1:642530/193‑1 (MQ=255)                                                                                                                 cTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGatat                                                                                                                         >  1:209070/1‑251 (MQ=255)                                                                                                                 cTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGatat                                                                                                                         >  1:552400/1‑251 (MQ=255)                                                                                                                 cTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTCAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGatat                                                                                                                         >  1:477058/1‑251 (MQ=255)                                                                                                                    cggcggCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCTTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATAttgt                                                                                                                      <  1:438949/251‑1 (MQ=255)                                                                                                                           gtgcgtgGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATAttg                                                                                                                       <  1:325803/243‑1 (MQ=255)                                                                                                                             gcgtgGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAatca                                                                                                                                                                                           <  1:355049/173‑1 (MQ=255)                                                                                                                                  gcgcCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACgg                                                                                                         <  1:407181/250‑1 (MQ=255)                                                                                                                                     ccgcGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATcc                                                                                                                                                                                   <  1:396939/173‑1 (MQ=255)                                                                                                                                          gcgTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAAcc                                                                                                                                                     >  1:346520/1‑198 (MQ=255)                                                                                                                                          gcgTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCATGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGa                                                                                                  >  1:455626/1‑249 (MQ=255)                                                                                                                                               gttTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGTCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATAttg                                                                                                                       <  1:429999/223‑1 (MQ=255)                                                                                                                                               gttTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAAcc                                                                                                                                                     <  1:664259/193‑1 (MQ=255)                                                                                                                                               gttTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATa                                                                                                                                    >  1:569890/1‑210 (MQ=255)                                                                                                                                               gttTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGAc                                                                                                 >  1:138355/1‑245 (MQ=255)                                                                                                                                               gttTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGaa                                                                                           >  1:374081/1‑251 (MQ=255)                                                                                                                                                tttttATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATaaa                                                                                                                                                       <  1:814340/190‑1 (MQ=255)                                                                                                                                                tttttATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAAt                                                                                          >  1:875970/1‑251 (MQ=255)                                                                                                                                                      tCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGccc                                                                                    <  1:633777/251‑1 (MQ=255)                                                                                                                                                         cGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCGTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGGTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTgg                                                                                                                   <  1:126417/217‑1 (MQ=255)                                                                                                                                                             aaTGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTAt                                                                             <  1:417483/251‑1 (MQ=255)                                                                                                                                                             aaTGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTAt                                                                             <  1:443886/251‑1 (MQ=255)                                                                                                                                                             aaTGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTAt                                                                             <  1:374494/251‑1 (MQ=255)                                                                                                                                                              aTGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAAtt                                                                                         >  1:796052/1‑238 (MQ=255)                                                                                                                                                                  ttAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTc                                                                        <  1:826730/251‑1 (MQ=255)                                                                                                                                                                        cAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTgg                                                                                                                               <  1:219978/190‑1 (MQ=255)                                                                                                                                                                        cAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCg                                                                   <  1:617928/250‑1 (MQ=255)                                                                                                                                                                          aTAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAatca                                                                                                                                                                                           <  1:631415/128‑1 (MQ=255)                                                                                                                                                                           aaaaaaaGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGc                                                                >  1:173408/2‑250 (MQ=255)                                                                                                                                                                            aaaaaaGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAAt                                                                                          >  1:3996/1‑223 (MQ=255)                                                                                                                                                                                aaGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGAcc                                                          <  1:38817/251‑1 (MQ=255)                                                                                                                                                                                       aaaccTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCgtag                                                   <  1:862262/248‑1 (MQ=255)                                                                                                                                                                                         tCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCgtagta                                                 <  1:728008/251‑1 (MQ=255)                                                                                                                                                                                          ccTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTAc                                                <  1:373835/251‑1 (MQ=255)                                                                                                                                                                                           cTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTAcg                                               <  1:746665/251‑1 (MQ=255)                                                                                                                                                                                            ttGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTAcg                                               >  1:357925/1‑250 (MQ=255)                                                                                                                                                                                                    ccTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATcc                                      <  1:834587/251‑1 (MQ=255)                                                                                                                                                                                                    ccTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATcc                                      <  1:514111/251‑1 (MQ=255)                                                                                                                                                                                                         aTTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCaa                                 >  1:322386/1‑251 (MQ=255)                                                                                                                                                                                                                           ggTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAatat                                        <  1:436139/226‑1 (MQ=255)                                                                                                                                                                                                                           ggTTTTAGCAAATTATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGGGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTAc                                                >  1:321579/1‑218 (MQ=255)                                                                                                                                                                                                                                       ttATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCAAATAAATCGCACCt                    <  1:12272/234‑1 (MQ=255)                                                                                                                                                                                                                                        tATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAat                                          <  1:734669/211‑1 (MQ=255)                                                                                                                                                                                                                                        tATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCAAATAAATCGCACCTTTTGGCGCATCAtaataa  <  1:2357/251‑1 (MQ=255)                                                                                                                                                                                                                                        tATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCAAATAAATCGCACCTTTTGGCGCATCAtaataa  <  1:657570/251‑1 (MQ=255)                                                                                                                                                                                                                                        tATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCAAATAAATCGCACCTTTTGGCGCATCAtaataa  <  1:833664/251‑1 (MQ=255)                                                                                                                                                                                                                                        tATCAAAGGTGCCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGCGGATTTTTTTCCTCCAGGGTATAAACCGGCGTGTTGCTCAGATAGGAGGCGCTATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCaaataaat                           >  1:631753/1‑226 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                GCGGCGGATTGTCCAATTTCGATTATTTATACGGCGCCTTACCGAAAGCTTTACCGCCGCATTTAATGCGCAGCGCCAAAGTGCCGGTGATTAAAAAAGCCAAATATGGCGACTCCGGCGGCGTGCGTGGCGCCGCGGCGTTGTTTTTATCCCGTTAATGTTTAACGCAATAAAAAAGGCTGAATCCTTGTTCAGCCTTTATTTTTGCCTGAAATTATGGTTTTAGCAAATTATCAAAGGTGTCCACTTCCAATCTGCCTTGTTCTACCCGAATGCCTTTAGCGGCTTGTTTAATCATCATATCCCGTGGATTTTTTTCGTCTAGGGTATAAACCGGCGTGTTGCTCAGATAGGTGGCGATATTGTTGGTGAGTAACGGCACCAGACTTTGAATTTGCCCCATGTATTGTTCAGGCGAGCCGGACCAGCGTAGTACGCGAATATCCTTCAAATAAATCGCACCTTTTGGCGCATCATAATAA  >  1/516048‑516529

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