Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,626,668	A→C	16.6%	T253P (ACC→CCC)	radA →	DNA repair protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,626,668	0	A	C	16.6%	199.8 / 3.1	79	T253P (ACC→CCC)	radA	DNA repair protein
Reads supporting (aligned to +/- strand):  ref base A (45/18);  new base C (8/8);  total (53/26)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCG  >  NC_000913/4626437‑4626894                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    atttgctggttCCGGGTGGTGGATGTGGGGGATGGCCGTTCCTCGGCTGGTGGCGGGGCGGGGGGGCTTGAAACGGCGGTTTATTTGACCCGCTTCGCCAAAACGGGCGGGGTGGCGATTGGCATGGTGGGGCCCGTAACCAAAGATGTTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTgcg                                                                                                                                                                                                                              <  2:796427/229‑1 (MQ=255) atgggctggATCCGGGTGACGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGGGCAGGTGCGTGAAACGGCGGCTTATTTGACGCGCTTCGCCAAACCGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTgcg                                                                                                                                                                                                                              <  1:949309/231‑1 (MQ=255)     aCTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACCGTTCGGGGC‑TTTGGATGGCGATG‑CCGG‑CTCCCGTTTTTGCACCTTGCCCAG‑cc                                                                                                                                                                                                                        >  1:383466/1‑239 (MQ=255)         gATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAATGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGC‑CTCCCGCTTTCGCCCCTTGCCCAG‑CCATa                                                                                                                                                                                                                     >  1:883404/1‑239 (MQ=255)          aTCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGGGCTGGAACACTGGATCGACTGTTCGGGGGTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTTGCCAG‑CCTTaa                                                                                                                                                                                                                    >  1:625986/1‑239 (MQ=255)              ggggaatacATTTGGGGGATGAACAGTCATCGCCTGGGAGCGGGGCGGAGGGGCGGGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAAc                                                                                                                                                                                                                <  1:722525/231‑1 (MQ=255)               ggTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGAAG‑CCGC‑CTCCCGTTTTCGACCCTTGCGCCG‑CCATAAAAAc                                                                                                                                                                                                                >  1:1007850/1‑238 (MQ=255)                gTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGGGTGGCGGTTGTCATGGGGGGGCACCTAACCAAAAAAGGTTTCCTGGCTGGGCCGAAAGGGGCGGGACCCGGGTTTGGCCGTTTGGGGGTTTTTGGTGGCGATG‑CCGA‑CTCCCCTTTTCCCCCCTTGCGCCC‑CCACaaaaaacc                                                                                                                                                                                                              >  2:65266/1‑236 (MQ=255)                   aTGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTGGGGCACGTAACCCAAAATGGTTTCCTGGCTGGCCCGAAAGGGCTGGAACACTGTATCGCCTGTTCGGGGCTTTTGGGGGGGGGTG‑ACCG‑CTCCCCTTTTTCCCCCCTTCGCAGCCCAAAAAAAcccg                                                                                                                                                                                                             >  2:259945/1‑236 (MQ=255)                    tGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGGGCTGGAACACTGTATCGCCTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGGGCAG‑CCCTAAAAACCGCCTc                                                                                                                                                                                                          >  1:413042/1‑239 (MQ=255)                     gCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGAGTGTCATGGTGGGGCGCGTAACCAAAGATGGTTCGCCGGCGGGGCCCAAAAGGCTGGAACACTGTTTCAACCGTTTGGTGCTTTTGGGTGGGGAAGCCCGC‑CCCCCGTTTT‑GCCCCTTGCCCAA‑CCATAAAAAACGCCTCg                                                                                                                                                                                                         >  2:757863/1‑239 (MQ=255)                     gCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTGGGGGACGTAACCCAAAAAGGTTCGCTGGCTGGCCCGAAAGTGCTGGGAACCTGGATACACTGTTCCGGGCTTTTGGAGGGCGGTG‑CCCACCCCCCGTTTTC‑CCCCTTGCGAAG‑CCATaaaaaaagcatcg                                                                                                                                                                                                         >  2:91860/1‑231 (MQ=255)                                           tCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTTTTCGGGGCTTTTGGATGGCGATG‑CCGA‑CCCCCGTTTTCCCACCTTGCGCAG‑CCATAAAAACCGCCTCGGCGCGGGGACTGaccgtgcg                                                                                                                                                                                     >  1:490160/1‑230 (MQ=255)                                                 gggtgggtggtgCGGGTGGGTGAACGGGGGGTTTTTTTGCCCCGCTGCGCAAAAGCGGGGGGTGTGGCGATGGCAATGGGGGGGCCCGTAACCAAAGAGGGTTCGCTGGCGGGCCCGAAAGTGCTGGAACACTGTATCGCCTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCCCCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTcgc                                                                                                                                                                             <  2:414977/233‑1 (MQ=255)                                                  gTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCCAAACGCGCGGGGTGGCGATTGTCATGTGGGGGCACGTAACCAACGATGGTTCGCTGGCTGGACCGAAATGGCTGGGCCACTCATTCGACCGTTCGGGCGTTTTGGATTGCGATG‑GCGA‑ACCCCGTTTTCGCCACCTGCGCAG‑CCATAAAAACCCCTTCGGGGCGGGGAAAGCGGGGGGCGGCTTcggc                                                                                                                                                                            >  2:573589/1‑237 (MQ=255)                                                  gTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGGTTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTCCGCAG‑CCATAAAAACCGCTTAGGCACGGTGAATGAGCCGGGCGTCTTcc                                                                                                                                                                              >  1:559538/1‑236 (MQ=255)                                                  gTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGGGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGGGCTTTTGGATGGCGATG‑CCGA‑GTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGCCGCGGTGGATGAGCTGGGGGTCTTccgg                                                                                                                                                                            >  1:951387/1‑236 (MQ=255)                                                       gTGGCGGGGGTGCGTGAAACGGGGGTTTATTTGACACGCTTCGCCAAAACGGGGGGTGGGGGGATTGTCAGGGGGGGGCACGTAACCAAAGAGGGTTCGCTGGCTGGCCAGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CCCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAc                                                                                                                                                                       <  2:850650/239‑1 (MQ=255)                                                        tgggggaGGTGGGTAAAAGGGCGGCTTATTTGACACGTTTCGCCAAAACGCGGGGTGGGGCGATTGTCAGGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGAc                                                                                                                                                                       <  1:17083/232‑1 (MQ=255)                                                        tGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGGGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCCCGGTGAATGAGCTGGGCGGCTTCGCGATGAcc                                                                                                                                                                      >  1:717211/1‑239 (MQ=255)                                                        tGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGGTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGGATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CGCCCGTTTTCGCCACGTGTGCCG‑CCACAAAAAACGCCTCGGCCCGGGGAATGGGCCGGGCGGCttgcggcagaac                                                                                                                                                                      >  2:342487/1‑229 (MQ=255)                                                           cgcAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCCCTTCGGCGCGGGGAATTAGCTGGGCGGCTTCCCGATACCCGaa                                                                                                                                                                   >  1:101517/1‑239 (MQ=255)                                                            gcAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCGGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGaa                                                                                                                                                                   <  1:13379/238‑1 (MQ=255)                                                             cAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGGGCCTTGCCGCTGCCCGAACa                                                                                                                                                                 >  1:411626/1‑239 (MQ=255)                                                             cAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGGGAATGAGCTGGGCGGCTTCGCGGTGACCGAACa                                                                                                                                                                 >  1:750608/1‑239 (MQ=255)                                                              aGGTGTGTGAACCGGCGGGTTATTTGCCACGCTTCGCAAAAACGGGCGGTGTGGCGATTGTCATGGTGGGGCCCGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTGGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACa                                                                                                                                                                 <  1:91860/238‑1 (MQ=255)                                                              aGGTGCGTGAAACGGCGGTTTATTTGACACGCTTCGCCAAACCGGGGGGTGTGGCGATTGTCAGGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACa                                                                                                                                                                 <  1:259945/238‑1 (MQ=255)                                                                        aaCGGGGGCTTATTTGCCACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCCCCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTg                                                                                                                                                      <  1:65266/239‑1 (MQ=255)                                                                                 tttttcgcaccGCTTCGCCAAGACGGGGGTGGTGGGGATGGTCGTGGGGGGGCTCGTAACCAAAGGTGGGTCGCTGGCTGGGCAGAAAGTGCTGGAACACTTTATAGACTGTTCGGTGCTTTTGGAGGGGGATG‑CCGA‑CTCCCGTTTTCGCCCCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc                                                                                                                                              <  2:274846/228‑1 (MQ=255)                                                                                                          cgcgGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCCGGAACACTGTATCGACTGTTCGGGGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTCCGCAG‑CCATAAAAACCGCTTCGGCGCGGGGAATGCGCTGGGCGCCTTCCCGATGACCGAAAAGGGGCGCCGTGAAGCCGGCAACCCTCCGCCAATTTTTTTaacccc                                                                                                                    >  2:686581/1‑235 (MQ=255)                                                                                                           ggggTGTGGGGATTGTCAGGGGGGGGCACGTAACCAAAGTGGGTTCGCTGGCTGCCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTcg                                                                                                                    <  1:566184/236‑1 (MQ=255)                                                                                                                      tgtgacagggtggGGAACGTAACAAAGGTTGGTTCGCTGGGTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTCCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAg                                                                                                         <  1:754466/230‑1 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGGCCGGACAGGGGCTGCGTGAAGTCCGCAACCCTTCGGCAATTTTCTTAATGCGCGGCG‑TGAAggggg                                                                                                     >  1:390994/1‑234 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGACCAGGGGCTGCGTGAGTTCAGCAACCCTTCGGCAATTTTCTTAATTCGCGGCGATGAAGGGa                                                                                                      >  1:831747/1‑238 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGa                                                                                                      >  1:343699/1‑238 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCCCGGCGATGAAGTGa                                                                                                      >  1:1066769/1‑238 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCGCCGGGAAGCCAGCAACCCTCCGGCAATTTCCTTAAGTTCCGCCGATGAAGGGa                                                                                                      >  1:751060/1‑238 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGATGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGa                                                                                                      >  1:260210/1‑238 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGAGGACCGAACAGGGGCTGCGTGAAGTCAGAAACCCTTCGGCAATTTTCTTAAGCCGCGGCGATGagttga                                                                                                      >  1:966508/1‑233 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGCCTTCGCGATGACCGAACAGGGGCTGCGTGAAGCAAGCAACCCTTCGGCAATTTTCTTAAGTCgcggcgtggaatggc                                                                                                      >  1:692758/1‑229 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGACTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTGGCGGCGATGAAGTGa                                                                                                      >  1:614236/1‑238 (MQ=255)                                                                                                                         gTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGCGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCCG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTTGGCGTCTTCGCGATGGCCGAACAGAGGCTGCGTGAAGCCAGCAACCCTTCTGCAATTTTCTTAAGTCCCGGCGTTGACGTGa                                                                                                      >  1:770335/1‑238 (MQ=255)                                                                                                                          tacgggtggGGCACGTACCAAAGGGGGGTTCGGTGGCTGGCCGGAAGGGGCGGGACCCTTGTTCCGACTGTTCGGGCCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTGGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGAcc                                                                                                    <  2:106998/235‑1 (MQ=255)                                                                                                                           cATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGAcc                                                                                                    >  1:541510/1‑238 (MQ=255)                                                                                                                                  gggCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCGGCGGGAAGTCAGCAACCCTTCGGCCATTTTCTTAAGACGCGGCGAGGAAGTGCACCCCGGTa                                                                                             >  2:99814/1‑238 (MQ=255)                                                                                                                                    gCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGACTGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGGCGCGGCGATGAAGTGAGCTCCGCTAGc                                                                                           >  1:996306/1‑238 (MQ=255)                                                                                                                                    gCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGATGAAGTGAACTCCGGTAGc                                                                                           >  1:516554/1‑238 (MQ=255)                                                                                                                                    gCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGCCGCGGAGTTGACCTCCGGTAGc                                                                                           >  1:606202/1‑238 (MQ=255)                                                                                                                                             cAAAGATGGTTCGCGGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGAtgg                                                                                 <  1:757863/239‑1 (MQ=255)                                                                                                                                                      ttCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGGAGCCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATggg                                                                        >  1:468206/1‑239 (MQ=255)                                                                                                                                                      ttCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTGGCCAAGGGTTGGGGGTTTggg                                                                        >  1:516063/1‑239 (MQ=255)                                                                                                                                                      ttCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCCCGGCGTTGAGAGGCCCTCCGGTGGCCCAGGGATGGGGGTGTggg                                                                        >  2:583318/1‑239 (MQ=255)                                                                                                                                                        cGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGGGAGGAAGTGCACCCCGGAAGCCCAGAGAGGggggaaggggaa                                                                      >  2:528225/1‑239 (MQ=255)                                                                                                                                                                ggcCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGGGGAATGGGAGGGAAcgc                                                               >  1:324504/1‑238 (MQ=255)                                                                                                                                                                    ggAAAGTGCGGGACCACTGTATCGACTGTTCGGTGCTTTTGGATGGCGAGG‑GCGA‑CCCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTcc                                                           <  1:835901/237‑1 (MQ=255)                                                                                                                                                                                    cTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGGGAAGTCAGCAACCCTTCGGAAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCCGTGATGGTGGTATGGGAAGGAACGCGTCAACTGCTGGGGGAGATTc                                          >  1:984281/1‑239 (MQ=255)                                                                                                                                                                                     tGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGACGTGACCTCCGGTAGCTCCGTGCTGGTGGTATGGGAAGGAACGCGGCCACTCCtggggggggtgtc                                         >  2:648970/1‑232 (MQ=255)                                                                                                                                                                                     tGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCCCTGCTGGTGGGGATTCa                                         >  1:931769/1‑239 (MQ=255)                                                                                                                                                                                      gTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATg                                                                                                            >  1:925379/1‑171 (MQ=255)                                                                                                                                                                                       tATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc                                                                                                                                              >  2:1040276/1‑136 (MQ=255)                                                                                                                                                                                       tATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGc                                                                                                                                              <  1:1040276/136‑1 (MQ=255)                                                                                                                                                                                               gtggtggtcTTTTGGCTGGCGACG‑GCGA‑CCCCCGTTTTCGCACCCTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTc                               <  1:793967/231‑1 (MQ=255)                                                                                                                                                                                               gTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGGGGCCCCTGCTGGGGGGGATTAAGGGGCTGGtg                               >  2:681539/1‑238 (MQ=255)                                                                                                                                                                                               gTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTc                               >  1:661783/1‑239 (MQ=255)                                                                                                                                                                                               gTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCGGGTc                               >  1:130934/1‑239 (MQ=255)                                                                                                                                                                                               gTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCCGGGCGGTGAGGTGACCTCCGGTGACTCAGGGATGGGGGGATGGGAAGGAAGCCGGCCACCGCTGGGGGGGATTCAGGGGCTgggg                               >  2:141559/1‑237 (MQ=255)                                                                                                                                                                                                          ttGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGGAGGAACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGAGCTGGTCGATCACTCgat                    >  1:943142/1‑239 (MQ=255)                                                                                                                                                                                                          ttGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGCAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGGCCACTGCTGGTGGAGATTCAGGGGCTGGTCGATCACTCgag                    >  1:499641/1‑238 (MQ=255)                                                                                                                                                                                                          ttGAATGGGGATG‑CCGA‑CCCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCgat                    <  1:583318/239‑1 (MQ=255)                                                                                                                                                                                                           tGGTTGGCGGGG‑CCGC‑CCCCGGTTTTCGCCCCTGGCGCAG‑GCATAAAAACCGGTTCGGGGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCgatg                   <  2:1007850/239‑1 (MQ=255)                                                                                                                                                                                                               gggggctG‑CCTA‑CCCCCTTTTTCGCCTCTGGGCCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATTACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATgg               <  2:883404/233‑1 (MQ=255)                                                                                                                                                                                                                    caaggcggccTCCCGTTTCCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCACCCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCg             <  1:573589/229‑1 (MQ=255)                                                                                                                                                                                                                        ccGA‑CTCCCGTTTTCGCCCCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGAt                                                                                                             <  2:222031/137‑1 (MQ=255)                                                                                                                                                                                                                        ccGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGAt                                                                                                             >  1:222031/1‑137 (MQ=255)                                                                                                                                                                                                                        ccGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGATGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAcccc        >  1:413468/1‑238 (MQ=255)                                                                                                                                                                                                                           cccccccGTTTTCCGCCTTTGCCCCG‑CCTAAAAAACCGCTTCGGCGCGGTGGATGAGCTGGGCGGCTTCGCGATGACCGAACCGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCAcg     <  2:657271/235‑1 (MQ=255)                                                                                                                                                                                                                             cTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCAcccc   >  1:56821/1‑239 (MQ=255)                                                                                                                                                                                                                              tCCTGTTTTCCCCCCTTGCCCAG‑CCATAAAAACCGCTTCGGGGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCAcgccg  <  2:130934/239‑1 (MQ=255)                                                                                                                                                                                                                              tCCCTTTTTCAGACCTGGGGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCAcgccg  <  1:686581/239‑1 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                    ATTGACTCGATCCAGGTGATGCATATGGCGGATGTACAGTCATCGCCTGGTAGCGTGGCGCAGGTGCGTGAAACGGCGGCTTATTTGACACGCTTCGCCAAAACGCGCGGTGTGGCGATTGTCATGGTGGGGCACGTAACCAAAGATGGTTCGCTGGCTGGCCCGAAAGTGCTGGAACACTGTATCGACTGTTCGGTGCTTTTGGATGGCGATG‑CCGA‑CTCCCGTTTTCGCACCTTGCGCAG‑CCATAAAAACCGCTTCGGCGCGGTGAATGAGCTGGGCGTCTTCGCGATGACCGAACAGGGGCTGCGTGAAGTCAGCAACCCTTCGGCAATTTTCTTAAGTCGCGGCGATGAAGTGACCTCCGGTAGCTCAGTGATGGTGGTATGGGAAGGAACGCGTCCACTGCTGGTGGAGATTCAGGCGCTGGTCGATCACTCGATGATGGCGAACCCACGCCG  >  NC_000913/4626437‑4626894

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