Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,476,150	G→A	S211S (AGC→AGT)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,476,150	0	G	A	98.6%	176.0 / ‑3.8	71	S211S (AGC→AGT)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (43/27);  minor base T (1/0);  total (44/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.43e-01

ATTTTTACTTTTGCTAAAAGGGAAATGAAGTTGAGCGGAAAAGTGCGGCAAAACACCGCCAATCCACACCGCACTTTTCAGGGCAGGAAAGTATTCCCATTGGCGCGCTAACGTGGCATCGTCGCCAACAAGGAGAATTAGCTGCCGCGCCATGCTTTACCAACCCCGCGAATACGGATCCATCCAGTCTACATTGATATTTTTATTGCTTTTATAGAGCGGTTCGCTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCGTCGTCTATCGTTTGTCGTTGCTTTTGCACTTGGTTTTCAACAGAACGACTGTCGCCCGCCAAGCCTTGCCGGCGTTTGCGTTCGTTGAGACGATCGCGAAAGTGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTGAAAATAATCG  >  1/1475925‑1476396                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       aTTTTTACTTTTGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCaa                                                                                                                                                                                                                               >  1:108532/1‑251 (MQ=255)            tGCTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGt                                                                                                                                                                                                                    <  1:198051/251‑1 (MQ=255)              cTAAAAGGGAAGTGAGGCTGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTg                                                                                                                                                                                                                  <  1:572754/251‑1 (MQ=255)                             ggctGAGCAGAAAAGTGCGGTGAAACGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCg                                                                                                                                                                                                   >  1:367699/4‑251 (MQ=255)                                           gtgcggtgaaaCGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGT‑TGAAAAccaaatgc                                                                                                                                                                                                >  1:123204/9‑232 (MQ=255)                                                gtgaaaCGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTg                                                                                                                                                                                                                  >  1:5892/4‑217 (MQ=255)                                                gtgaaaCGCCGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGt                                                                                                                                                                                 <  1:889965/247‑1 (MQ=255)                                                        ccGCCAATCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGAATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTg                                                                                                                                                                        >  1:214952/1‑251 (MQ=255)                                                              aTCCATACCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCtt                                                                                                                                                                     >  1:210647/1‑248 (MQ=255)                                                                  atacCGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCatt                                                                                                                                                              >  1:53669/3‑249 (MQ=255)                                                                      cGCACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGt                                                                                                                                                          <  1:573032/251‑1 (MQ=255)                                                                        cACTTTTCAGGGCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGtt                                                                                                                                                         <  1:246518/250‑1 (MQ=255)                                                                                   gCAGGAAAATGTTCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAAc                                                                                                                                              <  1:486308/250‑1 (MQ=255)                                                                                        aaaatgttCCCATTGGTGCGCTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGa                                                                                                                                            >  1:444568/7‑247 (MQ=255)                                                                                                          cgcTAACGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTgac                                                                                                                      <  1:709852/251‑3 (MQ=255)                                                                                                              aaCGTGGCATCGTCACCAATGAGAAGAATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCg                                                                                                                  >  1:42654/1‑251 (MQ=255)                                                                                                                      atcgtcaccaatgagaagaATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTAcgtt                                                                                                          <  1:259635/235‑1 (MQ=255)                                                                                                                          tcaccaatgagaagaATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGAtt                                                                                                                                                                                                                         <  1:164543/124‑1 (MQ=255)                                                                                                                          tcaccaatgagaagaATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTAcgttcgtt                                                                                                      >  1:585950/13‑251 (MQ=255)                                                                                                                                    gaagaATTAACTGCCGCGCCATAGTTTACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTcacc                                                                                                                                    >  1:88790/3‑208 (MQ=255)                                                                                                                                                       catagtttACCACCCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCCCCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAAt                                                                         <  1:61574/246‑1 (MQ=255)                                                                                                                                                       catagtttACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAAt                                                                         <  1:783882/246‑1 (MQ=255)                                                                                                                                                        atagtttACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCg                                                                                 >  1:285007/5‑242 (MQ=255)                                                                                                                                                         tagtttACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGtttacg                                                                                                            >  1:42119/4‑211 (MQ=255)                                                                                                                                                           gtttACCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATcgcg                                                                     >  1:795635/2‑251 (MQ=255)                                                                                                                                                               aCCAACCCCGTGAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCAt                                                                 >  1:898787/1‑251 (MQ=255)                                                                                                                                                                       cgtgAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAAtt                                                          >  1:106947/4‑250 (MQ=255)                                                                                                                                                                        gtgAATACGGATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCa                                                        >  1:508193/3‑251 (MQ=255)                                                                                                                                                                          gAATACGGATCCGTCCAATCTGCATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAAt                                                                         <  1:894766/232‑1 (MQ=255)                                                                                                                                                                                 gATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGt                                                                                                                                                                                 <  1:342901/121‑1 (MQ=255)                                                                                                                                                                                 gATCCGTCCAATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTa                                               >  1:863042/1‑251 (MQ=255)                                                                                                                                                                                      gtccaATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCTCTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTca                                           >  1:645200/2‑250 (MQ=255)                                                                                                                                                                                      gtccaATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTca                                           >  1:126249/2‑250 (MQ=255)                                                                                                                                                                                      gtccaATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTca                                           >  1:541652/2‑250 (MQ=255)                                                                                                                                                                                      gtccaATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTca                                           >  1:122577/2‑250 (MQ=255)                                                                                                                                                                                       tccaATCTACATTAATATTTTTATTGCTTTTATAGCGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTcac                                          <  1:158483/250‑1 (MQ=255)                                                                                                                                                                                       tccaATCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTcgc                                          <  1:321731/250‑3 (MQ=255)                                                                                                                                                                                        ccaatCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCacac                                        <  1:472438/247‑1 (MQ=255)                                                                                                                                                                                          aatCTACATTAATATTTTTATTGCTTTTATAGGGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCacac                                        <  1:474416/247‑1 (MQ=255)                                                                                                                                                                                          aatCTACATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCg                                      >  1:278285/3‑251 (MQ=255)                                                                                                                                                                                               aCATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAAc                                                                                                                                                   >  1:893543/1‑137 (MQ=255)                                                                                                                                                                                               aCATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACgg                                                                                                                    >  1:515434/1‑168 (MQ=255)                                                                                                                                                                                               aCATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCacac                                        >  1:714529/1‑244 (MQ=255)                                                                                                                                                                                                cATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCTCCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTg                                <  1:549811/251‑1 (MQ=255)                                                                                                                                                                                                cATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGAt                                  >  1:392045/1‑249 (MQ=255)                                                                                                                                                                                                cATTAATATTTTTATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTg                                >  1:716841/1‑251 (MQ=255)                                                                                                                                                                                                        tttttATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGttgttg                        <  1:115215/251‑1 (MQ=255)                                                                                                                                                                                                        tttttATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGttgttg                        >  1:172984/1‑251 (MQ=255)                                                                                                                                                                                                         ttttATTGCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTc                                            <  1:766268/230‑1 (MQ=255)                                                                                                                                                                                                                gCTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTcac                                          >  1:863329/1‑225 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAACCGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCTCCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACt                   <  1:662983/247‑1 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACTCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAgttg                           >  1:610864/1‑239 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCacac                                        >  1:748508/1‑226 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCacac                                        >  1:104400/1‑226 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTg                                <  1:759924/234‑1 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAgttg                           >  1:654266/1‑239 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTaaa               >  1:490626/1‑251 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTaaa               >  1:473484/1‑251 (MQ=255)                                                                                                                                                                                                                 cTTTTATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTaaa               >  1:114056/1‑251 (MQ=255)                                                                                                                                                                                                                  ttttATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAAt              <  1:403612/251‑1 (MQ=255)                                                                                                                                                                                                                   tttATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGaa                                                                                        >  1:833959/1‑176 (MQ=255)                                                                                                                                                                                                                   tttATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGc                                                                                   <  1:428050/181‑1 (MQ=255)                                                                                                                                                                                                                    ttATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTg            >  1:300540/1‑251 (MQ=255)                                                                                                                                                                                                                    ttATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTg            >  1:871342/1‑251 (MQ=255)                                                                                                                                                                                                                    ttATAGAGCGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTg            >  1:445568/1‑251 (MQ=255)                                                                                                                                                                                                                            cGGTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTGAAaataat    <  1:605750/251‑1 (MQ=255)                                                                                                                                                                                                                              gTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTcaccc                                                                                                                                   <  1:371048/122‑1 (MQ=255)                                                                                                                                                                                                                              gTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAg                              <  1:712332/223‑1 (MQ=255)                                                                                                                                                                                                                              gTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGAt                      >  1:412979/1‑231 (MQ=255)                                                                                                                                                                                                                              gTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTGAAAATAATCg  >  1:580052/1‑251 (MQ=255)                                                                                                                                                                                                                              gTTCACTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCATCGTCTATCGTTTGTCGTTGCTTCTGCATTTGGTTTTCAACGGAACGACTGTCACCCGCCAAGCCTTGACGGCGTTTACGTTCGTTGAGGCGATCGCGAAAATGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTGAAAATAATCg  >  1:65591/1‑251 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                       ATTTTTACTTTTGCTAAAAGGGAAATGAAGTTGAGCGGAAAAGTGCGGCAAAACACCGCCAATCCACACCGCACTTTTCAGGGCAGGAAAGTATTCCCATTGGCGCGCTAACGTGGCATCGTCGCCAACAAGGAGAATTAGCTGCCGCGCCATGCTTTACCAACCCCGCGAATACGGATCCATCCAGTCTACATTGATATTTTTATTGCTTTTATAGAGCGGTTCGCTTTCGCCGCGAATTTGGCGTTCAATAGATTCTCGTTGCGTGTGAAAACGGCGAAACGCGTCGTCTATCGTTTGTCGTTGCTTTTGCACTTGGTTTTCAACAGAACGACTGTCGCCCGCCAAGCCTTGCCGGCGTTTGCGTTCGTTGAGACGATCGCGAAAGTGCCGTTTCAAATCGCGGCATTCAAATTCACGAATGTTATTCACACCGTGATTGAGTTGTTGATACTTAAATTGAAAATAATCG  >  1/1475925‑1476396

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