BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,893,446	T→A	42.2%	K63* (AAA→TAA)	ygcW ←	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,893,446	0	T	A	42.2%	174.6 / inf	268	K63* (AAA→TAA)	ygcW	putative dehydrogenase
Reads supporting (aligned to +/- strand): ref base T (84/71); new base A (53/60); total (137/131)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 280 total aligned reads displayed.

TCCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAA  >  CP009273/2893298‑2893595

tCCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGtttt                                                                                                                                                       >  2:3205441/1‑150 (MQ=255)

   aGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGtttttaa                                                                                                                                                    >  2:3644584/1‑148 (MQ=255)

   aGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCa                                                                                                                                                    <  2:1203960/150‑1 (MQ=255)

     cACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAAt                                                                                                                                                  <  2:1674885/150‑1 (MQ=255)

       ccTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCa                                                                                                                                                >  1:2228789/1‑150 (MQ=255)

         ttATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAtt                                                                                                                                              <  1:846572/150‑1 (MQ=255)

          tATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCAttt                                                                                                                                             <  2:408072/150‑1 (MQ=255)

            ttCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTcc                                                                                                                                           >  1:1154822/1‑150 (MQ=255)

             tCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCt                                                                                                                                          <  2:3441306/150‑1 (MQ=255)

              cAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCtt                                                                                                                                         >  2:1935413/1‑150 (MQ=255)

                  ttACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGtt                                                                                                                                     <  2:2984788/150‑1 (MQ=255)

                   tACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGttt                                                                                                                                    >  1:206405/1‑150 (MQ=255)

                    aCAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTc                                                                                                                                   <  1:3402036/150‑1 (MQ=255)

                       aaTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGcc                                                                                                                                >  1:4119360/1‑150 (MQ=255)

                       aaTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGcc                                                                                                                                >  1:3317026/1‑150 (MQ=255)

                        aTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCg                                                                                                                               <  1:4273720/150‑1 (MQ=255)

                        aTACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCg                                                                                                                               >  1:2080058/1‑150 (MQ=255)

                          aCCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGtt                                                                                                                             >  2:2931246/1‑150 (MQ=255)

                          aCCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGtt                                                                                                                             <  1:677898/150‑1 (MQ=255)

                          aCCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGtt                                                                                                                             <  1:677993/150‑1 (MQ=255)

                            cGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTAt                                                                                                                           >  1:2034570/1‑150 (MQ=255)

                            cGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAAACATTTCCTTTGTTTCGCCGTTAt                                                                                                                           >  2:4249182/1‑150 (MQ=255)

                             ggCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATc                                                                                                                          >  2:2643979/1‑150 (MQ=255)

                                aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                       >  1:1344982/1‑150 (MQ=255)

                                aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                       <  2:2385566/150‑1 (MQ=255)

                                aTTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCttt                                                                                                                       <  2:630605/150‑1 (MQ=255)

                                 ttGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTg                                                                                                                      <  2:3357756/150‑1 (MQ=255)

                                   gTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGAc                                                                                                                    <  1:3490001/150‑1 (MQ=255)

                                      aaCCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATCTCCTTTGTTTCGCCGTTATCTTTGACGaa                                                                                                                 >  1:4359408/1‑150 (MQ=255)

                                          aGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTa                                                                                                             >  2:1358056/1‑150 (MQ=255)

                                          aGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTa                                                                                                             <  1:3860116/150‑1 (MQ=255)

                                                 cAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaa                                                                                                      >  2:1873324/1‑150 (MQ=255)

                                                  aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                     >  1:396157/1‑150 (MQ=255)

                                                  aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATaaa                                                                                                     >  2:2080027/1‑150 (MQ=255)

                                                  aaCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGACTaaa                                                                                                     >  2:451476/1‑150 (MQ=255)

                                                   aCTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAg                                                                                                    >  1:2442014/1‑150 (MQ=255)

                                                    cTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGa                                                                                                   <  2:4017976/150‑1 (MQ=255)

                                                      gTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTCTCTTTGACGAAACTAGGAATAAAGata                                                                                                 >  2:1583630/1‑150 (MQ=255)

                                                          cGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTg                                                                                             <  2:1810213/150‑1 (MQ=255)

                                                           gAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgc                                                                                            >  1:3136593/1‑150 (MQ=255)

                                                           gAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgc                                                                                            >  1:1233845/1‑150 (MQ=255)

                                                             aaCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGAATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTgcgc                                                                                          >  1:3240336/1‑150 (MQ=255)

                                                              aCGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGcc                                                                                         <  2:2495160/150‑1 (MQ=255)

                                                                gCTCACAGCAAGCAGCGATAAACTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAg                                                                                       >  2:2028414/1‑150 (MQ=255)

                                                                  tCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCt                                                                                     >  1:498579/1‑150 (MQ=255)

                                                                   cacaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCtt                                                                                    <  2:2129673/150‑1 (MQ=255)

                                                                    acaGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCttt                                                                                   >  1:2907419/1‑150 (MQ=255)

                                                                     caGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTg                                                                                  <  2:1238530/150‑1 (MQ=255)

                                                                      agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTgg                                                                                 >  2:87856/1‑150 (MQ=255)

                                                                      agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTgg                                                                                 <  2:1775274/150‑1 (MQ=255)

                                                                      agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTgg                                                                                 <  2:1516646/150‑1 (MQ=255)

                                                                      agcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTgg                                                                                 <  2:2172692/150‑1 (MQ=255)

                                                                       gcaagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGc                                                                                <  2:424082/150‑1 (MQ=255)

                                                                         aagcaGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCa                                                                              >  1:262295/1‑150 (MQ=255)

                                                                          agcagcGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCaa                                                                             >  1:2956741/1‑150 (MQ=255)

                                                                               cGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCa                                                                        <  1:1708723/150‑1 (MQ=255)

                                                                               cGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCa                                                                        >  1:969589/1‑150 (MQ=255)

                                                                               cGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCa                                                                        >  1:1783969/1‑150 (MQ=255)

                                                                                 aTAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATg                                                                      >  2:4205004/1‑150 (MQ=255)

                                                                                   aaTCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGc                                                                    >  2:591129/1‑150 (MQ=255)

                                                                                        tctGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATg                                                              <  2:3900233/150‑1 (MQ=255)

                                                                                               cgcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGcc                                                       <  1:2602293/150‑1 (MQ=255)

                                                                                                gcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCt                                                      <  1:182626/150‑1 (MQ=255)

                                                                                                gcgcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCt                                                      <  1:2357242/150‑1 (MQ=255)

                                                                                                  gcCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTaa                                                    >  2:93305/1‑150 (MQ=255)

                                                                                                    cTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAAc                                                  >  1:2533014/1‑150 (MQ=255)

                                                                                                        tGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCt                                              <  2:518347/150‑1 (MQ=255)

                                                                                                        tGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCt                                              >  2:4257043/1‑150 (MQ=255)

                                                                                                         gCGGTGATACCCACCTGCATGAAGTCCACCTCAACCCCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGCAAAATGCCTGGCCTAAACCGCt                                              <  1:1606932/150‑1 (MQ=255)

                                                                                                           ggTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTAtt                                           <  1:3555597/150‑1 (MQ=255)

                                                                                                           ggTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTAtt                                           >  1:1835851/1‑150 (MQ=255)

                                                                                                            gTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTc                                          >  1:4026688/1‑150 (MQ=255)

                                                                                                             tGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTcc                                         <  2:2814129/150‑1 (MQ=255)

                                                                                                                  cccACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCg                                    >  1:2308973/1‑150 (MQ=255)

                                                                                                                  cccACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCg                                    <  1:2633883/150‑1 (MQ=255)

                                                                                                                   ccACCTGCATGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCgg                                   <  1:215450/150‑1 (MQ=255)

                                                                                                                    cACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGt                                  <  1:660611/150‑1 (MQ=255)

                                                                                                                      ccTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTaa                                <  1:1956439/150‑1 (MQ=255)

                                                                                                                      ccTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTaa                                >  1:3708893/1‑150 (MQ=255)

                                                                                                                      catgaaTGAAGTCCAACTAAACACCCTGTTATTATATCATGTCATTTGTTTCTCCTTTATCTTTGACGAAACTAGGAATAAAGATATTTTAGCCATCTTTGGCCAACGCCATGGC‑AATTGCCTGGGCTAAACCGCTATTCACACAGGTaa                                <  2:2248430/145‑1 (MQ=255)

                                                                                                                          cATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAAt                            >  2:1904574/1‑150 (MQ=255)

                                                                                                                            tGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTg                          >  1:588038/1‑150 (MQ=255)

                                                                                                                            tGAAGTCCACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTg                          <  1:209994/150‑1 (MQ=255)

                                                                                                                             gAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGc                         <  2:3368417/150‑1 (MQ=255)

                                                                                                                               aGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCgg                       <  2:1546931/150‑1 (MQ=255)

                                                                                                                                 tCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGtt                     <  2:4324111/150‑1 (MQ=255)

                                                                                                                                  ccACCTCAACACCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGttt                    <  1:3339940/150‑1 (MQ=255)

                                                                                                                                   cACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGtttt                   <  1:102238/150‑1 (MQ=255)

                                                                                                                                     ccTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAc                 >  2:1155176/1‑150 (MQ=255)

                                                                                                                                      cTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTAcc                >  2:176220/1‑150 (MQ=255)

                                                                                                                                       tCAACACCCTTTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCt               <  1:2987313/150‑1 (MQ=255)

                                                                                                                                         aaCACCATGTTTTTCAATCATTTCCTTTGTTTCGCAGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTTGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCttt             <  2:3311124/150‑1 (MQ=255)

                                                                                                                                            acCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAg          >  1:4209633/1‑150 (MQ=255)

                                                                                                                                            acCCTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAg          <  2:2338439/150‑1 (MQ=255)

                                                                                                                                             cccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAgg         <  2:3701094/150‑1 (MQ=255)

                                                                                                                                              ccTGTTATTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAggg        >  1:691881/1‑150 (MQ=255)

                                                                                                                                               cTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGga       <  2:4245342/150‑1 (MQ=255)

                                                                                                                                               cTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGga       <  2:3136593/150‑1 (MQ=255)

                                                                                                                                                 gTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGgaga     >  2:164264/1‑150 (MQ=255)

                                                                                                                                                  tttttCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGaa    <  1:2325425/150‑1 (MQ=255)

TCCAGCACCTTATTCAGCTTACAAATACCGGCATTGTTAACCAGAATATCAACTGTACCGAAACGCTCACAGCAAGCAGCGATAATCTTCTGCGGCGCGCCTTCTGCGGTGATACCCACCTGCATGAAGTCCACCTCAACACCCTGTTTTTCAATCATTTCCTTTGTTTCGCCGTTATCTTTGACGAAACTAGGAATAAAGATATTTGCGCCAGCTTTGGCCAACGCCATGGC‑AAATGCCTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAA  >  CP009273/2893298‑2893595

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: `‑`