Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,620,830	G→A	G632D (GGC→GAC)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,620,830	0	G	A	100.0%	198.0 / NA	70	G632D (GGC→GAC)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/40);  total (30/40)

TGGTTAAACGTTTAACTACAGAAATCGCATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGGCAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTGAGTGGAAAAGATACCCTCACTAAAAT  >  1/1620587‑1621061                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        tgGTTAAACGTTTAACTACAGAAATCGCATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAAc                                                                                                                                                                                                                                   <  1:880952/250‑1 (MQ=255)  ggTTAAACGTTTAACTACAGAAATCGCATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACaa                                                                                                                                                                                                                                 >  1:413362/1‑251 (MQ=255)        aCGTTTAACTACAGAAATCGCATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGttt                                                                                                                                                                                                                             >  1:96807/1‑249 (MQ=255)                      aaaTCGCATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCActct                                                                                                                                                                                                                 <  1:851606/247‑1 (MQ=255)                            cATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTtg                                                                                                                                                                                                               <  1:575870/243‑1 (MQ=255)                            cATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTAAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGtt                                                                                                                                                                                                       >  1:437123/1‑251 (MQ=255)                               cGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTgg                                                                                                                                                                                                    >  1:97618/1‑251 (MQ=255)                               cGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTgg                                                                                                                                                                                                    >  1:572711/1‑251 (MQ=255)                               cGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTgg                                                                                                                                                                                                    >  1:352337/1‑251 (MQ=255)                                  cAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGt                                                                                                                                                                                                  <  1:477383/250‑1 (MQ=255)                                    aTACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGa                                                                                                                                                                                               >  1:100694/1‑251 (MQ=255)                                    aTACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGa                                                                                                                                                                                               >  1:395975/1‑251 (MQ=255)                                              aTCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTg                                                                                                                                                                                                <  1:556794/240‑1 (MQ=255)                                                cTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAg                                                                                                                                                                                   <  1:816184/251‑1 (MQ=255)                                                cTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAg                                                                                                                                                                                   <  1:188876/251‑1 (MQ=255)                                                cTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAg                                                                                                                                                                                   <  1:455498/251‑1 (MQ=255)                                                         aGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCACTTAGACAGAAAATTAGACAGAACAAGTTTCTAc                                                                                                                                                                                                                         >  1:430204/1‑204 (MQ=255)                                                             tgGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGcacgg                                                                                                                                                                       <  1:758756/250‑1 (MQ=255)                                                              ggCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGcacggt                                                                                                                                                                      >  1:490295/1‑250 (MQ=255)                                                                   aCGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGcacggtcacg                                                                                                                                                                  <  1:551514/249‑1 (MQ=255)                                                                                      ttCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTAc                                                                                                                                              >  1:586367/1‑250 (MQ=255)                                                                                                  ggAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATTCGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATAACCGCGCTAAAAg                                                                                                                                   <  1:571892/249‑1 (MQ=255)                                                                                                  ggAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAg                                                                                                                                   <  1:767363/249‑1 (MQ=255)                                                                                                     aTCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGcacggt                                                                                                                                                                      <  1:430800/211‑1 (MQ=255)                                                                                                            gAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAAt                                                                                                                                                                                             <  1:609229/181‑1 (MQ=255)                                                                                                              aaTAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGt                                                                                                                     <  1:597676/251‑1 (MQ=255)                                                                                                                           ccacAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACggtgg                                                                                                                                                              <  1:58172/197‑1 (MQ=255)                                                                                                                           ccacAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGa                                                                                                                                                           >  1:280959/1‑200 (MQ=255)                                                                                                                           ccacAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGt                                                                                                                <  1:231634/243‑1 (MQ=255)                                                                                                                                aGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATAc                                                                                                                                                        <  1:87015/198‑1 (MQ=255)                                                                                                                                       aCGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGc                                                                                                                                                                           <  1:267106/172‑1 (MQ=255)                                                                                                                                       aCGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATAcc                                                                                                                                                       <  1:282536/192‑1 (MQ=255)                                                                                                                                       aCGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGa                                                                                            >  1:801933/1‑251 (MQ=255)                                                                                                                                          atatatAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGt                                                                                                                <  1:585543/228‑1 (MQ=255)                                                                                                                                          atatatAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTg                                                                                              >  1:229918/1‑246 (MQ=255)                                                                                                                                           tatataACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAAc                                                                                        <  1:26064/251‑1 (MQ=255)                                                                                                                                            atataACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAAc                                                                                                                                                                                                                                   >  1:755727/1‑111 (MQ=255)                                                                                                                                            atataACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGtt                                                                                                                                                                                                                              <  1:183652/116‑1 (MQ=255)                                                                                                                                            atataACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTAc                                                                                                                       >  1:294879/1‑219 (MQ=255)                                                                                                                                            atataACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGAt                                                                                           >  1:330040/1‑247 (MQ=255)                                                                                                                                            atataACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACg                                                                                       >  1:482136/1‑251 (MQ=255)                                                                                                                                                                   ttGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACAtt                                                                 <  1:219092/250‑1 (MQ=255)                                                                                                                                                                    tGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTa                                                               <  1:675050/251‑1 (MQ=255)                                                                                                                                                                      ccAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTAc                                                              <  1:255207/250‑1 (MQ=255)                                                                                                                                                                       cAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATGCCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTc                                                            <  1:655136/251‑1 (MQ=255)                                                                                                                                                                       cAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTc                                                            <  1:865248/251‑1 (MQ=255)                                                                                                                                                                       cAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTc                                                            <  1:386647/251‑1 (MQ=255)                                                                                                                                                                         aTGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTc                                                            >  1:91764/1‑249 (MQ=255)                                                                                                                                                                         aTGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCAt                                                          >  1:573242/1‑251 (MQ=255)                                                                                                                                                                                  cTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAg                                                                                                                                                                                   >  1:17662/1‑121 (MQ=255)                                                                                                                                                                                  cTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAatgat                                                                                                          <  1:845478/194‑1 (MQ=255)                                                                                                                                                                                  cTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTa                                                               <  1:214859/237‑1 (MQ=255)                                                                                                                                                                                  cTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAAc                                                  >  1:834714/1‑250 (MQ=255)                                                                                                                                                                                         tGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAAcc                                          <  1:695081/251‑1 (MQ=255)                                                                                                                                                                                         tGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACGCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAAcc                                          <  1:427657/251‑1 (MQ=255)                                                                                                                                                                                              cAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATAc                                     <  1:569107/251‑1 (MQ=255)                                                                                                                                                                                                aTCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGcacggtcacggt                                                                                                                                                                <  1:628931/126‑1 (MQ=255)                                                                                                                                                                                                aTCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACgg                                   <  1:686565/251‑1 (MQ=255)                                                                                                                                                                                                  cTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACgg                                                                                      >  1:758988/1‑198 (MQ=255)                                                                                                                                                                                                  cTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAggc                                                                         <  1:246147/211‑1 (MQ=255)                                                                                                                                                                                                  cTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGt                                  >  1:311787/1‑250 (MQ=255)                                                                                                                                                                                                       aTACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTg                            >  1:884594/1‑251 (MQ=255)                                                                                                                                                                                                         aCATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTtggtgg                                                                                                                                                                                                           <  1:381909/74‑1 (MQ=255)                                                                                                                                                                                                          cATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTAc                                                                                                                                              <  1:157879/134‑1 (MQ=255)                                                                                                                                                                                                           aTTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTGAGTg                        >  1:342673/1‑251 (MQ=255)                                                                                                                                                                                                                aTTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAAtgatga                                                                                                        >  1:896926/1‑166 (MQ=255)                                                                                                                                                                                                                aTTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTGAGTGGaaaa                   >  1:879477/1‑251 (MQ=255)                                                                                                                                                                                                                       tCTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTGAGTGGAAAAGATAccc            <  1:735178/251‑1 (MQ=255)                                                                                                                                                                                                                        cTTTCCGGGCGCTTAGACAGAAAATTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACggtgg                                                                                                                                                              >  1:158578/1‑104 (MQ=255)                                                                                                                                                                                                                                             aaaTTAGACAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTGAGTGGAAAAGATACCCTCACTAAAAt  >  1:552726/1‑239 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        TGGTTAAACGTTTAACTACAGAAATCGCATCGGCAATACATGCCAATCTTGAACCTAGTGTGGCTGACGCTGAGGTTATCGCTAATTCGGAAGTTAGGGAATCTGAGGAAATAGAAGTAGCACCACAAGCCCCAACGATATATAACCCAAATGAAATTCTCGTTGCCAATGCGATCGCTTCCAATGTAGCAATCTTGAATACATTAAATTTAGATCTTTCCGGGCGCTTAGACAGAAAATTAGGCAGAACAAGTTTCTACCGCACTCTTGGTGGTGTTTGGGTTGAATACGTTAATAGCGATATGCACGGTCACGGTGGTGATACCAATAATTACCGCGCTAAAAGTAACCAAATTACGTTAGGTAATGATGAAGCCCAATTGGATAACGGCGTGGTTTTAGGCGGCACATTTACTCATGCCAAAACCGATAACCAATACGGTCCATTGAGTGGAAAAGATACCCTCACTAAAAT  >  1/1620587‑1621061

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