Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	648,741	A→G	intergenic (+151188/‑50602)	unknown → / → unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	648,741	0	A	G	100.0%	137.0 / NA	50	intergenic (+151188/‑50602)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (33/17);  total (33/17)

CATTGTAACTGTTTATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACAGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGTGCAACTTGTTGCACCAAAATATTTAAAATCGGTGCAACAAGTTGCACCCTACCA  >  1/648502‑648983                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   cATTGTAACTGTTTATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCaa                                                                                                                                                                                                                                         >  1:64470/1‑251 (MQ=255) cATTGTAACTGTTTATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGA‑GGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAg                                                                                                                                                                                                                                        >  1:137484/1‑251 (MQ=255)     gTAACTGTTTATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGata                                                                                                                                                                                                                                     >  1:59398/1‑251 (MQ=255)           gTTTATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTaa                                                                                                                                                                                                                               >  1:418965/1‑251 (MQ=255)             ttATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGa                                                                                                                                                                                                                             <  1:751922/251‑1 (MQ=255)                     aaaaaCTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATa                                                                                                                                                                                                                     >  1:106727/1‑251 (MQ=255)                         aCTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAAtt                                                                                                                                                                                                                 >  1:418420/1‑251 (MQ=255)                              gCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATaa                                                                                                                                                                                                             <  1:744363/250‑1 (MQ=255)                               cTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATg                                                                                                                                                                                                           <  1:17128/251‑1 (MQ=255)                                 ccccTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGtt                                                                                                                                                                                                         >  1:773345/1‑251 (MQ=255)                                      cGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGaa                                                                                                                                                                                                    >  1:429898/1‑251 (MQ=255)                                                       gggCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGa                                                                                                                                                                                   >  1:211182/1‑251 (MQ=255)                                                              gCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTg                                                                                                                                                                             >  1:465208/1‑250 (MQ=255)                                                               cccTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCt                                                                                                                                                                           <  1:658484/251‑1 (MQ=255)                                                                                 ctcAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTta                                                                                                                                                         <  1:180241/251‑1 (MQ=255)                                                                                      tttttAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAgg                                                                                                                                                    <  1:844594/251‑1 (MQ=255)                                                                                           aaGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGcaca                                                                                                                                                >  1:34102/1‑250 (MQ=255)                                                                                                    cTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACt                                                                                                                                                                                 <  1:315861/208‑1 (MQ=255)                                                                                                    cTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGc                                                                                                                                         >  1:464262/1‑248 (MQ=255)                                                                                                    cTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGAGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGa                                                                                                                                      >  1:475693/1‑251 (MQ=255)                                                                                                      tgtgAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTtat                                                                                                                                                        <  1:30540/231‑1 (MQ=255)                                                                                                      tgtgAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGa                                                                                                                                      >  1:389676/1‑249 (MQ=255)                                                                                                       gtgAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACt                                                                                                                                    >  1:344856/1‑250 (MQ=255)                                                                                                         gAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACggagga                                                                                                                                                                                                                                               >  1:31057/1‑141 (MQ=255)                                                                                                                    ggatcttgcctccTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGaa                                                                                                                      >  1:260431/12‑251 (MQ=255)                                                                                                                       aaaTTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCa                                                                                                                   >  1:83277/1‑251 (MQ=255)                                                                                                                       aaaTTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCa                                                                                                                   >  1:532777/1‑251 (MQ=255)                                                                                                                        aaTTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGt                                                                                                                                             >  1:575272/1‑224 (MQ=255)                                                                                                                         aTTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAg                                                                                                                 >  1:821476/1‑251 (MQ=255)                                                                                                                              cccTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCt                                                                                                                                                                           >  1:653669/1‑188 (MQ=255)                                                                                                                                        cTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAAt                                                                                                                                                                        <  1:707779/181‑1 (MQ=255)                                                                                                                                         ttCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAAt                                                                                                 >  1:851543/1‑251 (MQ=255)                                                                                                                                                                   aaaCTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGgcg                                                                        <  1:398401/250‑1 (MQ=255)                                                                                                                                                                            tactcaCCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATg                                                              >  1:211716/6‑251 (MQ=255)                                                                                                                                                                            tGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGaaa                                                                                            >  1:136431/1‑221 (MQ=255)                                                                                                                                                                             gCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGcac                                                                                                                                                 >  1:49057/1‑167 (MQ=255)                                                                                                                                                                             gCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACt                                                                                                                                    <  1:742038/180‑1 (MQ=255)                                                                                                                                                                             actcaCCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGt                                                                                                                                                                                                          <  1:734186/106‑1 (MQ=255)                                                                                                                                                                                 aCCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAggg                                                         <  1:44256/251‑1 (MQ=255)                                                                                                                                                                                  ccccTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGc                                                                                                                                                                                                                                           >  1:52328/1‑72 (MQ=255)                                                                                                                                                                                   cccTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGt                                                             <  1:106779/245‑1 (MQ=255)                                                                                                                                                                                    ccTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAAc                                                                                                    >  1:360850/1‑205 (MQ=255)                                                                                                                                                                                              cGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGt                                                        <  1:347698/239‑1 (MQ=255)                                                                                                                                                                                                  ggCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGTGCAACTTGTTGCACCa                                        >  1:637788/1‑251 (MQ=255)                                                                                                                                                                                                                     aCGCATTGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGt                                                             >  1:657567/1‑211 (MQ=255)                                                                                                                                                                                                                          ttGCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGTGCAACTTGTTGCACCAAAATATTTAAAATCGGTGcaacaa                <  1:850179/251‑1 (MQ=255)                                                                                                                                                                                                                            gCGTTTTGTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAg                                                           <  1:170311/206‑1 (MQ=255)                                                                                                                                                                                                                                   gTCCCCCGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGTGCAACTTGTTGCAc                                          >  1:410216/1‑216 (MQ=255)                                                                                                                                                                                                                                       cccGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAAt                                                                                                                     >  1:715573/1‑137 (MQ=255)                                                                                                                                                                                                                                        ccGTAAACGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGTGCAACTTGTTGCACCAAAATATTTAAAATCGGTGCAACAAGTTGCACCCTACCa  >  1:374398/1‑251 (MQ=255)                                                                                                                                                                                                                                            aaaCGGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGt                                                                                                                >  1:341011/1‑137 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   CATTGTAACTGTTTATATAAAAAAACTTGGCTCCCCTCGTTCGCGGGGGGAGCTGGGCGAAGCCCTGAGGGGGGCAAAATCTCAATTTTTAAGCATGATCTTGTGAATACATCGAGATAAATTCCCCCTTCCGCCCTTCGGGTACGTTATTCGCAAGCTCGGAAACTCACTTGCTTACCCCTTTGGGGTCGTTGGCAGAACCAACGCTCAAAACGCATTGCGTTTTGTCCCCCGTAAACAGAGGAAGGCAAGATATTTTAAGATATCAATAAATTATAATGTTTCGAACAGCTGTGGTTTCCCGACTCGTGCTAATGGCTCACTAAATTTATAAGGCACACGTCGGCAGACTCGCGGCGGTATGAATCAAGTGAATTAAAAAACAATTGAAAAAACGACCGCACTTTTGGCGCGTCGTTATGTAGGGTGCAACTTGTTGCACCAAAATATTTAAAATCGGTGCAACAAGTTGCACCCTACCA  >  1/648502‑648983

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