Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,098,391	0	A	C	29.3%	103.9 / 10.9	78	Y219S (TAC→TCC)	rhaS	DNA‑binding transcriptional activator RhaS
Reads supporting (aligned to +/- strand):  ref base A (37/16);  new base C (1/21);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTCTGGCCGAGGT  >  NC_000913/4098150‑4098638                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        gATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCCGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTTTCTTTCACTGCGTACGCT‑ACCTCGGCA‑GCTTAAGCAGCAAACGGGGCTGAC‑GCCTCAGCGATACCTGAACCg                                                                                                                                                                                                                                                >  2:432959/1‑251 (MQ=255) gATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAAcc                                                                                                                                                                                                                                                 <  1:429220/250‑1 (MQ=255)     tACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTGGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAACTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCt                                                                                                                                                                                                                                             <  1:193450/250‑1 (MQ=255)            gACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACtg                                                                                                                                                                                                                                      >  1:264875/1‑250 (MQ=255)               cGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAAACGCCTGCGACtgatg                                                                                                                                                                                                                                   >  1:599092/1‑250 (MQ=255)                                     aTGCAATTACCGCTCTCGCTGCGAAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCTCATCACGTCGCAACTTGCTGCTGGCCTGGCGGGAGGACCTTTTT‑GCCGATGAGGTGAATTGGGATGCCGGGGCGGATCAATTTTCTCTTTCCCTGCGTACGCT‑ACAGCGGCC‑GCTTAAGCAGCAAACGGGACTGCC‑GCCTCAGCGATCCCGAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGcc                                                                                                                                                                                                             <  2:568045/250‑1 (MQ=255)                                        cAATTACTGCTCTTGCTGCGTCAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTATGGCCTGGCTGGCGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCCGCGATACCTGAACCGCCCGCGACTGATGAAAGCCCGACATCTGCTACGCcac                                                                                                                                                                                                           >  2:686201/1‑249 (MQ=255)                                                     tgcgtcGTAAAAGCAGTTGGCAGGAGAACCTGGAAAACAGCGACTCCCTTCTCAACTTGCTTCTGGCCTGTCTGGAGGACCATTTT‑GCCGATGAGGGGAATTGGGAGGCGGTGGCGGATCAATTTTCTCTTTCCCTGCGTCCGCT‑CTTTCGGCA‑GCTTAAGCAGCAAACGGGGCTGAC‑GCCTGAGCGGTACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGc                                                                                                                                                                                               <  1:360922/243‑1 (MQ=255)                                                          cGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTg                                                                                                                                                                                        >  1:393288/1‑250 (MQ=255)                                                           gTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGAACATTTT‑GCCGCTGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGCTACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGAc                                                                                                                                                                                      >  1:657090/1‑251 (MQ=255)                                                             aaaaGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATTCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGCGGCCAGCGTTACTGACac                                                                                                                                                                                    >  1:230385/1‑250 (MQ=255)                                                             aaaaGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCC‑ACATCGGCA‑GCTTAAGAAACAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCCGCGACTGATGAAAACCCGCCATCTGCTACGCCCCAGCGAGGCCAGCGTTACTGAcct                                                                                                                                                                                    >  2:519943/1‑249 (MQ=255)                                                                    gTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAACGTTACTGACATCCCCTAt                                                                                                                                                                             >  2:514292/1‑251 (MQ=255)                                                                      ttGCAGGAGAACCTGGAAAACAGCGCATCACGTCCCAACTGGCTTCGGGCCTGGCTGGAGGCCCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACAGCGGAA‑GCTTAAGCAGCAAACGGGACTGAC‑CCCTCAGCGATACCGGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATc                                                                                                                                                                            <  1:15264/250‑1 (MQ=255)                                                                          agtagacccTGAAAACCAGCGCATCACGTCTCAACTTGCTTCGTCCCTGGCTGGAGGCCCATTTT‑CCCGATGAGGTGATTTGGATTGCCGTGGCGAATCATTTTTCCCTTTCCCTGCGTCCGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGAACTGCC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCTCCTATCGCtgt                                                                                                                                                                       <  1:532650/244‑1 (MQ=255)                                                                           ggAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑CCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGTTGAAAGCCCGACATCTGCTACGCCCCCGCGAGGCCCGCGTTACTGACACCGCCCATCGCtttg                                                                                                                                                                      >  2:556296/1‑248 (MQ=255)                                                                            gagaACCGGGAAAACAGCGCATCACGTCGCAACTTGCTTCGGGCCGGGCTGGAGGACCATTTT‑GCCGAGGAGGTGAATTGGGATGCCGTGGCGGATAAATTTTCTCTTTCACTGCGTACTCT‑ACTTGGGCA‑TTTTAAGCAGCAAACGGGACTGCC‑TCCTCAGCGTTACGGGAGCCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTGCTGACATCGCCTATCGCTGTgg                                                                                                                                                                     <  1:620918/251‑1 (MQ=255)                                                                               aaCCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGAtt                                                                                                                                                                  >  2:157350/1‑251 (MQ=255)                                                                                cccTGGAAAACAGCGCATCACTTCTCAACTTCCTTCTGGCCTGGCTGGAGCACCATTTG‑GCCGATGAGGGGAATTGGGTTGCCGTGGCGTTTCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGCC‑GCCTCAGCAATTCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTc                                                                                                                                                                 <  1:359245/250‑1 (MQ=255)                                                                                 ccTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCCCAGCGATACCTGAACCGCCTGCGACTTATGAAAGCCCGACATCTGCTACGCCACAACGAGGCCAGCGTTACTGACATCGCCCATCGCTGTGGATTcc                                                                                                                                                                >  2:502374/1‑250 (MQ=255)                                                                                 ccTGGAAAACAGCCCTCCCCGTCTCATCTGGCTTCTGGCCGGGCGGGAGG‑CCATTTTTGCCGATGGGGTGATTTGGGATGCCGTGGCGAAACAATTTTCTCTTCCGCTGGGTGCGTT‑ACATCGGCA‑GCTTAGGCAGCAAAGGGGACTGAC‑GCCTCGGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCa                                                                                                                                                                <  2:537076/251‑1 (MQ=255)                                                                                      aaaaCAGCCCATCACGTCTCAACTTGTTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGAATCAATTTTCTCTTTCACGGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATTCCGTAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGATGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGAc                                                                                                                                                           <  1:432959/251‑1 (MQ=255)                                                                                              gcaccACGTCTCAACTGGCTTCTGGCCTGGCGGGAGGACCATTTT‑GCCGATGGGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCCCTGCGTACGCT‑ACACCGGCA‑GCTTAAGCAGCAAAGGGGCCGGAC‑GCCTCAGCGATACCTAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCa                                                                                                                                                   <  2:139987/247‑1 (MQ=255)                                                                                                  cacttctcAACTTGCTTCGGGCCGGTCGGGAGGACCATTTT‑CGCGATGAGGTGAATTGGGATGCCGTGGCTGATAAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGAA‑GTTTAAGCAGCAAACGGAACTGAT‑GCCTCAGCTATCCCTGAACCGCCTGCGACGGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACttt                                                                                                                                               <  2:601083/247‑1 (MQ=255)                                                                                                   aCGTCTCAACTTGCTTCTGGCCTGGATGGAGGACCATCAT‑GCCGATGCGGAGAATTGGGATGCCGCGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCC‑GCTTACGCAGACAACGGGCATGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAACCCGACATCTGCTACGCCACAGCGAAGCCAGCGTTACTGCCATCGCCTAACGCTGTGGATTCAGCGACAGTAACCACtttt                                                                                                                                              >  1:282550/1‑251 (MQ=255)                                                                                                    cGTCTCAACTTGCTTCTGTCCTGGCTGGAGGACCATTTG‑GCCGAGGAGGTGAATTGGGATGCCGTGGCGGATTAATGGTCTCTTTCGCTGCTTACGCT‑ACAGCGGCA‑GCTTAAGCAGGAAGCGGTTCGGGC‑TCCTCAGCGATCCCTAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCTAGTGTGACTTATAGCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTc                                                                                                                                             <  1:606711/251‑1 (MQ=255)                                                                                                       ctcAACTTGCTCCTGGCCTGCGTGGAGGACCATGTT‑GCCGTTGGGGTGAATTGGGATGCGGGGGCGGATCAATTTTCTCTTTCACTGCTTTCGCT‑ACTTCGGCA‑GCTAAAGCAGCAAAGGGGACTGAC‑GCTTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACTCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCg                                                                                                                                            <  1:438028/249‑1 (MQ=255)                                                                                                             ttGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑AAATCGGCA‑GCTTAAGCCGCAAACGGGACTGAC‑GCCGCAGCGATACCTGAACCCCCTGCGAATGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAAACACTTTTCGACGCtttt                                                                                                                                    >  1:500276/1‑251 (MQ=255)                                                                                                              ttcTTCTGGCCGGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATGGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGAAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCttt                                                                                                                                     <  1:319297/247‑1 (MQ=255)                                                                                                               gctcctggcctggcTGGAGGCCCTTTTT‑GCCGATGAGGTAAATTGGAATGCCGTGGCGAATAAATTTTCTCTTTCCCGGCTTACGCT‑ACATCGGCA‑GCTTAACCAGCAAACGGGACGGAC‑GCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTc                                                                                                                                  <  1:145228/247‑1 (MQ=255)                                                                                                                         gtggttggaggcccttttgtccGATGAGTTGAATTGGGATGCCGGGGCGGAGAAATTTTCTTTTTCGCTGCGTAGGCG‑ACATCG‑CATGCTAAAGGAGAAACCGGGCCTGCC‑GCCTCAGCGATCCCTAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACGTTTCGACGCTGTTTCGCCGAGAGt                                                                                                                         <  2:301367/231‑1 (MQ=255)                                                                                                                              ggaggaCCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCCATTTTCTCTTTCACTGCGGACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGGTACCTGAACCGCCTGCGCATGATGAAAGCCCGACATCTGCTACGCCCCAGCg                                                                                                                                                                                                       >  2:168816/1‑167 (MQ=255)                                                                                                                              ggaggaCCATTTT‑GCCGATGAGGGAAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCTTACGCT‑ACACCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCg                                                                                                                                                                                                       <  1:168816/167‑1 (MQ=255)                                                                                                                               gaggaCCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTg                                                                                                                  >  1:151241/1‑251 (MQ=255)                                                                                                                                aggaCCTTTTT‑GCCGTGGAGGGGATTTGGGATGCCGGGGCGGATCAATTTTCTCTGTCCCTGCTTCCGCT‑ACAGCGGCA‑GCTTAAGAAGAAAACGTGACTGAC‑GCCTAAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACAGCGCCTATCGCTGTGGATGCAGCGACAGTAACCACTTTTCGACTCTTTTTCGCCGAGAGTTTAACTgg                                                                                                                 <  2:675504/251‑1 (MQ=255)                                                                                                                                aggaCCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCCATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGCGAGTTTAACTgg                                                                                                                 >  2:159612/1‑251 (MQ=255)                                                                                                                                  gaCCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTc                                                                                                               >  1:571847/1‑251 (MQ=255)                                                                                                                                   aCCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCCATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGAAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGCGAGTTTAACTTGTc                                                                                                               >  1:467849/1‑250 (MQ=255)                                                                                                                                          t‑GCCGAAGGGGTGAATGAGGATGCCGTGCAGGATCATTTTTCGCTTCCAGTGCGTACCCT‑CAATCGGCA‑GCTAAAGCAGCAAACGGACCGGCC‑GCCTAGGGGATACCTGACCCGCCTGCGACTGATGAAAGCCCGACATCGGCTACGCCACAGCGAGGCCAGCGTTTCTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTGAACTGGTCACCGCGt                                                                                                        <  2:121773/250‑1 (MQ=255)                                                                                                                                           ggcctatggggtgatgtGGGATGCGGTGGCGTATCAATTTTCTCTTGCACTGCGTCCGCT‑ACAGGGGCA‑GCTTAAGAAGCAAACGGGACTGAC‑GCCCCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTGTCTCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACGGGTCACCGCGt                                                                                                        <  2:334356/234‑1 (MQ=255)                                                                                                                                            gCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTCCACGGCGTCCGCT‑ACATCGGCA‑GCTTAAGCAGCAACCGGTACTGCC‑GCCTCAGCGATCCCTAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGa                                                                                                      <  1:330617/251‑1 (MQ=255)                                                                                                                                              cGATGAGGTAAATGGGGATGCCGTGGCGGATCATTTTTCTCTTCCACTGCGTCCGCT‑ACATCGGCA‑GCTAAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCGAAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGata                                                                                                    <  1:258345/251‑1 (MQ=255)                                                                                                                                               gATGAGGTAAATTGGGATGCGGTGGCGGATCATTTTTCTCTTTCACTGCTTCCGCT‑CCATCGCCA‑GCTTAAGCACCAACCGGGCCTTAC‑GCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCTCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGatat                                                                                                   <  1:686201/251‑1 (MQ=255)                                                                                                                                                     gTGAATTGGGATGCTGTGGCGGATAAATTTTCATTTTCACTGCGTCCGCT‑ACTTCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCCTTCCCGGAGCCGCCTGCGACTGATGAAAGCCCGACTTCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCa                                                                                             <  1:447624/251‑1 (MQ=255)                                                                                                                                                           tGGGATGCCGTGGCGGATCAATTTTCTCTTTCACGGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAACCGGGCCTGAC‑GCCTCGGCGATCCCGGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGAc                                                                                        <  1:519943/250‑1 (MQ=255)                                                                                                                                                             ggATGCCGTGGGGGATCAATTTTCTCTTTCACGGCGTACCCT‑ACACCGGCA‑GCTTAACCAGAAAGCGGGACTGAC‑GCCTCAGCGATCCCGGAGCCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACggg                                                                                     <  2:129417/251‑1 (MQ=255)                                                                                                                                                                     gggCGGAACAATTTTCTCTGTCACGGCGT‑CGGTGACTCCGCCA‑GCTGAACCAGCAAACGGGACTGAC‑CCCCCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATGCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCtt                                                                              <  2:590996/249‑1 (MQ=255)                                                                                                                                                                        gggATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTg                                                                          <  1:514292/250‑1 (MQ=255)                                                                                                                                                                          gATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGc                                                                         >  1:191847/1‑250 (MQ=255)                                                                                                                                                                            tCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTACCGCtgtg                                                                                                                                                                      >  2:425579/1‑155 (MQ=255)                                                                                                                                                                                 tttCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGCGTTTAACTGGTCACCGCGTGATATTCGCCCGGGACGGGATGGCTTTCTGCAATAAcgc                                                                 >  2:444383/1‑251 (MQ=255)                                                                                                                                                                                      ctTTCACTGCGTCCGCT‑ACACCGGCA‑CCTTAAGCAGAAAACGGGAC‑GACCGCCTCGGCGATCCCGGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCCTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGaa                                                              <  2:17291/249‑1 (MQ=255)                                                                                                                                                                                       tttCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAAtct                                                           >  2:503213/1‑251 (MQ=255)                                                                                                                                                                                        ttCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCATGCGACTGATGAAAGCCAGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGCAtctt                                                          >  2:640234/1‑251 (MQ=255)                                                                                                                                                                                            ccggggcacgtcacagcGGCA‑GCTAAAGCACGAATCGGGACGGAC‑GCCTCAGCGATCCCGGAGCCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTGCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGTTTTCTGCAGTAACTCGAATCTtctc                                                       <  1:571796/235‑1 (MQ=255)                                                                                                                                                                                               cGTACGCA‑CCATCGGCA‑GCTTAAGCAGCAAACGGGACTGCG‑GCCCCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACg                                                   <  2:264875/251‑1 (MQ=255)                                                                                                                                                                                                gtccGCT‑ACACCGGCA‑GCTTAAGCAGAAAACGGGACTGAC‑CCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATg                                                                                  <  1:368763/216‑1 (MQ=255)                                                                                                                                                                                                gTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATg                                                                                  >  2:368763/1‑219 (MQ=255)                                                                                                                                                                                                   cGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAAAGTAtt                                               >  2:486981/1‑251 (MQ=255)                                                                                                                                                                                                    gCT‑ACATCGCCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATCCCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTGCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGTTTTCTGCAATAACGCGAATCTTCTCAACGTAttt                                              <  1:423532/251‑1 (MQ=255)                                                                                                                                                                                                      t‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGt                                            >  2:21282/1‑251 (MQ=255)                                                                                                                                                                                                        aCATCGGCA‑GCTTAAGCAGCAAACGGGACTGAA‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGCCATCTGCTACGCCACAGCGCGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGCTATTCGCCAGGGACGGGATGGCTTTCTGCAATCACGCGAATCTTCTCAACGTATTTGt                                            >  2:407961/1‑250 (MQ=255)                                                                                                                                                                                                             ggCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGccc                                      >  2:471203/1‑250 (MQ=255)                                                                                                                                                                                                             ggCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGACCCCTGAACCGCATGCAACTGATGAAAGCCCGACATCTGCTACGACACAGCGAGGCCAGCGTTACTGACATCGCCTCTCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGCTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGcc                                       >  1:87441/1‑250 (MQ=255)                                                                                                                                                                                                                cggCTTAAGCAAAAAACGGGCCGGAC‑GCAGCAGGGATCCCGGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCata                                    <  2:571847/249‑1 (MQ=255)                                                                                                                                                                                                                  gCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATAtt                                  >  1:400400/1‑251 (MQ=255)                                                                                                                                                                                                                  ccTTAAGCAGCAACCGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATAtt                                  <  1:20357/250‑1 (MQ=255)                                                                                                                                                                                                                     tAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCg                               >  2:212133/1‑251 (MQ=255)                                                                                                                                                                                                                      aaGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCg                               >  1:471267/1‑250 (MQ=255)                                                                                                                                                                                                                          agcaAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGaata                           >  1:120249/1‑250 (MQ=255)                                                                                                                                                                                                                            caAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGCCGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATc                        >  1:72074/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaCGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGACTCTTCTCAACGTATTTGTACGCCATATTGCGAATAATc                        >  1:599338/1‑250 (MQ=255)                                                                                                                                                                                                                                 gggACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACtt                   >  1:488591/1‑251 (MQ=255)                                                                                                                                                                                                                                     cggAC‑GCCTCAGCGATACCTGACCCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGtt               <  1:682642/249‑1 (MQ=255)                                                                                                                                                                                                                                        gaccccccAGCGATACCCGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTtct             <  2:191847/244‑1 (MQ=255)                                                                                                                                                                                                                                              ccAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCGCTGGcc       <  1:207061/250‑1 (MQ=255)                                                                                                                                                                                                                                               cAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCa                                                      >  1:536583/1‑203 (MQ=255)                                                                                                                                                                                                                                                   gATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTCTGGCCGAGGt  >  1:389703/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        GATTTACCCTCGACCGCCAGTCGCGAGATCTTGTTTATGCAATTACTGCTCTTGCTGCGTAAAAGCAGTTTGCAGGAGAACCTGGAAAACAGCGCATCACGTCTCAACTTGCTTCTGGCCTGGCTGGAGGACCATTTT‑GCCGATGAGGTGAATTGGGATGCCGTGGCGGATCAATTTTCTCTTTCACTGCGTACGCT‑ACATCGGCA‑GCTTAAGCAGCAAACGGGACTGAC‑GCCTCAGCGATACCTGAACCGCCTGCGACTGATGAAAGCCCGACATCTGCTACGCCACAGCGAGGCCAGCGTTACTGACATCGCCTATCGCTGTGGATTCAGCGACAGTAACCACTTTTCGACGCTTTTTCGCCGAGAGTTTAACTGGTCACCGCGTGATATTCGCCAGGGACGGGATGGCTTTCTGCAATAACGCGAATCTTCTCAACGTATTTGTACGCCATATTGCGAATAATCAACTTCGTTCTCTGGCCGAGGT  >  NC_000913/4098150‑4098638

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