Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	433,649	T→G	25.8%	G65G (GGT→GGG)	ribD →	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,649	0	T	G	25.8%	95.5 / 16.6	93	G65G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase and 5‑amino‑6‑(5‑phosphoribosylamino) uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (19/46);  new base G (20/5);  total (42/51)
Fisher's exact test for biased strand distribution p-value = 2.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.79e-01

ATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCAC  >  NC_000913/433418‑433885                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               aTTGGGCGAAGAGATCGCGCGCCTGGAGGACTAGGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑Ac                                                                                                                                                                                                                                         <  1:874513/237‑1 (MQ=255)    tGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAAGGGC‑G‑AACGccc                                                                                                                                                                                                                                    >  1:930611/1‑238 (MQ=255)        gaaagaaTCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTGTTACATGGCGCGGGCGCTAAGGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGt                                                                                                                                                                                                                                <  1:468082/233‑1 (MQ=255)                cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAAGTGC‑G‑ACCGCCTATGTCCCACTCg                                                                                                                                                                                                                        >  1:506017/1‑239 (MQ=255)                 gcgcgcCTGGAGGACTAAGCCGTGCAGGACGGTTATTACAGGGCGCGGGCGCTAAAGCGGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGGCGGGGGCGTCATTGGCAAAGATGGCGAAATTGTCGGGGAAGGTTACCACCAA‑CGTGCGGGGGAACCACATGCCGAAGTACACGCGGTGCTTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGa                                                                                                                                                                                                                       <  2:770723/239‑1 (MQ=255)                 gccgggcgggaggaCAAAGCCGGGCAGGACAAGTTTTCACGGGCGCGGGCGCTAAAGCGGGCGCAACGAGGACGTTTTACCACGCACCCCAACCC‑GAATGTCGGGTGCGTCATTGTAAAAGATGGCAAAATTGTGGGTAAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGa                                                                                                                                                                                                                       <  2:506017/231‑1 (MQ=255)                      ccTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACCCGAACCCt                                                                                                                                                                                                                  >  1:663408/1‑239 (MQ=255)                        tGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGGGC‑G‑AACGCCTATGTTACACTCGGACCCTGt                                                                                                                                                                                                                >  1:307970/1‑239 (MQ=255)                         ggaggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                               <  1:728583/239‑1 (MQ=255)                          gaggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTa                                                                                                                                                                                                               <  1:1165723/238‑1 (MQ=255)                               ctagtcctggCAGGACGAGTTTTCAAGGCCGGGGGCGCAAAAGCTGGCGCAACGAGGACTTTTTACCACGCACCCCAACCG‑GAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAC‑CGTGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTTTGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCcat                                                                                                                                                                                                          <  2:223663/229‑1 (MQ=255)                                gaggCCGTGCGGGACGATTAATACAGGGCGCGGGCTGTAAAGCTGGCGCAACGAGGAGATTTTACCACGCATCCCACCCC‑GAATGGCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGGAAAGGTTCCCACCAA‑CGTGGGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCcatca                                                                                                                                                                                                        <  2:113017/236‑1 (MQ=255)                                     cGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATgg                                                                                                                                                                                                     <  1:1066040/237‑1 (MQ=255)                                       tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAA‑CCAAAGGGGC‑G‑GCCGCCCATGTCACACCCGCACCCCGTAGCCCTCAAGGTCGt                                                                                                                                                                                                 >  2:836279/1‑239 (MQ=255)                                       tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAA‑ACAAAGGGGC‑G‑ACCCCCTACGACACCCTCGAACCCGGTAAccgccactgtgc                                                                                                                                                                                                  >  2:149473/1‑228 (MQ=255)                                       tGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATTTCACACTCGAACCCTGTAGCCATCATGGTCGt                                                                                                                                                                                                 >  1:492185/1‑239 (MQ=255)                                              cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCACCCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTCCACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                           <  1:866675/238‑1 (MQ=255)                                              cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAA‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCa                                                                                                                                                                                           >  1:362530/1‑238 (MQ=255)                                               gAGTATTACGTGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTTTCGGTGAAGGTTACCACCAA‑CGGGCGGGTGAACCACATGCCGAAGTACCCGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCAcc                                                                                                                                                                                         <  2:239354/239‑1 (MQ=255)                                                 gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGGGAAAAAG‑CC‑AAGGGGC‑G‑GCCCCCGCCGTCCCACACCTACCCTGTAGCCACCACGGGCGTACGCCAcac                                                                                                                                                                                        >  2:1101172/1‑236 (MQ=255)                                                 gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGt                                                                                                                                                                                       >  1:838831/1‑239 (MQ=255)                                                 gTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGC‑G‑GCCGCCTATGTCCCACTCGAACCCTGTAGCCCACATGGTCGTACGccca                                                                                                                                                                                          >  1:888673/1‑234 (MQ=255)                                                  tATTACATGGCGGGGGCTCTGAAGCTGGGGCAACGAGGCCGTTTCGCCGCGCAGCCCACCCC‑GAAGTGCGGGTGCGACATTTTGAAAAAGGCCGAAATGTTCGGTGAGGTGTCCCACCAC‑CGGGCGGGTGAACCACATGCCGAAGTACCCGCTTTGCTTATGGC‑GGGGAAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCGCCGTg                                                                                                                                                                                      <  2:1135782/239‑1 (MQ=255)                                                      aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTAAGCCACCGGGCtgt                                                                                                                                                                                  >  1:799548/1‑239 (MQ=255)                                                      aCATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACAACCGAACCCTGTAGCCATCATGGGCGTACGCCACCGTGCtgt                                                                                                                                                                                  >  1:1183365/1‑239 (MQ=255)                                                       cATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCCCCGGGCtgtg                                                                                                                                                                                 >  1:1064158/1‑239 (MQ=255)                                                       cATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGCGGGAAAG‑ACAAAGGGGC‑G‑AACACCCCTGTCCCACTCGGACCCTGTAGCCACCAGGGGCggcagccacacggcttgg                                                                                                                                                                                 >  2:419630/1‑222 (MQ=255)                                                         tGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAA‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCCGTAGCCATCATGGTCGTACGCCCCCGGGCCGTGAc                                                                                                                                                                               >  1:943661/1‑239 (MQ=255)                                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAG‑CCAAAAGCGC‑G‑GACCCCCATGTTCCCATCGCACCCTGTAACCAACCTGGGCGTACGCCCCCGGGCTGAGACGc                                                                                                                                                                             >  2:199789/1‑238 (MQ=255)                                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAG‑GCCCAGGGGC‑G‑AGCGCCCATCTCACACAATCACCCAGTAGCCATCCTGGACGTAAGCccccctccgttggccc                                                                                                                                                                             >  2:1020116/1‑226 (MQ=255)                                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAG‑CCAAAAGGGC‑GAACCGCC‑ACGTCTCCCTCCAAACCCTCAGCCCTACTGGGGGTACGGCACCGTGCTCGGACGc                                                                                                                                                                             >  2:190287/1‑238 (MQ=255)                                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGAGAAAAAG‑CCAAAAGGGC‑G‑ACCGCCTATGTAACACTCGAACCCTGTAGCCATCCTAGTTGTACCCCACCCTGCTGTGagcg                                                                                                                                                                             >  1:30810/1‑235 (MQ=255)                                                            cgcgGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGCGCGGGTGAACCACATGCCGAAGTACACGAGTTGCGTATGGC‑GGGTGGCAAAG‑GCAAAAGTGG‑T‑ACCACCTATTGAACAATCCGACCCAGTAGCCCTCATGGGAGCACGCCAcccggcttgggagc                                                                                                                                                                             >  2:383864/1‑226 (MQ=255)                                                                 gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAG‑CCAAAGGGGG‑G‑ACCGCCCATTTCCCCCCCGACCCCTGTTGCCCTCATGGGTGTGACCCGCCCGGCTGTGACGCACtac                                                                                                                                                                        >  2:375549/1‑236 (MQ=255)                                                                 gcgcTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCCGTAACCATCATGGTCGTAAGCCCCCGTGCTGTGACGCAcccc                                                                                                                                                                        >  1:686502/1‑235 (MQ=255)                                                                                               aaccgggcaaccAACCCGGAAGGTGGGGTGCGTAATTGCAAAAGAGGGCGAAATGGGGGGTGAAGGTTACCCCCAA‑CGGGCGGAGGAACCACAGGCCGAAGTACACGCGTGGCGTAGGGC‑GGGGGAAAAAG‑CAAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCt                                                                                                                                          <  2:486480/229‑1 (MQ=255)                                                                                               aCCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAG‑CCAAAAGGGC‑G‑ACCGCCTATGTTACACTCGAACCCTGTAGCCCTCATGGTCGTAAGCCCCCGCGCTGTGCCCCACCCATCGCCCCCGGCCTAGCGCGCGCGGCTGccc                                                                                                                                          >  2:1218492/1‑237 (MQ=255)                                                                                                ccACGCATCCCACCCC‑GAATGTCGGGTGCGTCATTGTCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTGGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCg                                                                                                                                        <  1:168112/239‑1 (MQ=255)                                                                                                   cGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCCCGGGCGTAGGGCGCGCGGTTGCCTCgcg                                                                                                                                      >  1:867633/1‑236 (MQ=255)                                                                                                     tagcacAACCA‑GAAGGGCGGGGGCGTCATTGTAAAGGAGGGCAACATTGGCGG‑GAAGGGTTCCCCTAAGCGGGGGGGTAAACCACAGGCCGAAGTACCCGCTTTGGGTGTGGC‑GGGCAAAAACG‑CAAATGGGGC‑G‑GCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  2:960454/234‑1 (MQ=255)                                                                                                      atcccatccc‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATGGCCGGTGAAGGTTCCCCCCAA‑CGGGGGGGGGAACCACATGCCGAAGTCCCCGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  1:149473/238‑1 (MQ=255)                                                                                                      aTCCCAACCC‑GAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  1:221561/238‑1 (MQ=255)                                                                                                      aTCCCAACCC‑GAATGGCGGGTGCGTCATTGACAAAGATGGCGAAATTGTCGGTGAAGGTTCCCACCAA‑CGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAC‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                   <  1:865485/238‑1 (MQ=255)                                                                                                                 gAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCCCTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCCATGCAAGATCCTAAc                                                                                                                         >  1:829843/1‑238 (MQ=255)                                                                                                                   aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGA‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCCTCATGGTCGTACGCCACCGGGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGCGGGTGCCTCGGTGCCAGATCCCAAccc                                                                                                                       >  1:356264/1‑238 (MQ=255)                                                                                                                   aTGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGGAAAAG‑CCCAAAGTGC‑G‑ACCGCCTATTTGCCACTCGAACCCCGTTGCCATCATGGTCGGACGCCCCCGTGCTGTGACTCCCCCATCCCCCCGGGCGTAGCCCGCGCGGGTGTCCCGagtccaaaccttcaccc                                                                                                                       >  2:385548/1‑222 (MQ=255)                                                                                                                         ggTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCAAAGGTGC‑G‑CGCGCCTATTGCACACTCCAACCCCGGAGGCATCACGGTTGGACGCCACCGCCGCGAGTCGCACCCATCGCCCCTGCCGTCGCGCGCGCGGCTTCCCCGACGGAACAACCTAAccccccgct                                                                                                                 >  2:1059824/1‑232 (MQ=255)                                                                                                                            gCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAC‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                              <  1:971098/238‑1 (MQ=255)                                                                                                                            gCGTCATTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTACCACCAC‑CTTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGc                                                                                                              <  1:1020116/238‑1 (MQ=255)                                                                                                                                cATTGTAAAAGTAGGCTAAATGGACGGTGAGCGTGCCCGCAAC‑CGTGGGGGTGAACCAAATCCGGAAGTACGCGGGTTGGGTAGGGC‑GGGTTAAAAGC‑CCAAGGGGTC‑C‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                         <  2:861966/239‑1 (MQ=255)                                                                                                                                cATTGGAAAGGAGGGCGAAATTGTGGGTGAAGGTTACCCCCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGc                                                                                                         <  1:940914/239‑1 (MQ=255)                                                                                                                                  ttGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCTTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGt                                                                                                       <  1:987729/239‑1 (MQ=255)                                                                                                                                     tCAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAgcgcg                                                                                                                                                     <  1:230321/190‑1 (MQ=255)                                                                                                                                          cgtagcgaatgtgtgtggaaaGGATCCCCCCAG‑CGTGGGGGTGACCCCCACGCCGGAGTACACTCTTGGCG‑GTGGG‑GGGTGGAAAAGCCCAAAGGTGC‑G‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                               <  2:324680/220‑1 (MQ=255)                                                                                                                                            tGGCGAAATTGTGGGTGAAGGTTCCCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGG‑GGGTGAAAAAC‑CCAAGGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                               <  1:1218492/237‑1 (MQ=255)                                                                                                                                            tGGCGAAATGGGCGTTGAGGGTTACCACCAG‑CGTGAGGGGGAACCACATGCCAAAGTACACTGGTTGCTAAGGGC‑GGGTGAAAAAC‑CCAAAGGTCC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTGCGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                               <  2:356264/237‑1 (MQ=255)                                                                                                                                             ggCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGCGGTTGCCTCGATGCCAGATCCTAACCCGCAGGTCGCTGGGCGGGGACTTTAcc                                                                                             >  1:850421/1‑238 (MQ=255)                                                                                                                                               caaaaGTGTCGTTGAAGGTTACCCCCAA‑CGTGCGGGTGAACCACATGCCGAAGTACCCGTGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTAGGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGt                                                                                           <  1:1154660/236‑1 (MQ=255)                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGGAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTTACAATCGAACCCCGTCGCCATCACGGTCGTACGCCACCCTGCTGTGGCGCACTCATCGCCCCTGGCGTAGCGCGCGGGGGTGcc                                                                                                                                           >  2:230321/1‑190 (MQ=255)                                                                                                                                               cGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCAGAAGTACACGCGTTGCGTATGGC‑GGGGGAAAAAG‑CCCAAGGTGC‑G‑ACCCCCTCTGTCACCCTCGAAACCTGTAGCCATCAAGGTCGGACACCCCCCTGCTCTGACCCACACACCGCCCCGCGCGTCACGCCCGCGGTTGCCCCGGGGCAAGCCCCTAACCCGCCGGGCGCTGGGCggcgcgttttccgg                                                                                           >  2:40695/1‑225 (MQ=255)                                                                                                                                                       tCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCCGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAg                                                                                  >  1:318523/1‑239 (MQ=255)                                                                                                                                                         ggTGAAGGTTACCCCCAC‑CGTGCGGGTGACACCCAGGGCGAGTTACACGCGTGCCTTATGGC‑GGGTAAAAAAC‑CAAAGGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAgg                                                                                 <  2:214014/238‑1 (MQ=255)                                                                                                                                                                       ccaA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCa                                                                  >  1:383177/1‑239 (MQ=255)                                                                                                                                                                            cGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGG‑GGGTGAAAAAT‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                              <  1:1062616/239‑1 (MQ=255)                                                                                                                                                                            cGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                              <  1:1109934/239‑1 (MQ=255)                                                                                                                                                                             gggggggTGAACCAAAGGCCGAAGTACACGCTTTGCGTATGGC‑GGTGAAAAAAC‑CCAAAGGTCC‑G‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                              <  2:492185/234‑1 (MQ=255)                                                                                                                                                                             gggcgtgggaacccactGCGGAAGTACCCGTGTTGCGTACGGC‑GGGTAAAAAAC‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                              <  2:444818/222‑1 (MQ=255)                                                                                                                                                                             gTGCGGGTGACCCACATGCCAAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAGGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCa                                                              <  1:287940/238‑1 (MQ=255)                                                                                                                                                                               gCGGGTGAACCAAATGCCGAATTACCCGCGTTGCGAGTGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACg                                                            <  1:375549/238‑1 (MQ=255)                                                                                                                                                                                 gtgTGACCCACGTGCCGAGGTACCCGGGTCGCGTATGGC‑GGGTAAAAACG‑CCAAGGGTGC‑G‑GCCGCCTAGGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                                          <  1:1045159/236‑1 (MQ=255)                                                                                                                                                                                 gggTCAACCAAAGGCCGAAGTACACGCTTGGCGTATGGC‑GGGTGAAAAAG‑CAAAAGGTCG‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGc                                                          <  2:864458/238‑1 (MQ=255)                                                                                                                                                                                      acccACAGGCCGAAGTACACGCGTTGCGTATGGC‑GGGTAAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatg                                                    <  1:419630/237‑1 (MQ=255)                                                                                                                                                                                             tGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCCCCCTGATGTGACGCACCCATCGCCGCTGGCGTAGCGCGCGCGGTTGCCTCGGTGCCAGATCCTAAACCGCCGGGCGCTGGGCGTGGACTTTACCGTCTGCACCCGGCTGCCATTGGCGGCAGCCACGGCCTGCTGAtgtgtg                                             >  2:518765/1‑239 (MQ=255)                                                                                                                                                                                             tGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGGCGTACGCCCCCGTGCTGTGACGCACTCATCTCCGCTGGCGTAGCGCGCGCGGGTGCCTCGATGGCAGATCCTCACCCGCCGGCCGCTGGGGGGGGGCTTTTCCGTGTGCCACAGGCTGGCCTTGGAGTCGTCCACGGCCTGTTTATGAGTg                                             >  2:714232/1‑239 (MQ=255)                                                                                                                                                                                                cgaggtccccctgTTGCGTATGGC‑GGGTGAAAAAG‑CCAAGGGTGC‑G‑CCCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAg                                          <  2:903133/227‑1 (MQ=255)                                                                                                                                                                                                  aaGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGcc                                        >  1:903133/1‑239 (MQ=255)                                                                                                                                                                                                     aacagGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGACCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                     <  1:539408/238‑1 (MQ=255)                                                                                                                                                                                                      ccccgcgTGGCGTAGGGC‑GGGGAAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                     <  1:1142795/235‑1 (MQ=255)                                                                                                                                                                                                      acccggtctggtgaTGGC‑GGGTGAAAACG‑CCAACGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                     <  2:838831/225‑1 (MQ=255)                                                                                                                                                                                                      acacGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                     <  1:152593/238‑1 (MQ=255)                                                                                                                                                                                                        ccgtgTTGCTTATGGC‑GGGTGAAAAGT‑CAAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCaa                                  <  1:518765/235‑1 (MQ=255)                                                                                                                                                                                                         cgcgTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAt                                 <  1:1098218/239‑1 (MQ=255)                                                                                                                                                                                                             ttGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGGGGAATTTACCGGCTGGCACAGGCTGGCATTTGCGTCAGCCACGGCCTGATGATGAGGGAAGCCGAGCACTTGaa                             >  1:661555/1‑239 (MQ=255)                                                                                                                                                                                                              tGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGAGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAt                            >  1:809890/1‑239 (MQ=255)                                                                                                                                                                                                                    tGGC‑GGGTGAAAAAG‑CCAAAGGCGC‑G‑ACCGCCTATGTCCCCCTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGGCGCCCCCATCGCCGCTGGCGTACCGCGCGTGGTTGCCTCGCTGCAAGATCCTAACCCCCCGGGCTCTGGGCGCGGGCTTTCCCTTCTGCACCAGGCTGGCATTGGCGTCTGCCACGGGCTGCCGCTGACTGCAGCCGAGCAATTGAATAAAgg                       >  2:808838/1‑238 (MQ=255)                                                                                                                                                                                                                     ggC‑GGGGGAAAAAG‑CCAAAGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGGACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGTTGCCCCGATGCAAAATCCTAACCCCCAGGTCGCTGGGGGGGGAGTTTACCGGCTGCAACAGGACGGCCATGACGGCAGCCACGGCCTGTTGATGAGGGAAGACGGGGAAATGAATaactgt                      >  2:539408/1‑234 (MQ=255)                                                                                                                                                                                                                           gTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtct                 >  1:407916/1‑238 (MQ=255)                                                                                                                                                                                                                              aaaaaG‑CAAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg             <  1:190287/239‑1 (MQ=255)                                                                                                                                                                                                                               aaaaG‑CCAAAGGTGC‑G‑ACCGCCTATGTCCCACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg             <  1:1191877/238‑1 (MQ=255)                                                                                                                                                                                                                               aaaaG‑CAAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAg             <  1:793031/238‑1 (MQ=255)                                                                                                                                                                                                                                  aG‑CCCAAGGGGCGG‑ACCGCCTGAGTCTCCCTCGAACCCTGTGGCCATCATGGTCGTACGCCACCGTTCTGTGACGCAGTCATCGCCGCTGGAGTAGGGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAGCAGGCTGGCATTGACGTCAGCCGCGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCgg          <  2:712362/239‑1 (MQ=255)                                                                                                                                                                                                                                        aaGGGGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgcgc    <  1:1101172/239‑1 (MQ=255)                                                                                                                                                                                                                                           gTGG‑G‑ACCGCCTATGTCACACTCGACCCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCAc  <  2:686502/238‑1 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               ATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACGTTTTACCACGCATCCCAACCC‑GAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAA‑CGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGC‑GGGTGAAAAAG‑CCAAAGGTGC‑G‑ACCGCCTATGTCACACTCGAACCCTGTAGCCATCATGGTCGTACGCCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGCAC  >  NC_000913/433418‑433885

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