Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	612,537	T→G	A26A (GCT→GCG)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	612,537	0	T	G	93.4%	161.0 / 3.2	76	A26A (GCT→GCG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (39/32);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01

AATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAGATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAG  >  1/612292‑612721                                                                                                                                                                                                                                                      |                                                                                                                                                                                            aaTTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaaact                                                                                                                                                                                        <  1:211689/251‑3 (MQ=255)                  cacCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAGGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGtt                                                                                                                                                                       >  1:275174/1‑251 (MQ=255)                  cacCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGtt                                                                                                                                                                       >  1:215795/1‑251 (MQ=255)                                gCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTAt                                                                                                                                                         >  1:341581/1‑251 (MQ=255)                                    aaaTGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑a                                                                                                                                                   <  1:56730/251‑1 (MQ=255)                                            ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGcc                                                                                                                                             <  1:130955/249‑1 (MQ=255)                                            ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCtt                                                                                                                                           >  1:858432/1‑251 (MQ=255)                                                  cacaCCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGcc                                                                                                                                             >  1:68607/1‑243 (MQ=255)                                                   acacCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGaa                                                                                                                                    <  1:580042/251‑1 (MQ=255)                                                               tCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAg                                                                                                                                                                                 >  1:110595/1‑196 (MQ=255)                                                               tCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAAt                                                                                                                        >  1:26391/1‑251 (MQ=255)                                                                                        cTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                                                              <  1:254574/251‑1 (MQ=255)                                                                                        cTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                                                              <  1:427356/251‑1 (MQ=255)                                                                                                        ttCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTaaaat                                                                              <  1:828780/251‑4 (MQ=255)                                                                                                        ttCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTaaaat                                                                              <  1:185103/251‑4 (MQ=255)                                                                                                                gCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                                                                   >  1:503364/1‑222 (MQ=255)                                                                                                                 cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                                                             <  1:9422/227‑1 (MQ=255)                                                                                                                 cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                                                            >  1:516002/1‑244 (MQ=255)                                                                                                                 cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACAAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCAGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGCAAGAGGGCTTAAt                                                                                                                        >  1:190524/1‑201 (MQ=255)                                                                                                                  aCTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                                                              <  1:829796/225‑1 (MQ=255)                                                                                                                   cTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGc                                                                                                                             >  1:571814/1‑194 (MQ=255)                                                                                                                       gATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGAt                                                               >  1:342208/1‑251 (MQ=255)                                                                                                                            gCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                                                              <  1:411520/215‑1 (MQ=255)                                                                                                                               ccccAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAAt                                                                                                                        >  1:783098/1‑187 (MQ=255)                                                                                                                               ccccAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaag                                                        >  1:744744/1‑251 (MQ=255)                                                                                                                               ccccAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                                                            >  1:508810/1‑230 (MQ=255)                                                                                                                                           gacacggacacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGa                                                                >  1:629087/1‑230 (MQ=16)                                                                                                                                                ggacacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATc                                         >  1:879860/1‑250 (MQ=11)                                                                                                                                                 gacacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAc                                                                                               >  1:99360/1‑193 (MQ=21)                                                                                                                                                  acacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATccc                                       <  1:901112/250‑1 (MQ=11)                                                                                                                                                   cacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGc                                                                                                                             <  1:460274/162‑1 (MQ=25)                                                                                                                                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaaca                                 >  1:877292/1‑251 (MQ=9)                                                                                                                                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑aaaa                                                                                                              >  1:71407/1‑174 (MQ=21)                                                                                                                                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaa                                                      <  1:139465/230‑1 (MQ=9)                                                                                                                                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                                                                   >  1:578082/1‑187 (MQ=11)                                                                                                                                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaa                                   >  1:736451/1‑251 (MQ=11)                                                                                                                                                      cTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAAACTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaac                                  <  1:594780/251‑1 (MQ=11)                                                                                                                                                        aaTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAag                                                                                                                                  <  1:566356/152‑1 (MQ=255)                                                                                                                                                        aaTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                                                                   >  1:695174/1‑182 (MQ=21)                                                                                                                                                        aaTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaa                                >  1:528141/1‑251 (MQ=11)                                                                                                                                                        aaTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGAGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaaag                                                    <  1:469041/229‑1 (MQ=9)                                                                                                                                                         aTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaaaca                             >  1:851005/1‑251 (MQ=14)                                                                                                                                                         aTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaaaca                             >  1:829844/1‑251 (MQ=14)                                                                                                                                                                gTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTtcataaaatctcc                                                                                                                                                                           <  1:877259/105‑10 (MQ=14)                                                                                                                                                                gTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTaaaat                                                                              <  1:218990/197‑4 (MQ=11)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGGGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCa                 >  1:250326/1‑251 (MQ=9)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                                                                   <  1:762231/169‑1 (MQ=25)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                                                             >  1:372101/1‑175 (MQ=21)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAag                                                                                                                                  >  1:235025/1‑141 (MQ=2)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                                                             >  1:384126/1‑177 (MQ=11)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTaaaat                                                                              >  1:652554/1‑189 (MQ=11)                                                                                                                                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaaact                             <  1:380148/241‑2 (MQ=11)                                                                                                                                                                      aTTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTccc                   <  1:223326/248‑1 (MQ=14)                                                                                                                                                                         aTGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg             <  1:207233/251‑1 (MQ=14)                                                                                                                                                                         aTGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaag                                                        <  1:404710/210‑1 (MQ=11)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                                                                   <  1:161160/163‑1 (MQ=25)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                                                              >  1:68552/1‑168 (MQ=25)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGAt                                                               <  1:320903/199‑1 (MQ=12)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCg           >  1:450059/1‑251 (MQ=14)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAag                                                                                                                                  <  1:377876/135‑1 (MQ=2)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAcc                                                                                              >  1:880063/1‑170 (MQ=11)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAAc                                                                                        >  1:217527/1‑176 (MQ=11)                                                                                                                                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg             >  1:897128/1‑251 (MQ=11)                                                                                                                                                                                actcttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTaaaat                                                                              <  1:316459/179‑4 (MQ=12)                                                                                                                                                                                actcttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTaa      >  1:366174/1‑251 (MQ=14)                                                                                                                                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtatat                                                                                                   <  1:244382/155‑1 (MQ=35)                                                                                                                                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGAc                          <  1:656815/228‑1 (MQ=14)                                                                                                                                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTaaaaa   >  1:732945/1‑251 (MQ=14)                                                                                                                                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTCCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTaaaaa   >  1:297610/1‑251 (MQ=14)                                                                                                                                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg             >  1:249516/1‑241 (MQ=14)                                                                                                                                                                                   cttacCAATTGCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                                                            <  1:454698/180‑1 (MQ=11)                                                                                                                                                                                   cttacCAACTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGa                                                                         <  1:779677/181‑1 (MQ=17)                                                                                                                                                                                     tacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCAtt                                                                                                                                                            <  1:900396/99‑1 (MQ=255)                                                                                                                                                                                     tacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAg  >  1:80582/1‑250 (MQ=14)                                                                                                                                                                                       ccAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAAccc                                                                                             <  1:249728/157‑1 (MQ=16)                                                                                                                                                                                               gTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTtata                                                                                                    <  1:892280/144‑1 (MQ=11)                                                                                                                                                                                                                                                      |                                                                                                                                                                                            AATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAGATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAG  >  1/612292‑612721

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