Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	652,503	+TT	coding (31/201 nt)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	652,500	1	.	T	100.0%	163.4 / NA	58	G12S (GGC→AGC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (21/37);  total (21/37)
*	1	652,500	2	.	T	100.0%	161.7 / NA	58	G12S (GGC→AGC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base T (21/37);  total (21/37)

AACTACCGCTTCCCTGCGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑TTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAGACTTTAGATTGTTC  >  1/652250‑652678                                                                                                                                                                                                                                                            ||                                                                                                                                                                                   aaCTACCGCTTCCCTGCGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGcc                                                                                                                                                                                      <  1:471715/251‑1 (MQ=255) aaCTACCGCTTCCCTGCGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGcc                                                                                                                                                                                      <  1:566431/251‑1 (MQ=255)             ccTGCGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑tttttgctgtat                                                                                                                                                                        >  1:576530/1‑242 (MQ=255)                 cGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATc                                                                                                                                                                       <  1:192842/250‑1 (MQ=255)                 cGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCt                                                                                                                                                                      <  1:117860/251‑1 (MQ=255)                  ggTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTatat                                                                                                                                                                        >  1:641767/1‑248 (MQ=255)                    ttCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTa                                                                                                                                                                   <  1:258130/251‑1 (MQ=255)                    ttCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTa                                                                                                                                                                   <  1:45064/251‑1 (MQ=255)                    ttCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTGAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTa                                                                                                                                                                   <  1:449923/251‑1 (MQ=255)                       ttGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGa                                                                                                                                                                <  1:884325/251‑1 (MQ=255)                       ttGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGa                                                                                                                                                                >  1:279360/1‑251 (MQ=255)                               ttCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCaaa                                                                                                                                                        >  1:441536/1‑251 (MQ=255)                                       gagTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTa                                                                                                                                                                   >  1:320028/1‑232 (MQ=255)                                       gagTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAAc                                                                                                                                                       >  1:417022/1‑244 (MQ=255)                                       gagTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑tttttgc                                                                                                                                                                             <  1:801904/220‑5 (MQ=255)                                         gTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGa                                                                                                                                                                <  1:705058/233‑1 (MQ=255)                                                    cTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGAt                                                                                                                                   <  1:554352/251‑1 (MQ=255)                                                     tACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATc                                                                                                                                  <  1:63327/251‑1 (MQ=255)                                                         ctctTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAAccc                                                                                                                                                     <  1:511104/228‑1 (MQ=255)                                                           ctTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑tttttgc                                                                                                                                                                             <  1:583522/200‑5 (MQ=255)                                                           ctTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑tttttgc                                                                                                                                                                             >  1:617990/1‑196 (MQ=255)                                                                        ttAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑tttttgc                                                                                                                                                                             >  1:171676/1‑183 (MQ=255)                                                                               aTTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAAcc                                                                                                                                                      >  1:21330/1‑205 (MQ=255)                                                                               aTTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAAcc                                                                                                                                                      <  1:467669/205‑1 (MQ=255)                                                                               aTTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑tttttgctg                                                                                                                                                                           <  1:386573/182‑7 (MQ=255)                                                                                         gagaAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTaca                                                                                              >  1:499348/1‑251 (MQ=255)                                                                                         aagaAAACCTTTTAGATTTTTGGGTTTTTTGTT‑CTTAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTacac                                                                                             <  1:238228/250‑1 (MQ=25)                                                                                            aaaaaCCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTacac                                                                                             <  1:523780/249‑1 (MQ=255)                                                                                                 ccTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTa                                                                                       >  1:268703/1‑250 (MQ=255)                                                                                                  cTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTAcc                                                                                     <  1:555257/251‑1 (MQ=255)                                                                                                                 cTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGa                                                                       <  1:728314/250‑1 (MQ=255)                                                                                                                   tgttgtTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATc                                                                                                                                                                       <  1:605561/152‑1 (MQ=255)                                                                                                                   tgttgtTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAagag                                                                                                                                             >  1:165102/1‑178 (MQ=255)                                                                                                                       gttCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGa                                                                >  1:28098/1‑251 (MQ=255)                                                                                                                        ttCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGa                                                                >  1:784820/1‑250 (MQ=255)                                                                                                                               gaTTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAagag                                                                                                                                             <  1:740731/166‑1 (MQ=255)                                                                                                                               gaTTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACtt                                                        >  1:290116/1‑251 (MQ=255)                                                                                                                                           aTTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATtc                                            <  1:655185/251‑1 (MQ=255)                                                                                                                                                     cTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATtgaat                                  <  1:172036/251‑1 (MQ=255)                                                                                                                                                          gCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTacac                                                                                             <  1:104458/187‑1 (MQ=255)                                                                                                                                                          gCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATtc                                            <  1:149797/236‑1 (MQ=255)                                                                                                                                                          gCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTAt                                       >  1:757942/1‑241 (MQ=255)                                                                                                                                                            gTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGAt                                                                      <  1:98723/208‑1 (MQ=255)                                                                                                                                                              tAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGc                                                                                                                                                           >  1:496230/1‑121 (MQ=255)                                                                                                                                                               aaTTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGCTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTATATtg                                     >  1:352830/1‑238 (MQ=255)                                                                                                                                                                  tACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTAttt                     <  1:716630/251‑1 (MQ=255)                                                                                                                                                                       ttGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAg                <  1:873202/251‑1 (MQ=255)                                                                                                                                                                          cTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGcc                                                                                                    >  1:64728/1‑164 (MQ=255)                                                                                                                                                                          cTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGAt                                                                      >  1:496060/1‑194 (MQ=255)                                                                                                                                                                                ggCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAGACTCTAGa        <  1:635092/250‑1 (MQ=255)                                                                                                                                                                                 gCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAGACTTTAGAt       <  1:528182/250‑1 (MQ=255)                                                                                                                                                                                       gACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAGACTTTAGAt       <  1:80638/244‑1 (MQ=255)                                                                                                                                                                                        aCTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCTAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAg                                                                         <  1:255749/177‑1 (MQ=255)                                                                                                                                                                                                 tCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTacac                                                                                             <  1:415397/148‑1 (MQ=255)                                                                                                                                                                                                    ttATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTatat                                                                                                                                                                        <  1:98486/70‑1 (MQ=255)                                                                                                                                                                                                    ttATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAg                <  1:787843/222‑1 (MQ=255)                                                                                                                                                                                                     tATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGcc                                                                                                    >  1:375493/1‑137 (MQ=255)                                                                                                                                                                                                      aTCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTGTGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATtg                                     >  1:445615/1‑199 (MQ=255)                                                                                                                                                                                                      aTCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTAcc                                                                                     <  1:6973/151‑1 (MQ=255)                                                                                                                                                                                                      aTCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAg                                                                         >  1:757147/1‑163 (MQ=255)                                                                                                                                                                                                      aTCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATc                          <  1:189041/210‑1 (MQ=255)                                                                                                                                                                                                              gTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATc                                                                                                                                                              <  1:588889/70‑1 (MQ=255)                                                                                                                                                                                                              gTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATct                                                                                                                         <  1:543441/107‑1 (MQ=255)                                                                                                                                                                                                              gTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTacac                                                                                             <  1:773186/135‑1 (MQ=255)                                                                                                                                                                                                              gTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAg                <  1:191017/212‑1 (MQ=255)                                                                                                                                                                                                                                                    gTTATGCCTTTTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGGTAATTCTCTATTGAATGAAAAATCTATTTTATAGACTTTAGATTGTTc  >  1:312355/1‑188 (MQ=255)                                                                                                                                                                                                                                                            ||                                                                                                                                                                                   AACTACCGCTTCCCTGCGGTTCTTGCTGAATTCTCAAAGAGTAAAAACAGACTACGCTCTTTTTTTGCTCTTTAGACTATTTTTTATAGAGAAAAACCTTTAGATTTTTGAGCTTGTTGTTCCTGAGATTTAGGTTTAATTTTAGAACCTACGGCGTTAATTACTCTTGCTTCAGGGCGATTGACTTTTAATTCATTATCTAGCTGTTTAGGCTTGATAGCGCCAGATTTCAAGGCTTCTGCGGTTATGCC‑‑TTTGCTGTATATCTTTATGATCAGCAAACCCATCAAGAGCCTTTGAGATCTGGCTATCTTGCTTGATAATTTCATTAGCCTTTACACCGTTTACCAGAGCATCATAGGATTTTTGACCTAACTTTTCGGATAATTCTCTATTGAATGAAAAATCTATTTTATAGACTTTAGATTGTTC  >  1/652250‑652678

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: ‑