Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,021,098	T→C	14.8%	G42G (GGA→GGG)	Y75_RS19860 ←	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,021,098	0	T	C	14.8%	223.6 / 3.2	90	G42G (GGA→GGG)	Y75_RS19860	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base T (61/11);  new base C (2/14);  total (63/27)
Fisher's exact test for biased strand distribution p-value = 6.16e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGTTACA  >  NC_007779/4020906‑4021335                                                                                                                                                                                                   |                                                                                                                                                                                                                                               cAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATAc                                                                                                                                                                                                 >  2:1506155/1‑239 (MQ=11)        aaGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTAccc                                                                                                                                                                                           <  1:1506155/238‑1 (MQ=32)        aaGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTAccc                                                                                                                                                                                           <  2:2129055/238‑1 (MQ=32)        aaGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGGTTTTGTCCGGTAACTGCTGTGGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGTACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTAccc                                                                                                                                                                                           <  1:682313/238‑1 (MQ=17)             gATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCt                                                                                                                                                                                      <  1:1459953/238‑1 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACC‑CCGAACAGGCCCCGCCGGTGGGGGCACCAATAAGTACCCCGGCTGa                                                                                                                                                                                    >  2:903158/1‑239 (MQ=17)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGCTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGa                                                                                                                                                                                    >  1:1491602/1‑239 (MQ=17)                 tCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAc                                                                                                                                                                                  >  1:1779656/1‑238 (MQ=32)                 tCCGCAGGCATTTGCCAGCCGGGGCCGAGGCGCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTGCCCCGGCTGAAcc                                                                                                                                                                                 <  1:665656/239‑1 (MQ=17)                     cAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAcc                                                                                                                                                                                 >  2:457342/1‑235 (MQ=32)                     cAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAcc                                                                                                                                                                                 <  1:457342/235‑1 (MQ=32)                        gCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAAc                                                                                                                                                                           <  1:634547/238‑1 (MQ=32)                        gCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGCCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGCTGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACCCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACa                                                                                                                                                                          >  2:619431/1‑239 (MQ=17)                         cATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGaa                                                                                                                                                                            <  2:591737/236‑1 (MQ=32)                                        gCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCAt                                                                                                                                                             >  1:1409281/1‑236 (MQ=32)                                               ccccccACGGGCGCGGGCTTTTTTGGCGGGTACGTGGGGTGGGCGGGGGGGGGGGGGAAGGGCAGGGGGCGGGGCAGGGGGATGGCGGTTCCCCGGGGAAGGCCCTTGGCCCCGGGCAGGGGATCAG‑CGGGGTGGCCGGC‑CGTCACCCCGCCGGGGCACC‑CCAAACCGGCCACGCGGGGGGGGCCACCAATCCGTCCCCCGGCTAAACCCTGAACAATGCTACCGCCATACTGCGTCa                                                                                                                                                    <  2:872584/236‑1 (MQ=2)                                                cccccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACGGCGTCATg                                                                                                                                                  >  1:439857/2‑239 (MQ=17)                                                  cccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTc                                                                                                                                                >  1:176576/1‑239 (MQ=32)                                                  cccACCGGCGCGGGCGTTTTTGTCCGGAAACTGCTGTAGGCGCGGGAGAGGATGAACGGCAGCGCGCCGGGCAGGGCGAAGTCGGTTTCAACGGGCAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGGCGCCCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCCACAATACGTACCCCGGCTGAACCCCGTACAATGCTACCGCCAGACGGCGTCATGt                                                                                                                                                 >  2:1325528/1‑238 (MQ=11)                                                     aCCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGAAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCAc                                                                                                                                              >  1:247509/1‑238 (MQ=18)                                                       cGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCACTACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCAc                                                                                                                                              >  2:1193380/1‑236 (MQ=18)                                                          cgcgGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGAc                                                                                                                                        >  1:1086118/1‑239 (MQ=33)                                                            cgGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACACTGAACAATGCTACCGCCATACTGCCTCATGTCACCCTGCcg                                                                                                                                       >  2:635552/1‑238 (MQ=18)                                                              ggCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGAcgc                                                                                                                                      >  1:637647/1‑237 (MQ=14)                                                                   ttttGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGccgccg                                                                                                                               >  1:925389/1‑239 (MQ=35)                                                                        tCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGttt                                                                                                                           >  1:1390588/1‑238 (MQ=35)                                                                             tAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATGCTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGc                                                                                                                      >  1:635450/1‑238 (MQ=25)                                                                                ctgctgTGGGGGGGGGGGGGGAGGGACGGCGGGGGGCGGGGGGGGGGGAGGTCGGTTCCCCGGGGAGGGCCCTTTGCCCCGGGCGGGGGATTGA‑CGGGAGGGCCGGA‑CGTCACCCCGCCGGGGCCCA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCa                                                                                                                   <  2:925389/238‑1 (MQ=11)                                                                                          ggggggAGAGGATGGACGGCGGGGGGGCGGGCGGGGCGAGGTCGGTTCCACCGGGAAGGTCCTTGGCACCGAGCAGGGGTTTGC‑CCGGGTGGCCGGC‑CGCCCCCCCGCGGGCGCACA‑CCAACCAGGCCACGCCGGTGGGGGCACCAATACGTCCCCCGGCTGACCCCTCAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATc                                                                                                         <  2:627336/235‑1 (MQ=14)                                                                                           gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCt                                                                                                       >  1:650413/1‑239 (MQ=40)                                                                                           gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCt                                                                                                       >  2:634676/1‑239 (MQ=40)                                                                                             gggAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCACTACGTACCCCGGCTGAACCCTGAACAATGCAACCGCCATACTGCGTCATGTCACCCTGACGCGCCGACGGCTTTCCGCCCATCGTAAAATctgt                                                                                                      >  1:1170684/1‑235 (MQ=255)                                                                                             gggAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCttt                                                                                                     >  1:564987/1‑239 (MQ=255)                                                                                             gggAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGGTTTCCGCTCATCGTACAATCCttt                                                                                                     >  1:1838863/1‑239 (MQ=255)                                                                                             gggAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGATTTCCGCTCATCGTACAATCCttt                                                                                                     >  1:102911/1‑239 (MQ=255)                                                                                                    gATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTCAAATCCTTTTCtatta                                                                                              >  1:522238/1‑239 (MQ=255)                                                                                                    gATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCtatta                                                                                              >  1:315825/1‑239 (MQ=255)                                                                                                     aTGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGC‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTCCGTCATGCCACCCTGCCCCGCCGCCGGTTTTCCGCCCATCGTAAATACTTTTTCtattat                                                                                             >  2:825240/1‑239 (MQ=14)                                                                                                         aCGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                         >  1:1166693/1‑239 (MQ=255)                                                                                                         aCGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                         >  1:1440700/1‑239 (MQ=255)                                                                                                         aCGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                         >  1:1104865/1‑239 (MQ=255)                                                                                                            cgggggggcggggcggggcgTGGTCGGTTCCCCCGGGAGGGCCCTTGGCCCCGAGGGGGGGATTGA‑CGGGAGGGCTGGC‑CGTCACCCCGCGGGGCCACC‑CGAAACAGGCCCCGCCGGTGGGGCCACCAATACGTCCCCCGGCTGAACCCTGAACAATTCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAAt                                                                                       <  2:744124/233‑1 (MQ=17)                                                                                                             cAGCGGGCCGGGCAAGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATCTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAATATCCTTTTCTATTATTTTAATg                                                                                      >  1:354823/1‑238 (MQ=255)                                                                                                              aGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGAc                                                                                    >  1:684385/1‑239 (MQ=255)                                                                                                              aGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGAc                                                                                    >  1:1416591/1‑239 (MQ=255)                                                                                                              aGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTATGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGTGCACA‑CCGAACAGGCCACGCCGGTGAGGGCGCCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCCGTTTTCCGATCCTCGTAAAAACCTTTTCTATTATTTTAATGa                                                                                     >  2:1646936/1‑238 (MQ=14)                                                                                                                 gggcggggcgGGGAGAGTTGGGTTCCGCCGGGAAGGACCTTGGCCCCGAGCAGGGGGTGGA‑CCGGATGGCGGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACGGGCCACGCCGGTGGGGGCACCAATCCGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAg                                                                                  <  2:662679/238‑1 (MQ=17)                                                                                                                  ggccgggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTACTGACAGTa                                                                                >  1:1977483/1‑239 (MQ=255)                                                                                                                  ggccgggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGCCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTa                                                                                >  1:1594819/1‑239 (MQ=255)                                                                                                                       ggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTc                                                                           >  2:898703/1‑239 (MQ=255)                                                                                                                            ggCGATGTCGGTTTCACCGGGAAGGACCTTTGCACAGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACg                                                                      >  1:1517023/1‑239 (MQ=255)                                                                                                                             gCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGa                                                                     >  1:91010/1‑239 (MQ=255)                                                                                                                                 tGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGaaaaa                                                                 >  1:1876258/1‑239 (MQ=255)                                                                                                                                  gTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCCGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGGACCAATATGTACCCCGGCTGAACCCTGAACAATTCTACCGCCATCCCGCGCCAGGTCCGCCTGACGCGCCGGCGGTTTCCCGCTCATAGGAACATTCTTCTCTATTATTCTCACGACAGTATCTTcctcgcaaacag                                                                >  2:409432/1‑228 (MQ=14)                                                                                                                                   tgggctccgccGGGAGGCACTTTTGC‑CCGGGGCGGGGGTTGACCCGGAGGGCGGGA‑CACCCCGCGGCCGGGCCACA‑CCGACCAGGCCACGCCGGGGGGGCCACCAATCCGTACCCCGGCGAAACCCTGAACAATCCTCCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAg                                                                <  2:1838863/229‑1 (MQ=18)                                                                                                                                       tttCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGT‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAACATCCTTTTCTATTATTTTAATGAAAGTATCTCCATACGAAAACGCgacg                                                           >  2:1900246/1‑237 (MQ=255)                                                                                                                                           aCCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGaaaa                                                       >  1:706199/1‑239 (MQ=255)                                                                                                                                           aCCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGaaaa                                                       >  1:932926/1‑239 (MQ=255)                                                                                                                                             cGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGaaaa                                                       >  1:535217/1‑237 (MQ=255)                                                                                                                                               ggAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACtt                                                   >  1:28099/1‑239 (MQ=255)                                                                                                                                                    gtccTTTGCCCAGAGCGGGGGATTG‑‑CCGGGCGGCCGGAGCCGCACCCCGCCGGGGACCC‑CCGAAAAGCCCACGCGGGTGGGCGCCCCAATCCGTCCCCGGGCTGACCCCAGAACAAGGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTgg                                               <  2:1379318/236‑1 (MQ=21)                                                                                                                                                    gccttgggacccGGGCGGGGGATAGC‑CCGGGTGGCGGGG‑CGCCCCCCCGCGGGAGCACC‑CCCAACAGGCCACGCGGGGGGGGCCCCCAATCCGTCCCCCGGCTAAACCCTAAACATTCTTCCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTgg                                               <  2:684385/228‑1 (MQ=21)                                                                                                                                                     ccctgtccccggAGCGGGGGCTTGC‑CCGGGTGGCGGGA‑CGTCCCCCCG‑CGGGGCCCACCCGACCGGGCCCCGCGGGTGGGGCAACCAATCCTCCCCCCGGCTGACCCCGAAAAAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAATATCTTCATACGAAAAAGCGAAGAAAAACTTGTgg                                               <  2:559689/226‑1 (MQ=18)                                                                                                                                                     aCCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTc                                             >  1:1906906/1‑239 (MQ=255)                                                                                                                                                               cGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAAt                                                                                       >  1:513716/1‑187 (MQ=255)                                                                                                                                                                gAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATCCGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTg                                  >  1:155700/1‑239 (MQ=255)                                                                                                                                                                gAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTg                                  >  1:2088522/1‑239 (MQ=255)                                                                                                                                                                 aGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                 >  1:129924/1‑239 (MQ=255)                                                                                                                                                                 aGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                 >  2:2094974/1‑239 (MQ=255)                                                                                                                                                                 aGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                 >  2:1598250/1‑239 (MQ=255)                                                                                                                                                                         aTTGC‑CGGGTTGGCCGGC‑CGTCACCCCGCCGGGACACA‑CAGAACGGGCCCCGCGGGTGGGGGCACAAATACTTACCCCGGCTCACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGAtttt                          <  2:176576/238‑1 (MQ=25)                                                                                                                                                                         aTTGA‑CGGGAGGGCGGGA‑CGTCACCCCCGCGGGGCCCA‑CCGAACAGCCCACGCGGGGGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGAtttt                          <  2:1152037/238‑1 (MQ=25)                                                                                                                                                                           tGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTAccc                                                                                                                                                                                           >  2:141634/1‑75 (MQ=39)                                                                                                                                                                           tGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTAccc                                                                                                                                                                                           <  1:141634/75‑1 (MQ=39)                                                                                                                                                                           tGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCAt                                                                         >  1:715082/1‑189 (MQ=255)                                                                                                                                                                            gA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACCCGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAgg                      >  1:1780259/1‑239 (MQ=255)                                                                                                                                                                             a‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGa                     >  1:852527/1‑239 (MQ=255)                                                                                                                                                                                 gggtGGCGGGA‑CTTCCCGCCGCCGGGGCCCA‑CCGAACACGCCACGCCGGTGGGGGCACCAATCCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat                   <  2:2088522/235‑1 (MQ=25)                                                                                                                                                                                 ggATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAata                  >  1:1983720/1‑239 (MQ=255)                                                                                                                                                                                   aTGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAatata                >  1:536178/1‑239 (MQ=255)                                                                                                                                                                                   aTGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAatata                >  1:374491/1‑239 (MQ=255)                                                                                                                                                                                   aTGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAatata                >  1:1107453/1‑239 (MQ=255)                                                                                                                                                                                    agggccggggctCACCCGGCGGGGCCACC‑CCGACCCGGCCCCGCGGGGGGGGCCACAAATACTGCCCCCGGCTGACCCCGGAACAATGCAACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat                   <  2:373529/225‑1 (MQ=25)                                                                                                                                                                                         cgggcatCACCCCCCGGGCGCCCA‑CCAAACAGGCCCCGCGGGGGGGGCCACCAATCCTCACCCCGCCTGAACCCCAAACAATTCTCCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAg           <  2:129924/233‑1 (MQ=25)                                                                                                                                                                                           accGTCACCCCGCCGGGCCACC‑CCGCAAAGCCCGCGCCGGGGGGGCCCCCATTACATCCCCCGGCTGAACCTTGACTACTGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa         <  2:91010/237‑1 (MQ=37)                                                                                                                                                                                           a‑CGCCACTCCGCCGGGGAACA‑CCGACCAGGCCCCGCCGGTGGGGCAACAATTACGTCCCCCGTTTGAACCCTGACCAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa         <  2:155700/238‑1 (MQ=37)                                                                                                                                                                                            ccGTCACCCCCCCGGGACACA‑CGCAACGGGCCACGCCGGTGGGGGCACCAATCCGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa         <  1:1754159/237‑1 (MQ=255)                                                                                                                                                                                             cgccACCCCGCCGGGGCACC‑CCGACCGGCCCACGCGGGTGGGGGCACAAATCCGTACCCCGGCTGACCCTTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACGATGAATTGAAGATTTTTAGGAATATATTTAGAACg       <  2:536178/236‑1 (MQ=25)                                                                                                                                                                                                cACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGtt     >  1:701740/1‑238 (MQ=255)                                                                                                                                                                                                  cTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGTTACa  >  1:2049413/1‑239 (MQ=255)                                                                                                                                                                                                   |                                                                                                                                                                                                                                               CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTGCACCGAGCAGGGGATTGA‑CCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACA‑CCGAACAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGTTACA  >  NC_007779/4020906‑4021335

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