Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,582,479	A→C	14.1%	S191A (TCG→GCG)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,582,479	0	A	C	14.1%	179.0 / 3.5	81	S191A (TCG→GCG)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand):  ref base A (52/12);  new base C (2/10);  total (55/26)
Fisher's exact test for biased strand distribution p-value = 3.11e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGA  >  NC_000913/4582251‑4582716                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              ggTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑AACCAGAGTGTTGCCCAGACGGATTGCGCCGCCCTGGGCGGGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGGGCCAGGAGGCGGGTGCCGGGCACCCGTTCGAGGCCGATAAA‑CGCGCGGCGGGGCTGGAAATCC‑TGCGTGTCGCCATCAAGGGCGTCCAGGTCAGTGGTTTG‑GGCCCTAACCTAg                                                                                                                                                                                                                                     >  2:1187335/1‑239 (MQ=255)            aCAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTtc                                                                                                                                                                                                                          >  1:324157/1‑239 (MQ=255)                   gcgcCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCaaa                                                                                                                                                                                                                   >  2:233464/1‑239 (MQ=255)                     gcCTTCGGGA‑GTTTTCCCCGCCGGTGACCCGGGGGTTTGCCCGGGCGGAGTCCACCCCGGGGGCCGAGTTTGCCTTAAATACCGTGCGGCAGGCAGTTCATCAGTGCCAGACGACGGGTCCCGGGCGCCAGTCCGAGGCGGATAAA‑CGCGCGGTGGTTCTGGAAATCC‑TGCGTGTCCCCATCAAGGTCGCCCAGATCATTGGTTTG‑CGCCTTAACATGGCGGTCGGTTTCAATCaaaaa                                                                                                                                                                                                                 <  2:967234/239‑1 (MQ=255)                       cTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCCATCAAAAAg                                                                                                                                                                                                                >  1:438782/1‑238 (MQ=255)                        ttCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCTAGGCAGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGAGTGTAGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGc                                                                                                                                                                                                               >  1:435751/1‑238 (MQ=255)                        ttCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGc                                                                                                                                                                                                               >  2:209439/1‑238 (MQ=255)                          cGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCG‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGATTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGccc                                                                                                                                                                                                             >  2:940228/1‑238 (MQ=255)                          cGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGGCCAGATCATTGGTTTG‑CGACTTCACATAGCGGTCGGCTTCACCCAAAAAGccc                                                                                                                                                                                                             >  2:4963/1‑238 (MQ=255)                          cGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGCTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACCACGGGCGCCGGGCACCAGGTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCACCAAGGTCGTCACGGTCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCGATCAAAAAGccc                                                                                                                                                                                                             >  2:957837/1‑238 (MQ=255)                          cGGCAGGTTTTCACCGTCCCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGACGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGCCGGCTTCAATCACAAAGccc                                                                                                                                                                                                             >  2:297299/1‑238 (MQ=255)                                      gccgccgCT‑CCCCAGAGTGTGGCCCGGACGGTTTGCGCGGCGGTGGGCGAGGTTGCCTTCAAGAGCGGGCAGCAGGCGTTTAAGCAGTGCCAGAGGACGGTGGCGGGGCCCCGGTTCGAGGCCGATAAC‑CGCGGGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGCTTTACCAGGGCGGTCGGCTCCAATAAAAAAGCCCGCCGTACCTgccg                                                                                                                                                                                                <  2:1100383/238‑1 (MQ=255)                                          tCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGCCTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                                                                                                                                                                                            >  1:222693/1‑239 (MQ=255)                                          tCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCgg                                                                                                                                                                                            >  1:207798/1‑239 (MQ=255)                                                aCCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATGAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCAGTACCTGCCGCCGGGt                                                                                                                                                                                          >  1:562834/1‑236 (MQ=255)                                                 cccAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTACCATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGCCCt                                                                                                                                                                                       >  2:1009603/1‑238 (MQ=255)                                                   cAGAGGGTAGCCCAGAGGGATGTCCCCACCGTGGTGGAGGTCGCCTCCAATCTGGTGCAGCAGGCAGTTAATCAGGGCCGGAGGGCGGGTGCCGGGCCCCAGTGGGAGGAATAAAAC‑CGGGGGGTGGTTTTGAAACCCT‑TGCGGTTCGCCTTAAAGGTCGTCCAGTTCTTTGGTTGG‑CGCTTTAACAGCGCGGTGGGCTTCAATAAAAAACCCCGCCGTACCTGCCGCCGGGTCCTGca                                                                                                                                                                                    <  2:223953/239‑1 (MQ=255)                                                      agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGcacc                                                                                                                                                                                  >  1:925305/1‑238 (MQ=255)                                                      agTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGca                                                                                                                                                                                    >  1:611577/1‑236 (MQ=255)                                                           tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGGCTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTc                                                                                                                                                                             >  1:132191/1‑238 (MQ=255)                                                           tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTc                                                                                                                                                                             >  1:150387/1‑238 (MQ=255)                                                           tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGGGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGGGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGATGGATCTGGAAATCC‑TGCGGGTCGCCATCAAGGTCGTCCCGTTCATTGGTTTG‑CGACTTAACATATCGGTCGGCTTCAGTCACAAAGGCCGCCGAACCGGCCGCCGGGCCCTGCACCACTTc                                                                                                                                                                             >  2:1033658/1‑238 (MQ=255)                                                           tGCCCAGACGGATTGCGCCGCCGTGGTCGAGGCTGCCTTCAATATCGTGCAGCAGGCAGATCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCACGGTCGTCCAGATAATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAACCCCGCCGACCCTGCCGACAGGTCCTGCACCACTTc                                                                                                                                                                             >  2:932492/1‑238 (MQ=255)                                                                   gggagcgcgcCGCCGTGTGCAGGGTGGCCTGCAATAGCGTGCGGCAGGCGGTTCATCAGTGCGAGAGGAGGGTTGCCGGGCCCCAGTTAGAGGCCGACAAA‑CGCGGGGTGGTTCTGGAAACCT‑TGCGTTTCCCAAGCAAGGTCGCCCAAATCTTTGGTTTG‑CGCCTAAACATAGCGGCGGGCTAAAACCAAAAGCCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCt                                                                                                                                                                     <  2:1128621/232‑1 (MQ=255)                                                                            ccgccgTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGCTTcag                                                                                                                                                           >  2:672441/1‑239 (MQ=255)                                                                                tgTGGTCGGGGTGGCCTGCAAGAGCGTGCCGCAGGCAGTACATCAGGGCGAGACGAGGGGTGCCGGGCGCCGGTGCGGGGCGGATAAA‑CGCGGGGTGGATCTGAAAATCC‑GGCGTGTCGCAACCAAGGTCATCCAGATCGTTGGTTTG‑CGCTTTAACATGGCGGTCGGCTTCAAACAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTcagcag                                                                                                                                                        <  2:595876/237‑1 (MQ=255)                                                                                cggggTCGGGGGTGCCTTCAATATGGCGCAGCAGGCGGTAAACAAGTGCCAGCGGAGGGGTGCCGGGCACCAGTTCGAGGCGGAAAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAAATCATTGGTTTG‑CGATTTAACATAGCGGGCGGCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTcagcag                                                                                                                                                        <  2:925305/235‑1 (MQ=255)                                                                                cgTGGTCGGGTTTGCCTACAATATCGTGCAGCAGGCAGTTCATCAGGGCAAGCGGCCGGGTGCCGGGCACCAGTTCGAGGCAGATAAA‑CGCGGGGTGGATCGGGAAATCC‑TGCGTGTCGCCATAAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTcagcag                                                                                                                                                        <  1:1046574/238‑1 (MQ=255)                                                                                            tGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGtt                                                                                                                                           >  2:345422/1‑239 (MQ=255)                                                                                              ccTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGACTAATGGtttt                                                                                                                                         >  2:768657/1‑239 (MQ=255)                                                                                                   aaTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCa                                                                                                                                    >  1:360817/1‑239 (MQ=255)                                                                                                   aaTATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCa                                                                                                                                    >  1:446189/1‑239 (MQ=255)                                                                                                          tGCAGCAGGCAGTTCCTCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAACCCGCCGTACCTGCCGCCGGGTCCTGCCCCACTTCACGCGGCGGCGGCTTCAGCAGATGAGTAATGGTTTTAATCAGCggca                                                                                                                              >  1:614125/1‑236 (MQ=255)                                                                                                          tGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGTGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCTCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTAAGCAGATGAATAATGGGTTTAATCAGCGGAc                                                                                                                              >  2:145302/1‑238 (MQ=255)                                                                                                          tGCAGCAGGCAGTTCATCAGTGCCAGAAGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGTTCGGCTTCAATCAAAACGCCCCCCGTCCCTGCCGCCGGGTCCTGCACAACTACACGCGGCTGCGGTTTCAGCAGATGACTCCTGGCTTTAATCAGCGGAc                                                                                                                              >  2:1014387/1‑238 (MQ=255)                                                                                                          cacagcagGCCGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TCCGTGTCGCCATCAGGGGCGTCCAGATCATTGGTTTG‑CGACTTAACTTAGCGCTCGTCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCCGCGGATTCAGAAGTTGAACAATGGTTTTAATGATCGGAc                                                                                                                              >  2:716042/3‑238 (MQ=255)                                                                                                           gcaggagagagTGCATCGGGGCCAGAGGACGGGGGACGGGCGCGAGGCCGA‑GCCGAAAAAGCGCGCGTTGGTTCTGAAACTCC‑TGGTTGTCGCCATCAGGGTCGTCAAAATCATGGGTTGG‑CCACTAAACATGGCGGTCGGCTACAACAAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGAcg                                                                                                                             <  2:484379/229‑1 (MQ=255)                                                                                                                  gcgtttCAGCAGGGCCAGACGGGGGGGGCGGGACCCCGGTTGGAGGCCGATAAC‑CGGGGGGTGGATCGGGAAACCC‑TGCGTGCCCCCATAAAGGCCGTCCATATCATTGGTTTG‑CGACTTAACATAGCGGTCGGTTTAAACAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTg                                                                                                                      <  1:1057285/234‑1 (MQ=255)                                                                                                                  gcagTTCATCAGTGCCAGACGACGTGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATGC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGGTTG‑CGACTTAACATAGCGGTCGGCTTCAGTCAAAAAGCCGGCCGTACCTGCCGCCGGGCCATGCCCCACTTCACGCCGCCGCGGTTTCGGCAGATGAATAATGGTTTTAATGAGCGGACTCGGGGTgt                                                                                                                     >  2:870489/1‑238 (MQ=255)                                                                                                                      ttCGTCAGTGCCAGGCGGCGGGGGCCGGGCCCCCGCGCGAGGCCGATAAC‑CGGGGGTTGGAT‑TGGAAACCCGCGCGTTTCGCCACCAGGGTCTACCAGATCGTTGGTTGG‑GGTCCAAACATAGCGGTCGTCTTCAACCAAAAACCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTa                                                                                                                 <  2:72715/239‑1 (MQ=255)                                                                                                                                ccAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGc                                                                                                        >  1:228893/1‑238 (MQ=255)                                                                                                                                ccAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGAGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGATTCAGCAGATGCATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGc                                                                                                        >  1:789800/1‑238 (MQ=255)                                                                                                                                      gaggggTGCCGGGCCCCAGTTCGAGGCGTATACA‑GCGGCGGTGGGTCT‑GAAATCC‑TGGGTGGCGCCATCAAGGCCGCCAAGATCATGGGTTTG‑CGTCTTAACATGGCGGTCGGCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGAc                                                                                                 <  2:1156683/235‑1 (MQ=255)                                                                                                                                        cgGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACtt                                                                                               >  1:523881/1‑239 (MQ=255)                                                                                                                                        cgGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGAACTGGCCTGCACCCGACtt                                                                                               >  2:255645/1‑239 (MQ=255)                                                                                                                                        cgGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGACGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAGAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGGGAAGTACTGGCCTGCACCAGACtt                                                                                               >  2:622679/1‑239 (MQ=255)                                                                                                                                              ccGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTc                                                                                         >  1:1129178/1‑239 (MQ=255)                                                                                                                                                     ccAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCTGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTcgcg                                                                                  <  1:304458/239‑1 (MQ=255)                                                                                                                                                     ccAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTcgcg                                                                                  >  1:76868/1‑239 (MQ=255)                                                                                                                                                                  gATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACa                                                                     <  1:209439/239‑1 (MQ=255)                                                                                                                                                                   aTAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACa                                                                     >  1:1037038/1‑238 (MQ=255)                                                                                                                                                                   aTAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACa                                                                     >  2:708899/1‑238 (MQ=255)                                                                                                                                                                     taacggcggggggtcttgggAAACCG‑TGCGTGCCCCCAGCAAGGTCCCCAAGACCATGGGTTTG‑CGATTTAACAGGGGGGTCGTCTCAAACCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGc                                                                   <  2:573113/221‑1 (MQ=255)                                                                                                                                                                         tttaaagaaagcaggaaattgc‑TGCGGGCCGCGATCAGGGTCATCCACACTGTTGTTTTG‑CGCTGTAACAGATCGGTGGTCTCCGACAAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCt                                                                <  2:507720/217‑1 (MQ=255)                                                                                                                                                                                tGGATTTGAAAATCC‑GGCGTGTCCCCAACACGGCCGTACAGATTTTTGGTTCG‑CGACTTAACATGGGGGCCGTCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACa                                                         <  1:1187335/238‑1 (MQ=255)                                                                                                                                                                                          aaTCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATACTGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAg                                              >  2:63762/1‑239 (MQ=255)                                                                                                                                                                                             cc‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGCTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATCTCGCCGAAGTCa                                           >  2:874486/1‑239 (MQ=255)                                                                                                                                                                                              c‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCa                                           >  1:718341/1‑238 (MQ=255)                                                                                                                                                                                                   gtgtCGCCATCAAGGTTCTCCAGAACTTTGGTTTG‑CGTCTTAACATAGCGGTCGGCTACAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATcgc                                       <  1:552812/238‑1 (MQ=255)                                                                                                                                                                                                   gtgtCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACatat                                                      >  1:507720/1‑223 (MQ=255)                                                                                                                                                                                                        gCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACTCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGCCt                                 >  1:642095/1‑239 (MQ=255)                                                                                                                                                                                                                ggTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgc                         >  1:461502/1‑239 (MQ=255)                                                                                                                                                                                                                 agagcaaggATCATGGGTTTG‑CGACTAAACATGGGGGTCGCTTTAAACAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  2:835038/231‑1 (MQ=255)                                                                                                                                                                                                                  tggtcCAGATCTTTGGTTTG‑GGCTTAAACCTAGGGGTGGGCTACAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:1163259/236‑1 (MQ=255)                                                                                                                                                                                                                  tctcgcagcattttGGTTTG‑CGTCTTAATATGGCGGTCGGCTTCAATAAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:310327/226‑1 (MQ=255)                                                                                                                                                                                                                  tcgccCATATCATTGGTTGG‑CACTTTAACATAGCGGTCGGCTACAACCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:811559/238‑1 (MQ=255)                                                                                                                                                                                                                  tcgccCAGTTCA‑TGGTTTGCCGCCTTAACATGGCGTCCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcg                        <  1:833559/238‑1 (MQ=255)                                                                                                                                                                                                                   cctcCATATCTTTGGTTCG‑CGATTTAACATAGGGGTCGTCTTAAATCAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTgcgc                       <  1:230073/236‑1 (MQ=255)                                                                                                                                                                                                                      ccAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATg                                                                                                                                              >  1:1039108/1‑116 (MQ=255)                                                                                                                                                                                                                           tCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTCCATATCGCCGAAGTCCTCGCCCGACTTACCGTGCGCGCCGTTga               >  2:310327/1‑237 (MQ=255)                                                                                                                                                                                                                           tCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGt               >  1:493208/1‑238 (MQ=255)                                                                                                                                                                                                                            ctttGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAAACAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTAc             <  1:4963/237‑1 (MQ=255)                                                                                                                                                                                                                                gttttG‑CGCCTAAACGGAGCGGGGGGCTTAAACAAAAAGGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGt         <  1:957837/237‑1 (MQ=255)                                                                                                                                                                                                                                ggTTTG‑CGACTTAACATGGCGGTCTGCTTCAAACAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCGTTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGt         <  1:940228/239‑1 (MQ=255)                                                                                                                                                                                                                                    tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTa              >  1:835038/1‑230 (MQ=255)                                                                                                                                                                                                                                    tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  1:941535/1‑238 (MQ=255)                                                                                                                                                                                                                                    tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  1:774718/1‑238 (MQ=255)                                                                                                                                                                                                                                    tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  1:762383/1‑238 (MQ=255)                                                                                                                                                                                                                                    tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCa      >  2:949551/1‑238 (MQ=255)                                                                                                                                                                                                                                    tG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGCGCGCGCCGTTGTACCAGTCCa      >  2:14154/1‑238 (MQ=255)                                                                                                                                                                                                                                       cGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCg   >  1:464472/1‑239 (MQ=255)                                                                                                                                                                                                                                        caCTTACCAGTGCGGTCGGCTCAAAACAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGa  <  1:815903/238‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                              GGTTAGTGGCGACAATATGCGCCTTCGGCAGGTTTTCACCGTCGCT‑ACCCAGAGTGTTGCCCAGACGGATTGCGCCGCCGTGGTCGAGGTTGCCTTCAATATCGTGCAGCAGGCAGTTCATCAGTGCCAGACGACGGGTGCCGGGCACCAGTTCGAGGCCGATAAA‑CGCGCGGTGGATCTGGAAATCC‑TGCGTGTCGCCATCAAGGTCGTCCAGATCATTGGTTTG‑CGACTTAACATAGCGGTCGGCTTCAATCAAAAAGCCCGCCGTACCTGCCGCCGGGTCCTGCACCACTTCACGCGGCTGCGGTTTCAGCAGATGAATAATGGTTTTAATCAGCGGACGCGGGGTGAAGTACTGGCCTGCACCAGACTTGGTTTCATTCGCGTTCTTCTGCAACAGCCCTTCGTACATATCGCCGAAGTCATCGCGCGACTTACCGTGCGCGCCGTTGTACCAGTCCAGCGA  >  NC_000913/4582251‑4582716

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