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%	222.9 / 4.4	98	G42G (GGA→GGG)	Y75_RS19860	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base T (57/24);  new base C (3/14);  total (60/38)
Fisher's exact test for biased strand distribution p-value = 7.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGTTAC  >  NC_007779/4020906‑4021334                                                                                                                                                                                                  |                                                                                                                                                                                                                                              cAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCCCTCCGCCGGGGCACACCGAA‑CAGGCCCCGCCGGTGGGGGCACCAATAc                                                                                                                                                                                                >  2:922048/1‑239 (MQ=11)   gCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGt                                                                                                                                                                                              >  1:909808/1‑239 (MQ=31)   gCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCAAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGt                                                                                                                                                                                              >  1:1821008/1‑239 (MQ=14)     tGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACg                                                                                                                                                                                               <  1:1675648/236‑1 (MQ=31)     tGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTAc                                                                                                                                                                                            <  1:23044/239‑1 (MQ=32)            gggtgtCCGAAGGCATGTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGGCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCt                                                                                                                                                                                     <  2:1530256/234‑1 (MQ=12)            ggATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCt                                                                                                                                                                                     <  1:1253738/239‑1 (MQ=32)            ggATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCt                                                                                                                                                                                     <  1:1125944/239‑1 (MQ=32)            ggATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCt                                                                                                                                                                                     <  1:1612988/239‑1 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCCCACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGa                                                                                                                                                                                   >  2:1021777/1‑239 (MQ=17)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGGTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGa                                                                                                                                                                                   >  1:1526622/1‑239 (MQ=17)                atCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAc                                                                                                                                                                                 <  1:352434/239‑1 (MQ=32)                atCCGCAGGCATTTTCCAGCCGCGGCCGAGGCACCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGACAGCGGGCCGGGCAGGGCGATGACGGTTTCACCGGGAAGGCCCTATG‑CACCGAGCAGGGGATCGACCGGATGGCCGGACGTCACTCAGCGTGGGGCCACCGAA‑CAGGCCCCTCCGGCGGGGGCCCCAATACGTGCCCCGGCCGaa                                                                                                                                                                                  >  2:84458/1‑238 (MQ=11)                 tCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAcc                                                                                                                                                                                <  1:135674/239‑1 (MQ=32)                 tCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAc                                                                                                                                                                                 >  1:919828/1‑238 (MQ=32)                  ccGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAAccc                                                                                                                                                                               <  1:3564/239‑1 (MQ=32)                         cATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGaa                                                                                                                                                                           <  1:1518657/236‑1 (MQ=32)                         cATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACa                                                                                                                                                                         <  1:758146/238‑1 (MQ=32)                         agtagtCCCGCCGGGCCCGCGCCGCCCCACCGGCGCGGGCGTTTTGGTCCGGTAACTGCGGTAGGTGCGGGAGAGGATGACCGGCCGCGGGCCGGTCAGGGCGATGTCGGTTTCACCGGGAAGGGCCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACa                                                                                                                                                                         <  2:1066915/233‑1 (MQ=11)                                      gggCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTCCCCCGGCCGAACCCTGCACAATGGTACCGCCAt                                                                                                                                                            >  2:1405105/1‑238 (MQ=12)                                            ggtgcccccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CAGCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCg                                                                                                                                                      >  1:1414086/6‑238 (MQ=17)                                            aGGCCCCCCACCGGCGCGGGCGTTTTTGTCCGGTACCGGCTGTAGGTGCGGGAGAGGATGAACGGCGGGGGGCCGGGCAGGGCGATGTCGGTTTCCCCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCCCTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCg                                                                                                                                                      <  1:552806/238‑1 (MQ=255)                                                cccccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATg                                                                                                                                                 >  1:339328/2‑239 (MQ=17)                                                        ggCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTg                                                                                                                                         >  1:926825/1‑239 (MQ=33)                                                         gcgcGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCGGAACAATGCTCCCGCCATACTGCGTCCTGTCCCCCTg                                                                                                                                         >  2:1146859/1‑238 (MQ=12)                                                          cgcgGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGAc                                                                                                                                       >  1:929209/1‑239 (MQ=33)                                                            cgGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGAcgc                                                                                                                                     >  1:851665/1‑239 (MQ=33)                                                              ggCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCAGGTCACCCTGAcgcc                                                                                                                                    >  2:952584/1‑237 (MQ=18)                                                                                  gctgTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCg                                                                                                               >  2:1284484/1‑239 (MQ=37)                                                                                    tgTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTa                                                                                                             <  1:965497/239‑1 (MQ=37)                                                                                        ggTCCGGGAGAGAATGAACGCCAGCGGGCCGGGCAGGGCGATGTCGGTTTCGCCGCGGCGGACCTTTG‑CACCGAGCACGGGATTGACCGGACGGCCGGACGTCCCCCCGCCGGGGCACACCGAA‑CAGGCCACGCCCGCTCGGGCACCCACACGTAACCCGGCTGAAACCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCACCGACGGTTTTCCGCCCCTCGTaaca                                                                                                          >  2:1002237/1‑236 (MQ=11)                                                                                           gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCCTACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCt                                                                                                      >  2:1184388/1‑239 (MQ=255)                                                                                           gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCt                                                                                                      >  1:167027/1‑239 (MQ=40)                                                                                            cGGGAGAGGATGCACGGCAGCGGGCGGGGCAGGGCGATGGAGGTTTGACCGGGAAGAACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGTCGGACGTCAGCGCGCCGGGGCACACCGAA‑CAGGCCCGGCCGGTGGGCGCACGAATCCGTACCCCGGCTGAACCACGAACCACGCCACGGCCCTAACGCGTCATGTCACACTGACGCGCCGCCGGTTTTCCGACCATCGTACACTCCtt                                                                                                     >  2:687343/1‑239 (MQ=11)                                                                                             gggAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCttt                                                                                                    >  1:1844583/1‑239 (MQ=255)                                                                                                    gATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCtatta                                                                                             >  1:15139/1‑239 (MQ=255)                                                                                                    gATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATAATGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCtatta                                                                                             >  1:1440106/1‑239 (MQ=255)                                                                                                     aTGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCtattat                                                                                            >  2:1753333/1‑239 (MQ=255)                                                                                                      tGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTT‑CACCGAGCAGGGGGTTGACCGGATGGCCGGCCGGCCCTCCCCCGGGGCACACCGGA‑AAGGCCACGCCGGTGGGGGCCCCAATACGTACCCCGGCTGAACCCCGAACAATGCTACCGCCATACCGCGCCATGTCACCCGTACGCGCCGCCGGTTTTTCCCTCATCGTAAAATCCTTTTCCattatt                                                                                           >  2:938441/1‑239 (MQ=14)                                                                                                      tGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTattatt                                                                                           >  1:914364/1‑239 (MQ=255)                                                                                                         aCGGCGGCGGGCGGGGCAGGGCGATGTCGGTTTCCCCGGGAAGGCCCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCCCTCCGCGGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                        <  1:1894325/239‑1 (MQ=255)                                                                                                         aCGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGGTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                        >  2:550847/1‑239 (MQ=255)                                                                                                          cGGCTGGGGGCGGGGCAGGGGGATG‑CGGGTTCCACGGGAAGGCCCCTCGCCCCCGCGCGGGGGATTTGCGTGATGGCGGGACTCCCCCCCGCGGGGCCACACCAAA‑CAGGCCCCGCCGGTGGGGACCCCATTGCTTCCCCCGGCTGACCCCTGAACAATTCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                        <  2:1742944/238‑1 (MQ=9)                                                                                                           ggCAGGGGGCGGGGCAGGGCGATGTCGGTTTCACCGGGAAGGCCCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACCCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAAt                                                                                      <  1:1267807/239‑1 (MQ=255)                                                                                                           ggCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGCCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAAt                                                                                      <  1:1021777/239‑1 (MQ=255)                                                                                                               gCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACACTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACa                                                                                  >  1:1038477/1‑239 (MQ=255)                                                                                                                ggggcggggcggggcgAGGCGGGTTTCGCCGGGAAGGCCCTTTG‑CCCCGGGAGGGGGATTGACCGGAGGGCGGGACG‑CCCCCCCCCGGGGCCCCCCCAACCAGCCCACGCGGGTGGGGGCACCAATACTTCCCCCGGCTCAACCCTGAACATTGCTGCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACa                                                                                  <  2:914364/237‑1 (MQ=17)                                                                                                                       ggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTc                                                                          >  1:687083/1‑239 (MQ=255)                                                                                                                            ggCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAACTCCTTTTCTATTATTTTAATGACAGTACCCTCATACg                                                                     >  2:1768206/1‑239 (MQ=255)                                                                                                                             gCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGa                                                                    >  2:1027487/1‑239 (MQ=255)                                                                                                                                 tGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGa                                                                                                                                                                            <  2:1455365/131‑1 (MQ=34)                                                                                                                                 tGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGa                                                                                                                                                                            >  1:1455365/1‑131 (MQ=34)                                                                                                                                 tGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCt                                                                            >  1:179019/1‑227 (MQ=255)                                                                                                                                 tGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGaaaaa                                                                >  1:190081/1‑239 (MQ=255)                                                                                                                                      gTTTCACCGGGCTGGCACTTTG‑AACCGAGCAGGGGATTGACCGGAATGACGCACGTCACTCAGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCACTACGTACCCCGGCTGAACACTGAACCATTCTACCGCCATACTGGGTCCTGTCACCCTGACGCGCCGCCCGGTATCCTCTCCTCGTAACATCCTTTTCTAGTATTTTAATGAAAGGATCTTCTTACGAAAAAGCgaa                                                           >  2:1129143/1‑239 (MQ=17)                                                                                                                                      gTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCgaa                                                           >  2:225383/1‑239 (MQ=255)                                                                                                                                        ttCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTTTTCTTTTAATGACATTATCTTCATACGAAAAAGCgaaaa                                                         >  2:52277/1‑237 (MQ=255)                                                                                                                                                gAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTtg                                                 >  2:1341862/1‑239 (MQ=255)                                                                                                                                                 aaGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTtgt                                                >  1:377530/1‑239 (MQ=255)                                                                                                                                                    gccttttG‑CCCCGAGAAGGGGATTGCCGGGATGGCCGGACGTCACCCCGCGGGGGCACCCCAAA‑CAGGCCACGCCGGTGGGGGCCCCAATACGTCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTgg                                              <  1:952584/234‑1 (MQ=21)                                                                                                                                                    gACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTgg                                              >  1:1509906/1‑238 (MQ=255)                                                                                                                                                     aCCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTc                                            >  2:1002010/1‑239 (MQ=255)                                                                                                                                                     aCCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTc                                            >  2:1321791/1‑239 (MQ=255)                                                                                                                                                      cctttgaccgcGAGCGGGGAGTGGACCGGGGGCCCGGGCGCCCCCCCCCGGGGACACCCCAAA‑CAGGCCCCGCCGGTGGGGCCACAAATACGTCCCCCGCCTGACCCTTAAACAGTGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGt                                             <  2:906093/228‑1 (MQ=21)                                                                                                                                                        tttG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTcaca                                         <  1:1284484/239‑1 (MQ=255)                                                                                                                                                         ttggccccGAGCAGGGGATTGACCGGAGGGCCGACCTTCCCCCCGCGGGGGCACCCCGAA‑CAGGCCACGCGGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTcaca                                         <  1:922048/233‑1 (MQ=21)                                                                                                                                                         gtcgcccggggCAGGGGTTGGCCCGGCGGGCGGGACGTCCCCCCGCGGGGGCCCCCCGAA‑CGGGCCACGCCGGTGGGGGCACCAATCCTTCCCCCGGCTGACCCCTGAACAATGTTCCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTcaca                                         <  2:1553261/229‑1 (MQ=21)                                                                                                                                                          tG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAAt                                       >  1:1553261/1‑239 (MQ=255)                                                                                                                                                          tG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAAt                                       >  1:1224053/1‑239 (MQ=255)                                                                                                                                                                cGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAAtt                                  >  2:1566447/1‑239 (MQ=255)                                                                                                                                                                  aGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                >  1:1112235/1‑239 (MQ=255)                                                                                                                                                                  aGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                >  1:567965/1‑239 (MQ=255)                                                                                                                                                                  aGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                >  2:622188/1‑239 (MQ=255)                                                                                                                                                                  aGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                                >  2:728307/1‑239 (MQ=255)                                                                                                                                                                  aGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGCTTTCCGCTCATCGTAAAATCCTGTTCTATTATTTTACTGACAGTATCTTCATACGAAGAAGCGAAGAAAAACTTGTGGTCACAAGGAATTGa                                >  2:507030/1‑239 (MQ=255)                                                                                                                                                                          aTTGCCCGGAGGGCCGGACGCCCCCCCGCGGGGGCACCCCAAA‑CAGGCCACGCCGGTGGGGGCACAAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGAtttt                         <  1:1004912/238‑1 (MQ=25)                                                                                                                                                                           ttGCCCGGAGGGCCGGACGTCACCCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTCCCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTa                       <  1:1693806/239‑1 (MQ=255)                                                                                                                                                                           ttGACCGGATGGCCGGACGTCACTCCGCCGGGGCACCCCGAA‑CAGCCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGAttttt                        <  1:938441/238‑1 (MQ=255)                                                                                                                                                                            tgacaagttGGCCGGAAGTCTCTCCGCCGACGCACACAGGA‑CAGGCCACGCCGGACGTAGCGCCGCTACGTACCCCTGCTGCGCCCAGTACAATTCTTCAGCCATACTGCGTCATGTCACCCACACCCGACGCCGCTGTTCAGCGCATCGTACAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAACGCGAAGACAACTTTGTGGTTACACTGAATTGAAGATTTTTAg                      >  2:1921208/9‑239 (MQ=255)                                                                                                                                                                            tGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAg                      >  1:666555/1‑239 (MQ=255)                                                                                                                                                                            tGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAg                      >  1:40565/1‑239 (MQ=255)                                                                                                                                                                            tGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAg                      >  2:1530513/1‑239 (MQ=255)                                                                                                                                                                             gACCGGATGGCCGGACGTCCCTCCGCCGGGGCACACCACA‑CAGGCCACGCCGGTGGGGGCACCAATACTTCCCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAgg                     <  1:831243/239‑1 (MQ=255)                                                                                                                                                                               cgggAGGGCCGGACGTCCCTCCGCCGGGGCACACCCAA‑CAGTCCCCGCCGGTGGGGGCACCATTACGTCCCCCGGCTAACCCCTGAACAGTGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGaa                   <  1:348493/237‑1 (MQ=25)                                                                                                                                                                                cGGGTGGCCGGCCGCCACTCCGCGGGGGCCCCCCGAC‑CAGCCCACGCCGGGGGGGGCACCAATCCGTACCCCGGCTGACCCCGGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat                  <  2:94653/239‑1 (MQ=25)                                                                                                                                                                                cGGGTGCCCGCACGTCACCCCCCGGGGACACCCCGAA‑CAGGCCACGCCGGTGGGGGAACCATTCCGTCCCCCGGCGGAACCCTAAACAATGTTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat                  <  2:1006436/239‑1 (MQ=25)                                                                                                                                                                                cGGATGGCCGGACGTCCCCCCGCGGGGGAACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTCCCCCGGCTGACCCCTGAACAATGTTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat                  <  1:289964/239‑1 (MQ=255)                                                                                                                                                                                cGGAGGGCCGGCCGTCACTCCGCCGGGGCACCCCGAA‑CAGCCCACGCGGGTGGGGGCACCAATACGTCCCCCGGCTGAACCCTAAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat                  <  1:1504716/239‑1 (MQ=255)                                                                                                                                                                                   atcgCCGGACGTCACTCTGCCGGTGGACACCAAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAAGGTTACCGCTGTCCTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAACACTTTTCTATTATTTTCATGACAGTATCTTCATACGAAAAAGCACAGAAAAACTTGCGGTCACAATGAATTGAAGATTTTTAGGAatata               >  2:1274052/4‑239 (MQ=21)                                                                                                                                                                                   aTGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTaaaa                                                                                                          >  2:258783/1‑148 (MQ=255)                                                                                                                                                                                   aTGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTaaaa                                                                                                          <  1:258783/148‑1 (MQ=255)                                                                                                                                                                                   aTGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTCTTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTTTGTTCAACATGAATTGAAGATTTTTAGGAatata               >  2:70611/1‑239 (MQ=255)                                                                                                                                                                                           ccgccACCCCGCCGGGGCCACCCGAA‑GTGGCCCCGCGGGGGGGGCCACCAATACGTACCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa        <  2:1414086/234‑1 (MQ=37)                                                                                                                                                                                           ccGTCCCCCCGCCGGGCCACCCCGAC‑AGGGCCACGCCGGGGGGGGCACAAATACTTACCCCGGCTGAACCCTGAACAATGTTTCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa        <  2:279566/237‑1 (MQ=37)                                                                                                                                                                                           aCGTCACCCCGCCGGGGCCCACCGAA‑CGGGCCACGCGGGTGGGGGCCCCAATCCGTACCCCGCCTGAACCCTGAACAATGCTCCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa        <  2:1112235/238‑1 (MQ=25)                                                                                                                                                                                            cGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACg      >  1:1330980/1‑239 (MQ=255)                                                                                                                                                                                            cGTCACCCCGCCGGGCCACACCGCA‑CGGGCGACCTCGGTGGGGGCACCAATCCGGACCCCGGCTGCACCATGAACAATACTACCGCCAAACAGCGTCATGTTACCCTGACGCGCCAACGGCTTTCCGCTCACTGTAAAATCCTTTTCTATTATGTTAATGACAGTATCTTTATACGAACAAACGAAGAAAAACTTGTGATCACAAGGAATTGAAGATTTTTAGGAATATATTTAGaatg      >  2:1708311/1‑237 (MQ=255)                                                                                                                                                                                                aCTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGTTAc  >  2:1348962/1‑239 (MQ=255)                                                                                                                                                                                                  |                                                                                                                                                                                                                                              CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTTG‑CACCGAGCAGGGGATTGACCGGATGGCCGGACGTCACTCCGCCGGGGCACACCGAA‑CAGGCCACGCCGGTGGGGGCACCAATACGTACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAACGTTAC  >  NC_007779/4020906‑4021334

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