Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,446,592	C→A	P294P (CCG→CCT)	aroC ←	chorismate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,446,592	0	C	A	98.8%	229.7 / NA	80	P294P (CCG→CCT)	aroC	chorismate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (45/34);  minor base G (0/1);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACG  >  NC_000913/2446357‑2446835                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ccatcagaccAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGAGCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCatc                                                                                                                                                                                                                                      <  2:603631/242‑2 (MQ=255)     cagcGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATggg                                                                                                                                                                                                                                  >  1:28843/1‑251 (MQ=255)      agcGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGc                                                                                                                                                                                                                                 >  1:430712/1‑251 (MQ=255)         gcAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAAt                                                                                                                                                                                                                              >  1:222034/1‑251 (MQ=255)            aTCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGGAGGTTTCATCGCCATATggggcatgat                                                                                                                                                                                                                           >  2:691842/1‑244 (MQ=255)               gcgGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTtg                                                                                                                                                                                                                        <  1:181123/251‑1 (MQ=255)                               tttACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGAt                                                                                                                                                                                                         >  1:182186/1‑250 (MQ=255)                                 tACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAcc                                                                                                                                                                                                      >  2:167714/1‑251 (MQ=255)                                  accagctgggAATATCAGTCTTCACATCGGCATTTTGCGCCCTTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAcc                                                                                                                                                                                                      <  2:605568/242‑1 (MQ=255)                                  aCCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccg                                                                                                                                                                                                     <  1:140768/251‑1 (MQ=255)                                    cAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCTCGCTGCTGATAccgcc                                                                                                                                                                                                   >  1:176006/1‑251 (MQ=255)                                       cGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGGGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATAccgccg                                                                                                                                                                                                  >  2:338386/1‑249 (MQ=255)                                          ggAATATCAGTCTTGACATCGGCATTTGGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCGGCGATCGGCTTTGCGGGGGTC‑CCGAC‑ACAGGGATCGTGGCG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTCCGGCCCGCTGCTGATACCGCCGAGAAt                                                                                                                                                                                             <  1:585015/251‑1 (MQ=255)                                          ggAATATCAGTCTCCACATCGGCATTTTGCTCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                                                                                                                                             <  1:782601/251‑1 (MQ=255)                                          gaaaTATCAGTCTTCACTCCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATAAAAACGATCGCCAGCAGCGCTTCTGCGATCGGCACGGCGCGGACC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTCCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGTAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAAt                                                                                                                                                                                             <  2:371886/249‑1 (MQ=255)                                               atCAGTCTACACATCGGCATGTTGCGCCCGTTGCCGTAACAGGTGACCCATTAAAACGACCCCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTCCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTTCTGCCCGCTGCTGATACCGCCGAGAATgccg                                                                                                                                                                                         <  2:739515/250‑1 (MQ=255)                                                acagccTTCACATCGGCTTTTTGCCCCCGTTGCCGTAACAGGGGATCCATTAAACCGATCGCCAGCACCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCTTTTGGTGATCATCTCATCTTCTTCGCCAAAGCGGTTAATGGTACGCCCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATgccgcc                                                                                                                                                                                       <  2:691210/246‑1 (MQ=255)                                                   gTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGCC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGCGGTAGGTTTCAGCGCCATATGGGCAATGTTTTGCTGCCCGCTGCTGATACCGCCGCGAATGccgcccgc                                                                                                                                                                                    >  2:257409/1‑251 (MQ=255)                                                          caTCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTg                                                                                                                                                                             >  1:25285/1‑251 (MQ=255)                                                           aTCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGc                                                                                                                                                                            >  1:94632/1‑251 (MQ=255)                                                             gggCATTTCGCGCCCGTTGCCGTAAGAGGTGACCCATTAAACCGATCCCCCGAATCGCTTCTGCGACCGGCCCTGCGCGGTCC‑CCGAC‑ACAGGGATCGGGACG‑GCCTTGGGTGATCTTCTCAACTTCGCCGCCAAAGCGGTTAAGGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTGTCAGCGCCATAGGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGctc                                                                                                                                                                          <  2:791468/250‑1 (MQ=255)                                                             cggctttggcaccccGTTGCGGTACCAGGGGATCAATAAAAACATTCGCCACCTCCGCTTCTGCGATCGGCCCGGC‑CGG‑TCACCCACCACGGGGATCGTGACG‑GCCTTTGTGGATCATCTCAACTTCGTCGCCAAAGCGTTTAATGGTACGCCCCGGCACGGTAAGGCTGGAGGTAGGTCTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGct                                                                                                                                                                           <  2:139454/238‑1 (MQ=255)                                                                cATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGAGGATC‑CCGAC‑ACAGGGATCGTGACG‑GACTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGAGCCATATGGGCAATGAATTGCTGCCCGCTGCTGATACCGCAGAGAATGCCGCCCGCATGGTTGCTCTgg                                                                                                                                                                       >  1:792965/1‑251 (MQ=255)                                                                 aTTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTgg                                                                                                                                                                       >  1:552085/1‑250 (MQ=255)                                                                   ttggcccccGTTGCCGTAACAGGTGATCCAAAAAAACGATCGCCAGCACCGCTTTGGCTATCGGCACTGCGCGTATC‑CCGAC‑ACAGGGATCGTGACG‑GCTTTGGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTGTCAGCGCCATAGGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGaa                                                                                                                                                                     <  2:430712/244‑1 (MQ=255)                                                                       cgcccgccgCCGTACCAGCTGCTCCATTAAAACGACCCACAGCATCGCTTCTGCGATCGGCACTGCGCGGCTC‑CCGAC‑CCAGGGCTCGTGACG‑GCCTTTGGTGATCATCTCACCTTCTTCGCCCAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGCGAATGCCGCCCGCATGGTTGCACTGGaacacgt                                                                                                                                                                >  2:106933/9‑246 (MQ=255)                                                                       ccccctttGCCGTAACAGGTGTTCCATTAAACCGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGTGGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGt                                                                                                                                                                <  1:452553/245‑1 (MQ=255)                                                                            gTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTg                                                                                                                                                           >  1:351549/1‑251 (MQ=255)                                                                                  gTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCCATGATTTGCTGCCCGCGGCTGATACCGCCGAGCATGCCGCCCGCATGGTTGCTCTGGTAACCGTCCTTggggctt                                                                                                                                                     >  2:268881/1‑246 (MQ=255)                                                                                      cAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCCCGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCCATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTcct                                                                                                                                                 >  2:498418/1‑249 (MQ=255)                                                                                           gATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGCAACCGTCTTTGGTGATTTCATCGCgg                                                                                                                                            >  1:324353/1‑251 (MQ=255)                                                                                           gATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCgg                                                                                                                                            >  1:169126/1‑251 (MQ=255)                                                                                            aTCCATAAAAATGATCGCCAGCACCGCTTCTGCGACCGGCTCTGCGCGGATC‑CCGAC‑ACAGGGATCGTGATG‑GCCTTTGTTGATCATCTCAACTTTTTCGCCAAAGCGGTTAATGGTACGCCCCGGCACGGTAATGCGGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGt                                                                                                                                           <  1:691842/251‑1 (MQ=255)                                                                                                aTTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTg                                                                                                                                       >  1:605110/1‑251 (MQ=255)                                                                                                     aaCGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGc                                                                                                                                  >  2:785205/1‑251 (MQ=255)                                                                                                        gATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgc                                                                                                                               >  1:104867/1‑251 (MQ=255)                                                                                                         aTCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCcgcg                                                                                                                              >  1:239007/1‑251 (MQ=255)                                                                                                                  aTCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGccac                                                                                                                     >  1:188953/1‑251 (MQ=255)                                                                                                                    cgtttcgcggatgtccACTGCGGGGATC‑CCGAC‑TGAGGGATGGTGCCG‑GCCTTTGGGGACCACCCAAACGCCTCCGCCAAGGCGTTTAAGGGGACGACCCGGCACGTGAAGGCGGGAGGTAGGTTTCAGCGCAATATGGGCAATGATTTGCTGCCAGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGccacc                                                                                                                    <  2:28843/235‑1 (MQ=255)                                                                                                                         cTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTc                                                                                                               >  1:734547/1‑250 (MQ=255)                                                                                                                            cGATCGGCACTGCGCGGATC‑CCGAC‑ACCGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTCATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGACCGCATGGTTGCTCTGGAACCCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCCCGTCaa                                                                                                             >  2:521573/1‑249 (MQ=255)                                                                                                                            cGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCaa                                                                                                             >  1:293999/1‑249 (MQ=255)                                                                                                                             gATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCg                                                                                                                                                                 >  1:682321/1‑196 (MQ=255)                                                                                                                              aTCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGcc                                                                                                         >  1:167159/1‑251 (MQ=255)                                                                                                                                ccgcccTGCGCGGCTC‑CCGAC‑ACAGGGTTTGTG‑CGTGCCTTTGGTTATCATCTCACCTTCTTCGCCAATGCGTTTAAGGTTACGACCCGCCCCTGTAATGCTGGAGGTAGGTTTCAGCCCCTTATGGCCAATCATTTGCTGCCCGCGTCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGt                                                                                                       <  2:182186/246‑1 (MQ=255)                                                                                                                                 gcaccTGCGCGTGTC‑CCGAC‑CCAGGGATCGTGCCG‑GCTTTTGGTGATCATCTCAACTGCTTCGCCAAAGCGGTTAATGGTACGACCCGGCCCGGTAATGCTGGGGGTAGGTTTCAGCGCCATAGGGGCATTGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTc                                                                                                      <  2:316387/247‑1 (MQ=255)                                                                                                                                   cACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGAGATCATCTCAACTTATTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAGTAATGTGGGAAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCACGCTGCTGATACCGACGAGAATTCCGCCCGCATGGTTGCTCTGTAAAACAACTTTGGTGATTTCAACCCGTTACTGGCGTCCGCGCAGCGCCACCACGTCAAAGACGTCGc                                                                                                    >  1:423773/1‑251 (MQ=255)                                                                                                                                        cgggtgttcccacc‑acaGGGACCGTGACG‑GCTCTTGGTGATCACCTCACCTCCCTCCCCAAAGCGGTTAAGGGTACGACCCGGCACGAGAGTGCGGGAGGTAGGTTTCATCCCCATAGGGGCAAGGATTTGCTGCCCGGTGCTGACCCCGCCGAGAAGGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCa                                                                                                  <  2:473002/236‑1 (MQ=255)                                                                                                                                         gcgGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAACCCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCCGCGCCACCCCGTCAAAGCCGTCGCCACttt                                                                                              >  2:327122/1‑251 (MQ=255)                                                                                                                                          aggATC‑CCGCA‑ACAGGGATCGTGCCG‑GCCTTTGGTGATCATCTCACCTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCCCGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGACACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCCAAGCCGTCGcacatttc                                                                                             >  2:734872/2‑244 (MQ=255)                                                                                                                                                    ac‑acAGGGAGCGTGCCG‑GCTTTTGTTGATCATCGCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTGCAGCGCCATATGGGCAATGATTTGCTCCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGTAACCCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCttt                                                                                    <  1:31446/251‑1 (MQ=255)                                                                                                                                                         aGGGATCGTGACG‑GCCTTTGGTGATCCTCCCAACTTTTTCGCCAAAGCGGTTAATGGGACGACCCGGCACGGTAAGGCTGGAGGTAGGTTACAGCGCCATATGGGCAGTGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTTCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCAcc                                                                                <  2:667090/251‑1 (MQ=255)                                                                                                                                                            gATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGTAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACcgcg                                                                             <  1:585067/251‑1 (MQ=255)                                                                                                                                                               cGTGACG‑GCTTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGtt                                                                           <  1:227736/250‑1 (MQ=255)                                                                                                                                                                 tGACT‑GCCTTTGGTGATCATCTCAACTTCTTCGCCCAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGCTACCGCCGAGCATGCCGCCCGCATGGTTGCTCTGGCAACCGACTTTGGAGATTTCATCCCGGTTCTGGCTGCCGCGCCGCGCCACCACGTCAAAGCCCGCCCCAATTTGAACGCCTTTCACCGCGtttat                                                                        >  2:333330/1‑248 (MQ=255)                                                                                                                                                                  gACG‑GCTTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGAt                                                                        <  1:336900/250‑1 (MQ=255)                                                                                                                                                                              ggtatcatcTCAACTTCTTCGCGAAAGCGGTTAAGGGTACGCCCCGGCACGGTAATGCTGGAGGTAGGTTGCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAgcgc                                                            <  2:25285/248‑1 (MQ=255)                                                                                                                                                                                gATCACCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAAGGCTGGAGGTAGGTTGCAGCGCCATAGGGGCAATGATTTGCTGCCCGCTGCTTATCCCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCAt                                                          <  1:366766/251‑1 (MQ=255)                                                                                                                                                                                   catcTCAATGTCTTAGCCAAAGCGGTTATTGGCACGCCCGGGCACGGTAATGCTGGAGGTAGGTTTAAGCGCCATAGGGGTAATGATTTGCTGCCCTCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATgtg                                                       <  1:816722/251‑3 (MQ=255)                                                                                                                                                                                   catcTCAACTGCTCCGCCAAAGCGGTTAATGGTCCGACCCGGCAGGGTAATGCGGGAGGGAGGTTTCAGGGCCATATGGGCAATGATTTGCGGCCCGCTGCTGGTACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATggg                                                       <  2:297768/251‑1 (MQ=255)                                                                                                                                                                                     tctcAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCg                                                     >  1:431206/1‑251 (MQ=255)                                                                                                                                                                                        catcttcttCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGACACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGc                                                      <  1:32242/244‑1 (MQ=255)                                                                                                                                                                                           cttcttCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTc                                                >  1:48863/1‑250 (MQ=255)                                                                                                                                                                                           cttcttCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTc                                                >  2:634473/1‑250 (MQ=255)                                                                                                                                                                                             tctgccccAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCa                                               <  1:338386/243‑1 (MQ=255)                                                                                                                                                                                              cttcGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGCGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCCCGTCAAAGCCGTCGCCAATTTCCCCGCCTTTCACCCCGTTGATTCTCATCAGCGCATGGGCTCTGTCAGc                                             >  2:667375/1‑250 (MQ=255)                                                                                                                                                                                                ccGCCAAAGGGGTTAATGGTACCACCCGGCACGGTAATGGTGGGGGTAGGTTGCCGCGCCGTATGGGCAATGATTGGCGGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGGTGATGCTCATCAGCGCATGGGCGATGTCAGCATc                                          <  2:94632/250‑1 (MQ=255)                                                                                                                                                                                                 cGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATcc                                         >  2:382759/1‑251 (MQ=255)                                                                                                                                                                                                       aGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCCTGGGCGATGTCAGCATCCAGGCgg                                   >  2:744719/1‑251 (MQ=255)                                                                                                                                                                                                           gTTAATGGTACGACCCGCCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCaa                               <  1:175301/251‑1 (MQ=255)                                                                                                                                                                                                                tGGTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGAcc                          >  2:99602/1‑251 (MQ=255)                                                                                                                                                                                                                  gTACGACCCGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCgg                        >  1:814279/1‑251 (MQ=255)                                                                                                                                                                                                                     cgccccGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATAGGGGCAATGATTGGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTc                     <  1:167714/248‑1 (MQ=255)                                                                                                                                                                                                                          cGGCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCa                 >  1:134088/1‑250 (MQ=255)                                                                                                                                                                                                                           ggCACGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGccc                 >  1:164535/1‑248 (MQ=255)                                                                                                                                                                                                                              ccGGTAATGCGGGGGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCGGCTCGCTGCTGATCCCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGAGTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTcc            <  2:176006/250‑1 (MQ=255)                                                                                                                                                                                                                              aCGGTAATGCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCAGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGCTGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTcc            >  1:456419/1‑251 (MQ=255)                                                                                                                                                                                                                                    aTGCTGGCGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGACAAAGACCGGCTctcccaatcccgccgg      >  2:337628/1‑236 (MQ=255)                                                                                                                                                                                                                                      gCTGGAGGTAGGTTTCAGCGCCATATGGGCAATGATTTGCGGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGa     <  1:498418/250‑1 (MQ=255)                                                                                                                                                                                                                                        gggAGGTGGGTTCCAGCGCCATATGGGCAATGTTTTGCTGCCCGCTGCGGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACg  <  2:605110/250‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    CCAGCAGCGCAATCGCGGTTTTATTCATTTTTTACCAGCGTGGAATATCAGTCTTCACATCGGCATTTTGCGCCCGTTGCCGTAACAGGTGATCCATTAAAACGATCGCCAGCATCGCTTCTGCGATCGGCACTGCGCGGATC‑CCGAC‑ACAGGGATCGTGACG‑GCCTTTGGTGATCATCTCAACTTCTTCGCCAAAGCGGTTAATGGTACGACCCGGCACGGTAATGCTGGAGGTCGGTTTCAGCGCCATATGGGCAATGATTTGCTGCCCGCTGCTGATACCGCCGAGAATGCCGCCCGCATGGTTGCTCTGGAAACCGTCTTTGGTGATTTCATCGCGGTTCTGGCTGCCGCGCAGCGCCACCACGTCAAAGCCGTCGCCAATTTCCACGCCTTTCACCGCGTTGATGCTCATCAGCGCATGGGCGATGTCAGCATCCAGGCGGTCAAAGACCGGCTCGCCAAGTCCGGCAGGAACG  >  NC_000913/2446357‑2446835

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