Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	612,541	T→C	F28L (TTT→CTT)	unknown →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	612,541	0	T	C	93.3%	148.4 / 1.3	75	F28L (TTT→CTT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (39/31);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-03

CACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAGATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAG  >  1/612309‑612721                                                                                                                                                                                                                                         |                                                                                                                                                                                        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)                                                                                                                                                                                                                                         |                                                                                                                                                                                        CACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAGATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAG  >  1/612309‑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: ‑