Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,336,601	A→C	12.3%	G14G (GGT→GGG)	yfaA ←	DUF2138 family protein, putative host defense protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,336,601	0	A	C	12.3%	188.6 / 2.4	86	G14G (GGT→GGG)	yfaA	DUF2138 family protein, putative host defense protein
Reads supporting (aligned to +/- strand):  ref base A (40/32);  new base C (2/10);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGG  >  NC_000913/2336363‑2336829                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     tCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGAGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGAGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCGAGCGGAAAAC‑AGTGGTATTGAGCCCAAACCTA‑CAAGa                                                                                                                                                                                                                                      >  1:51507/1‑239 (MQ=255)   tGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTCCGCCAAAACCTA‑CAAGAcc                                                                                                                                                                                                                                    >  1:218352/1‑239 (MQ=255)     ataataGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCAt                                                                                                                                                                                                                                  >  2:932944/1‑239 (MQ=255)     ataataGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACATA‑CAAGACCAt                                                                                                                                                                                                                                  >  2:523796/1‑239 (MQ=255)         taGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACAACCAGACGCGACAAAGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGa                                                                                                                                                                                                                               >  1:527779/1‑238 (MQ=255)          aGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGaa                                                                                                                                                                                                                              >  2:948920/1‑238 (MQ=255)          aGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGaa                                                                                                                                                                                                                              >  2:807149/1‑238 (MQ=255)          aGAAAACGAAATCTTCGCTCAGAACATAGTGGAGAAACGGACTGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGGCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCAGCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑CGTGCTATTGCACCAAAACATA‑CACGCCCATAGaaa                                                                                                                                                                                                                             >  1:836076/1‑238 (MQ=255)                   aaTCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCA‑CACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGG‑AAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCtt                                                                                                                                                                                                                   >  2:773586/1‑238 (MQ=255)                         cGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCg                                                                                                                                                                                                                >  1:1027267/1‑237 (MQ=255)                               gAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCtttttt                                                                                                                                                                                                        >  2:730326/1‑239 (MQ=255)                                   aTCGTGGAGAACCGGAATGGTGAGAGGTTCTTTGGGAAGCTGTCTCAGAGAGTCGCTGTCAAGAAAGAGTCCCGCCTCATGGAGAT‑CGATTTGCAGATTGTGGTGCACCACCAGCGGCGACAACGTTTTTTATGACCCACTGCAAGCATATTGCAACCCCCAGCCCCCGGCGAAAAGC‑AGTGCTATTGCGCAAAACCCAA‑CAAGTCCATAGAACCGCCAGCCTTTCGCCTTTTTTTcac                                                                                                                                                                                                    <  1:730326/239‑1 (MQ=255)                                     cGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTcacc                                                                                                                                                                                                   >  2:987546/1‑238 (MQ=255)                                      gTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTcacca                                                                                                                                                                                                  >  1:662429/1‑238 (MQ=255)                                      gTGGAGAAACGGAATGGTGAGGAGATCTTGGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTGGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCAAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTcacc                                                                                                                                                                                                   <  2:6118/237‑1 (MQ=255)                                      gTGGAGAAACGGAATGGTGAGGAGACCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTcacca                                                                                                                                                                                                  >  2:176731/1‑238 (MQ=255)                                          agaAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCAACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATCGAACGCCCAGCCTTTCGCCTTTTTTTCACCACTCAt                                                                                                                                                                                             >  1:155268/1‑239 (MQ=255)                                          aaaaaCGGAATGGTGAGGAGATTTTTGGGAAGCTGGCTCGGAGAGTCGCTGTCGAGAAAGAGTTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTGGTGCACCACCAGCGGCGAAAACGTTTTTTCTGGCCCCCTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCa                                                                                                                                                                                              <  1:638067/236‑1 (MQ=255)                                             aaCGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTg                                                                                                                                                                                           >  1:35143/1‑238 (MQ=255)                                                   agggtggggatacgtgtggtatgCTCGCACAGAGAGACGCTGTGGATAAAGAGGGCCGGCAAATGGAGAT‑GGATTTGCAGAGTGTTGCCCATCCCCAGCGCCGACAACGTTGTTTCTGCCCTCTTGCCAGCATCTTGCAACGCCCGTACCCCA‑CGGGAAGCAAGTGTGATGGCGCAAAACCCGA‑AAAGCCCATGGACCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACAt                                                                                                                                                                                     <  1:444064/216‑1 (MQ=255)                                                     ggTGAGGAGATCTTGGGGAAGCGGGCTCAGAGAGCCGCTGTCGAAAAAGAGGTCCGGCTCATTGAGAT‑CGATTGGCAGATTGTTGTGCACCCCCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACCCCAAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCGGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTcc                                                                                                                                                                                  <  2:722035/239‑1 (MQ=255)                                                      gTGAGGGGATCTTTGGGAAGCTGGGCAGGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATGGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCACCGCCCAGACGCAAGCGGAAAGC‑AGTGCTTTTGCGCCAAACCCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTcc                                                                                                                                                                                  <  1:317681/238‑1 (MQ=255)                                                        gaggagATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCCAAACCTA‑CAAGACCATAGAACCGCCCGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtg                                                                                                                                                                               >  2:325931/1‑239 (MQ=255)                                                         aggagATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTtgt                                                                                                                                                                              >  2:52547/1‑239 (MQ=255)                                                            agaTCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGta                                                                                                                                                                           >  1:347761/1‑239 (MQ=255)                                                                     ggAAGCTGGCTCAGAGAGTGGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCACCGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATagccagc                                                                                                                                                                   <  2:851402/238‑1 (MQ=255)                                                                          cTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGCTTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCTAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCCTTGCCACATTCTTTGTGTATAGCCAGCCAtttt                                                                                                                                                             >  2:705967/1‑239 (MQ=255)                                                                          cTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCGGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt                                                                                                                                                             >  2:615511/1‑239 (MQ=255)                                                                          cTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttt                                                                                                                                                             >  2:991787/1‑239 (MQ=255)                                                                             gCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCAtttttt                                                                                                                                                           >  2:43215/1‑238 (MQ=255)                                                                              cTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTa                                                                                                                                                          <  2:347761/238‑1 (MQ=255)                                                                                   agagtcgtgtgggaaaaaGGGGGCCGACTCATTAAGAG‑CTATTTGCAGATTGTTGGGCACCACCAGCGGCGAAAACGTTTTTTCTGGCCCACCCCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTTCTATTGCGCCAAACCCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGc                                                                                                                                                     <  1:352790/223‑1 (MQ=255)                                                                                       tCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTAGGGGCACAg                                                                                                                                                 >  2:915215/1‑238 (MQ=255)                                                                                       tCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGAGTATAGCCAGCCATTTTTTACGGGCACAg                                                                                                                                                 >  2:186139/1‑238 (MQ=255)                                                                                                  aagagggtcaggtgcATGGAGGT‑GGATTGGGAGATTTTTGTGCACCCCAGGCGGCGACAACGTTTTTCCGGGCCCACTGCCAGCATATTCCACCGCCCGGCCGCCAGCGGAAAGC‑AGTGCTTTTGCGCCAAAACCAA‑CAAGACCAGAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTAcc                                                                                                                                     <  1:888594/225‑1 (MQ=255)                                                                                                  aaaGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTTGTGTGCACCACCGGCGTCAACAACGTTTTTTCTGGCCCACTGCCATCATATTGCAACGCCTAGACGCCGGCGGAAAGC‑AGTGCTTTTGCCCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTAcc                                                                                                                                     <  1:1037243/239‑1 (MQ=255)                                                                                                   aaGAGGTCCGGCTCTGTGAGTT‑CGATTGGCAGTTTGTTGTGCCCCACCAGCGGAGACAACTTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCg                                                                                                                                    <  1:1175766/239‑1 (MQ=255)                                                                                                   aaGAGGTCCGGCTCATTGAGAT‑CGTTTTGCAAATTGTTGTCACCCACAGGCGGCTAAACCTTTTTTTCGGGCCCGCTGCCAGCATTTTGCAACGCCCAGCCGCCAGCGAAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCg                                                                                                                                    <  1:932944/239‑1 (MQ=255)                                                                                                   aaGAGGGCAGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCCCCACCAGCGGCGACAACGTTTTTTCTGGCCCTCTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCAAAACCCTA‑CAAGTCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTAc                                                                                                                                      <  2:773353/237‑1 (MQ=255)                                                                                                   aaGAGGCCCGGTTAATTG‑GATGCAATTTGCGGATAGTTGCGGACCACCAGCGGCAAAAACATTTTGTTTGGCCGCCGCACCGCATGCGGCAGCGACAAGCCGCCGGCGGAAAGT‑ATTTCTATTGGAATAAACCCAA‑CAAACCAATAGATCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCg                                                                                                                                    <  1:52547/239‑1 (MQ=255)                                                                                                   aaGAGCTCCGGCTCATTGAGTT‑CGTTTTGCAGATTGTTGTGCA‑CACCAGCGGCGAAAACGTTTTTTCTGCCCCCCTGCCAGCATATTGAAACGCCCAGACGCCAGCGG‑AAGC‑AGTGCTATTGCGCAAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCg                                                                                                                                    <  1:773586/237‑1 (MQ=255)                                                                                                       ggTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGc                                                                                                                                 >  1:546341/1‑238 (MQ=255)                                                                                                       ggTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGc                                                                                                                                 >  1:211449/1‑238 (MQ=255)                                                                                                         tCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTTCCCt                                                                                                                              >  2:215852/1‑239 (MQ=255)                                                                                                         tCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCt                                                                                                                              >  1:734712/1‑239 (MQ=255)                                                                                                               tCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATGGCGCAAAACCCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATAc                                                                                                                         <  2:187418/238‑1 (MQ=255)                                                                                                                  ttGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACACCGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACa                                                                                                                     <  2:35143/239‑1 (MQ=255)                                                                                                                  ttGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCGGGCCCACTGCCAGCATATTGCAGCCCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACa                                                                                                                     <  2:857079/239‑1 (MQ=255)                                                                                                                  ttGAGAG‑CGATTTGCAGATTGTGCTGCCCGACCAGCGGCAACAGCGTTTTTGCTGGCCCCTGGCTGGCATATTGCAAGACCCAGCGGCCGGCGGAAGGA‑AGTGCTCTTGCGCCAAACCCTC‑AAAGCCCAAAGAACCGCCAGCCTTTCGCCTTTTTTTCGCCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACa                                                                                                                     <  1:742655/239‑1 (MQ=255)                                                                                                                  gtGATAT‑TGATTTGCAGATTGTTGTGCACCACCAGAGGCGACATCGTTTTTTTTGCCGCACTGCCAGCATATTCCACCGCCCAGCCGCCAGCGGGAAGC‑CTTGCTGTTGCACCAATACCTA‑CGATCTCATACACCCACCCGCGCTTCGCCTTCTTTTCACCACTCATTCCCACCTTCCTTGTGTATAGCCAGCCATTTTTTACGGGCCCAGACAAACTTCACCGCGCCCTAACACGACa                                                                                                                     >  1:498192/2‑239 (MQ=255)                                                                                                                   ttgaatccGATTTGAAGTTTGTTGTGCACCACCAGCGGCGACAACTTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACa                                                                                                                     <  1:991787/232‑1 (MQ=255)                                                                                                                                  ggttcgttgtGCACCGCAAGCGGCGAAAACGTGTTTCCTGGCCAACCACCAACACATTAAAACGCCCAGACACCAGCGCAAAGA‑AGTGCAATTGCGCCAAAACCAA‑CAAGACCAGAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGTGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACttt                                                                                                       <  1:591305/233‑1 (MQ=255)                                                                                                                                  aGATTGTTGTGCCCCCCCAGCGGCGACAACGTTTTTGCTGCCCCACTGCCAGCATATTGCACCGCCCAGACGCCAGCGGAAAGC‑ATTGCTATGGCCCAAAAACCTA‑CAAGACCAAAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACttt                                                                                                       <  1:710981/238‑1 (MQ=255)                                                                                                                                  aGATTGTGGTGCACCACCAGCGGCGACAACGTTTTTCCTGGCCAACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAACCCTA‑AAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACttt                                                                                                       <  2:218352/238‑1 (MQ=255)                                                                                                                                     ttgttgTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTg                                                                                                      <  2:662429/236‑1 (MQ=255)                                                                                                                                              accaccAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACtt                                                                                          >  1:665454/1‑239 (MQ=255)                                                                                                                                                caccaGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACTACAAAAGCCCAGACTTTGCAGCCTGGACttt                                                                                         >  1:43692/1‑238 (MQ=255)                                                                                                                                                caccaGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACttt                                                                                         >  2:1100539/1‑238 (MQ=255)                                                                                                                                                     gcggcgACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttcaat                                                                                    >  1:642897/1‑238 (MQ=255)                                                                                                                                                       ggcgACACGGTTTTTTCTGCCCCACTGCCAGCATATTGCACCGCCCAGACGCCAGCGAAAAGC‑ATTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttcaat                                                                                    <  2:51507/236‑1 (MQ=255)                                                                                                                                                       ggccACAACGTTTTTTCTGGCCCACGGCCAGCATATTGCACCGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCCCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTttcaat                                                                                    <  2:473829/236‑1 (MQ=255)                                                                                                                                                       agaaacggcGTTTTTCGTGGCCCTCTGCAAGAAAATTGCACCGCCCAGGCGACAGGTTAAAGC‑TGGTTTATTGCAAAAAACCCTC‑CAAGCCCATAAAGCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGAGCACAGCCAAACTTTACCGAGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTtacaattca                                                                                 <  1:438365/231‑1 (MQ=255)                                                                                                                                                          gACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAAc                                                                              >  1:719315/1‑239 (MQ=255)                                                                                                                                                          caCAACGTTTTTTCTGGCCCAATGCCAACAAATTGCAACGCCCAGACGCCGGTGGAAAGC‑AGTCATTTTGCGAAAACACCTA‑AAAAACCATGAAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAAc                                                                              <  1:915215/238‑1 (MQ=255)                                                                                                                                                           aCAACGTTTTTTCTGTCCCACTCCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACa                                                                             <  2:908320/239‑1 (MQ=255)                                                                                                                                                              ccGTTTTTTCTGGCCCCCTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAgg                                                                          <  2:989437/238‑1 (MQ=255)                                                                                                                                                                       cTGCCCCCCTGACAGCATATTGCAACGCCCAGACGCCAGCGAAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCt                                                                 <  1:390668/239‑1 (MQ=255)                                                                                                                                                                       cTGCCCAACTGCAAGCATATTGCACCGCCCAGACCGCAGCGGAAAGC‑AGTTCTATTGCGCCAAACCCTA‑CAACACCATATACCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCt                                                                 <  2:479196/239‑1 (MQ=255)                                                                                                                                                                          gCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGgag                                                                       >  2:438365/1‑230 (MQ=255)                                                                                                                                                                                    aGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATg                                                    >  2:430923/1‑239 (MQ=255)                                                                                                                                                                                       atatTGCAACGCCCGCGCGCCAGCGGAAAGC‑AGTGTTATTGCGCAAAAACCAA‑CACGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaa                                                  <  2:208464/238‑1 (MQ=255)                                                                                                                                                                                       atatTGCAACGCCCAGCCGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAACCCTC‑CAAGCCCATAGATCCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATga                                                   <  2:527779/237‑1 (MQ=255)                                                                                                                                                                                       atatTGCAACGCCCAGACGCCAGCGAAAAGC‑AGTGCTATTGCGCCAAACCCTA‑CAAGCCCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATgaag                                                 <  2:670872/239‑1 (MQ=255)                                                                                                                                                                                                cGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAACAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAat                                        >  1:929305/1‑239 (MQ=255)                                                                                                                                                                                                 gCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAatt                                       >  1:901465/1‑239 (MQ=255)                                                                                                                                                                                                 gCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAatt                                       >  1:1187643/1‑239 (MQ=255)                                                                                                                                                                                                           cAGGGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAGGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttctt                              <  2:1027267/238‑1 (MQ=255)                                                                                                                                                                                                            aGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTAttcttct                               >  1:971623/1‑236 (MQ=255)                                                                                                                                                                                                               ggAAAGC‑AGTGTTATTGCGCAAAAACCTA‑CAAGACCAAAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATtcttcttct                            <  2:1096244/236‑1 (MQ=255)                                                                                                                                                                                                                     c‑AGTGCTGTTG‑GCCAACGCCTCCCAAGACCTGGGAGCCGCCAGCCGTTCGCCTTTCTTTCACCACTCATGGCCGCATTCCTTGTGTATAGCCAGCCATTATTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCtcgtcg                   <  1:615511/239‑1 (MQ=255)                                                                                                                                                                                                                          gCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCa                <  2:665454/238‑1 (MQ=255)                                                                                                                                                                                                                              ttGCGCCAAACCCTA‑AAAGACCAAAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTc           <  1:43215/239‑1 (MQ=255)                                                                                                                                                                                                                              ttGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTc           <  1:1183769/239‑1 (MQ=255)                                                                                                                                                                                                                                 cgcCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCAc        >  1:446963/1‑239 (MQ=255)                                                                                                                                                                                                                                   ccAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACt       >  2:651573/1‑238 (MQ=255)                                                                                                                                                                                                                                       aaCCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCgg  <  2:423117/239‑1 (MQ=255)                                                                                                                                                                                                                                       aaCCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCgg  <  1:1100539/239‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                     TCTGATAATAGAAAACGAAATCTTCGCTCAGAACATCGTGGAGAAACGGAATGGTGAGGAGATCTTTGGGAAGCTGGCTCAGAGAGTCGCTGTCGAGAAAGAGGTCCGGCTCATTGAGAT‑CGATTTGCAGATTGTTGTGCACCACCAGCGGCGACAACGTTTTTTCTGGCCCACTGCCAGCATATTGCAACGCCCAGACGCCAGCGGAAAGC‑AGTGCTATTGCGCCAAAACCTA‑CAAGACCATAGAACCGCCAGCCTTTCGCCTTTTTTTCACCACTCATTGCCACATTCCTTGTGTATAGCCAGCCATTTTTTACGGGCACAGCCAAACTTTACCGTGCCCTAATACGACAAAAGCCCAGACTTTGCAGCCTGGACTTTTCAATTCAAACAAGGGAGATAGCTCCCTTTTGGCATGAAGAAGTAAAATTATTCTTCTTCTGGCTCGTCGTCAACGTCCACTTCCGG  >  NC_000913/2336363‑2336829

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