BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,269,897	T→G	3.7%	noncoding (40/66 nt)	rdlB →	sRNA antisense regulator affects LdrB translation, proposed addiction module in LDR‑B repeat, with toxic peptide LdrB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,269,897	0	T	G	3.7%	485.0 / 5.9	147	noncoding (40/66 nt)	rdlB	sRNA antisense regulator affects LdrB translation, proposed addiction module in LDR‑B repeat, with toxic peptide LdrB
Reads supporting (aligned to +/- strand): ref base T (59/80); new base G (4/3); total (63/84)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 148 total aligned reads displayed.

AAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCT  >  NC_000913/1269623‑1270182

aaaaTGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTg                                                                                                                                                                                                                                                                                   <  2:1369932/289‑1 (MQ=34)

       ccGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTTCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGt                                                                                                                                                                                                                                                                            >  1:1166128/1‑289 (MQ=21)

               tgtttagaagagacagggTTTCGGGCGCATATTGCCATCACGGCACGCTGACGCCCGTTTTCACCGGACTTCCGGTTAAGCCACCAGCCGACAATCGCTGCGATAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCATAAAGTCATGGCAAACTGCGCGAGTGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTctct                                                                                                                                                                                                                                                                     <  1:445214/272‑1 (MQ=11)

                           ggCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                                                                                                                                                                                                                                                         >  2:875436/1‑272 (MQ=34)

                           ggCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGtttt                                                                                                                                                                                                                                                                         <  1:875436/272‑1 (MQ=34)

                             ccGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                                                                                                                                                                                                                                                       <  1:1433177/288‑1 (MQ=34)

                                gTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcc                                                                                                                                                                                                                                                    >  1:922802/1‑288 (MQ=34)

                                   tcaacctcttacagCCTTCACGCCGTCCCGACTCCCGTTTTCACCTTACTTCCGGTTACGCAACAAGCCGAAAATCGCTGCGGTAATAATCCCCGCAAGGATCGGTGCCGCCACGCCGTGCCAGAAATTCATGGCAAAATGCGCGAGCTTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTT‑TCTGGTTTCAAGATTTGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca‑cca                                                                                                                                                                                                                                               <  1:1094700/276‑1 (MQ=11)

                                   tCGTGCGCATATAGCCATCCCGGCAGCCGGACGCCCGTTTTCACCTTACTTCCGGGTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGGGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca‑cca                                                                                                                                                                                                                                               <  1:1354599/289‑1 (MQ=255)

                                   tCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca‑cca                                                                                                                                                                                                                                               <  2:26063/289‑1 (MQ=34)

                                          cATATTGCCATCACGGCAGCCTGCCGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATCATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCACCGTTCGTGGGTTTTCTCTTTCCAGCAAc                                                                                                                                                                                                                                                           >  2:445214/1‑271 (MQ=19)

                                           atatTGCTATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCCGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATaaa                                                                                                                                                                                                                                       >  2:11347/1‑289 (MQ=21)

                                           atatTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACCATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAACCTGCGCGAGCGTCATATAGCCGCCTGTTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTCCACGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCACCGTGCGGGGTTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGCTaaa                                                                                                                                                                                                                                       >  1:1472734/1‑289 (MQ=255)

                                             atTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGc                                                                                                                                                                                                                                     <  2:261959/289‑1 (MQ=34)

                                             atTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGc                                                                                                                                                                                                                                     <  2:307432/289‑1 (MQ=34)

                                               tGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAAAATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTCACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAGCCCA‑CAGGGGATAAAGct                                                                                                                                                                                                                                    >  1:123066/1‑287 (MQ=255)

                                                       tctgaGCCTAGCGCCCGGTGTCCCCGTACTTCCTGTTACGCCCCCAGCCGACAAACGCAGCGGTAAGCATTCCCGCCGGGAGCGGTGAAGCCAGGGTGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTAATTTAGCTGCCTTGTTGTAATGACAACATTTGGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGGTCTGTTGACAGGTTCTACCTCTCTATGGTCGGGGGTGTTAGCTTTCTAGAAGAGACTGCCA‑CCAGGGATAAAGCCCCCGCAACa                                                                                                                                                                                                                           <  1:1112226/285‑1 (MQ=2)

                                                          cAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGCTAAAGCCCCCGCAACAtt                                                                                                                                                                                                                         >  2:277928/1‑288 (MQ=37)

                                                           aGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCGGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTg                                                                                                                                                                                                                        <  2:1384777/288‑1 (MQ=37)

                                                            cccTGCCGCCCGTTTTCCCCTTACTTCCGGTTCCGCCACCGGCCGAAAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATTGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTgcg                                                                                                                                                                                                                      <  1:1434988/288‑1 (MQ=255)

                                                                  cGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                                                                                                                                                                                                                 <  2:572896/288‑1 (MQ=37)

                                                                  cGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCa                                                                                                                                                                                                                 <  1:392013/288‑1 (MQ=37)

                                                                   gCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCTGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCAcc                                                                                                                                                                                                               <  2:642140/289‑1 (MQ=17)

                                                                    cccGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCg                                                                                                                                                                                                              <  1:596586/289‑1 (MQ=37)

                                                                      cGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGa                                                                                                                                                                                                            <  2:591820/289‑1 (MQ=37)

                                                                               cTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGGTTTCTCTTTCCAGCAACCAATCCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCg                                                                                                                                                                                                    >  2:1254180/1‑288 (MQ=255)

                                                                                  aCTTCCGGTTACGCCACCAGCCGAAAATCGCTGCGGGAAAAATTCCCGCCAGGATCGGGGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                                                                                                                                                                                <  2:123066/289‑1 (MQ=255)

                                                                                     tCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                                                                             <  1:728051/289‑1 (MQ=39)

                                                                                     tCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                                                                             <  1:920168/289‑1 (MQ=39)

                                                                                     tCCGGTTACGACACCTGACGACAATGGGTGCGGTAATAATTCCCGCCTGGAGCGTGGACTGCAGGTGGTGCCAGAAAGGCATTGCAAACTGCGCGGGAGTCATATATCCGGCTTTTTTTAATGTCAAGTTGTTGCGGTTATTCTTTATGTG‑TCTTGTTTGAAGATTCGCCTCCGTTTTGTTGTCAGGTTGTAGCTCGCAACGAGCGGGAGTCTTCTCTTTCCAGCAA‑CAGCGTCAGATAGGTATAAAAGCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                                                                             <  1:1166808/289‑1 (MQ=255)

                                                                                      ccGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgtg                                                                                                                                                                                            <  2:169401/289‑1 (MQ=39)

                                                                                         gTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                                                                             <  2:246577/285‑1 (MQ=39)

                                                                                         gTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGgt                                                                                                                                                                                             >  1:246577/1‑285 (MQ=39)

                                                                                         gTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTgg                                                                                                                                                                                         >  2:240646/1‑289 (MQ=39)

                                                                                         gTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTgg                                                                                                                                                                                         >  1:1551585/1‑289 (MQ=39)

                                                                                              gCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAAAATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGTAACCAATGCCA‑CCAGGTATAAAGCCCCCGCCACATTGCGCCTAACCGGAAAATTCCGGCTTGGTGTGGAtac                                                                                                                                                                                     >  1:1256580/1‑288 (MQ=255)

                                                                                                caccaGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGAtacta                                                                                                                                                                                   <  1:258033/288‑1 (MQ=39)

                                                                                                aaccaGCGGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCTGCGAGGTCTGGCCAGAAAGTCATTGCAAAATGCTCGAGCGTCATATAGCCGCCTTGTTTTAATGACAACATTTTGCGGCTATTCTTGAATTG‑TCTGGTTTCAAGATTAGCCCACGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGAtacta                                                                                                                                                                                   <  1:775978/287‑1 (MQ=11)

                                                                                                         acactgccTGGGGTAATAATTCCCGCCCGGATCGGGGGCGACATGTCGTGCCAGAAAGTCATGGCACACTGCGCGAGCGTCATATAGCCGCCTTTTGCGAATGACAACATTTTGCGTCTATTTTTTAAGGT‑TCTGGTTTCAAAATTAGCCCCCGTTCCGTTGTCAGGTGTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGCGATAAATCCCCCGATACATTGCGCCTCATCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAAt                                                                                                                                                                         <  1:355485/281‑1 (MQ=14)

                                                                                                             tCGCTGAGGTAAAAATTCCCGCCAGGATCGGGGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTGTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATc                                                                                                                                                                    <  1:506572/289‑1 (MQ=255)

                                                                                                                cTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGCTCTGTTGTGAGGTTTTACCTCTCAACGTGCGGGGGTTTTTTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGc                                                                                                                                                                                                                                     >  2:1085221/1‑222 (MQ=21)

                                                                                                                cTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATAGCAAACTGCGCGAGCGTCCTATAGCCGCCTTGTTGTAATGACCAACTTTTGCGGCTCTTCTTGCCTTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTCGTCAGGTTGTACCCATCAA‑GTGCGGGGGTTTTCTCTTTCCAGCAACCAATGACA‑CCAGGGATCACGCCCCCACAACCTTGCGTCTCACCGGATAATTCAGGATTGGTGTGGCTACTCC‑GTCTCAATTCCActtca                                                                                                                                                                >  1:26343/1‑289 (MQ=11)

                                                                                                                       aataatTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATcttcacttca                                                                                                                                                           <  2:1140743/288‑1 (MQ=39)

                                                                                                                              cccGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGaa                                                                                                                                                   >  1:1176309/1‑289 (MQ=39)

                                                                                                                               ccGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCg                                                                                                                                                                                                                                                                                 >  2:1424293/1‑164 (MQ=35)

                                                                                                                                cGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTTTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAAt                                                                                                                                                 >  2:1307537/1‑289 (MQ=39)

                                                                                                                                   cAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGt                                                                                                                                               <  2:1176309/288‑1 (MQ=39)

                                                                                                                                   cAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGt                                                                                                                                               <  2:922802/288‑1 (MQ=39)

                                                                                                                                        tCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                                                                          <  1:1339160/288‑1 (MQ=39)

                                                                                                                                        tCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATGTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                                                                          <  2:294925/288‑1 (MQ=39)

                                                                                                                                           gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGa                                                                                                                                                                                        >  1:34791/1‑240 (MQ=40)

                                                                                                                                           gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATaa                                                                                                                                      >  2:866680/1‑289 (MQ=40)

                                                                                                                                           gTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTAGTGTGGa                                                                                                                                                                                        <  2:34791/240‑1 (MQ=255)

                                                                                                                                            tGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCa                                                                                                                                                                                                                                                                <  1:673415/168‑1 (MQ=35)

                                                                                                                                                gccAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTAATCTTCCCTTCATCCCTGAAATGCTTGCAATAAAGAGt                                                                                                                                 >  1:117771/1‑289 (MQ=255)

                                                                                                                                                                        gcaaatgtctccggCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTTGTGTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCTGTATTCAAAGCCA‑CCAGGTATAAATCCACAGCAACATACCGCCTCACCTGGGAAGCCTGGGATGGGTAGAATAGTCC‑GACTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTg                                                                                                          <  1:1063365/275‑1 (MQ=11)

                                                                                                                                                                         cAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCACGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAACCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGtt                                                                                                        >  2:959871/1‑289 (MQ=12)

                                                                                                                                                                            aCTGCGCGAGCGTCATATAGCTGCCTTGTTGTAATGACAACATTTTGCGGATATTCTAGAAGTG‑TCTGGTTTCAAGATTAGCCCCAGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTATTTCCAGCAACCAATg                                                                                                                                                                                                                                                      >  1:409328/1‑146 (MQ=11)

                                                                                                                                                                             cTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATcttcacttca                                                                                                                                                           >  2:760598/1‑234 (MQ=38)

                                                                                                                                                                                  cgAGCGTCATATGGCGGCCTTGTTGTAATGGCAACATTTTGCGGCTATTCTTGAAGTG‑TCTGTTTGCAAGAGTAGCCCCCGTTCTGTTGTAAGGTTTTACCTCTCCACGTGCGGGGGTTTTCTCTTTCCTGCAACCAATGCCA‑CCAGTGATAAGGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCGATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGt                                                                                                <  1:1254180/288‑1 (MQ=14)

                                                                                                                                                                                         cATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGTTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGTGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTGTCAACAGCTGTTGCAGTCGTTTCATGa                                                                                         <  2:1472734/288‑1 (MQ=255)

                                                                                                                                                                                         cATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGa                                                                                         <  2:1166128/288‑1 (MQ=40)

                                                                                                                                                                                         cATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGa                                                                                         <  2:60964/288‑1 (MQ=40)

                                                                                                                                                                                         cATATAGCCGCATTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGa                                                                                         <  2:1019352/288‑1 (MQ=255)

                                                                                                                                                                                             taGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                                                                    <  2:1297455/289‑1 (MQ=255)

                                                                                                                                                                                             taGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                                                                    <  2:1551585/289‑1 (MQ=255)

                                                                                                                                                                                             taGCCGCCTTGTTGTAATGAAAACATTTTTCGGCTATTCTTGAAGTG‑TCTGTGTTTAATATTGTCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTCCGGGGGTTTTCTCTTTCCAGCAACCAGTGCCA‑CTAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGct                                                                                    <  2:1482846/289‑1 (MQ=14)

                                                                                                                                                                                                 cgccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCCGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGGGTGCTCTgg                                                                                <  2:736491/289‑1 (MQ=255)

                                                                                                                                                                                                   ccTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgat                                                                              >  2:322134/1‑289 (MQ=255)

                                                                                                                                                                                                     ttgttgTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGgatga                                                                            <  1:240646/289‑1 (MQ=255)

                                                                                                                                                                                                        ttgtAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TTTTGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGcca‑cc                                                                                                                                                                                                                                                <  1:1021022/123‑1 (MQ=17)

                                                                                                                                                                                                             aTGACAACATTTTGCGGCTATTCTTGAATTGTTCTGG‑TTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCa                                                                    >  2:1113805/1‑289 (MQ=16)

                                                                                                                                                                                                              tGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCATTAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTgg                                                                                >  2:62598/1‑276 (MQ=255)

                                                                                                                                                                                                              tGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAg                                                                   <  2:1338605/289‑1 (MQ=255)

                                                                                                                                                                                                                acaacaTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTAGCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCACTTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTTAACAGCTGTTGCAGTCGTTTCATGCGTGCCCTGGATGATGCTTCCAgt                                                                  >  1:322681/1‑287 (MQ=255)

                                                                                                                                                                                                                 aacatttttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTTTGGATACTACTGTC‑CAATTCAACTTCACTCCATCCCTGAAATGCTTGCAAGAAAGAGTACATTCTGGCTTTTCAACAACTATTGAAGTCGTTTCATGATTGCTCTTGATGATGCTTCCAGCTc                                                                >  2:147050/6‑289 (MQ=14)

                                                                                                                                                                                                                    caTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCggg                                                             <  2:1002304/289‑1 (MQ=255)

                                                                                                                                                                                                                     aTTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTTCAGTCGTTTCATGAGTCCTCTGGATGATGCTTACCGCTCGGGt                                                            >  1:1537586/1‑289 (MQ=255)

                                                                                                                                                                                                                       tttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtt                                                                                                                     >  2:1279112/1‑230 (MQ=40)

                                                                                                                                                                                                                       tttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCtt                                                                                                                     <  1:1279112/230‑1 (MQ=40)

                                                                                                                                                                                                                        ttGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGc                                                         <  1:322134/289‑1 (MQ=255)

                                                                                                                                                                                                                          gCGGCTATTCGTGAATTG‑TTTTGGTTCAAGATTAGCCCCCGTGCGGTTGTAAGGTTTTACCTCTCAACGTACGGGGGTTTTCTCTTTCCAGCAACCAATGACA‑CCAGGGATAAAACCCCCGCAACATTGCGCATCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTgaaa                                                       <  1:959871/289‑4 (MQ=11)

                                                                                                                                                                                                                                           g‑tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGaa                                       <  2:1537586/288‑1 (MQ=255)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAg                                                                                                                                                                                                                                      >  1:28590/1‑97 (MQ=33)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                                                                          <  2:1543140/188‑1 (MQ=40)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCa                                                                                                                                          >  1:1543140/1‑188 (MQ=40)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAagag                                                                                                                                  >  2:1111368/1‑196 (MQ=255)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAagag                                                                                                                                  <  1:1111368/196‑1 (MQ=255)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTAc                                               <  1:133810/279‑1 (MQ=255)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAaga                                     >  1:1302735/1‑289 (MQ=255)

                                                                                                                                                                                                                                             tCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAaga                                     >  2:1337711/1‑289 (MQ=255)

                                                                                                                                                                                                                                              cTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTAGTTTCATGAGTGCTTTGGATGATGCTTCCATCTCGGGTTGCCAATATTTCCTTGTGGAagag                                    >  1:438129/1‑289 (MQ=255)

                                                                                                                                                                                                                                                       aGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATaa                                                                                                                                      <  2:699496/182‑1 (MQ=21)

                                                                                                                                                                                                                                                       aGATTAGCCCCCGTTCTGTTGTCAGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAACTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATaa                                                                                                                                      >  1:699496/1‑182 (MQ=14)

                                                                                                                                                                                                                                                                  cGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCg                 >  1:737726/1‑288 (MQ=255)

                                                                                                                                                                                                                                                                  cGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTAATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGCATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGc                >  2:72148/1‑289 (MQ=255)

                                                                                                                                                                                                                                                                   gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                                                                                                                                                                                                                                                       >  2:525146/1‑59 (MQ=35)

                                                                                                                                                                                                                                                                   gTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAAt                                                                                                                                                                                                                                                       <  1:525146/59‑1 (MQ=35)

                                                                                                                                                                                                                                                                    ttCTGTTTTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAg              <  1:620173/289‑1 (MQ=255)

                                                                                                                                                                                                                                                                              aGGTTGTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCtt                                                                                                                                                                                                >  1:836013/1‑102 (MQ=17)

AAAATGGCCGCAGGCTGTTACCCCTGCGGCCGGTTTCGGGCGCATATTGCCATCACGGCAGCCTGACGCCCGTTTTCACCTTACTTCCGGTTACGCCACCAGCCGACAATCGCTGCGGTAATAATTCCCGCCAGGATCGGTGCCGCCAGGTCGTGCCAGAAAGTCATGGCAAACTGCGCGAGCGTCATATAGCCGCCTTGTTGTAATGACAACATTTTGCGGCTATTCTTGAAGTG‑TCTGGTTTCAAGATTAGCCCCCGTTCTGTTGTCAGGTTTTACCTCTCAACGTGCGGGGGTTTTCTCTTTCCAGCAACCAATGCCA‑CCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGATAATTCCGGCTTGGTGTGGATACTAC‑GTCTCAATTCATCTTCACTTCATCCCTGAAATGTTTGCAATAAAGAGTACATTCCGGCTTTTCAACAGCTGTTGCAGTCGTTTCATGAGTGCTCTGGATGATGCTTCCAGCTCGGGTTGCCAATATTTACTTGTGGAAGAGATAAAGACAAAAATGGCCGCAGGCTGTTACCCCT  >  NC_000913/1269623‑1270182

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