Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	624,544	Δ1 bp	4.4%	intergenic (‑34/‑341)	fepB ← / → entC	ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	624,544	0	C	.	4.4%	277.3 / 6.1	90	intergenic (‑34/‑341)	fepB/entC	ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base C (41/42);  new base . (2/2);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

GTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCG  >  NC_000913/624265‑624827                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                            gTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACa                                                                                                                                                                                                                                                                                    >  2:793664/1‑289 (MQ=255)         cgACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGCCCTGTTAATAGAAGGGCGTTGCGGCAGAGAGGGGCGAGTCTCACAAATCAGCTTTCTGTTATTAATAAGGTTAAGGGCGTAATGACAAAATCGag                                                                                                                                                                                                                                                                            >  1:791942/1‑288 (MQ=255)            cgcgGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTCATTAAAAAGGTTAAGGGGGTAATGGCAAAATCGACaca                                                                                                                                                                                                                                                                         >  1:369165/1‑287 (MQ=255)                  gtttCGGCGTGGTCGCGCGGACGGGCAGTACCGGAGCATCACTCGCCAGCAGTGGGCCGGTCAGGGTGCCGCTGGTGGAAACAATACGCTGCGGCTAGCGTTCCAGT‑GTATGTGTGCCAGGGCTGTCAGTAATCGGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAACGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGca                                                                                                                                                                                                                                                                    <  1:638660/286‑1 (MQ=255)                        gctttgTCGGGCCGCTGCCGATGACCGGAGCATGAATCGCCAGCAGTGGTCCGGTCAGGGTGACGCTGGTTGAAAGAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGTCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTTTTAATAGAAGTGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGTGCGTAATGACAAATTCGACAAAGCGCACAATcc                                                                                                                                                                                                                                                              <  1:1143771/283‑1 (MQ=255)                             gTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCTCGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCATGAAAGCACAAGTCCTGTTAATAGAAGGGCGTTCCGGTAGCGCGGGGCGAGTCTCACACACCAGCTTCCTGTTATTTATAAGGTTAAGGGCGTACTGACAAATTCGCCAAAGCGCACAATCTGTccc                                                                                                                                                                                                                                                         >  1:273668/1‑288 (MQ=255)                                gcgcCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACTCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCg                                                                                                                                                                                                                                                     <  2:976341/289‑1 (MQ=255)                                gcgcCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCg                                                                                                                                                                                                                                                     <  2:323663/289‑1 (MQ=255)                                gcgcCGCTGGCGATCACCGGAGCAGCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCg                                                                                                                                                                                                                                                     <  1:416468/289‑1 (MQ=255)                                   ccgcTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGATATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGcc                                                                                                                                                                                                                                                   >  2:64938/1‑288 (MQ=255)                                      cTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCtt                                                                                                                                                                                                                                               >  1:1068102/1‑289 (MQ=255)                                            aTCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCCCAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         >  1:670368/1‑289 (MQ=255)                                             tCACCGGAGCATCAATCGCCAGCAGTGAGCCGTTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  2:331810/288‑1 (MQ=255)                                             tCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  2:1068102/288‑1 (MQ=255)                                             tCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  1:328377/288‑1 (MQ=255)                                             tCACCGGAGCATCAATCGCAAGCAGTGAGCCGGTCAGTGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTAGGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  1:320181/288‑1 (MQ=255)                                                         cAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         >  1:246585/1‑276 (MQ=255)                                                         cAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  2:246585/276‑1 (MQ=255)                                                                cagcagTGGGCCGGTCAGGGTGCCGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCCCGGCTGTCAGTCATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTTTTAATAGAAGGGCGTTGGGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTATGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                                                                                     <  2:369165/289‑1 (MQ=255)                                                                   cgtgaaGCCGGTCGGGGGGAGGGGGGTGGAACCATGACGCGGCGGCTGGCGTT‑CAGTGGTATGTATGCCCGGGCTGCCAGTAAGCTGACGGGGACG‑TGCAGCGACCTTAACGTGC‑GCTATCCCTGAAAGAAAAAGTCCTGTTAATAGACGGGCTTTGCGGTAGGGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATGCGACAAAGCGCAGAATCCGTCCCCTGGCCCCTTTGGGGAGCGGGTTAGGTTGAGGGGAGca                                                                                                                                                                                                                   <  1:1390574/283‑1 (MQ=255)                                                                             gTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCCGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGCCGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTTAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGACCAGCCAGCCCTgg                                                                                                                                                                                                        >  1:435623/1‑289 (MQ=255)                                                                             gTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTAGTCCTGACAGCAAAAGTCCTGGTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAAACAGCTTCCTGATATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAACGCCCAACCCGCCCCCTCGCCGCCTTGGGTAGCGGGTGAGGGTGAGGGGACcagccagccatgg                                                                                                                                                                                                        >  1:1159118/1‑284 (MQ=255)                                                                                  ggTGACGCTGGGGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GGAGGTGGGCCACGGCTTTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGGGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                                                                                                                    <  1:88751/288‑1 (MQ=255)                                                                                       cGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCC‑GTCTGCGGCTTGAAC‑TGCGGCTATTCCTGAAGGCAAAAGTCCCGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGCGTCTCACAAATCAGCATCCCGTTGTTAATAAGGTTAAGGGCCGAATGACAAATTCGACAAAGAGCAACACCCCTGCCCTAGCCGCTTTTGGGAAGGGGTTAGGGTGTGGGGAACAGCCCGCCCTTGCGCGAACAtt                                                                                                                                                                                              >  1:676061/1‑289 (MQ=255)                                                                                            gtggAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAAcc                                                                                                                                                                                          >  2:597414/1‑288 (MQ=255)                                                                                              ggAAACAATCCCGTGCTGCCGGTCTTCCAGT‑TGATGTGTGCCAGGGTGTTCAGTAATCAGAGCCGGCAA‑GTCAGCGGCCTGAAC‑TGCGGCTTATCCTGAAAACAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCAAAAATCAGCTTCCTGTTAGTAATAAGGTTAAGGGCGTAATGACAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                                                                                                                        <  1:793664/288‑1 (MQ=255)                                                                                                 aaCAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACaa                                                                                                                                                                                                                                                                          >  2:173736/1‑203 (MQ=255)                                                                                                 aaCAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACaa                                                                                                                                                                                                                                                                          <  1:173736/203‑1 (MQ=255)                                                                                                            gCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGACTGGCAGTAATCTGGCGGTGCCA‑GTCAGCGTCCTAAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCAAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCA‑CCCAGCCCTCAcc                                                                                                                                                                        <  2:37997/289‑1 (MQ=255)                                                                                                              ggctggctTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTCTTCCTGAAAGCAAACGTCCTGTTAATAGAAGGGCGTTGCGGTACAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTCTTAATAAGGGTAAGGGCGTAATGACAAATTCGACAAAGCGAACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATGAACCCTCACCCCAGcccccccc                                                                                                                                                                        >  2:111672/1‑283 (MQ=255)                                                                                                                              gtATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAgggg                                                                                                                                                        >  1:442581/1‑289 (MQ=255)                                                                                                                               tATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGCAGGGAGAgggg                                                                                                                                                        >  2:572163/1‑288 (MQ=255)                                                                                                                               tATGTGTGCAACGGCTGTCAGTAACCCGACGCGGACA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGTAAAAGTCCTGTTAATAGAAGGGCGTTGAGGGAGAGCGGGGCGAGTCTAACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAACCCGTCCCCTCCCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGAACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggggg                                                                                                                                                       <  2:943750/289‑1 (MQ=39)                                                                                                                                    gtgCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTACTAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGa                                                                                                                                                   >  2:627209/1‑288 (MQ=255)                                                                                                                                       ccACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTCATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCAAGCACTGGTTCGAACATTAACCCTCACCCCGGCCTCAACCCTGGAAGGGCGAGGGGGGAGAACgg                                                                                                                                               >  1:943750/1‑289 (MQ=255)                                                                                                                                       ccACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATAAGCTTCATGTTATTAATACGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCATTGGGGAGAGGGTTAGGGCGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCCCCCTCAAGCCTTACCCCTGGCAGGGTGAGGGGGCAGACCg                                                                                                                                                >  1:332745/1‑288 (MQ=255)                                                                                                                                         aCGGGTTTCAGTAATCTGACGCGGCCA‑TGAAGGGGCCTGAAC‑TGCGCGTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCAAAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGGGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                                                                                                             <  1:597414/289‑1 (MQ=39)                                                                                                                                         aCGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                                                                                                               >  2:751345/1‑207 (MQ=255)                                                                                                                                         aCGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAgg                                                                                                                                                                                                                               <  1:751345/207‑1 (MQ=255)                                                                                                                                           ggCTGTCAGTAATCTGACGGGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCTTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGTTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                                                                                           <  2:1346404/289‑1 (MQ=255)                                                                                                                                            gcggtagggaaTGTGAAGCGGCCA‑AGCAGAGGCCTGAAG‑GCGGGCTATGCCTGAAAGAAAGAGTCCTTTTAATAGAAGGGCGTTGGTGGAGAGCGGGGCAGGTCTAAGGATGCAGCTTGTTGGGATTAATAAGGTTAAGGGGGTAATGAAAAAGTCGACGAAGCCCACAAGCCGACCCCTGTACCCTTTGGAGAGAGGGTTAGGGGGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCTCCCCAGCCCTCACCCTTGTAGGGTGAGGGGGCAGAACGGCGCa                                                                                                                                           <  1:39369/279‑1 (MQ=255)                                                                                                                                                      tgtctaaggggCCAGGTCAGCGGCCGGAAC‑GGCGGGTATTGCGTAAAGAAAAGGCCCGGTTGGTGGGAGGGCGTGGCGGGAGAGCGGGGTGGGTCTAAGAAATACGTTTCCGGTTATAAATAGGGTTAAGGGCGTAAGCACAAGTGCAAAAAAGCGCAAACTCCGCCGCCTCGCCCCTTGGGGGGGAGGGTTAGGGTGGGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                                                                                                                 <  1:111672/280‑1 (MQ=255)                                                                                                                                                           aCGCGGCCA‑GTCAGGGGCCTGAAC‑TGCGGCTATGCCTAAAAGCAAAAGTCCTGTTAATAGAAGGGGGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgcg                                                                                                                            <  2:791942/288‑1 (MQ=255)                                                                                                                                                            cgcgGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCt                                                                                                                          >  1:768370/1‑289 (MQ=255)                                                                                                                                                            cgcgGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCt                                                                                                                          >  1:381763/1‑289 (MQ=255)                                                                                                                                                                gCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATg                                                                                                                      >  1:580996/1‑289 (MQ=255)                                                                                                                                                                  cA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCg                                                                                                                    <  2:442581/289‑1 (MQ=255)                                                                                                                                                                      tCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  1:148452/169‑1 (MQ=255)                                                                                                                                                                      tCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         >  2:148452/1‑169 (MQ=255)                                                                                                                                                                      tCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                                                                                                                  <  2:381763/288‑1 (MQ=255)                                                                                                                                                                       cAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATc                                                                                                                <  2:580996/289‑1 (MQ=255)                                                                                                                                                                        aGCGGCCTGAAC‑GGCGGGTATTCCTGAAAGCAAAATGCCTGTTAATAGAAGGGCTTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                                                                                                               <  2:73254/289‑1 (MQ=255)                                                                                                                                                                                gAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                                                                                                                    >  1:936039/1‑196 (MQ=255)                                                                                                                                                                                gAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCg                                                                                                                                                                                                    <  2:936039/196‑1 (MQ=255)                                                                                                                                                                                                     cccAAAAGTTCTTTTAATGTAGGGGCGTTGCGGTTGAGCGGGGCGAGTCTCACAAAGCAACTTCCTTGTATTAATAAGGTTCAGGGCGCAATTACAACTTCGACAAAGCGCACAAGCCGTCCCCTCCCCCCTTTGGGGAGAGGGTTAGGGTGAGGAGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTa                                                                                   <  2:676061/287‑1 (MQ=255)                                                                                                                                                                                                     aGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         <  2:577473/139‑1 (MQ=255)                                                                                                                                                                                                     aGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                                                                                                                         >  1:577473/1‑139 (MQ=255)                                                                                                                                                                                                                  ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGG‑‑GTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                                                                                                                  >  1:865623/1‑211 (MQ=255)                                                                                                                                                                                                                  ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGG‑‑GTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                                                                                                                  <  2:865623/211‑1 (MQ=255)                                                                                                                                                                                                                  aaatagaaggggtgtggggggggggtggggCG‑GTCTCAAAAATCAGATTCATTGTATTAAAAAGGTAAAGGGCGAAATGGCAAATGCGAAAAAGCGCAGAAGCCATCCCCGCACCCCGTTGGGGGGAGGGTTAGGGTGAGGGGAGCAGCCAGCACGGGTGCGAACATTAACCCTCGCCCCAGCCCTCACCCGGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGAt                                                                      <  1:1168306/262‑1 (MQ=255)                                                                                                                                                                                                                   tAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATg                                                                     <  2:768370/289‑1 (MQ=255)                                                                                                                                                                                                                         aaGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGca                                                                                                                                                                                                                                                                    >  2:325734/1‑92 (MQ=255)                                                                                                                                                                                                                         aaGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGca                                                                                                                                                                                                                                                                    <  1:325734/92‑1 (MQ=255)                                                                                                                                                                                                                           gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAgcg                                                                                                                                                                                                                                                                      <  1:76401/88‑1 (MQ=255)                                                                                                                                                                                                                           gggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAgcg                                                                                                                                                                                                                                                                      >  2:76401/1‑88 (MQ=255)                                                                                                                                                                                                                            ggCGTTGCGGTGGAGCGGGGCGAGTCTAACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGAGTTAGGGGGAGGGGAACAGCGAGCACTGGTGCGAACAATAACCGTCACACCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGATCGGCGCAGGACAGCACATTGCGCTTATGCGAATCCATCAATAATGCTTCCAATTTTCATTGTAAGCACAACCAGATGCACCCCCGa                                                            >  1:676747/1‑289 (MQ=255)                                                                                                                                                                                                                            ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGCAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGa                                                            >  2:1009717/1‑289 (MQ=255)                                                                                                                                                                                                                            ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTCTGCGAATCCATCAATAATGCTTCTCCTTTTCATTGTAACCACAACCAGATGCAACCCCGa                                                            >  2:953336/1‑289 (MQ=255)                                                                                                                                                                                                                            ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGa                                                            >  1:551434/1‑289 (MQ=255)                                                                                                                                                                                                                              cGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                                                                                                            <  2:217413/175‑1 (MQ=255)                                                                                                                                                                                                                              cGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                                                                                                            >  1:217413/1‑175 (MQ=255)                                                                                                                                                                                                                              cGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGt                                                          >  2:1204847/1‑289 (MQ=255)                                                                                                                                                                                                                              cGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGt                                                          >  2:1009207/1‑289 (MQ=255)                                                                                                                                                                                                                                          ggggggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAACCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                                                                                                                      <  2:1300733/199‑1 (MQ=255)                                                                                                                                                                                                                                           cgggtggattgcccaaaacaggCTTCCTTTTATTAAAATGGTTAAGGGCTTAAGGAAAAATCCAAAAAGGCGCCCAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTCAGGGTGGGTGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGt                                              <  1:627209/267‑1 (MQ=255)                                                                                                                                                                                                                                           cggggcgtagcccacaATTCAGTTTCCTTTTTTTAATCATGTTCAGGGAGTAATGACAAATGCGACAAAGCGCACAAGCCGTCCCCTCGCCCTTTTGGGGAGAGGGTTAGGGTTAGGGAAACAGCCCGCACTGGTGCGAACATTAACCCTCACCCCAGTCCTCACCCTTGAAGGGAGAGGGCGCAGAACGGCTCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGt                                              <  1:64938/276‑1 (MQ=255)                                                                                                                                                                                                                                           cGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                                                                                           <  2:717079/195‑1 (MQ=255)                                                                                                                                                                                                                                           cGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                                                                                                           >  1:717079/1‑195 (MQ=255)                                                                                                                                                                                                                                                   gTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGAGAGGGGAACAGCCAGCACTGGAGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaata                                                                                                        >  2:485942/1‑222 (MQ=255)                                                                                                                                                                                                                                                   gTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGG‑GTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAaga                                    >  2:356242/1‑289 (MQ=255)                                                                                                                                                                                                                                                    tctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCACACCCTGGCAGGGAGAGGGGGCAGAACGCCGCAGGACAGCACATTgcg                                                                                                                            >  1:1300733/1‑201 (MQ=255)                                                                                                                                                                                                                                                     ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAag                                     <  2:551434/287‑1 (MQ=255)                                                                                                                                                                                                                                                          aaaTCAGCTTCCTGTTATTAATAAGGTTAAGGG‑GTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACa                             <  1:356242/289‑1 (MQ=255)                                                                                                                                                                                                                                                                gCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                                                                                                                  >  2:1176501/1‑87 (MQ=255)                                                                                                                                                                                                                                                                gCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                                                                                                                  <  1:1176501/87‑1 (MQ=255)                                                                                                                                                                                                                                                                     cTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGt                   <  2:670368/289‑1 (MQ=255)                                                                                                                                                                                                                                                                     aagTTATTAATAAGGTTTAGGGCGTAAGGAAAAAGTCGAAAAAGCGAACAAGCCGCCCCCCCGGCCTGTTGGGGGTGGGGTTAGGGCGGGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGt                   <  1:1222479/287‑1 (MQ=255)                                                                                                                                                                                                                                                                                aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAgg        >  1:1365287/1‑289 (MQ=255)                                                                                                                                                                                                                                                                                  ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAg         >  2:1222479/1‑286 (MQ=255)                                                                                                                                                                                                                                                                                       aGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCTTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCg  >  1:501329/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                            GTCATCCGCGACGCGGTTATTCGGCGTGGTCGCGCCGCTGGCGATCACCGGAGCATCAATCGCCAGCAGTGAGCCGGTCAGGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGT‑GTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCA‑GTCAGCGGCCTGAAC‑TGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCG  >  NC_000913/624265‑624827

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