Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	899,038	T→C	T20T (ACT→ACC)	unknown →	16S rRNA (cytosine(1402)‑N(4))‑methyltransferase (EC 2.1.1.199)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	899,038	0	T	C	100.0%	214.5 / NA	73	T20T (ACT→ACC)	unknown	16S rRNA (cytosine(1402)‑N(4))‑methyltransferase (EC 2.1.1.199)
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (47/26);  total (47/26)

GAACAAGACATGGCATTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACTGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCGGCAAGGTTG  >  1/898798‑899279                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  gAACAAGACATGGCATTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACgg                                                                                                                                                                                                                                           >  1:842327/1‑249 (MQ=255) gAACAAGACATGGCATTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACgg                                                                                                                                                                                                                                           >  1:541408/1‑249 (MQ=255) gAACAAGACATGGCATTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAA‑TTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCtt                                                                                                                                                                                                                                        >  1:486733/1‑251 (MQ=255)             gCATTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTgaaag                                                                                                                                                                                                                             >  1:241057/1‑251 (MQ=255)               aTTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTgaaagaa                                                                                                                                                                                                                           >  1:131973/1‑251 (MQ=255)                                cTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGa                                                                                                                                                                                                          >  1:729799/1‑251 (MQ=255)                                cTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGa                                                                                                                                                                                                          >  1:44734/1‑251 (MQ=255)                                cTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGa                                                                                                                                                                                                          >  1:256313/1‑251 (MQ=255)                                   ttGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATgg                                                                                                                                                                                                       <  1:284297/251‑1 (MQ=255)                                           tCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTg                                                                                                                                                                                               <  1:128166/251‑1 (MQ=255)                                                       aaaaCGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGa                                                                                                                                                                                    >  1:313027/1‑250 (MQ=255)                                                               gTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGt                                                                                                                                                                           >  1:578615/1‑251 (MQ=255)                                                                 cATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGttt                                                                                                                                                                         <  1:73071/251‑1 (MQ=255)                                                                     ataAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATc                                                                                                                                                                     <  1:828293/251‑1 (MQ=255)                                                                      taAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATcc                                                                                                                                                                    <  1:754034/251‑1 (MQ=255)                                                                      atagagGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATcc                                                                                                                                                                    <  1:462711/249‑1 (MQ=255)                                                                            gCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTcc                                                                                                                                                               >  1:863499/1‑250 (MQ=255)                                                                                  tCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTAt                                                                                                                                                        >  1:504239/1‑251 (MQ=255)                                                                                  tCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTAt                                                                                                                                                        >  1:469598/1‑251 (MQ=255)                                                                                  tCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTAt                                                                                                                                                        >  1:415816/1‑251 (MQ=255)                                                                                     caGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCgg                                                                                                                                                     >  1:182479/1‑251 (MQ=255)                                                                                        aaCAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCgg                                                                                                                                                     <  1:685214/248‑1 (MQ=255)                                                                                          cAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGacac                                                                                                                                                                                 >  1:459467/1‑218 (MQ=255)                                                                                             cgacACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGac                                                                                                                                                                                   <  1:509332/210‑1 (MQ=255)                                                                                                    aGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACg                                                                                                                                                                                                                       >  1:857245/1‑170 (MQ=255)                                                                                                           gATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTgg                                                                                                                                                                                                                                               >  1:148096/1‑139 (MQ=255)                                                                                                                    cAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGa                                                                                                                      >  1:369100/1‑251 (MQ=255)                                                                                                                         aCGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGc                                                                                                                                  >  1:15365/1‑234 (MQ=255)                                                                                                                              gcTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTcc                                                                                                                                                                              >  1:470040/1‑185 (MQ=255)                                                                                                                              gcTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATcc                                                                                                                                                                    >  1:78740/1‑195 (MQ=255)                                                                                                                                 ctATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGAt                                                                                                                                                                                                         <  1:640467/155‑1 (MQ=255)                                                                                                                                          ctTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGac                                                                                                                                                                                   <  1:388875/168‑1 (MQ=255)                                                                                                                                              aTTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTATTCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACaaa                                                                                            >  1:750956/1‑251 (MQ=255)                                                                                                                                                        gTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGa                                                                                  >  1:333474/1‑251 (MQ=255)                                                                                                                                                          aaTTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGAt                                                                                                                                                                                                         >  1:447067/1‑130 (MQ=255)                                                                                                                                                          aaTTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATc                                                                                <  1:694862/251‑1 (MQ=255)                                                                                                                                                              tttatttaCCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGTCGGTAGCTGAAGCACAAAAAATTCAAGATCcgcg                                                                            >  1:249564/1‑251 (MQ=255)                                                                                                                                                                    taCCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTAt                                                                                                                                                        >  1:626248/1‑169 (MQ=255)                                                                                                                                                                           ggTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGaa                                                               >  1:682270/1‑251 (MQ=255)                                                                                                                                                                               aTGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAAcaca                                                           >  1:100596/1‑251 (MQ=255)                                                                                                                                                                                 ggATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCatat                                                                   <  1:813741/241‑1 (MQ=255)                                                                                                                                                                                 ggATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAac                                                              <  1:710069/246‑1 (MQ=255)                                                                                                                                                                                 ggATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAAt                                                         >  1:236139/1‑251 (MQ=255)                                                                                                                                                                                  gATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATa                                                        >  1:456045/1‑251 (MQ=255)                                                                                                                                                                                      aTGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACttt                                                    <  1:644573/251‑1 (MQ=255)                                                                                                                                                                                             aTTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGaa                                             >  1:666230/1‑251 (MQ=255)                                                                                                                                                                                              tttAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGc                                                                                                                                                                                                                                     >  1:543705/1‑66 (MQ=255)                                                                                                                                                                                              tttAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATcc                                                                                                                                                                    >  1:85594/1‑131 (MQ=255)                                                                                                                                                                                              tttAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATcgc                                                                                                                        >  1:807975/1‑175 (MQ=255)                                                                                                                                                                                                      ccTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCg                                                                                                            <  1:317089/179‑1 (MQ=255)                                                                                                                                                                                                      ccTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGAt                                                                                 <  1:871690/206‑1 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTatat                                                                                                                                                                                                             <  1:716670/78‑1 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGtt                                                                                                                                                                                                 >  1:153346/1‑90 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGtt                                                                                                                                        >  1:635450/1‑147 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCatat                                                                   <  1:681837/216‑1 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTg                                                                 <  1:206850/218‑1 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACttt                                                    <  1:755620/231‑1 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACttt                                                    <  1:305285/231‑1 (MQ=255)                                                                                                                                                                                                          ttcttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGAtt                                 >  1:328549/1‑250 (MQ=255)                                                                                                                                                                                                            cttcAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTg                                                                                                     <  1:135657/180‑1 (MQ=255)                                                                                                                                                                                                               cAACTGAGCATTTGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATT‑‑TATTGATGGTACGTTTGGGCGTGGTGGacac                                                                                                                                                                                 >  1:642043/1‑99 (MQ=255)                                                                                                                                                                                                                         tttGACGGTTTTATTACATGAAACCGTGGACGGCTTGGCGTCGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACTCTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACttt                                                    <  1:249048/216‑1 (MQ=255)                                                                                                                                                                                                                                gTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCgg          >  1:410001/1‑251 (MQ=255)                                                                                                                                                                                                                                gTTTTATTACATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCgg          >  1:358534/1‑251 (MQ=255)                                                                                                                                                                                                                                     attaCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCGGCAAgg     <  1:327557/251‑1 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTcc                                                                                                                                                               >  1:418632/1‑94 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCat                                                                     <  1:814305/184‑1 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACttt                                                    >  1:653618/1‑201 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTtg                               >  1:641692/1‑222 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGt                  >  1:582277/1‑235 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCGGCAAGGTTg  <  1:212878/251‑1 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCGGCAAGGTTg  >  1:271437/1‑251 (MQ=255)                                                                                                                                                                                                                                        aCATGAAACCGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCGGCAAGGTTg  >  1:778426/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                  GAACAAGACATGGCATTAGGTGCAACGGATACTTTTGCATTATCGGAAAAATTAAAAACGCTGTCATTATAAGAGGCAAACTCACAGAACAGCGGCACTAGCAAAAGATAAGCGGCAGTTACGCCGCTCTATCTTCTCTTTATTTTACATTGTAATTTTTATTTACCATTGGTAATGGATTATGACGCATTTAAATACCTTTTCTTCAACTGAGCATTTGACGGTTTTATTACATGAAACTGTGGACGGCTTGGCGTTGAAAGAAAACGGCATTTATATTGATGGTACGTTTGGGCGTGGTGGACACTCCCGTTTAATCCTTTCCCAATTATCGGTGAACGGCAAGTTGATTGCCATCGATCGCGATCCGCAGGCGGTAGCTGAAGCACAAAAAATTCAAGATCCGCGTTTTCATATTGAACACAATACTTTCTCGGAAATCCTGCCGATTTGTGAAAAATTCGGTTTAGTCGGCAAGGTTG  >  1/898798‑899279

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