BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,495,063	A→G	8.6%	K397E (AAA→GAA)	yjhB →	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,495,063	0	A	G	8.6%	909.1 / 55.8	278	K397E (AAA→GAA)	yjhB	putative transporter
Reads supporting (aligned to +/- strand): ref base A (131/123); new base G (8/16); total (139/139)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 283 total aligned reads displayed.

TTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATG  >  CP009273/4494915‑4495212

ttatAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTaa                                                                                                                                                      <  2:1238920/150‑1 (MQ=255)

     aCTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAAttt                                                                                                                                                 <  2:4147087/150‑1 (MQ=255)

        tAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTtata                                                                                                                                              <  2:4192638/150‑1 (MQ=255)

        tAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTtata                                                                                                                                              <  2:4192015/150‑1 (MQ=255)

        tAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTtata                                                                                                                                              <  2:3109333/150‑1 (MQ=255)

         aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTTATaa                                                                                                                                             >  1:2571721/1‑150 (MQ=255)

         aGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTTATaa                                                                                                                                             >  2:2212720/1‑150 (MQ=255)

           gggCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAAc                                                                                                                                           <  2:2649932/150‑1 (MQ=255)

                aCTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATcc                                                                                                                                      >  2:1991452/1‑150 (MQ=255)

                aCTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATcc                                                                                                                                      <  2:2804359/150‑1 (MQ=255)

                 cTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAAATATCCg                                                                                                                                     >  1:4015048/1‑150 (MQ=255)

                   ggaggaATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTg                                                                                                                                   <  2:2035519/150‑1 (MQ=255)

                   ggaggaATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTg                                                                                                                                   <  2:4347178/150‑1 (MQ=255)

                   ggaggaATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTg                                                                                                                                   <  2:4348320/150‑1 (MQ=255)

                   ggaggaATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTg                                                                                                                                   >  2:1429137/1‑150 (MQ=255)

                     aggaATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGc                                                                                                                                 >  2:564987/1‑150 (MQ=255)

                           ggCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATg                                                                                                                           <  2:1222562/150‑1 (MQ=255)

                             ccGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGAt                                                                                                                         >  2:2776724/1‑150 (MQ=255)

                              cGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGAtt                                                                                                                        <  2:1607701/150‑1 (MQ=255)

                                cACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTgg                                                                                                                      <  2:2102636/150‑1 (MQ=255)

                                cACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTgg                                                                                                                      <  2:2101445/150‑1 (MQ=255)

                                 aCCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTTATAAACTATCCGTGGCTAAATGATTGGa                                                                                                                     <  2:876741/150‑1 (MQ=255)

                                    tGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGc                                                                                                                  <  1:4460888/150‑1 (MQ=255)

                                        ttAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTa                                                                                                              <  1:3430433/150‑1 (MQ=255)

                                          aGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGTGGCTTTATg                                                                                                            >  2:2572514/1‑150 (MQ=255)

                                           gCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGa                                                                                                           >  2:3967889/1‑150 (MQ=255)

                                            cTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGat                                                                                                          <  1:4245323/150‑1 (MQ=255)

                                            cTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGat                                                                                                          <  1:4244766/150‑1 (MQ=255)

                                             tacatacaTTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGatt                                                                                                         <  1:761002/150‑1 (MQ=255)

                                             tacatacaTTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGatt                                                                                                         >  2:4378486/1‑150 (MQ=255)

                                               catacatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaa                                                                                                       >  2:3090721/1‑150 (MQ=255)

                                               catacatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaa                                                                                                       >  2:3090689/1‑150 (MQ=255)

                                                atacatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaat                                                                                                      >  1:1006440/1‑150 (MQ=255)

                                                atacatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaat                                                                                                      >  1:2816159/1‑150 (MQ=255)

                                                atacatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGATACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaat                                                                                                      >  1:1006570/1‑150 (MQ=255)

                                                 tacatTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTTTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGattaatt                                                                                                     <  1:3704036/150‑1 (MQ=255)

                                                    atTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATg                                                                                                  >  1:3762408/1‑150 (MQ=255)

                                                    atTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTATGTTAGTTGGTTATGATATACCAGTTAAAATTTATAAACTATCATTGGCAAAATGATAGGAGGCTTAAAGATTAATTAtt                                                                                                  >  2:3919756/1‑149 (MQ=255)

                                                    atTTCAGGATATTAGGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATg                                                                                                  <  2:1358192/150‑1 (MQ=255)

                                                     tttCAGGATATTCTGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATgg                                                                                                 <  1:2706485/150‑1 (MQ=255)

                                                           gATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCgttgt                                                                                           >  1:3983351/1‑150 (MQ=255)

                                                            atatTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCgttgtt                                                                                          >  2:1295216/1‑150 (MQ=255)

                                                             tattatGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGttgttg                                                                                         >  1:1990087/1‑150 (MQ=255)

                                                               ttatGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTtgttggt                                                                                       >  1:4114766/1‑150 (MQ=255)

                                                                tatGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTtgttggtg                                                                                      <  2:1845929/150‑1 (MQ=255)

                                                                   ggCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTtgttggtgttg                                                                                   >  1:2650749/1‑150 (MQ=255)

                                                                        aGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATAc                                                                              >  1:810272/1‑150 (MQ=255)

                                                                         ggTGTTTCATTATTCATTGTTAAGGTTGCATTCTATGCCTTATTAATTTTGTTAGTTGGGTTTGATATTCAATGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGAGTTGTTGGTGTTGGATaat                                                                             >  2:505376/1‑148 (MQ=255)

                                                                            gtTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTg                                                                          >  2:1004711/1‑150 (MQ=255)

                                                                            gtTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTg                                                                          <  2:2944033/150‑1 (MQ=255)

                                                                            gtTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTg                                                                          <  2:3811682/150‑1 (MQ=255)

                                                                            gtTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTg                                                                          >  1:3293842/1‑150 (MQ=255)

                                                                              ttCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGgc                                                                        <  1:3973515/150‑1 (MQ=255)

                                                                              ttCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGgc                                                                        <  1:3219779/150‑1 (MQ=255)

                                                                              ttCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGgc                                                                        >  1:3432098/1‑150 (MQ=255)

                                                                               tCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTTTTGGTGTTGGATACTTTGgcg                                                                       >  2:2084059/1‑150 (MQ=255)

                                                                                  ttattCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTg                                                                    >  2:2289248/1‑150 (MQ=255)

                                                                                   tattCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGa                                                                   >  1:3642157/1‑150 (MQ=255)

                                                                                         ttGTTACGGTTGCATTCTCTTCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGc                                                             <  2:2851018/150‑1 (MQ=255)

                                                                                          tGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGTCGCTGAATTAGCt                                                            >  2:3825421/1‑150 (MQ=255)

                                                                                           gTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTc                                                           >  1:2973671/1‑150 (MQ=255)

                                                                                            ttACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCg                                                          <  2:4316771/150‑1 (MQ=255)

                                                                                            ttACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCg                                                          >  2:3368285/1‑150 (MQ=255)

                                                                                            ttACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCg                                                          >  1:3829457/1‑150 (MQ=255)

                                                                                            ttACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCg                                                          >  1:2692667/1‑150 (MQ=255)

                                                                                              aCGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGtt                                                        >  1:1373307/1‑150 (MQ=255)

                                                                                              aCGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGtt                                                        >  1:4377008/1‑150 (MQ=255)

                                                                                                   tGCATTCTCTGCCTTATTAATTTTGTTTGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATg                                                   <  2:1730426/150‑1 (MQ=255)

                                                                                                    gCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTGAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGa                                                  <  2:3770149/150‑1 (MQ=255)

                                                                                                    gCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGa                                                  <  1:3574029/150‑1 (MQ=255)

                                                                                                    gCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAAATATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGa                                                  >  2:3672088/1‑150 (MQ=255)

                                                                                                        tctcTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAatat                                              <  1:2440511/150‑1 (MQ=255)

                                                                                                            tGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCAt                                          <  1:831094/150‑1 (MQ=255)

                                                                                                            tGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCAt                                          <  1:2245361/150‑1 (MQ=255)

                                                                                                            tGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCAt                                          <  2:3413651/150‑1 (MQ=255)

                                                                                                             gCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATg                                         >  1:2893160/1‑150 (MQ=255)

                                                                                                             gCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATg                                         <  2:1885284/150‑1 (MQ=255)

                                                                                                               cTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGAt                                       >  2:1619732/1‑150 (MQ=255)

                                                                                                                   ttaATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATg                                   >  2:355794/1‑150 (MQ=255)

                                                                                                                   ttaATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATg                                   <  2:194256/150‑1 (MQ=255)

                                                                                                                     aaTTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCa                                 <  1:1931255/150‑1 (MQ=255)

                                                                                                                       ttttGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCaaa                               <  1:1746142/150‑1 (MQ=255)

                                                                                                                           gTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAAtt                           >  1:1125711/1‑150 (MQ=255)

                                                                                                                            ttAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTa                          >  1:1173348/1‑150 (MQ=255)

                                                                                                                             tAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTAc                         <  1:134521/150‑1 (MQ=255)

                                                                                                                              aGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACa                        >  2:3206287/1‑150 (MQ=255)

                                                                                                                              aGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACa                        <  1:536812/150‑1 (MQ=255)

                                                                                                                                ttGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACAtg                      <  1:564983/150‑1 (MQ=255)

                                                                                                                                 tGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACAtgt                     <  2:1127896/150‑1 (MQ=255)

                                                                                                                                  ggTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACAtgtg                    <  1:2854244/150‑1 (MQ=255)

                                                                                                                                   gTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATgtgt                   >  2:2770872/1‑150 (MQ=255)

                                                                                                                                      ttGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGtata                >  1:1130932/1‑150 (MQ=255)

                                                                                                                                       tGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATAc               <  2:1125711/150‑1 (MQ=255)

                                                                                                                                           atTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATc           <  1:1801224/150‑1 (MQ=255)

                                                                                                                                           atTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATc           >  2:1546314/1‑150 (MQ=255)

                                                                                                                                               cAGGTAAAATTTATAATCTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGa       <  2:2358595/150‑1 (MQ=255)

                                                                                                                                                aGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGaa      >  2:3477536/1‑150 (MQ=255)

                                                                                                                                                  gTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGaaaa    >  1:2044445/1‑150 (MQ=255)

                                                                                                                                                  gTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGaaaa    >  1:2045059/1‑150 (MQ=255)

                                                                                                                                                   tAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAAt   >  1:4329936/1‑150 (MQ=255)

TTATAACTTAGGGGCAACTGGAGGAATGGCCGCACCTGTATTAGCTACATACATTTCAGGATATTATGGCTTAGGTGTTTCATTATTCATTGTTACGGTTGCATTCTCTGCCTTATTAATTTTGTTAGTTGGTTTTGATATTCCAGGTAAAATTTATAAACTATCCGTGGCTAAATGATTGGAGGCTTTATGATTAATTATGGCGTTGTTGGTGTTGGATACTTTGGCGCTGAATTAGCTCGTTTTATGAATATGCATGATAATGCAAAAATTACATGTGTATACGATCCTGAAAATG  >  CP009273/4494915‑4495212

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