BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,558,907	T→A	K207M (AAG→ATG)	iadA ←	isoaspartyl dipeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,558,907	0	T	A	95.0%	277.5 / ‑4.1	101	K207M (AAG→ATG)	iadA	isoaspartyl dipeptidase
Reads supporting (aligned to +/- strand): ref base T (0/5); new base A (42/52); total (42/59)
Fisher's exact test for biased strand distribution p-value = 7.04e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 104 total aligned reads displayed.

CTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558661‑4559152

ctACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACATTCGGCAGGGGCGACCGGTTCGTCCATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCATGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCACCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                                                                                                                                                                                                     >  1:303201/1‑251 (MQ=255)

 atccGTTGCGGGCGGGGCTGAGGGGGACGCGTCCCAGCTGAATACCTGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTAAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑TCGTGGGTCG‑GCAGCAGCTTTCTGATCGGCACAGC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTATCGCCatttt                                                                                                                                                                                                                                                     <  2:161187/248‑1 (MQ=255)

  aCCGTTGCCGTCGGAGCTGAGGGTTACGCGGGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCatttt                                                                                                                                                                                                                                                     <  1:146103/249‑1 (MQ=255)

        cccGTCGGGGCGGGGGGTGACGGCTGCCAGCGAAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCGGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGACCTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑Gc                                                                                                                                                                                                                                            <  2:260367/250‑1 (MQ=255)

         ccGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCTGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTTCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTTGGTCG‑GCAGCAGCTTGCTGATCGGCACGTC‑GCA‑GTT‑TG‑CCAGCAGGTCATAAATAGGCTTTAACGCCATTTTACTGTC‑Gcc                                                                                                                                                                                                                                           <  2:436401/251‑1 (MQ=255)

         ccGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑Gcc                                                                                                                                                                                                                                           >  2:261893/1‑251 (MQ=255)

              ggAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCACTACCTTCGGCAGGGGCGACCGGTTCGTCAATGATGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAAAA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCAtgt                                                                                                                                                                                                                                      >  2:248644/1‑251 (MQ=255)

              ggAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTGCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCGTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGGTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCAtgt                                                                                                                                                                                                                                      <  2:171039/251‑1 (MQ=255)

                         tGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGACCTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTg                                                                                                                                                                                                                           >  2:157805/1‑251 (MQ=255)

                          gACGCGTGCCAGCGGAATCCCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGATGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGa                                                                                                                                                                                                                          <  1:412272/251‑1 (MQ=255)

                           gcgcgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATTCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGAc                                                                                                                                                                                                                         <  1:466128/250‑1 (MQ=255)

                           ccgcgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGAc                                                                                                                                                                                                                         <  2:449581/250‑1 (MQ=255)

                           aCGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGAc                                                                                                                                                                                                                         >  1:106500/1‑251 (MQ=255)

                            cgcgTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAAGAGCGACCGGTTCGTCCATGATGCTGGTGCTCTCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCCGGTCCTAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACg                                                                                                                                                                                                                        >  1:312431/1‑251 (MQ=255)

                               gTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGcc                                                                                                                                                                                                                      >  1:104435/1‑250 (MQ=255)

                                       ggAATACCCGCCTGCACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTAcc                                                                                                                                                                                                             >  1:122176/1‑251 (MQ=255)

                                          aTACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTAccgcc                                                                                                                                                                                                          <  1:417262/251‑1 (MQ=255)

                                            aCCCGCCTGAACGGGCCGGCCAATACATCCGGCAGGGGCGGCCGGTTCGTCAATGCTTCTGGTGAGATCGATGGTGCCGCCTTTGTGCGGAAACTCCAGCGCCTGCTCGAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑CCAGCAGCTTGCTGAGCGGCACATC‑GCA‑GTT‑TT‑CCAGCTGGACATAATTTGGCTGCAACGCCCTATTACTGTC‑GCCAATCTGGAACACGGTGAGGCCGGGTTTAccgacg                                                                                                                                                                                                         <  2:354990/250‑4 (MQ=255)

                                             cccGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAg                                                                                                                                                                                                       <  1:154226/251‑1 (MQ=255)

                                               cGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCa                                                                                                                                                                                                     <  1:155628/251‑1 (MQ=255)

                                                gCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCaa                                                                                                                                                                                                    <  1:97797/251‑1 (MQ=255)

                                                  cTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc                                                                                                                                                                                                  >  1:1642/1‑251 (MQ=255)

                                                    gAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCGGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGGGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAAc                                                                                                                                                                                                  <  2:33945/249‑1 (MQ=255)

                                                      aCGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATGTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑Gcc                                                                                                                                                                                             <  2:309132/251‑1 (MQ=255)

                                                       cGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑Gcc                                                                                                                                                                                             >  1:166714/1‑250 (MQ=255)

                                                         gcgcGGGAAATCCCTGCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGTGGCCGCCTTTGCGCGCGAACCCCAGCGCCGGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAAc                                                                                                                                                                                          <  2:188120/251‑1 (MQ=255)

                                                           gcgGGCAATCCCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAAcgc                                                                                                                                                                                        <  1:52403/251‑1 (MQ=255)

                                                              ggCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCggg                                                                                                                                                                                     <  1:158799/251‑1 (MQ=255)

                                                                 aatcccTTCGGCAGGGGCCACCGGTCCGTCAATGCTGCGGGTGATATCGATGGTGCCGCCTGTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGGACGTTGCGGTTA‑ACGTGGGCCG‑GCAGCAGCTTGCTGATCGGCACATC‑CCA‑GTT‑TT‑CCAGCATGTCATAAATAGGCTGTAGCGCCATTTTACGGTC‑GCCCAGGTCGAACACGTTTACGCCAGTTTTACCGCCCGGCAAACC‑GCCAACGCGGGAtt                                                                                                                                                                                  <  2:444460/247‑1 (MQ=255)

                                                                 aaTACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGCGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGAtt                                                                                                                                                                                  >  2:249379/1‑251 (MQ=255)

                                                                        tCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTgcc                                                                                                                                                                             <  1:157805/249‑1 (MQ=255)

                                                                         cGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCCAACC‑GCCAACGCGGGATTCTgccgcc                                                                                                                                                                          >  2:274824/1‑251 (MQ=255)

                                                                          ggCAGGGGCGACCGGTCCGTCAATGCTGCTGGCGATATCGCTGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGCCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGGCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTgccgccc                                                                                                                                                                         >  2:281791/1‑250 (MQ=255)

                                                                           gCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCat                                                                                                                                                                        >  2:385368/1‑251 (MQ=255)

                                                                            cAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAaca                                                                                                                                                                                                                                >  1:51781/1‑195 (MQ=255)

                                                                            cAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTGTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCATCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACTCCATTTTACTGTC‑GCCCATGTGGAaca                                                                                                                                                                                                                                <  2:51781/195‑1 (MQ=255)

                                                                                  cGACCGGTTCGTCAATGCTGCTGGTGATCTCGATGGTGCCGCCTTTGCGCGCGACCTCCAGCGCCTGCTCAAACC‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGcc                                                                                                                                                                 >  2:229519/1‑251 (MQ=255)

                                                                                     ccGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGCTTCTGCCGCCATATTGGCCAgg                                                                                                                                                              >  2:225626/1‑251 (MQ=255)

                                                                                      cGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGcccggg                                                                                                                                                             >  2:268706/1‑247 (MQ=255)

                                                                                            gTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCACACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATaa                                                                                                                                                        >  2:227639/1‑250 (MQ=255)

                                                                                             tCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAAc                                                                                                                                                      >  1:241632/1‑251 (MQ=255)

                                                                                                      ctgGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTgg                                                                                                                                             >  2:68787/1‑250 (MQ=255)

                                                                                                        gggaatatggatgggcccgccgTTGTGCGCGAACTCCAGCCCGTGTTCAGAAA‑ACGTTCCGTGGCGGTTC‑CCGGGGGTCG‑CAAGCAGCTTTCTGCGCGGCCCATC‑‑CACGTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTTTC‑GCCAATGTGGAACCCGGTGACGCCAGGTTTACCGCCGAGAAAACC‑GCCAACGCGGGATTCTTCCGCCATATTGGCCAGGTGATAAACGTCCGGTGc                                                                                                                                             <  2:278732/233‑1 (MQ=255)

                                                                                                        ggTGATACCGATGGTGCCGCCTTCGCGCGCGAACTCCAGCGCCTGCTAAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGGCG‑GCAGCAGCTGGCTGATCGGCACAGC‑GCA‑GTT‑TT‑CCAGCAGCGCATAAATAGGCTGGAACGCCATTTTACTGTC‑GCCCATGTGGACCACGGTGACGCCAGGTTTACCGCGGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCgg                                                                                                                                           <  2:130036/251‑1 (MQ=255)

                                                                                                         gTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGCGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGc                                                                                                                                          >  2:38219/1‑251 (MQ=255)

                                                                                                             aatCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTCCGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGaa                                                                                                                                      <  1:415605/250‑1 (MQ=255)

                                                                                                              agcgcgggTGCCGTCTTTGCGTGCGAACCCCAGCGCC‑GGTCAAAAACACGGTACGTTGCG‑TTACACGTGGGTCG‑GAAGCTGCTTGCTGTTCGCTACGTC‑GCA‑GTT‑TT‑CCAGCAGGTCATGAATGGGCTGTAACGCCATTTTCCTGTC‑GCCCATGTGGAACACGGTGACGCGTGGTTTACGGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAAc                                                                                                                                     <  2:355332/245‑1 (MQ=255)

                                                                                                               ccgacggTGCCCCCTTTGCGCGCGAACGCCACCTCCTGCGCAAACA‑ACGTCACGTTTCGTTTA‑ACGTGGGTCG‑GCAACAGTGTGCCGATCGGCTCATC‑GCA‑TTT‑TG‑CCAGCAGGTCATAAATAGGCGTTAACCCCATTTTACTGTC‑GCGCAGGTGGAGCGCGGTGACGCGAGGTTTACCTCCGAGCAATCC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACg                                                                                                                                    <  1:69765/246‑1 (MQ=255)

                                                                                                                 gagggTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGGTGCGGTTA‑CCGTGGGTCG‑GCAGCAGCTTGCTGGTCGGCCCATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGTTGTAACGCCATTTTACTGTC‑GCCCATGTGGACCAGGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGt                                                                                                                                  <  1:469392/248‑1 (MQ=255)

                                                                                                                     gTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑CCGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGAAACGCCATTTGTCTTTC‑GCCCAGGTGGAACACGGGGACGCCAGGTGTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATc                                                                                                                              <  1:229519/251‑1 (MQ=255)

                                                                                                                       gcctccTTGGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCag                                                                                                                            <  1:257970/247‑1 (MQ=255)

                                                                                                                       gccgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCTATGTGGAACTCGGTGACGCCAGGTTGCCCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCag                                                                                                                            <  1:38219/251‑1 (MQ=255)

                                                                                                                         cgccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCag                                                                                                                            >  1:140333/1‑249 (MQ=255)

                                                                                                                          gccTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTCCGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAgaga                                                                                                                         <  1:385368/251‑1 (MQ=255)

                                                                                                                              ttGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATgg                                                                                                                      <  1:112253/250‑1 (MQ=255)

                                                                                                                                gcgcgcGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGgcg                                                                                                                    <  1:172880/250‑1 (MQ=255)

                                                                                                                                     cgAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTc                                                                                                              >  1:217974/1‑251 (MQ=255)

                                                                                                                                        aCTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACg                                                                                                           >  1:277706/1‑251 (MQ=255)

                                                                                                                                        aCTCCAGCACCTGCTAAACCA‑ACGGTACGTTGCGGTTA‑CCGGGGGTCG‑GCAGCAGCTTGCTGATCGGCACTTC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACg                                                                                                           <  2:167367/251‑1 (MQ=255)

                                                                                                                                              gctccggCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACTTC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCGGTATCGCCATTTTACTGTC‑GCCCATTGGGAACACGGTGAGGCCAGGTTTACCGCCGAGCAAGCC‑GCCAACGCGGGATTCTGCCGCCATATTGGCGAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                     <  1:292593/245‑1 (MQ=255)

                                                                                                                                               cgcCTGCTCCAACA‑ATGGAACGATGAGGTTA‑CCGTGGGTTG‑GCAGCAGCTTGCTGATAGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTACCGCCATTTTACTGTC‑GCCCCTGTGGAACACGGTGACGCCAGGTTTACCGCAGAGCAAAAC‑GCAAACGCGGGATTCTGCCGCCATATTGGCCAGGTGCTAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                     >  2:331560/1‑250 (MQ=255)

                                                                                                                                               cgcCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATc                                                                                                     >  2:144920/1‑250 (MQ=255)

                                                                                                                                                    gCTCAAACA‑ACGGTACGTTGCGGATA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcg                                                                                                >  2:53641/1‑250 (MQ=255)

                                                                                                                                                    gCTCAAACA‑ACGGTACGGTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACAcga                                                                                               <  1:249379/251‑1 (MQ=255)

                                                                                                                                                     cctaaaaa‑acGGTGCGTTGCGTTTA‑CCGTGGGTCG‑GCAGCAGCTGCCTGATCGGCACAGC‑GC‑‑GTTGTC‑CCAGCAGGTCATAAATGGCCGGTAACGCCATTTTACGGCC‑GCCCTTGGGGAACACGGTGCCCCCGGGTTTCCCGCGGAGAAACCC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACTCCAATCACAcgat                                                                                              <  2:279231/244‑1 (MQ=255)

                                                                                                                                                                 gTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGcc                                                                                   >  1:369782/1‑251 (MQ=255)

                                                                                                                                                                    cGTTGCGGTTA‑ACGGGGGTCG‑GACGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAAGAGGCTGTAACGCCATTTTCCTGTC‑GCCCATGGGGAACAGGGTGCCGCCAGGTTTACCGCCGAGAAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACg                                                                                <  2:30355/251‑1 (MQ=255)

                                                                                                                                                                     gTTGCGGTTA‑ACGTGGGTCG‑GTAGTAGCTTGCTGACCGGACCATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAAAAGGCTGTAAGGCCATTTTACTTTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGAGCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                               <  1:278982/251‑1 (MQ=255)

                                                                                                                                                                     gTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                               >  2:143439/1‑251 (MQ=255)

                                                                                                                                                                       ttccggtgaccggtggTCG‑GCAGCAGCTTGCGGAGCGGCCCTTC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTACCGCCATTTTACTGTC‑GCCCAGTTGGAACGCGGGGACGCCAGGTTTACCGCCGAGCAACCC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGt                                                                               <  2:280333/236‑1 (MQ=255)

                                                                                                                                                                       tGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTc                                                                              <  1:168690/250‑1 (MQ=255)

                                                                                                                                                                       tGCGGTAA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTc                                                                              <  1:29183/250‑1 (MQ=255)

                                                                                                                                                                       ggCGGTTA‑ACGGGGGGAG‑GCAGGTGCTTCCTGATCGGCA‑ATCCGCA‑GTT‑TG‑CCAGAAGGTAATAAAGAGGCGGTACCGCTTTTTTACGGTC‑CCCCATGGGAAACCCGTGGCCGCCGGGTTTACCGCCGAGAAAACC‑GCCAACTGGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCt                                                                             <  2:317881/250‑1 (MQ=255)

                                                                                                                                                                           gTTA‑ACGTGGGTCGCGC‑GCACCTTGCTGATCGGCCCAGC‑GCA‑GTT‑TT‑CCAGCATGTCATAATAAGTCGGTATCGCCATTTTAC‑GTCCGCCCATGTGGAACCCGGTGACGCCAGGTTTCCCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGATATGGCGCATTTCACGCCAAGCACACGATCGATAATCGCCACGTCttttt                                                                         <  1:254582/251‑1 (MQ=255)

                                                                                                                                                                                aacggggtccggcgccggcTTGCTGATCGGCACATG‑TCA‑GTT‑TT‑CCTGCGGGGCATAAAAAGGCTGAAACCCCATTTTACTGTC‑GCCTTTGTGGAACCCGGGGACGCCAGGTTTGCCGCCGAGCAAACC‑GCCAACGCGGAATTCTGCCGCCATATTGGCCAGGTGAGAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACGCTATCGATAATCGTCGCGTCTTTTTcc                                                                       <  2:78026/233‑1 (MQ=255)

                                                                                                                                                                                  gTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACgg                                                                   >  1:209130/1‑251 (MQ=255)

                                                                                                                                                                                  gTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCAGGTGGAACCCGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACgg                                                                   <  1:150434/251‑1 (MQ=255)

                                                                                                                                                                                    gggTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGACCACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGaa                                                                 <  1:260050/251‑1 (MQ=255)

                                                                                                                                                                                       ccG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TG‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGGGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAAccc                                                              <  1:8770/250‑1 (MQ=255)

                                                                                                                                                                                             agcagcTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAAt                                                         <  1:7833/251‑1 (MQ=255)

                                                                                                                                                                                              gcagcTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAAGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATg                                                        >  2:96315/1‑251 (MQ=255)

                                                                                                                                                                                                agcTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGt                                                      >  1:34679/1‑251 (MQ=255)

                                                                                                                                                                                                     gCTGATCGGCACATA‑GCC‑GTT‑TT‑ACAGCCGGTCATAAATAGGCTGTAACGCCATTTTACTCTG‑GCCCATGTGGCACAAGGAGACGCCAGGTTTACCGACGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATCTTGGCCAGGAGCTAAACGTCCGGTGCGGCCGAACGGTGATCCGAGATGGAGCATTTCCAGACAATCCCACGATCGATAATCGCCACGTCTTTTTCAACGGATCACGTAAGGGTGCggg                                                 >  2:396266/1‑251 (MQ=255)

                                                                                                                                                                                                       gccggtgggcccaccga‑GTT‑TTACCAGCAGGTAAAAAAAGGGCTGTAACGCATTTGTACGGTC‑GCCAGTGTGAAACCCGTTGGCGCAAGGTTTACCGCGGAGAAACCC‑GCCAACGCGGGATGCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAAGCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCggg                                                 <  2:304504/235‑1 (MQ=255)

                                                                                                                                                                                                        gATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGa                                                >  2:78979/1‑249 (MQ=255)

                                                                                                                                                                                                            ggCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGTGAc                                          >  2:90523/1‑251 (MQ=255)

                                                                                                                                                                                                             gcccaTC‑GCA‑GTT‑TT‑CCAGCAGGTCAAAAATGGGCTGTAACGCCTTTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGTTTTCCCGCCGAGAAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACa                                         <  2:382674/248‑1 (MQ=255)

                                                                                                                                                                                                              cacaTC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCCCGGAACCCGTAATGGTGCGGGAAGGGACAt                                        >  2:300254/1‑251 (MQ=255)

                                                                                                                                                                                                               ccttC‑GCA‑TTT‑TC‑CCAGCAGGTCATAAATAGGCTGTAACGCCATTTTACTGCC‑GCCCATGTGGACCACGGTGCGGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACTGCTTTTTCCACGGAGCCCGTAAGGGTGCGGGAAGGGACATg                                       <  1:386893/248‑1 (MQ=255)

                                                                                                                                                                                                                caTC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAAGAGGCTGAAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGa                                      <  1:227639/251‑1 (MQ=255)

                                                                                                                                                                                                                              ccccAGCGGGTAATAAATGGGCTGAAACCCCTTTTGCCTGTC‑GCCCGTGGGGCGCACGGTGAGGCCGGGTTTCCCGCCG‑GCAAACCCGCCAACGAGGGATTCGTCCGCCATATTGGCCAGGTGATAAACTTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACCCGATCGATAATCGCCACGTCTTTGTCCACGGAACCCGTAATGGTCGGGGAACGGACATGATTAGCGCCg                             <  2:293189/248‑1 (MQ=255)

                                                                                                                                                                                                                                 cagcagGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                          >  1:84265/1‑250 (MQ=255)

                                                                                                                                                                                                                                  agcagGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTc                          >  2:359301/1‑249 (MQ=255)

                                                                                                                                                                                                                                     agGTCATAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAAGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCAt                     <  1:96315/251‑1 (MQ=255)

                                                                                                                                                                                                                                           tAAATAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgc               >  2:148643/1‑251 (MQ=255)

                                                                                                                                                                                                                                           aaaaTGGGCTGGAACGCCATTTTACTGTC‑CCCCAGGGGGAACACGGGGACCCCTGTTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgc               <  2:369782/250‑1 (MQ=255)

                                                                                                                                                                                                                                            aaagaGGCTGTAACGCCATTTTACTGTC‑GCCCAGGGGGAACCCGGGGAGGCCGGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGgcg              <  2:277706/247‑1 (MQ=255)

                                                                                                                                                                                                                                               tAGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTg           >  1:257977/1‑251 (MQ=255)

                                                                                                                                                                                                                                                aGGCTGTAACGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGATTACCGCCGCGCAAACC‑GCCAACGCGGGATTTTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTg           >  1:436260/1‑250 (MQ=255)

                                                                                                                                                                                                                                                       aaCGCCATTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCcttc   >  2:123186/1‑251 (MQ=255)

CTACCGTTGCCGTCGGAGCTGAGGGTGACGCGTGCCAGCGGAATACCCGCCTGAACGGCGCGGGCAATACCTTCGGCAGGGGCGACCGGTTCGTCAATGCTGCTGGTGATATCGATGGTGCCGCCTTTGCGCGCGAACTCCAGCGCCTGCTCAAACA‑ACGGTACGTTGCGGTTA‑ACGTGGGTCG‑GCAGCAGCTTGCTGATCGGCACATC‑GCA‑GTT‑TT‑CCAGCAGGTCATAAATAGGCTGTAACGCCTTTTTACTGTC‑GCCCATGTGGAACACGGTGACGCCAGGTTTACCGCCGAGCAAACC‑GCCAACGCGGGATTCTGCCGCCATATTGGCCAGGTGATAAACGTCCGGTGCGGCAGAACGGTGATCAGAGATGGCGCATTTCACGCCAATCACACGATCGATAATCGCCACGTCTTTTTCCACGGAACCCGTAATGGTGCGGGAAGGGACATGATAAGCGCCGGTCAGCATCCAGGCGCTGATGCCTTCT  >  NC_000913/4558661‑4559152

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