Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,531,873	C→T	T216T (ACG→ACA)	unknown ←	Putative DNA‑binding protein in cluster with Type I restriction‑modification system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,531,873	0	C	T	100.0%	214.8 / NA	71	T216T (ACG→ACA)	unknown	Putative DNA‑binding protein in cluster with Type I restriction‑modification system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/44);  total (27/44)

ATTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCCGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTATCAC  >  1/1531627‑1532112                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                aTTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCg                                                                                                                                                                                                                                             <  1:879739/251‑1 (MQ=255)        ttAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTcccc                                                                                                                                                                                                                                      <  1:901150/251‑1 (MQ=255)             cAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCttttt                                                                                                                                                                                                                                 >  1:596048/1‑251 (MQ=255)                  gTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCa                                                                                                                                                                                                                              >  1:587726/1‑249 (MQ=255)                      tACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACg                                                                                                                                                                                                                        <  1:201599/251‑1 (MQ=255)                             cTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGACGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGttt                                                                                                                                                                                                                 <  1:456341/251‑1 (MQ=255)                                gCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTaa                                                                                                                                                                                                              >  1:810608/1‑251 (MQ=255)                                     tAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGTCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTc                                                                                                                                                                                                                                <  1:451803/228‑1 (MQ=255)                                        gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGTCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAAt                                                                                                                                                                                                             <  1:344599/244‑1 (MQ=255)                                        gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                      >  1:78679/1‑251 (MQ=255)                                        gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                      >  1:623770/1‑251 (MQ=255)                                        gtgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                      >  1:335900/1‑251 (MQ=255)                                         tgtTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAAt                                                                                                                                                                                                     <  1:249776/251‑1 (MQ=255)                                                    ccGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATa                                                                                                                                                                                          >  1:698017/1‑251 (MQ=255)                                                       cAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGt                                                                                                                                                                                       <  1:3103/251‑1 (MQ=255)                                                        aaaTACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTc                                                                                                                                                                                      <  1:832347/251‑1 (MQ=255)                                                        aaaTACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGCCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTc                                                                                                                                                                                      <  1:672193/251‑1 (MQ=255)                                                            aCCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCt                                                                                                                                                                                  <  1:612487/251‑1 (MQ=255)                                                             ccATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTa                                                                                                                                                                                 >  1:456425/1‑251 (MQ=255)                                                                 aaaaaCAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAacac                                                                                                                                                                             <  1:552502/251‑1 (MQ=255)                                                                 aaaaaCAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAacac                                                                                                                                                                             <  1:213860/251‑1 (MQ=255)                                                                 aaaaaCAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAacac                                                                                                                                                                             <  1:417930/251‑1 (MQ=255)                                                                      cAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACa                                                                                                                                                                        >  1:652429/1‑251 (MQ=255)                                                                      cAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACa                                                                                                                                                                        >  1:489334/1‑251 (MQ=255)                                                                         gAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAg                                                                                                                                                                       <  1:463568/249‑1 (MQ=255)                                                                          aaaGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAAc                                                                                                                                                                    <  1:767028/251‑1 (MQ=255)                                                                                  cGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCAtattat                                                                                                                                                            <  1:654778/251‑1 (MQ=255)                                                                                      cGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                          >  1:23182/1‑249 (MQ=255)                                                                                      cGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTg                                                                                                                                                        >  1:891714/1‑251 (MQ=255)                                                                                       gTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTg                                                                                                                                                        >  1:313659/1‑250 (MQ=255)                                                                                             tGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGACGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGc                                                                                                                                                 <  1:497704/251‑1 (MQ=255)                                                                                               ccATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGc                                                                                                                                                                                                                       <  1:406308/179‑1 (MQ=255)                                                                                               ccATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAAc                                                                                                                                                                    >  1:466187/1‑230 (MQ=255)                                                                                                         aCGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTGATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGtt                                                                                                                                     >  1:569589/1‑251 (MQ=255)                                                                                                                          ataaAATACAATAAATTCGCAGCTCGGCTAGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTGCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCGCTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                          <  1:185487/213‑1 (MQ=255)                                                                                                                              aaTACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGTAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGcc                                                                                                                <  1:453119/251‑1 (MQ=255)                                                                                                                              aaTACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGcc                                                                                                                <  1:634262/251‑1 (MQ=255)                                                                                                                                  cAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAAc                                                                                                                                                                         <  1:215737/190‑1 (MQ=255)                                                                                                                                   aaTAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGaaa                                                                                                                                                                                                      >  1:706422/1‑160 (MQ=255)                                                                                                                                        aTTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAaac                                                                                                      <  1:877711/251‑1 (MQ=255)                                                                                                                                             ccgcTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCg                                                                                                                                       <  1:467632/213‑1 (MQ=255)                                                                                                                                                cTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAaactccaactc                                                                                              <  1:796490/251‑1 (MQ=255)                                                                                                                                                          aaCATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAAt                                                                                    <  1:24877/251‑1 (MQ=255)                                                                                                                                                            cATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATc                                                                                   >  1:768504/1‑250 (MQ=255)                                                                                                                                                                  gATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGcac                                                                             <  1:848744/250‑1 (MQ=255)                                                                                                                                                                     aaaaaaGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCacaaca                                                                         >  1:274131/1‑251 (MQ=255)                                                                                                                                                                         aaGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAAc                                                                                                                                                                         <  1:703540/151‑1 (MQ=255)                                                                                                                                                                             gtgtTACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGAtt                                                                                                                                            <  1:21824/176‑1 (MQ=255)                                                                                                                                                                                tCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGc                                                              <  1:269498/251‑1 (MQ=255)                                                                                                                                                                                       aTCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATcc                                                                                  >  1:543255/1‑224 (MQ=255)                                                                                                                                                                                        tCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAGCTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCtt                                                      >  1:527664/1‑251 (MQ=255)                                                                                                                                                                                        tCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCacaac                                                                          <  1:890863/231‑1 (MQ=255)                                                                                                                                                                                           ttGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTg                                                                                                                                                                                                           >  1:898018/1‑99 (MQ=255)                                                                                                                                                                                           ttGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCg                                                                                                                                       >  1:840134/1‑167 (MQ=255)                                                                                                                                                                                           ttGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGt                                                   >  1:792670/1‑251 (MQ=255)                                                                                                                                                                                            tGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTc                                                  <  1:272966/251‑1 (MQ=255)                                                                                                                                                                                                  aTCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAAc                                                                                                          <  1:569253/189‑1 (MQ=255)                                                                                                                                                                                                       ggCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGc                                                              <  1:210457/228‑1 (MQ=255)                                                                                                                                                                                                        gCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAg                                       <  1:710472/250‑1 (MQ=255)                                                                                                                                                                                                           aTGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAgattg                                   <  1:763339/251‑1 (MQ=255)                                                                                                                                                                                                             ggCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTccc                                                                                                                                                                                                                                       <  1:547567/53‑1 (MQ=255)                                                                                                                                                                                                                   ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCTCTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGAc                                                                                                                                 >  1:119044/1‑149 (MQ=255)                                                                                                                                                                                                                   ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                          <  1:289670/124‑1 (MQ=255)                                                                                                                                                                                                                   ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCa                           >  1:153483/1‑251 (MQ=255)                                                                                                                                                                                                                   ggCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCa                           >  1:318235/1‑251 (MQ=255)                                                                                                                                                                                                                        aaCCGGAATTTCGTTGTTGTCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTc                      <  1:277146/251‑1 (MQ=255)                                                                                                                                                                                                                                 ttCGTTGTTGGCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTTATGTAACACCAACAGAACCATAGTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACACCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGc             <  1:431641/251‑1 (MQ=255)                                                                                                                                                                                                                                          gtCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGt                                                                                                                                                                                       <  1:363127/72‑1 (MQ=255)                                                                                                                                                                                                                                          gtCCAAACAACGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTAt     >  1:74566/1‑250 (MQ=255)                                                                                                                                                                                                                                                acaacGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTATCAc  <  1:172635/247‑1 (MQ=255)                                                                                                                                                                                                                                                 caacGCTGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTAttt                                                                                                                                                          <  1:801063/94‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                ATTGTCATTAATCAATTGTTCTACCGGCCTTGCCAATAAGTGTTGATTCAACCGCAAATACCATAAAAACAGGAAAGATCCCGTGCGTTTATTGCCATCGGTCAACGGATGATTTTTGATAATAAAATACAATAAATTCGCCGCTCGGCTGGCAACATTCGGATAAAAAAGCGTGTCACCGAATCCTTGTTCAATCGTGGCAATGGCAGAGGCTAAACCGGAATTTCGTTGTTGTCCAAACAACGCCGTCGCTTCCCCTTTTTCCACTAACGCCTGTTTTAATTGAGAAATTGCCTGTAATACGTCATCTAACACCAACAGAACCATATTATTTTGCTTATGCTGATTAGCCGTTAGACTTTGCTCATCGTAAGCCTGTAACAAACTCCAACTCCTAGCATAATCCTGCACAACAGACAATACGGCTTTGCCTTCGTCAGATAACAAAGATTGATTACTCAAGGTCTGCCCAAGCAAACCTATCAC  >  1/1531627‑1532112

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