Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	330,906	T→C	E147E (GAA→GAG)	unknown ←	Type III restriction‑modification system methylation subunit (EC 2.1.1.72)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	330,906	0	T	C	100.0%	200.5 / NA	60	E147E (GAA→GAG)	unknown	Type III restriction‑modification system methylation subunit (EC 2.1.1.72)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (37/23);  total (37/23)

GGCTTGTCCCTGCCGTTCGTATAAACGGTTTTTGATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCTTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGATTTTCTTTAAGGAAAGA‑TTGGGGAAACGTTCCAACGCCTTTTCCGCCGTGAT  >  1/330664‑331151                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         ggCTTGTCCCTGCCGTTCGTATAAACGGTTTTTGATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCGACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTccc                                                                                                                                                                                                                                                   >  1:840871/1‑250 (MQ=255)               gTTCGTATAAACGGTTTTTGATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCCAACGCCATTTGTGCACCAAGCTGTTGAACATTTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGAATCCACCAATCCTCCTCCCCGCCAAACAATcc                                                                                                                                                                                                                                      >  1:768976/1‑250 (MQ=255)                   gTATAAACGGTTTTTGATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGcc                                                                                                                                                                                                                                   >  1:24859/1‑248 (MQ=255)                    tataAACGGTTTTTGATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGa                                                                                                                                                                                                                               >  1:224291/1‑251 (MQ=255)                                  gATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCa                                                                                                                                                                                                                   >  1:649576/1‑249 (MQ=255)                                        gATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATCTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTggg                                                                                                                                                                                                            <  1:597712/250‑1 (MQ=255)                                          tttCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTccc                                                                                                                                                                                                                                                   >  1:414390/1‑209 (MQ=255)                                          tttCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCtt                                                                                                                                                                                                         >  1:224904/1‑251 (MQ=255)                                             ccAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCttttt                                                                                                                                                                                                      >  1:144635/1‑251 (MQ=255)                                                  gtgGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGc                                                                                                                                                                                                  >  1:789421/1‑250 (MQ=255)                                                             ccAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTggc                                                                                                                                                                                      <  1:79561/251‑1 (MQ=255)                                                             ccAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTggc                                                                                                                                                                                      <  1:325711/251‑1 (MQ=255)                                                                               cTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAga                                                                                                                                                                    <  1:374191/251‑1 (MQ=255)                                                                                 tCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTggc                                                                                                                                                                                      <  1:42602/231‑1 (MQ=255)                                                                                  cTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAgacga                                                                                                                                                                 >  1:674953/1‑251 (MQ=255)                                                                                    gtgtATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGcc                                                                                                                                                                                                 >  1:717625/1‑217 (MQ=255)                                                                                      gtATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGGTTGTTTTTAGACGACCtt                                                                                                                                                             >  1:158177/1‑251 (MQ=255)                                                                                            aaCAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCg                                                                                                                                                         <  1:752023/249‑1 (MQ=255)                                                                                            aaCAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCg                                                                                                                                                         <  1:483834/249‑1 (MQ=255)                                                                                                         cttctGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCgg                                                                                                                                          >  1:577896/1‑251 (MQ=255)                                                                                                         cttctGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCgg                                                                                                                                          >  1:314028/1‑251 (MQ=255)                                                                                                         cttctGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCgg                                                                                                                                          >  1:305406/1‑251 (MQ=255)                                                                                                                                aacaaGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGGTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCCGCCGGGTGGCg                                                                                                                   >  1:829378/1‑251 (MQ=255)                                                                                                                                 acaaGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGAc                                                                                                                                                                                                                       >  1:610710/1‑150 (MQ=255)                                                                                                                                          tATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGc                                                                                                        <  1:655051/251‑1 (MQ=255)                                                                                                                                          tATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGc                                                                                                        <  1:364912/251‑1 (MQ=255)                                                                                                                                           aTGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCt                                                                                                       <  1:776965/251‑1 (MQ=255)                                                                                                                                                ccc‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAAt                                                                                                  <  1:76891/251‑1 (MQ=255)                                                                                                                                                 cc‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAAtc                                                                                                 <  1:831797/251‑1 (MQ=255)                                                                                                                                                     cGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAAtcgtc                                                                                              <  1:844635/251‑1 (MQ=255)                                                                                                                                                     cGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAAGCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCGCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAAtcgtc                                                                                              <  1:343335/251‑1 (MQ=255)                                                                                                                                                                caacaaGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCAC‑TTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCAt                                                                                    <  1:369223/249‑1 (MQ=255)                                                                                                                                                                  acagcttGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGGTTCCGCCATCGGCAGGGTGGCGG‑CGTTTTGGCTGTAATCGGCGGGTCCCCCTTc                                                                                  >  1:85164/7‑250 (MQ=255)                                                                                                                                                                       cTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGcc                                                                                                                                  <  1:728917/198‑1 (MQ=255)                                                                                                                                                                       cTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGc                                                                                                        <  1:800769/223‑1 (MQ=255)                                                                                                                                                                            gCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGGGG‑CGTTTTGGCTGTAATCGTCGTGTCCCCCTTCGCAACGCGcc                                                                        >  1:366118/1‑250 (MQ=255)                                                                                                                                                                               caATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAAtcgtcg                                                                                             <  1:110189/226‑1 (MQ=255)                                                                                                                                                                                 aTTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGG‑CGTTTTGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAAt                                                                   >  1:377879/1‑250 (MQ=255)                                                                                                                                                                                 aTTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGG‑CGTTTTGGCTGTAATCGTCGTGTCCCCCTTCGCAACGCGCCAAAAt                                                                   >  1:703930/1‑250 (MQ=255)                                                                                                                                                                                   ttCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGA‑CGGTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCCAAATCA‑t                                                               >  1:481542/1‑251 (MQ=255)                                                                                                                                                                                            cAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGAt                                                      >  1:588852/1‑251 (MQ=255)                                                                                                                                                                                             aGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAga                                                                                                                                                                    <  1:218614/142‑1 (MQ=255)                                                                                                                                                                                              gCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGG‑CGTTTTGGc                                                                                                        >  1:138380/1‑200 (MQ=255)                                                                                                                                                                                                   gCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGGTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGA‑CGTTTAGGCTGTAATCGTCGGGTCCCCATTCGCAACGCGCCCAAATCA‑TTTTCGGGGTTTTCttt                                               >  1:748305/1‑251 (MQ=255)                                                                                                                                                                                                    ccTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGATTTTCtt                                                <  1:533654/249‑1 (MQ=255)                                                                                                                                                                                                                  cGCATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGcc                                                                                                                                  <  1:54224/155‑1 (MQ=255)                                                                                                                                                                                                                    cATATATTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGATTTTCTTTAAGGAAAGA‑TTggggga                             >  1:271974/1‑249 (MQ=255)                                                                                                                                                                                                                       atatTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTTAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGcc                                                                                                                                  >  1:352424/1‑150 (MQ=255)                                                                                                                                                                                                                       atatTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGA‑CGTTTTGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTTGGGGTTTTCTTTAAGGAAAGA‑TTGGGGGAACg                          >  1:658380/1‑251 (MQ=255)                                                                                                                                                                                                                         atTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTTTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGGGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCCCAACGCGCCCAAATCA‑TTTTTGGGGTTTTCTTTTAGGAAAAA‑TTGGGGAAAAGtt                        >  1:541265/1‑251 (MQ=255)                                                                                                                                                                                                                         atTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGATTTTCTTTAAGGAAAGA‑t                                    <  1:235846/239‑1 (MQ=255)                                                                                                                                                                                                                         atTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGATTTTCTTTAAGGAAAGA‑t                                    <  1:267422/239‑1 (MQ=255)                                                                                                                                                                                                                         atTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGG‑CGGTTAGGCTGTAATTGTCGTGTCCCCATTCGCAACGCGCCCAAATCA‑TTTTCTGGATTTTTTTTAAGGAAAGA‑TTGGGGGAACggt                        >  1:633551/1‑249 (MQ=255)                                                                                                                                                                                                                         atTTTAGATTGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTCGTTTTTAGACGACCTTGGCGCTTTTTCAAGCTCGGGTTCCGCCATCGGCAGGGTGGCGAGCG‑TTAGGCTGGAATTGTCGTGTCCCCCTTCCCAACGCGCCCAAATCA‑TTTTCCGGGTTTTCTTTAAGGAAAGA‑TTtggggctgttt                        >  1:413153/1‑244 (MQ=255)                                                                                                                                                                                                                                  tGCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGGTTGTTTTCAGACGACCCTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTTGCGA‑CGTTTTGGCTGTAATCGTCGTGTCCCCCTTTGCCAAGCGCCAAAAACA‑TTTTTGGGGTTTTCTTTAAGGGAAGA‑TTgggggaaaggtcccacgcc                >  1:402272/1‑235 (MQ=255)                                                                                                                                                                                                                                   gCGGGTATCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGGTTGTTTTCAGACGACCTTTGCGTTTTTTTAAGCTCGGGTTTCGCCCTCCGCAGGTTTGGGA‑CGGTTTGGGTGTTATTCTTGTGGCCCCCTTCGCCACGCGCCCAAAACATTTTTGGGGTTTTTTTTTAGGGAAAAA‑TTGGGGaaagcttcaacgcc                 >  1:565090/1‑238 (MQ=255)                                                                                                                                                                                                                                         aTCCATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGGGA‑CGGTTAGGCTGTAATCGTCGTGTCCCCATTCGCCACGCGCCAAAATCA‑TTTTCGGGGTTTTTTTTAAGGGAAGATTTGGGGAAA‑GTTCCCACGCCCTTTc           >  1:494820/1‑248 (MQ=255)                                                                                                                                                                                                                                           ccATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGG‑CGTTTTGGCTGTAATCGTCGTGTCCCCCTTCGCAACGCGCCCAAATCA‑TTTTTGGGATTTTCTTTTAGGGAAGA‑TTGGGGAAAAGTTTCAACGCCTTTTccccc       >  1:512584/1‑247 (MQ=255)                                                                                                                                                                                                                                           ccATCAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGGTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCCTTCGCAACGCGCCCAAATCA‑TTTTCGGGGTTTTCTTTTAGGAAAGA‑TTGGGGAAACGTTCCAACGCCTTTTccgccg      >  1:783595/1‑251 (MQ=255)                                                                                                                                                                                                                                               cAATCCTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTTGGCTGTAATCGTCGTGTCCCCCTTCGCAACGCGCCAAAATCA‑TTTTCGGGGTTTTCTTTAAGGGAAGA‑TTGGGGAAAAGTTCCAACGCCTTTTCCGCCGTGAt  >  1:428170/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                         GGCTTGTCCCTGCCGTTCGTATAAACGGTTTTTGATGTGGATTTCCAATGTGGTACGCGACCAACTTTCAGCGAGTGCCTTCTGTGTATAGAACAAACGGGCTTCTTCTGTTTTAATCCGGTTTAGCAACAAGACGTTATGCCCCC‑ACGGCAATTGTGCAACAAGCTGTTGCACAATTTCAAAATCCAGCCAAGCCTCAGCAAAAGCACGCATATATTTTAGATTGCGGGTATCCATCAATCTTCCTCCCCGTCAAACAATCCGCCCTGATTTTCGACCGCACTTTGGGCTTTTTTTGCCTTGGTTTTGGCGGCGGGTTTGTTTTCAGACGACCTTGGCGTTTTTTCAAGCTCGGTTTCCGCCATCGGCAGGTTGGCGA‑CGTTTAGGCTGTAATCGTCGTGTCCCCATTCGCAACGCGCCAAAATCA‑TTTTCGGGATTTTCTTTAAGGAAAGA‑TTGGGGAAACGTTCCAACGCCTTTTCCGCCGTGAT  >  1/330664‑331151

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