Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,145,868	A→G	14.8%	E868G (GAG→GGG)	yegE →	putative diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,145,868	0	A	G	14.8%	229.3 / 7.0	95	E868G (GAG→GGG)	yegE	putative diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base A (33/47);  new base G (12/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 9.58e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCC  >  NC_000913/2145632‑2146103                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               tGACTTTCGTTTTAAATGGGAAGGCCGTGTGCAACGGGGAGGTGCCAGTGCCGGGATTACCTTGATTGATGCCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAgaga                                                                                                                                                                                                                                              <  2:721255/239‑1 (MQ=255)  gACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGaa                                                                                                                                                                                                                                             <  1:982746/239‑1 (MQ=255)    ggtataagaggcaggggAAGGCCGTGTGCACCGGGTAGGTGCCAGTTGCGGGATTCCCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑c                                                                                                                                                                                                                                          <  1:1005372/225‑1 (MQ=255)               tatGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCTTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAAGGAAT‑CCGTTgaggat                                                                                                                                                                                                                                >  1:626983/1‑238 (MQ=255)               tatGGGAAGGCCGTGGGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGatgat‑g                                                                                                                                                                                                                              <  1:916059/239‑1 (MQ=255)                     aaGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGGGCCGGTTTACGAACCGCAGCAACCTGCCGCACCTAGCGAACGGGCGGCGATGTCGCTTGA‑GTAAAGTGGCGGGGTGATTCAAAAGGAC‑AAGTtgatgtggacgggg                                                                                                                                                                                                                         >  2:957028/1‑230 (MQ=255)                       ggCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGCTGAACAGTGGC‑GGATGGTTAAAAGAAAT‑CAGTTGTGGGA‑GCTCGccc                                                                                                                                                                                                                       >  1:1031789/1‑238 (MQ=255)                         ccGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGCACCGCAGCAAGCTGCCGCACATAGCGAACGGCGGGCGGGTGCGCTTGATGAAACGGGCC‑GGGGGATTACAGAGAAA‑CAgaggttggtgccaccccacg                                                                                                                                                                                                                     >  2:199057/1‑220 (MQ=255)                            tgtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGGTGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAAACGCAGCAAGCTGCCGCACCTAGCGAACGGGCGGCCGACTCCCTTGGTGAACCGGGGC‑GGGCGACTCAAAAGAAT‑Cacattatggaacgcgcccccaggg                                                                                                                                                                                                                  >  2:430830/1‑216 (MQ=255)                            tgtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGGTTAAAGAGAAT‑ACGTTGGTGAT‑GTTCGCCCACgggg                                                                                                                                                                                                                 >  1:693686/1‑237 (MQ=255)                             gtgCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACATTGGC‑GGATGATTAAAGGGAAT‑CAGTTGGTGAT‑GCTGGCCCACgggg                                                                                                                                                                                                                 >  1:1075549/1‑236 (MQ=255)                                        aGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCa                                                                                                                                                                                                     <  1:810942/238‑1 (MQ=255)                                         tggtccagggcgcggATTACCTTGATTGGTGACAACAATCATCAAGCGGCAGAAGGGATGCCGCCGGGTGATATCCGCTGTTATTCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCggctcg                                                                                                                                                                                                    <  2:216196/226‑6 (MQ=255)                                          gggCCAGTGCCCGGATTTCCTTGATTGATGACAACAATCAGCAGGCGGCAGAAGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACg                                                                                                                                                                                                   <  1:674036/236‑1 (MQ=255)                                          gTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGTTGATTAAAGGGAAT‑CAGTTGATGAT‑GCTCCCCCACGGTGGCCCTTCG‑CCCCgg                                                                                                                                                                                                  >  1:688703/1‑239 (MQ=255)                                             ccAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGCTGATTTAAAGGAAT‑CGAGTGATGAT‑GCTGGCCCCCGGTTGCGCTTCG‑CCCCgggt                                                                                                                                                                                                >  1:1082035/1‑236 (MQ=255)                                               aGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGCTGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGGGTCGCTTTG‑CCACGGAcccc                                                                                                                                                                                             >  1:125787/1‑239 (MQ=255)                                                  ggtgggTTTCCCTTGATTGATGACAACAATCATCAGTGGGCAGAAGGGAGGTCGCAGGCTGATATCGCTTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGa                                                                                                                                                                                          <  2:626983/236‑1 (MQ=255)                                                         tttccTTGATGTATGAAAAAAATCTTCAGGCGGCAGAAGGGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGcgcg                                                                                                                                                                                    <  1:1016010/235‑1 (MQ=255)                                                                 attgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAACTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCCTAGCGAACGGGCGGCGATGTCGCTTGATGAACGGGGGC‑GGGTGATTAAAGAGGAA‑CAGTTGGTGATCGCT‑GCCCAACGGGTGGCGTAG‑CCCGGGACCCCGGAAgcaggaaatctgag                                                                                                                                                                           >  2:416132/1‑227 (MQ=255)                                                                 attgatGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGGGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTGGATGAACAGTGGC‑GGATGATTTAAGGGAAT‑CCGTTGGTGAT‑GCCCGCCCACGGGGGCGCTTCC‑CCCCCGGACCCGGAAAcgcggatattttg                                                                                                                                                                           >  2:676471/1‑230 (MQ=255)                                                                     atgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCCAGCTGCCGCACATAGCGAACGGGGGGCGATGTCGCTTGATGAACCGTGGC‑GGCAGATTAAAGGAGAT‑CAGGTGGTGAT‑TCACCCCCCCGGGGTCGCTCTG‑CCACGGGTCCCGGAAGGGCGAAAATTGTGGGTg                                                                                                                                                                       >  2:82601/1‑239 (MQ=255)                                                                      tgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAAGGAAT‑CAGTTGATGAT‑GCTCCGCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGGGGGTAATTTTTGGCTGa                                                                                                                                                                      >  1:681949/1‑239 (MQ=255)                                                                      tgaCAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGGCGCTTGATGAACAGTGGC‑GGACGATTAAAGAGAAC‑TCGTTGATGAT‑GCTCGCCCACGGGTGCGCTTCG‑CCACCGGGCCCgggaaccggagattttgggccgg                                                                                                                                                                      >  2:958277/1‑218 (MQ=255)                                                                       gaCAACCATCATCAGGCGGCAGAAGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  1:540025/239‑1 (MQ=255)                                                                       gaCAACAATCATCGGGGGGGAGAAGTGGTGGCGCTGGCTGATATCCGCCGTTATGCCTCCAAAAATGGTGGCCGGGGGCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCAAATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCTCGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  2:894777/239‑1 (MQ=255)                                                                       gaCAACAATCATCAGGCGGCAGAAGGGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGAt                                                                                                                                                                     <  1:512085/239‑1 (MQ=255)                                                                            caaTCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCCCATAGCGAACGGGCGGGGATGTACCTTGCGGTAAAGTGGC‑GGATGGTTACAAGGAATACAGTTGATG‑T‑GCTCCCCCACGGGGGTCCTTCGCCCAGGGACCCCGGAA‑CGCGAGATTTTTGGcgagttttcc                                                                                                                                                                >  2:1093676/1‑230 (MQ=255)                                                                              atcatcAGGCGGCAGAAGTGATGTCGCAGGCTGATTTCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                                                                                                                              <  1:1010859/239‑1 (MQ=255)                                                                                    cggcggcAGAAGTGGTGTCGCAGGCTGGTTTCGCCTGTTATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGGAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGctc                                                                                                                                                        <  2:499532/238‑1 (MQ=255)                                                                                              aGTGATGGCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTGATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                                                                                                                              <  1:543967/239‑1 (MQ=255)                                                                                                tGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGAGTAAAGAGAAT‑CAGTAGAAGAG‑TCTCGCCCACGGTGTGGCCTCG‑CCACCGGTGCCGGAAACGCGAGCTTTGGTGGTGTTTTCACCTaaagcctggcggtgggga                                                                                                                                            >  2:949956/1‑223 (MQ=255)                                                                                                  gtGTCGCGGGGTGATATCGCCTGGTCTGCCTCCAAAAATGGGGGACGGGGCCGGGGGCAGGTTTACGAACCGCAGCAAGCTGCCGGACATAGCGAACGGGCGGCGATGGCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAgg                                                                                                                                          <  2:1027661/238‑1 (MQ=255)                                                                                                   tgtcggagcgggttaccGCTTGTTAGCCCTCAAAAAATGGTGGGCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGACGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACGTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:1016533/223‑1 (MQ=255)                                                                                                   tggcgcctgCTTATTTCGCCGGTTCTGCCGCAAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGAACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:688703/231‑1 (MQ=255)                                                                                                   tGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAAGGGAT‑ACGGTGATGGT‑GGCCGCCACCCGGGTCCGTTCG‑CCACGAGTCCCGGAAGCGCGTAAATTTTGGGTGATATCCATTaaaaccaggggggtggaaagg                                                                                                                                         >  2:526773/1‑220 (MQ=255)                                                                                                   gatttgtataagggcccGCCTGTTATGCCTCCAAAAAGGGTGGCCGTGGCCGGGTGGCGGGTTACGAACCGCAGCAAGCTGCCGCATATAGGGAACTGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGc                                                                                                                                         <  2:337400/223‑1 (MQ=255)                                                                                                         aGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGat                                                                                                                                    <  1:380713/238‑1 (MQ=255)                                                                                                          ggCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGGAATTTGTGGCTGATTTCACTTAAGCTCTGGGGTTGCGAAGGCGGGatta                                                                                                                                  >  1:60469/1‑239 (MQ=255)                                                                                                             tgaatacGCCTTTTAGGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               <  1:430830/233‑1 (MQ=255)                                                                                                               atatCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGTGAAGGCGGGATTATtgat                                                                                                                             >  1:227406/1‑239 (MQ=255)                                                                                                                 atCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTTACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACtt                                                                                                                                                              <  1:626775/204‑1 (MQ=255)                                                                                                                  tCGCCTGTTATGCCTACAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGA‑GATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaac                                                                                                                         <  1:812685/239‑1 (MQ=255)                                                                                                                    gCCTGTTATGCCTCCAAAAAGGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGaaca                                                                                                                        <  1:82601/239‑1 (MQ=255)                                                                                                                     ccTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGTATCCCGGAAGCGCGGAATTTGTGGCTGATTTCACTTAAGCCCTGGAGTTGCCAAGGCGAGATTATTGATGaaaaa                                                                                                                       >  1:226648/1‑239 (MQ=255)                                                                                                                         ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTTAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTTGCTTCG‑CCCCGGATCCCGGAAGCCCGGAATTTTTGGCTGATTTCCCTTTAGCCCTGGGGTTTGGCAGGGGAGGTTTTtgagaaaaaaaaaa                                                                                                                   >  2:839367/1‑228 (MQ=255)                                                                                                                         ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTTACTTTAGCTCTGGGGTTGCGAAGGCGAGATTATTGATGAACAAACAt                                                                                                                   >  1:483333/1‑239 (MQ=255)                                                                                                                         ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAt                                                                                                                   >  1:245847/1‑239 (MQ=255)                                                                                                                         ttATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGAGTCGCTTCG‑CCACGGCTCCCGGAAGCGCGTAATTTTGGGCTGATTTTCTTTAAGCCCTGGAGTTGCGAAGGGCAGATTATTGATGGACAAACAt                                                                                                                   >  1:592521/1‑239 (MQ=255)                                                                                                                          tATGCCTCCAAAAATGGGGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACAtt                                                                                                                  <  1:358459/239‑1 (MQ=255)                                                                                                                                tCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTTTtgatggttccttttcaacatt                                                                                                            >  1:337400/1‑223 (MQ=255)                                                                                                                                   aaaaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGGGGTA‑CCGTTGAGGAT‑CTTCGCCCACGGTGTCGCTTCG‑CCCCGGCTCCCCGCAGCCCGTAATTTGTGGGCGATTTTCCTTAAGCTCTGGAGGTGGCGAGCGGGGAGTATTGGTGAAAAAACATTtctgccggg                                                                                                         >  2:848283/1‑232 (MQ=255)                                                                                                                                   aaaaaTGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGGGTCGCTTCG‑CCACGGAACCCGGAAGCGCGTAATTTGTGGCTGATTTTACTTAAGCTCTGGGGTGGGCAAGGGGGGATTCTTGAGGaacaaaaaattttggggag                                                                                                         >  2:521304/1‑226 (MQ=255)                                                                                                                                            ggccggggccggGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTTGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTCGCTTCAGc                                                                                                  >  1:831083/1‑237 (MQ=255)                                                                                                                                                ggggccgggGGACGGTTTACGAACCGCAGCAAGCTGCCGCGCATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCg                                                                                            <  1:1036290/239‑1 (MQ=255)                                                                                                                                                  ggccgggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGCAGGCGGGATTATTGATGAACAAACATTTCTTCTTAGCTTCAGCGTTCCgc                                                                                           >  1:574734/1‑237 (MQ=255)                                                                                                                                                       ggTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGCTGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGGGCtt                                                                                      >  1:1028535/1‑238 (MQ=255)                                                                                                                                                           aCGGTTTGCGAACCGCAGCAAGCTGCCGGACATAGAGAACGGGCGGGGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                  <  1:354871/238‑1 (MQ=255)                                                                                                                                                           aCGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGGTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGcc                                                                                  <  1:1093676/238‑1 (MQ=255)                                                                                                                                                            cGGTTTACGAACCGCAGCAAGCTGCCGCAAATAGCGGACGGGCGGCGATGTCGCTTGATGAACAGGGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCa                                                                                 <  1:101428/238‑1 (MQ=255)                                                                                                                                                                  tcggaCCGCAGCAAGCTGCCGCAAATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCtt                                                                          <  1:866734/235‑1 (MQ=255)                                                                                                                                                                       ccGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAcc                                                                      >  1:383787/1‑238 (MQ=255)                                                                                                                                                                        cGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGaa                                                                                                                                                                                         <  2:779425/122‑1 (MQ=255)                                                                                                                                                                        cGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGaa                                                                                                                                                                                         >  1:779425/1‑122 (MQ=255)                                                                                                                                                                         gcagcaAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGGGAGATTATTGATGAACAAACATTTGTTCGTAGCTTCAGCGATCCGGGGCTTGGCCATACTCTTGAccgcc                                                                   >  2:504615/1‑239 (MQ=255)                                                                                                                                                                          gagcaAGTCGCCGTACATAGCGAACCGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGAccgccg                                                                  <  1:506539/238‑1 (MQ=255)                                                                                                                                                                                  tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTTATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATGTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               <  2:54534/170‑1 (MQ=255)                                                                                                                                                                                  tGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATtg                                                                                                                               >  1:54534/1‑170 (MQ=255)                                                                                                                                                                                    ccGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGATCTTGACCGCCGGGGATTCCAc                                                        >  1:971651/1‑239 (MQ=255)                                                                                                                                                                                     cGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGGGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAACGCGGAATTTGTGGCTGATTTAACTTAAGCTCTGGAGTTGCGAAGGGGAGACTATTGTTGAACCAACATTTCCTCGTTACTTTCGCGTCACCGCGCTTAGCCATGTTCTTGACCGCCGGGGGTTCCAcc                                                       >  2:1068945/1‑238 (MQ=255)                                                                                                                                                                                       cccgttcgctatgtgcggcgATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGaa                                                     >  1:4302/15‑239 (MQ=255)                                                                                                                                                                                            aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCg                                                                                                                                              <  1:526160/145‑1 (MQ=255)                                                                                                                                                                                            aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                       <  2:857595/200‑1 (MQ=255)                                                                                                                                                                                            aGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCt                                                                                       >  1:857595/1‑200 (MQ=255)                                                                                                                                                                                              cGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGGTTTTTGATGAACAAACCTTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATTCTCTTGACCGCCGGGCATTCCACGAATTTTTcc                                              >  1:655822/1‑239 (MQ=255)                                                                                                                                                                                                  cgggcgggGATTTTGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  2:245847/239‑1 (MQ=255)                                                                                                                                                                                                  cGGGCGGCGGTGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  1:199057/239‑1 (MQ=255)                                                                                                                                                                                                  cGGGCGGCGATGTCGGTTGATGAACAGTGGC‑GGATGATTAAAGGGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGGGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCcagca                                          <  1:785819/239‑1 (MQ=255)                                                                                                                                                                                                      cggcgATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTATCCATGCACTTGACCGCCGGGTATTTCCCGCAATTTTCCAGCAGGc                                       >  1:1104268/1‑238 (MQ=255)                                                                                                                                                                                                                  ttGGTGAAAAGTGTC‑GGAGAATTAAAGAGAAT‑CAGGTGAGGAG‑GCGCGCCCACGGGGTCGTTTAG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                          <  1:242691/239‑1 (MQ=255)                                                                                                                                                                                                                  ttGATGAACGTTGGC‑GGTTGATTAAAGAGAAT‑CAGTTGATTAT‑GCTCGCCCACGGTGTCGCGTCG‑CCCCGGATCCCGGAAGCGCGTAATTTGGGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGtt                          <  2:1028535/239‑1 (MQ=255)                                                                                                                                                                                                                   gtgtgaAGAGTGGG‑GGATGAGTAAAGAGAAT‑CAGGTGGTGAT‑GCGCGCCCGCGGGGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGGGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGGGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTg                         <  2:592521/236‑1 (MQ=255)                                                                                                                                                                                                                    ttggaACAGTTGG‑GGGTGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCCCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGGGGGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                        <  2:1031789/236‑1 (MQ=255)                                                                                                                                                                                                                        aaCAGTGGG‑GGGTGATTAAAGGGAAT‑CAGTTGATTAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGc                        <  2:254934/235‑1 (MQ=255)                                                                                                                                                                                                                         aCAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                   <  1:416132/239‑1 (MQ=255)                                                                                                                                                                                                                         aCAGGGGG‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTa                   <  1:1068945/239‑1 (MQ=255)                                                                                                                                                                                                                               gC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAACCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCCCGGATTTTTCCCACAGGCCCCAAAAGCGGTTGCCAGTAAAGGCa             >  1:439542/1‑239 (MQ=255)                                                                                                                                                                                                                                c‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCa             <  1:958277/238‑1 (MQ=255)                                                                                                                                                                                                                                  ggctgtttAAAGAGATC‑CAGTTGATGAT‑GCGGGTCCGCGGTGTCGCTTCG‑CCACGGATCCAGGAAGGGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGGTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt            <  2:1104268/232‑1 (MQ=255)                                                                                                                                                                                                                                  ggATGATTAAAGGGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCAt            <  1:938700/238‑1 (MQ=255)                                                                                                                                                                                                                                   gatgatTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGc              >  1:254934/1‑235 (MQ=255)                                                                                                                                                                                                                                   gatgatTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCACCAGGCCGCAAAACAGGTTGCCAGTAAAGGCATaa          >  1:72498/1‑239 (MQ=255)                                                                                                                                                                                                                                     tgatTAAAGAGAAT‑CAGTTTATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAg         <  1:949956/238‑1 (MQ=255)                                                                                                                                                                                                                                       atTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCTCGGTTTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCAt      <  1:264319/239‑1 (MQ=255)                                                                                                                                                                                                                                          aaaGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCCGGAGTTGCGAAGGCGAGATTATTGATGCACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGCATTCCTCGAATTTTTCCAGTAGGCCGCAAAATCTGTTGCCAGGTAAGGCATAAGCAtcc    >  1:75227/1‑237 (MQ=255)                                                                                                                                                                                                                                           aaGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCTTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATctcc  <  1:782678/239‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               TGACTATCACTTTATATGGGAAGGCCGTGTGCATCGGGTAGGTGCCAGTGCCGGGATTACCTTGATTGATGACAACAATCATCAGGCGGCAGAAGTGATGTCGCAGGCTGATATCGCCTGTTATGCCTCCAAAAATGGTGGCCGGGGCCGGGTGACGGTTTACGAACCGCAGCAAGCTGCCGCACATAGCGAACGGGCGGCGATGTCGCTTGATGAACAGTGGC‑GGATGATTAAAGAGAAT‑CAGTTGATGAT‑GCTCGCCCACGGTGTCGCTTCG‑CCACGGATCCCGGAAGCGCGTAATTTGTGGCTGATTTCACTTAAGCTCTGGAGTTGCGAAGGCGAGATTATTGATGAACAAACATTTCGTCGTAGCTTCAGCGATCCGGCGCTTAGCCATGCTCTTGACCGCCGGGTATTCCACGAATTTTTCCAGCAGGCCGCAAAAGCGGTTGCCAGTAAAGGCATAAGCATCTCC  >  NC_000913/2145632‑2146103

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