BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,407,295	T→A	6.0%	intergenic (+29/‑80)	accC → / → yhdT	acetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,407,295	0	T	A	6.0%	685.7 / 26.0	200	intergenic (+29/‑80)	accC/yhdT	acetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein
Reads supporting (aligned to +/- strand): ref base T (94/94); new base A (4/8); total (98/102)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-02

Only 100 of 203 total aligned reads displayed.

TCTACGCGGGCTACACCGTACCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGTTTGAGATGGCCTGCATCCTGACGCCGCTGCTGTTTA  >  NC_000913/3407028‑3407546

tCTACGCGGGCTACACCGTACCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGcc                                                                                                                                                                                                                                                       <  4:182627/275‑1 (MQ=255)

    cgcgGGCTACACCGTACCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCtttt                                                                                                                                                                                                                                                   >  4:570103/1‑275 (MQ=255)

                cGTACCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgggg                                                                                                                                                                                                                                       <  2:439057/275‑1 (MQ=255)

                   aCCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTTATCTGCAGCTCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACCTCCACTATCTGGAGACAANACTCGGACTTCAGGAAACATAAGACTGCAAAAG‑CGTCACAAGGCCAGATTTTCCAGACCTTTttattattacgtagc                                                                                                                                                                                                                                    >  3:160010/1‑266 (MQ=255)

                   aCCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCCTGGTGGCACTACCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATc                                                                                                                                                                                                                                    >  2:230653/1‑275 (MQ=255)

                   aCCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAGGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATc                                                                                                                                                                                                                                    >  2:387044/1‑275 (MQ=255)

                                 aCTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCATTTTTTATTACTGGGGATCG‑ACAACACCCAt                                                                                                                                                                                                                       >  1:17827/1‑273 (MQ=255)

                                     aaTGATCGGTAAGCTGATTTGCTACGGTGACAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAAATCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCAGGATTTTCCGGCCTTTTTTATTACTGGGGAACGCACAA‑CCCCATAAgg                                                                                                                                                                                                                   >  2:527070/1‑273 (MQ=255)

                                     aaTGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAgg                                                                                                                                                                                                                   >  2:259396/1‑273 (MQ=255)

                                                            aCGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                                                                                                                                                                                     >  4:645677/1‑233 (MQ=255)

                                                                                 cGATTGCACGCATGAAGAATGCGCTGCAGGAGCTGATCATCGTCGGTAGCAAAACCAACGTTGATCTGCAGAGCAGCATCATGAATGACGAGAACTGCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGacac                                                                                                                                                                      <  2:543706/274‑1 (MQ=255)

                                                                                                                        tcGACTGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAtttt                                                                                                                                                                                                                                                             <  6:426822/149‑1 (MQ=255)

                                                                                                                        tcGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAtttt                                                                                                                                                                                                                                                             >  6:1055154/1‑149 (MQ=255)

                                                                                                                        tcGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTttat                                                                                                                                                                                                                                               >  4:736610/1‑163 (MQ=255)

                                                                                                                        acGACGGTATTAAAACCAAAGTTGATCTGCAGATCCGAATACTGAATGACGAGAACTTACAGCATGGGGGCACGAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgct                                                                                                                               <  4:259396/273‑1 (MQ=255)

                                                                                                                          gACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTcc                                                                                                                                                                                                                                                           <  5:301445/149‑1 (MQ=255)

                                                                                                                              gTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCa                                                                                                                                                                                                                        <  2:638270/179‑1 (MQ=255)

                                                                                                                               tATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCt                                                                                                                                                                                                                                                      <  6:979915/149‑1 (MQ=255)

                                                                                                                               tATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACCCGTTTTGTTCAGGCCCATCAAAAGGCGGGCTGGGCGCTGGGGCGGa                                                                                                                       >  4:327843/1‑275 (MQ=255)

                                                                                                                               tATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAAAATCCGGCCt                                                                                                                                                                                                                                                      <  6:76401/149‑1 (MQ=255)

                                                                                                                                aTCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCtt                                                                                                                                                                                                                                                     <  6:98522/149‑1 (MQ=255)

                                                                                                                                aTCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAACATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTGCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGAGCTGGGGCTGa                                                                                                                       <  4:230653/274‑1 (MQ=255)

                                                                                                                                 tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                                                                                                                                                                                     >  6:1077660/1‑148 (MQ=255)

                                                                                                                                  cAAAACCAACGATGATCTGCAGATCCGCATGATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGCAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAAAATCCGGCCtttt                                                                                                                                                                                                                                                   <  8:903502/149‑1 (MQ=255)

                                                                                                                                   aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGCGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCttttt                                                                                                                                                                                                                                                  >  8:510238/1‑149 (MQ=255)

                                                                                                                                   aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCttttt                                                                                                                                                                                                                                                  >  5:926732/1‑149 (MQ=255)

                                                                                                                                   aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCttttt                                                                                                                                                                                                                                                  >  5:597279/1‑149 (MQ=255)

                                                                                                                                   aaaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCttttt                                                                                                                                                                                                                                                  >  8:605058/1‑149 (MQ=255)

                                                                                                                                    aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCtttttt                                                                                                                                                                                                                                                 >  6:270543/1‑149 (MQ=255)

                                                                                                                                    aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCtttttt                                                                                                                                                                                                                                                 >  8:454131/1‑149 (MQ=255)

                                                                                                                                    aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCtttttt                                                                                                                                                                                                                                                 >  7:11988/1‑149 (MQ=255)

                                                                                                                                         aaCGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAANTTCCGGCCNTTTTTATTAc                                                                                                                                                                                                                                            >  3:43682/1‑149 (MQ=255)

                                                                                                                                             ttGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTggg                                                                                                                                                                                                                                        <  8:165965/149‑1 (MQ=255)

                                                                                                                                                aTCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                                                                                                                                                                                     >  5:556881/1‑149 (MQ=255)

                                                                                                                                                 tCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTAGCACTAACATCCACTATCTGGAGAAAAAACTCGGTTTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                                                                                                                                                <  7:512643/137‑1 (MQ=255)

                                                                                                                                                     cAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGCTTTTCCGGCCTTTTTTATCACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGCACAAAAAATGTACCCCCGTTTTGTTCAGGCCCCTTAACAGGCCCACTTGGGCCTCGGGCTGACCCTTTTGGTTCTTGCAGttgg                                                                                                 >  4:484773/1‑273 (MQ=255)

                                                                                                                                                        aTCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAAcc                                                                                                                                                                                                                            >  6:163299/1‑149 (MQ=255)

                                                                                                                                                           cGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAAccccc                                                                                                                                                                                                                         >  7:1020326/1‑149 (MQ=255)

                                                                                                                                                           cGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAAccccc                                                                                                                                                                                                                         <  7:605058/149‑1 (MQ=255)

                                                                                                                                                           cGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGGCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAAccccc                                                                                                                                                                                                                         <  5:852722/149‑1 (MQ=255)

                                                                                                                                                            gCATCATGAATGACGAGAACGTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTtagt                                                                                           <  2:559375/274‑1 (MQ=255)

                                                                                                                                                             catcatGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTg                                                                                                                                                              >  1:88185/1‑206 (MQ=255)

                                                                                                                                                              atcatGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATa                                                                                                                                                                                                                      >  7:965268/1‑149 (MQ=255)

                                                                                                                                                               tcatGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGAAGTTTGGTTAGTAGcc                                                                                       >  2:437367/1‑275 (MQ=255)

                                                                                                                                                                 atgaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAgg                                                                                                                                                                                                                   <  7:1025659/149‑1 (MQ=255)

                                                                                                                                                                   gaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTa                                                                                                                                                                                                                 >  7:332143/1‑149 (MQ=255)

                                                                                                                                                                    aatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACACNCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGAACACTCGTTTTGTTCAGGCCCATAAANAGGCGCGCAGGGCGGGGGGGCTGACCCTTTNNTACCGGGCGGTTTGGTtagtgacgag                                                                                     >  3:77711/1‑266 (MQ=255)

                                                                                                                                                                       gaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGctta                                                                                  <  2:327843/272‑1 (MQ=255)

                                                                                                                                                                        aCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATc                                                                                                                                                                                                            <  5:1040215/149‑1 (MQ=255)

                                                                                                                                                                             aaCTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGc                                                                                                                                                                                                       <  5:1055154/149‑1 (MQ=255)

                                                                                                                                                                               cTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCtt                                                                                                                                                                                                     >  7:867255/1‑149 (MQ=255)

                                                                                                                                                                                      cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAc                                                                                                                                                                                              <  6:639400/149‑1 (MQ=255)

                                                                                                                                                                                           tggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATc                                                                                                                                                                                         >  7:940642/1‑149 (MQ=255)

                                                                                                                                                                                            ggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGAGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCa                                                                                                                                                                                        <  8:699044/149‑1 (MQ=255)

                                                                                                                                                                                               aCTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAggg                                                                                                                                                                                     >  7:651499/1‑149 (MQ=255)

                                                                                                                                                                                                      tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaa                                                                                                                                                                              >  5:135001/1‑149 (MQ=255)

                                                                                                                                                                                                            agggggAAAAAAACCTCGGTCNNCAGGANAAAT‑AGGCTTCTGAAGTCTTCGAAAGGCCCGATTTTCCGGGCTTCTGTCTTTCTGTGGGTCG‑ACAACCCCCATAAGGTGCAATCCCGGCTTTCTTAACCCCTCAGGGACAAAAAGTGGACACTCGTTTTGTTCAGGCCCATAAAGGGGCGCGCTGGGCGCTGAGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCcgcg                                           <  3:17827/270‑1 (MQ=255)

                                                                                                                                                                                                             tCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGCTCG‑ACAACCCCCATAAGGTACAACACCCGCTTTCTTCACCCCTCCGGGACCAAAAATGCACACTCGTTTTGTTCAGACCCATAAAGTGGCGCACTTGGAGATGGGGCGGAc                                                                                                                      >  1:117623/1‑198 (MQ=255)

                                                                                                                                                                                                             tCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGaca                                                                                                                                                                       >  5:266721/1‑149 (MQ=255)

                                                                                                                                                                                                               tGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACt                                                                                                                                                                     >  2:727465/1‑149 (MQ=255)

                                                                                                                                                                                                               tGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACt                                                                                                                                                                     >  6:961384/1‑149 (MQ=255)

                                                                                                                                                                                                               tGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgc                                                                                                                                <  4:402593/186‑1 (MQ=255)

                                                                                                                                                                                                                ggAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTc                                                                                                                                                                    <  7:696277/149‑1 (MQ=255)

                                                                                                                                                                                                                  agaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGtt                                     >  2:566590/1‑274 (MQ=255)

                                                                                                                                                                                                                   gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAGG‑CTTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTACGTCACCCATCAGGGACAAAAAATGGACACACGtt                                                                                                                                                                 >  6:534652/1‑149 (MQ=255)

                                                                                                                                                                                                                    aaaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGttt                                                                                                                                                                <  6:549308/149‑1 (MQ=255)

                                                                                                                                                                                                                    aaaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                                                                                                                                                                                     >  5:746088/1‑65 (MQ=38)

                                                                                                                                                                                                                    aaaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                                                                                                                                                                                     <  6:746088/65‑1 (MQ=37)

                                                                                                                                                                                                                     aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtttt                                                                                                                                                               <  6:235320/149‑1 (MQ=255)

                                                                                                                                                                                                                      aaaaCGCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTg                                                                                                                        <  1:69634/187‑1 (MQ=255)

                                                                                                                                                                                                                        aaCTGGGTCTTCAGAAAAAATAAGACTGCTAAAG‑CGTAAAAAGGCCGGATTTTCCGGCCTTTTTTATTACT‑GGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTc                                                                                                                                                           <  6:926732/149‑1 (MQ=255)

                                                                                                                                                                                                                                cTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCa                                                                                                                                                    <  8:639883/149‑1 (MQ=255)

                                                                                                                                                                                                                                 ttCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCTTTTTGTTCAGGCCCAt                                                                                                                                                   >  8:363441/1‑149 (MQ=255)

                                                                                                                                                                                                                                    aGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAggg                                                                                                                                                                                     >  4:547995/1‑112 (MQ=255)

                                                                                                                                                                                                                                    aGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATaaa                                                                                                                                                <  6:884456/149‑1 (MQ=255)

                                                                                                                                                                                                                                      gAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaa                                                                                                                                                                              <  5:544230/117‑1 (MQ=255)

                                                                                                                                                                                                                                      gAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAaga                                                                                                                                              >  6:565653/1‑149 (MQ=255)

                                                                                                                                                                                                                                      gAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAaga                                                                                                                                              >  6:307978/1‑149 (MQ=255)

                                                                                                                                                                                                                                        aaaaTAAGACTGCTAAAG‑CGGCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGTTTGAGATGGCCTGCATCCTGAcgc              <  1:133964/275‑1 (MQ=255)

                                                                                                                                                                                                                                          aaTAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTCCAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGTTTGAGATGGCCTGCATCCTGAcgccg            <  2:15722/275‑1 (MQ=255)

                                                                                                                                                                                                                                           aTAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGgcgc                                                                                                                                         >  7:177584/1‑149 (MQ=255)

                                                                                                                                                                                                                                                aCTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTGCGCGCTGGTTTGAGATGGCCTGCATCCTGACGCCGcggctg      >  2:3137/1‑275 (MQ=255)

                                                                                                                                                                                                                                                  tGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTCTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggc                                                                                                                                  >  8:1022766/1‑149 (MQ=255)

                                                                                                                                                                                                                                                    cTAAAG‑CGTCAAAATGCCGGATTTTCCAGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATACTGTACAATCCCCGATTTCTTCACGCATCAGGCACAAAAAATGGACCCTCGTGTTGGTCAGGCTCATAAAGAGGCgcgctgggcgc                                                                                                                                >  8:923177/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     tAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTGTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCAATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGTTTGAGATGGCCTGCATCCTGACGCCGCTGCTGTTTa  <  2:92628/274‑1 (MQ=255)

                                                                                                                                                                                                                                                        aG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGtt                                                                                                                                                            <  2:180/117‑1 (MQ=255)

                                                                                                                                                                                                                                                            gTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCt                                                                                                                         >  8:562707/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               aaaaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGTTTGAGATGGCCTGCATCCTGACGCCGctgctg      >  4:163875/1‑261 (MQ=255)

                                                                                                                                                                                                                                                                aaaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGAcc                                                                                                                     <  7:768664/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt                                                                                                                                                                 <  3:28922/104‑1 (MQ=255)

                                                                                                                                                                                                                                                                   ggCCGGATTTACCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATa                                                                                                                                                                                                                      <  1:41337/49‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCTCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt                                                                                                                <  5:435864/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGGCAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt                                                                                                                <  8:463581/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt                                                                                                                <  8:177584/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt                                                                                                                <  6:673617/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTgg                                                                                                >  4:148876/1‑165 (MQ=255)

                                                                                                                                                                                                                                                                      cGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTg                                                                                                               >  5:66992/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       ggATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGt                                                                                                              >  7:816033/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                         aTTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTtag                                                                                            <  2:148876/165‑1 (MQ=255)

                                                                                                                                                                                                                                                                          ttttCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATc                                                                                                           <  8:12681/149‑1 (MQ=255)

TCTACGCGGGCTACACCGTACCGCCGTACTATGACTCAATGATCGGTAAGCTGATTTGCTACGGTGAAAACCGTGACGTGGCGATTGCCCGCATGAAGAATGCGCTGCAGGAGCTGATCATCGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAG‑CGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCG‑ACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGTATCTGGCAGTTTGGTTAGTAGCCGCTTACTTATCTGGCGTTGCCCCCGGTTTTACCGGCTTTCCGCGCTGGTTTGAGATGGCCTGCATCCTGACGCCGCTGCTGTTTA  >  NC_000913/3407028‑3407546

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 11 ≤` `ATCG/ATCG` `< 19 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`