Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	925,658	T→C	I670T (ATC→ACC)	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	925,658	0	T	C	100.0%	117.8 / NA	44	I670T (ATC→ACC)	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 (24/20);  total (24/20)

TTCCCCATCTTATGAAAAAGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAATCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGAGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAATCTTGAATTACGCA  >  1/925425‑925900                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   ttCCCCATCTTATGAAAAAGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTa                                                                                                                                                                                                                                   <  1:122161/251‑1 (MQ=255) ttCCCCATCTTATGAAAAAGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTa                                                                                                                                                                                                                                   >  1:297098/1‑251 (MQ=255)    cccATCTTATGAAAAAGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAAcc                                                                                                                                                                                                                                <  1:204035/251‑1 (MQ=255)    cccATCTTATGAAAAAGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAAcc                                                                                                                                                                                                                                <  1:534752/251‑1 (MQ=255)               aaaaaGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAAc                                                                                                                                                                                                                           >  1:560416/1‑245 (MQ=255)                              aGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCg                                                                                                                                                                                                       >  1:309718/1‑250 (MQ=255)                                         ccTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATc                                                                                                                                                                                             <  1:847917/249‑1 (MQ=255)                                               aaTATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTgg                                                                                                                                                                                      <  1:239573/250‑1 (MQ=255)                                                                               tAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAata                                                                                                                                                                                                 <  1:160168/207‑1 (MQ=255)                                                                                aaTGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAAc                                                                                                                                                     >  1:221239/1‑250 (MQ=255)                                                                                        tGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAAt                                                                                                                                                                                                                                     >  1:743408/1‑162 (MQ=255)                                                                                        tGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAAc                                                                                                                                                                                   <  1:348550/212‑1 (MQ=255)                                                                                        tGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATAc                                                                                                                                                                <  1:279173/231‑1 (MQ=255)                                                                                             gTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACaaaaa                                                                                                                                       >  1:711949/1‑251 (MQ=255)                                                                                             gTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCAGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAAc                                                                                                                                                     <  1:4207/237‑1 (MQ=255)                                                                                              tGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAAcc                                                                                                                                                                                                                                <  1:694570/161‑1 (MQ=255)                                                                                                  gAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGgt                                                                                                                                                             >  1:211853/1‑224 (MQ=255)                                                                                                  gAATTAGATCATACGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAAc                                                                                                                                                                                                                                 >  1:268100/1‑156 (MQ=255)                                                                                                                                   gATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGa                                                                                                 >  1:628631/1‑251 (MQ=255)                                                                                                                                                  aaaCCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTtg                                                                                                                                                                                         <  1:505899/148‑1 (MQ=255)                                                                                                                                                            cTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTgg                                                                                                                                                                                      >  1:416261/1‑141 (MQ=255)                                                                                                                                                                   gAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAg                                                                                               >  1:90875/1‑221 (MQ=255)                                                                                                                                                                          atTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAAt                                                                           >  1:141002/1‑234 (MQ=255)                                                                                                                                                                                               gTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTAttt                                                                                                                                                                      <  1:261791/122‑1 (MQ=255)                                                                                                                                                                                               gTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAg                                     >  1:516200/1‑251 (MQ=255)                                                                                                                                                                                               gTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGAGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAg                                     >  1:229103/1‑251 (MQ=255)                                                                                                                                                                                                      gTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCt                                                                                                                                                                                           <  1:603598/94‑1 (MQ=255)                                                                                                                                                                                                      gTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTACAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGaa                              >  1:97241/1‑251 (MQ=255)                                                                                                                                                                                                      gTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAg                                                                                               >  1:53319/1‑186 (MQ=255)                                                                                                                                                                                                      gTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGaa                              >  1:65354/1‑251 (MQ=255)                                                                                                                                                                                                       tCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAg                             <  1:711145/251‑1 (MQ=255)                                                                                                                                                                                                       tCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAg                             <  1:741276/251‑1 (MQ=255)                                                                                                                                                                                                       tCTTATTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAg                             <  1:329331/251‑1 (MQ=255)                                                                                                                                                                                                           aTTCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGaa                                                       >  1:632876/1‑221 (MQ=255)                                                                                                                                                                                                            ttCTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATAc                         >  1:498317/1‑250 (MQ=255)                                                                                                                                                                                                              cTTGACTACTTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAg                             >  1:793200/1‑244 (MQ=255)                                                                                                                                                                                                                     actTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAAt               >  1:706546/1‑251 (MQ=255)                                                                                                                                                                                                                     actTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAAt               >  1:723507/1‑251 (MQ=255)                                                                                                                                                                                                                     actTCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAAt               >  1:817720/1‑251 (MQ=255)                                                                                                                                                                                                                        tCGCAGGTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTgg                                                                                                                                                                                      <  1:440410/81‑1 (MQ=255)                                                                                                                                                                                                                              gTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAATCTTGAAt        <  1:380936/249‑1 (MQ=255)                                                                                                                                                                                                                              gTAGCGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAATCTTGAAt        <  1:426130/249‑1 (MQ=255)                                                                                                                                                                                                                                  cGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGt                                                                      <  1:418323/183‑1 (MQ=255)                                                                                                                                                                                                                                  cGGCACAACCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGGGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAATCTTGAATTACGCa  >  1:899077/1‑251 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                   TTCCCCATCTTATGAAAAAGAAGAATATGAGCCAGAGAAGCCTTGGAATATCATTAATTTAAATTTTAATGTTGGCACTAATGAAAATGCTAGTGCAGAATTAGATCATATGTTTTATCAGAATTTTTCTGATAAGCGTTATCCCAAACCAGTTTCTTTAATGAAATATATTTTGAATTTTAATTTAAATGTAGAAAGTCTTATTCTTGACTACTTCGCAGGTAGCGGCACAATCGCCCACGCTGTAATTAACCTCAACCTTGAAGATAACGGCAATCGCAAATATATCCTTGTGGAACAAGGGGAGTATTTCGATACGGTGTTAAAACCTCGCGTACAAAAAGTGATTTTTAGCGAAAAATGGAAAGAGGGGAAGCCAGAAGCGGGAGAAAATGGTTTGAATGGCGTATCGCAAATTGTGAAAGTGTTGAAGTTGGAAAGCTATGAAGATACCTTGAACAATCTTGAATTACGCA  >  1/925425‑925900

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