BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,108,219:1	+A	100%	coding (1159/1395 nt)	wcaM ←	colanic acid biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,108,219	1	.	A	88.2%	704.7 / 92.8	211	coding (1159/1395 nt)	wcaM	colanic acid biosynthesis protein
Reads supporting (aligned to +/- strand): ref base . (11/14); new base A (96/90); total (107/104)
Fisher's exact test for biased strand distribution p-value = 5.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 217 total aligned reads displayed.

CTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAAT  >  CP009273/2108072‑2108363

cTTTTTCATTCACCGCATGAACATTGGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGGCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑a‑a                                                                                                                                                 <  2:901492/150‑1 (MQ=255)

 tttttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑a‑aa                                                                                                                                                >  2:3139354/1‑150 (MQ=255)

 tttttCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑a‑aa                                                                                                                                                <  2:3252245/150‑1 (MQ=255)

     tcattcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑a‑aaaa‑aa                                                                                                                                           >  2:291258/1‑150 (MQ=255)

       attcaCCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑a‑aaaa‑aaat                                                                                                                                         <  1:2882290/150‑3 (MQ=255)

                 tGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGTTATTTCGCATAA‑AAAA‑AATGTTGCGGttt                                                                                                                                >  2:4759431/1‑150 (MQ=255)

                    aCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTttgtt                                                                                                                             <  1:4360224/150‑1 (MQ=255)

                    aCATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTttgtt                                                                                                                             <  1:2190842/150‑1 (MQ=255)

                      aTTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTtgttgt                                                                                                                           >  1:1864621/1‑150 (MQ=255)

                        ttGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTgtgt                                                                                                                         <  2:2919836/150‑1 (MQ=255)

                             agagACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTc                                                                                                                    <  1:197129/150‑1 (MQ=255)

                              gagaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTcc                                                                                                                   <  2:3348973/150‑1 (MQ=255)

                              gagaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTcc                                                                                                                   >  1:4567072/1‑150 (MQ=255)

                              gagaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCa                                                                                                                  <  1:3858659/150‑1 (MQ=255)

                                gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                                 >  1:2645982/1‑150 (MQ=255)

                                gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                                 <  1:3470468/150‑1 (MQ=255)

                                gaCAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAg                                                                                                                 >  1:1967872/1‑150 (MQ=255)

                                   agcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                             >  1:2214492/1‑150 (MQ=255)

                                    gcagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                             <  1:3900306/150‑1 (MQ=255)

                                    gcagTGTTTCTTATGTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGac                                                                                                             <  2:2661173/150‑1 (MQ=255)

                                     cagTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacg                                                                                                            >  2:4892470/1‑150 (MQ=255)

                                        tgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcac                                                                                                         >  1:4689592/1‑150 (MQ=255)

                                        tgtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcac                                                                                                         <  1:4134735/150‑1 (MQ=255)

                                         gtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacg                                                                                                        >  2:570539/1‑150 (MQ=255)

                                         gtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacg                                                                                                        <  1:3511393/150‑1 (MQ=255)

                                         gtTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacg                                                                                                        >  2:566015/1‑150 (MQ=255)

                                          tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacgc                                                                                                       >  2:2145119/1‑150 (MQ=255)

                                          tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGacgcacgc                                                                                                       >  2:4410799/1‑150 (MQ=255)

                                          tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCt                                                                                                      >  2:3798046/1‑150 (MQ=255)

                                          tttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCt                                                                                                      >  2:4533446/1‑150 (MQ=255)

                                           ttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCt                                                                                                      <  1:3526903/150‑1 (MQ=255)

                                           ttCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCt                                                                                                      >  1:2970709/1‑150 (MQ=255)

                                              ttttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTGGTTGGGTAACTCCAGCGACGCACGCTTCa                                                                                                   >  2:499594/1‑150 (MQ=255)

                                               tttttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAt                                                                                                  >  2:286187/1‑150 (MQ=255)

                                                  ttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtctc                                                                                               >  2:5054745/1‑150 (MQ=255)

                                                  ttCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCAtctc                                                                                               >  2:3996496/1‑150 (MQ=255)

                                                       cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                          >  1:1355568/1‑150 (MQ=255)

                                                       cATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGGTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgt                                                                                          >  2:1726409/1‑150 (MQ=255)

                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtt                                                                                         >  1:3191166/1‑150 (MQ=255)

                                                        aTAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgtt                                                                                         <  2:4180538/150‑1 (MQ=255)

                                                             aaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtta                                                                                    >  2:1113607/1‑150 (MQ=255)

                                                             aaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATgttagtta                                                                                    >  2:472741/1‑150 (MQ=255)

                                                             aaCACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGttagttag                                                                                   >  2:3890876/1‑150 (MQ=255)

                                                              acacCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGGTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGgtagttag                                                                                   >  2:1748236/1‑150 (MQ=255)

                                                               caccacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgc                                                                                  >  1:2111414/1‑150 (MQ=255)

                                                               caccacGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAgc                                                                                  >  1:4697757/1‑150 (MQ=255)

                                                                    cGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcac                                                                             >  2:4280899/1‑150 (MQ=255)

                                                                     gAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcaca                                                                            <  1:738922/150‑1 (MQ=255)

                                                                     gAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCcaca                                                                            >  2:5082546/1‑150 (MQ=255)

                                                                      aaCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaa                                                                           >  2:4028673/1‑150 (MQ=255)

                                                                       aCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaaa                                                                          <  2:3630543/150‑1 (MQ=255)

                                                                       aCGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACaaa                                                                          <  1:3419691/150‑1 (MQ=255)

                                                                         gTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAgg                                                                        <  1:3329971/150‑1 (MQ=255)

                                                                            tttGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGaaa                                                                     <  1:5115540/150‑1 (MQ=255)

                                                                                cgcATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGc                                                                 <  2:2341874/150‑1 (MQ=255)

                                                                                 gcATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCa                                                                <  2:712182/150‑1 (MQ=255)

                                                                                     gTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTcc                                                            >  2:4637932/1‑150 (MQ=255)

                                                                                     gTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTcc                                                            >  1:139199/1‑150 (MQ=255)

                                                                                     gTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTcc                                                            <  1:2327817/150‑1 (MQ=255)

                                                                                     gTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTcc                                                            >  1:267996/1‑150 (MQ=255)

                                                                                       cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          >  1:1818261/1‑150 (MQ=255)

                                                                                       cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          >  1:1103909/1‑150 (MQ=255)

                                                                                       cAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCcg                                                          <  1:4869290/150‑1 (MQ=255)

                                                                                           gTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGga                                                     >  2:1089951/1‑150 (MQ=255)

                                                                                           gTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGga                                                     <  1:2076979/150‑1 (MQ=255)

                                                                                           gTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGga                                                     <  1:4663715/150‑1 (MQ=255)

                                                                                             tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    >  1:3741154/1‑150 (MQ=255)

                                                                                             tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    >  1:4683081/1‑150 (MQ=255)

                                                                                             tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgag                                                    <  1:478574/150‑1 (MQ=255)

                                                                                             tCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGgaga                                                   <  1:4871578/150‑1 (MQ=255)

                                                                                                 gCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGAttt                                                <  1:1285587/150‑1 (MQ=255)

                                                                                                  cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                               <  1:1182210/150‑1 (MQ=255)

                                                                                                  cTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTg                                               <  2:5136929/150‑1 (MQ=255)

                                                                                                    cAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            <  1:1894716/150‑1 (MQ=255)

                                                                                                     aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            <  2:365742/150‑1 (MQ=255)

                                                                                                     aaTGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAAt                                            >  2:4017538/1‑150 (MQ=255)

                                                                                                              ccAACTGAGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCaa                                   >  2:807337/1‑150 (MQ=255)

                                                                                                                    gAGGATTCCTGCATCACAGTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTtata                             >  2:1267793/1‑150 (MQ=255)

                                                                                                                     aGGATTCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATaa                            <  2:3593625/150‑1 (MQ=255)

                                                                                                                          tCCTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAg                       <  2:2658345/150‑1 (MQ=255)

                                                                                                                            cTGCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                     <  1:3931637/150‑1 (MQ=255)

                                                                                                                            cTGCATCACATTGATATTTCGCAT‑A‑AAAACAATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGAt                     <  2:1311042/150‑1 (MQ=255)

                                                                                                                              gCATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGa                   >  1:3235997/1‑150 (MQ=255)

                                                                                                                               cATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAg                  <  2:1818261/150‑1 (MQ=255)

                                                                                                                               cATCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAg                  <  1:4914440/150‑1 (MQ=255)

                                                                                                                                aTCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 <  1:2548950/150‑1 (MQ=255)

                                                                                                                                aTCACATTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 <  1:2439402/150‑1 (MQ=255)

                                                                                                                                aTCACATTGATATTTCGCAT‑ATAAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGt                 <  1:4469897/150‑1 (MQ=255)

                                                                                                                                 tCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtat               <  2:3598121/150‑1 (MQ=255)

                                                                                                                                  cacaTTGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGtat               >  1:651882/1‑150 (MQ=255)

                                                                                                                                       tGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc          <  2:3102363/150‑1 (MQ=255)

                                                                                                                                       tGATATTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATcc          >  2:2568599/1‑150 (MQ=255)

                                                                                                                                         atatTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAg        <  2:2002932/150‑1 (MQ=255)

                                                                                                                                          tatTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGt       >  1:3089811/1‑150 (MQ=255)

                                                                                                                                           atTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt      <  2:4962259/150‑1 (MQ=255)

                                                                                                                                           atTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt      >  2:2370312/1‑150 (MQ=255)

                                                                                                                                           atTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt      >  1:3903986/1‑150 (MQ=255)

                                                                                                                                           atTTCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGtt      >  2:1384490/1‑150 (MQ=255)

                                                                                                                                            tttCGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTg     <  2:4765670/150‑1 (MQ=255)

                                                                                                                                               cGCATAA‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAat  <  1:1624347/150‑1 (MQ=255)

CTTTTTCATTCACCGCATGAACATTTGCAAGAGACAGCAGTGTTTCTTTTTTCGCCATAAAAACACCACGAACGTCTTTGCGCATGTCAAAGTTCATGCTCAATGCGGGTCCAACTGAGGATTCCTGCATCACATTGATATTTCGCAT‑A‑AAAA‑AATGTTGCGGTTTGTTGTGTAACTCCAGCGACGCACGCTTCATCTCAATGTTAGTTAGCGCCACAAAGGAAACAGCATTCCCGGCGGAGATTTGAATGCCGCGCAATTTATAAGCAAGATGAGTATTATCCAGTTGAAT  >  CP009273/2108072‑2108363

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