Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	433,634	0	T	G	32.4%	73.8 / 11.7	75	G60G (GGT→GGG)	ribD	fused diaminohydroxyphosphoribosylaminopyrimidine deaminase/5‑amino‑6‑(5‑phosphoribosylamino)uracil reductase
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base G (22/2);  total (49/26)
Fisher's exact test for biased strand distribution p-value = 6.84e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGC  >  NC_000913/433392‑433883                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                tACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGGGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑Gc                                                                                                                                                                                                                                                       >  1:457767/1‑251 (MQ=255)          ttCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTg                                                                                                                                                                                                                                              <  1:81909/250‑1 (MQ=255)           tCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑GAAAA‑GCC‑AAAgaagcg                                                                                                                                                                                                                                            >  2:403795/1‑246 (MQ=255)               aGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCg                                                                                                                                                                                                                                        <  2:521521/251‑1 (MQ=255)               aGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGGGCGACCg                                                                                                                                                                                                                                        >  1:172412/1‑251 (MQ=255)                          aaTTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCACCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTcac‑a                                                                                                                                                                                                                            <  1:45203/251‑1 (MQ=255)                                    agagACCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACACGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTCCCCCGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCt                                                                                                                                                                                                                  <  1:390180/251‑1 (MQ=255)                                          cgcgcgCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAACTTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGGGCGAACGCCTATGTCACCAC‑CGACCCCTGTAGcc                                                                                                                                                                                                            >  1:394413/1‑251 (MQ=255)                                           gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑ACAGGGGCGAACGCCTATGTCAC‑ACTCGAACCCAGTAGcc                                                                                                                                                                                                            >  1:324320/1‑250 (MQ=255)                                           gcgcgcCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑ACC‑AAAGGGGCGACCGCCTATTTCAC‑CCTCGAACCCTGTAGcc                                                                                                                                                                                                            >  1:539255/1‑250 (MQ=255)                                                      ggaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑CAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAAACCTTTAGCCCTC‑ATGGGCGt                                                                                                                                                                                                >  1:400471/1‑250 (MQ=255)                                                       gaCTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAACGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTTCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGAG‑AGAAA‑CAC‑CAAGGAAAGCCCGCCTATGTCCA‑ACTCGAAACCTGTAACCAAC‑ATGGTCGTa                                                                                                                                                                                               >  1:389527/1‑250 (MQ=255)                                                                        cGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCCGTAGCCATC‑ATGGTCGTACG‑CCACCGGGCTGTGACg                                                                                                                                                                            >  1:255832/1‑251 (MQ=255)                                                                         gagtatcacaaagcgagggCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCCTCCCAACCCGAATGTCGGGCGCGTCATTGTAAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGCGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCCCAAAGG‑GCAGCCCCCTATTTCCC‑ACTCCAACCCCGTAGTCATC‑ATGGCCGTAAG‑CCACAGTGCttgtacgc                                                                                                                                                                           >  2:371770/17‑244 (MQ=255)                                                                             aTTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTCTGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGGCGTACG‑CCACCGTGCTGTGACGCCCTc                                                                                                                                                                       >  1:3893/1‑251 (MQ=255)                                                                                  aTGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGGGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGATACCACCAACGTGCGGGTGAACCACATGCCGAAGAAAAAGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGGGCGGCCGCCTATATCAC‑ACTCGGACACTGGAGCCAAC‑ATGGACGTACG‑CAACAGAGCTGTGACgaaatcaacgc                                                                                                                                                                  >  1:46481/1‑241 (MQ=255)                                                                                    ggCGCGGGCGCTAAAGCTGGCGCAACGAGGACT‑TTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGTTTACCACCAACGTGCGGGTAAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATcgcc                                                                                                                                                                 <  1:194544/250‑1 (MQ=255)                                                                                         gggCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑ACC‑AAAGGTGCGACCGCCCATGTCAC‑ACTCGAACCCTGTAGCCATC‑CTGGTCGACCG‑CCACCGTGCTTTGACGCACTCATCGCCGCTggg                                                                                                                                                           >  1:519265/1‑250 (MQ=255)                                                                                         gggCGATAAAGCGGGCGCAGCGAGG‑CGCTTTTGCCACGCATCCCAACCCGAAGGTCGGTGGCGGCATTGAAAAAGAGGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAGCCCGAAGGACACGCGGTGCGTAGGGGGGGGG‑AAAAA‑GCC‑AAAGGTTCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGc                                                                                                                                                           <  2:351862/251‑1 (MQ=255)                                                                                           gcgcTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAAGTGCGACCGCCTATGTTCC‑ACTCCAACCCTTTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTTGCGCACCCATCCCCCCTGGCGt                                                                                                                                                         >  2:414359/1‑251 (MQ=255)                                                                                              cTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑ACAAA‑GAC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCCGTAGACATC‑CTGGGCGTACG‑CCACCCGGCTTTGACGCACTCacctccgccgcgggatg                                                                                                                                                      >  2:186640/1‑235 (MQ=255)                                                                                               tAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTTACGCACTCCTCGCCGCTGGCGTAgcg                                                                                                                                                     >  1:435499/1‑251 (MQ=255)                                                                                                          cAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑CAAGATGCGACCGCCTATGGCAC‑ACTCGAACCCTGTAACCATC‑ATGGTCGTACG‑CCACCGAGCTGTGACGCCATCATCGCCGCTGGCGTAGAGCGCGTGGTTGc                                                                                                                                          >  1:371891/1‑251 (MQ=255)                                                                                                           aaCGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGCACCCTGTAGCCATC‑ATGGTCGTACG‑CCACAGTGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGTTGcc                                                                                                                                         >  1:511009/1‑251 (MQ=255)                                                                                                               aGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑GAAAA‑GCC‑CAAGGTGCGACCCCCTATGTCAC‑ACCCGAAACCTGGAGCCCTC‑ATGGGCGGATC‑CCCCCGTGCTGGGACGCCCTCATCGCCGCTTGCGGAgcgcccgcgtgtgcctccg                                                                                                                                     >  2:152952/1‑236 (MQ=255)                                                                                                                  gactttttACCACGCATCCCAACCCGAATGCCGGGTGCGTCAGTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCTCGCGTGGTTGCCTCGATg                                                                                                                                   <  1:174927/247‑1 (MQ=255)                                                                                                                   cG‑TTTTACCACGCATCCCACCCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCa                                                                                                                                 <  1:209687/251‑1 (MQ=255)                                                                                                                          aCCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCCTAAG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGGGGTTGCCTCGATGCAAGATcc                                                                                                                           >  1:380383/1‑251 (MQ=255)                                                                                                                            cACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGAAGCACTCATCGCCGCTGGCGGAGCGCGCGTGGTTGCCTCGGTGCAAGATCCt                                                                                                                          >  1:390919/1‑250 (MQ=255)                                                                                                                                  tCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCCGTAGCCCTC‑CTGGGCGTACG‑CCCCCGTGGTGTGACGCACTCCTCTCCGCTGGCGGAGCGCGCGTGGTTGCCTCGATGCAAAATCCTACCCCGc                                                                                                                   >  2:310169/1‑251 (MQ=255)                                                                                                                                               gTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑ACAAA‑GCC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCCTC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGGTGCCCCGCTGCAAGATCCTAACCCGCAGGGCGCTGGGCg                                                                                                      >  2:95121/1‑251 (MQ=255)                                                                                                                                                  gggTGCGTCATTGTCAAAGGGGGCGAAATGGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAg                                                                                                                               <  1:307392/223‑1 (MQ=255)                                                                                                                                                   ggtgcgcaaTTGTCAAAGATGGCGAAATGGTCGGTGAAGGTTACCACCAAGGGGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGAATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGa                                                                                                  <  2:307927/243‑1 (MQ=255)                                                                                                                                                     tGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGGTGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACt                                                                                                >  1:115866/1‑251 (MQ=255)                                                                                                                                                        gTCATTATCAAAGATGGCGAAATTTGCGGTGAAGGTTACCACCACCGGGCGGGTAAACCACATGCCTAAGTACACGCGTTGGGTGTGGCGGGTC‑AAAAA‑GAC‑AAAGTTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTa                                                                                             <  1:102948/251‑1 (MQ=255)                                                                                                                                                          cATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGGGTGGACTTTAcc                                                                                           >  1:369152/1‑251 (MQ=255)                                                                                                                                                           aTTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACAACCAACGCGCTGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATA‑CTGGTCGTACG‑CCACCGTGCTGTGAAGCAATCATCCCCGCTGGCGCAGCGCGCGCGGTTTCCCAGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCg                                                                                          >  1:44241/1‑251 (MQ=255)                                                                                                                                                                aaaaGATGTCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACAGGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGc                                                                                     <  1:307887/250‑1 (MQ=255)                                                                                                                                                                     aTGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑ACC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATTGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTTGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAAc                                                                                  >  1:42022/1‑249 (MQ=255)                                                                                                                                                                        gCGAAATTGTCGGTGAAGGTGCCCACAAACGTGCGGGTGACCCACATGCCGAAGTACACGCGGTGCGTAGGGCGGGTG‑AAAAA‑GCC‑AAGGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCt                                                                             <  2:115866/251‑1 (MQ=255)                                                                                                                                                                         cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCCAAACCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCACAGCCGCTGGCGGAGCGCGCGTGGTGGCCTCGGTGCGAGGGCCTTAACCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTg                                                                            >  2:427877/1‑251 (MQ=255)                                                                                                                                                                         cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGAAGCACTCCTAGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTg                                                                            >  1:75397/1‑251 (MQ=255)                                                                                                                                                                         cGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAAACCC‑AAAGGT‑CGACCCCCTATGTCAC‑CCCCGCACCCCGTCGCCATCCAGGGCCG‑ACG‑CCACCCTGCTGGGACGCCCTCATGCCCGCCGGCCTAACCCGCGTCGTTGCCTCGATGCACGATCCTCACCCCCCGGTCCATGCGCGTGGACGTTTCCGGCTGCCACAGGCTg                                                                            >  2:218022/1‑251 (MQ=255)                                                                                                                                                                           caaTTGTCGGTGAAGGTTACCACCAACGTGCGGGAGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAAGGGAGACCGCCTATTTTCC‑ACTCGAACCCTGTTACCATC‑ATGGTTGTACG‑CCACCGGTCTGTGACGCACTCATCGCCTCTGGCGTAGCGAGAGTGTTAGCCTCGATGCCAGATCCTAACCCGCAGGGCGCTTGGCGTGGACTTTACCGttatacacaggctggg                                                                          >  1:474463/2‑236 (MQ=255)                                                                                                                                                                               tGTCGGTGAAGGATACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAAGTGCGCCCGCCCATGTCAC‑CCTCGCACCCTGTAGCCATC‑ATGGGCGTACG‑CCCCCGTGCTGTGCCGCACTCATCGCCGCTGGCGTCGCGCGCGTGGTTGCCTCGATGCAAGATACTAACCCGCCGGTCGCTGGGCGTAGACCTTACCGGCTGCAACCGGCTGGCATTg                                                                      >  2:269232/1‑251 (MQ=255)                                                                                                                                                                                gTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCCGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTTCCGTCCGCAACAAGCTGGCATTgt                                                                     >  2:193168/1‑250 (MQ=255)                                                                                                                                                                                    gTGAAGTTTACCCCCAACTGGCGGGTGAACCACAGGCGGAAGTACACGCGTTGCGTAGGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTc                                                                 <  1:6674/251‑1 (MQ=255)                                                                                                                                                                                     ggagggTTACCCCCAACGGGCGGGTGACCCACATGCCGAAGTACACGCGTTGCTTATGGCGGGTA‑AAAAA‑GCC‑AAGGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGGCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCa                                                                <  2:205585/247‑1 (MQ=255)                                                                                                                                                                                      gAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑ACC‑AAAAGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGAAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTc                                                                 >  1:20336/1‑249 (MQ=255)                                                                                                                                                                                        aGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCCGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGcc                                                             >  1:94264/1‑251 (MQ=255)                                                                                                                                                                                            tACCCCCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGGTCCTAACCCGCAGGTCGCTGGGCGTGGACt                                                                                                <  1:496855/212‑1 (MQ=255)                                                                                                                                                                                            tACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGCCCCCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGGGCTGTGGCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGGTCCTAACCCGCCGGTCGCTGGGCGTgggct                                                                                                >  2:496855/1‑209 (MQ=255)                                                                                                                                                                                            tACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGAGTATGGCGGGTG‑AAACA‑GCC‑AAAGCTGACACCGCATCTCTCAC‑ACTCGCACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCAATCTTCGCCGCTAGAGAAGCGCGCGTGGTTGCCTCGCTCCACCGTCCTCACCCGCAGGTCGCGGGGAGTGGACTTTACCGGCTGCAACAGACTGCGATTGACGACAGcccagg                                                         >  1:247309/1‑247 (MQ=255)                                                                                                                                                                                                    aCGTGCGGGTGACCCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTA‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCtgatga                                                 <  2:255832/251‑1 (MQ=255)                                                                                                                                                                                                    aCGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑CCTCCAACCCTGTAGCCCTA‑ATGGTCGTCCG‑CCCCCGGGCTGTGACGCACTCATCGCCGCTGGCGTATCGCGCGTGGGTGCCTCGATGCAAGATCCCAACCCCCCGGCCGCTGGGCGTGGAATTTACCCTCTGCAAACGGCTGGCATTGACGTCCGCCACGGCCcgatgg                                                 >  2:329220/1‑250 (MQ=255)                                                                                                                                                                                                      gTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGTGAGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACGACCAC‑GTGCTGTGACGCACTCATCGCCGCTGGCTTAGCGCGCGTGGTTGCCTCAATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTAACGTCTGCAACAGGCTGGCATTGACGTCAACCACGGCCtggtgatg                                               >  1:308178/1‑251 (MQ=255)                                                                                                                                                                                                          gggTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGATGGAAAAG‑GCC‑AAAGG‑GCGCCCCCCCACGTCAC‑CCTCGGACCCTGTAGACATC‑ATGGTCGGAAC‑CCACCGTGCTGGGACGCACTCATAACCCCTGGCGTAGCGCCCGTGGGTTCCCCGATGCCAGATCCTAAACCGCAGGGCCCTGGGCGTGGACGTTACCGTCTGCACCACGCGCGCATTGACggacgcccaggcctgttggtgcgtg                                           >  2:94771/1‑227 (MQ=255)                                                                                                                                                                                                           gggaaaCCACATGCCGAAGTACCCGCGTTGCGTATGGGGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGa                                          <  2:130957/247‑1 (MQ=255)                                                                                                                                                                                                           ggTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGa                                          <  1:403795/251‑1 (MQ=255)                                                                                                                                                                                                                  cacaTGCGGAAGTCCCCGCGTTGCGTATGGCGGGTA‑AAAAC‑GCC‑AATGGTGCTACCGCCTATGTCAC‑TCTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCACGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAg                                   <  2:135833/251‑1 (MQ=255)                                                                                                                                                                                                                    caTGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑ACC‑AAAGGTGCGACCGCCTATTTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCCGTGCCAGATCCTAACCCGCAGGGCGCTGGGCGGGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCCCGGCCTGATGATGAGTGAAGCCgggcc                                 >  2:122934/1‑247 (MQ=255)                                                                                                                                                                                                                       ccggatgcccccgcattGCTTATGGCGGGTG‑AAAGA‑GCC‑AAAGGTGCGACCGCATATGGCGC‑TCTCGAACCCTGTAGTCATC‑ATGGACGTACG‑CCACCGTGCGGTGACGCACTCATCGCCGCTTGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTGCCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAAtt                              <  1:535742/236‑1 (MQ=255)                                                                                                                                                                                                                           aaGTACACGGGTTGCGTATGGCGGGTG‑AAAAA‑GCA‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAAt                          <  2:3893/251‑1 (MQ=255)                                                                                                                                                                                                                           aaGTACACGCGTTGCGTATGGCGGGGG‑AAAAA‑GCC‑AAAAGGGCGCCCGCCTATCTCAC‑ACCCGAACCCTGTAGCCCTC‑CTGGTCGTTCG‑CCACCGTGCTGTGACGCACTCATTGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCCAACCCGCAGGGCGCCGGGCGTGGACTTTACCCTATGCAACAGGCTGGGATTGACGTCAGCCCCGGCCTGATGTGGAGGGGAGCCGAGCCATTGAAt                          >  2:520951/1‑251 (MQ=255)                                                                                                                                                                                                                            agaacacGCGTTGCGTATGGCGGGTA‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGGCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGGAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGGCAGCCAccgc                                                        <  1:359259/217‑4 (MQ=255)                                                                                                                                                                                                                                  cgcgTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGACATC‑ATGGTCGTACG‑CCCCCGTGCTGTGCCGCACTCATCGCCGCTGGCGTAGCGCGCGTGGGTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCt                   >  1:341766/1‑251 (MQ=255)                                                                                                                                                                                                                                   gtgTTGCGTAGGGGGGGTG‑AAAAA‑CCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCtt                  <  2:270433/249‑1 (MQ=255)                                                                                                                                                                                                                                    ggTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGGGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCttt                 <  2:61296/250‑1 (MQ=255)                                                                                                                                                                                                                                       gggggagggggggTG‑AAAAA‑GCC‑AAAGGTGCGCCCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTtct               <  2:525153/240‑1 (MQ=255)                                                                                                                                                                                                                                          gTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCa             >  2:118555/1‑249 (MQ=255)                                                                                                                                                                                                                                             tGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCCTTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGCGCAATTGAATAAATGCTTTCTCAAGCgg        >  2:3494/1‑251 (MQ=255)                                                                                                                                                                                                                                              ggCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGGGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCgg        >  1:268688/1‑250 (MQ=255)                                                                                                                                                                                                                                                cGGGTG‑AAAAA‑GCC‑AAAGGTGCGACAGCCTATGTCAC‑CCTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGGGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATg     >  2:544224/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGCGTGAAGCCGAGCAATTGCATAAAGGCTTTCTCAAGCGGATgcgc  >  1:371695/1‑251 (MQ=255)                                                                                                                                                                                                                                                   gTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATgcgc  >  1:426633/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                TACCGCAGTTTCGAAGATATCAAAGAATTTGGCGAAGAGATCGCGCGCCTGGAGGACTAAGCCGTGCAGGACGAGTATTACATGGCGCGGGCGCTAAAGCTGGCGCAACGAGGACG‑TTTTACCACGCATCCCAACCCGAATGTCGGGTGCGTCATTGTCAAAGATGGCGAAATTGTCGGTGAAGGTTACCACCAACGTGCGGGTGAACCACATGCCGAAGTACACGCGTTGCGTATGGCGGGTG‑AAAAA‑GCC‑AAAGGTGCGACCGCCTATGTCAC‑ACTCGAACCCTGTAGCCATC‑ATGGTCGTACG‑CCACCGTGCTGTGACGCACTCATCGCCGCTGGCGTAGCGCGCGTGGTTGCCTCGATGCAAGATCCTAACCCGCAGGTCGCTGGGCGTGGACTTTACCGTCTGCAACAGGCTGGCATTGACGTCAGCCACGGCCTGATGATGAGTGAAGCCGAGCAATTGAATAAAGGCTTTCTCAAGCGGATGCGC  >  NC_000913/433392‑433883

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