BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,072,028	Δ1 bp	7.3%	intergenic (+523/+50)	ykgB → / ← USA300HOU_1927	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,072,028	0	T	.	7.3%	376.5 / 25.5	124	intergenic (+523/+50)	ykgB/USA300HOU_1927	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (53/60); new base . (6/3); total (60/64)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 127 total aligned reads displayed.

GTAGTATTAATGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCAAATGC  >  CP000730/2071888‑2072167

 tagtaTTAATGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAAtt                                                                                                                                            <  1:253421/141‑1 (MQ=32)

 tagtaTAAATGGCGAGGCTTCTGAGGGAGCACTGCCCGTCGAAGACCAAGGCGGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAGGATAATTTTCAATTA‑GAAAAt                                                                                                                                             <  2:903543/141‑1 (MQ=11)

     aTTAATGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCt                                                                                                                                        >  2:419475/1‑141 (MQ=32)

        aaTGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATGATTTTC‑AATTAGAAAA‑TATCttgc                                                                                                                                     <  1:915579/141‑3 (MQ=11)

        aaTGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATAATTTTC‑AATTAGAAAA‑TATCTTAc                                                                                                                                     <  2:55381/141‑1 (MQ=11)

        aaTGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTAc                                                                                                                                     <  2:352637/141‑1 (MQ=32)

        aaTGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTAc                                                                                                                                     <  1:363510/141‑1 (MQ=32)

        aaTGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTAc                                                                                                                                     <  2:1593981/141‑1 (MQ=32)

          tGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActg                                                                                                                                   >  1:836693/1‑141 (MQ=32)

           ggCGAGACTCGTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgc                                                                                                                                  <  1:1417581/141‑1 (MQ=17)

            gCGGGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgct                                                                                                                                 <  1:764380/141‑1 (MQ=17)

            gCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATGATTTTC‑AATT‑AGAAAATATCTTgctgct                                                                                                                                 <  1:109572/141‑1 (MQ=11)

            gCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgct                                                                                                                                 <  1:468148/141‑1 (MQ=32)

            gCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgct                                                                                                                                 <  1:343893/141‑1 (MQ=32)

            gCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgct                                                                                                                                 <  1:766297/141‑1 (MQ=32)

            gCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgct                                                                                                                                 <  2:1300467/141‑1 (MQ=32)

            gCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTActgct                                                                                                                                 <  1:1089041/141‑1 (MQ=32)

             cGAGACTCCAGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGATGATTTTCAATTA‑GAAAA‑TATCTTGctgctg                                                                                                                                <  1:542009/141‑1 (MQ=11)

               agaCTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGtt                                                                                                                              >  1:1478089/1‑141 (MQ=32)

                 aCTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTAc                                                                                                                                     <  1:1304287/132‑1 (MQ=33)

                     cTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAg                                                                                                                        >  1:1292179/1‑141 (MQ=32)

                      tGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAgg                                                                                                                       >  1:866912/1‑141 (MQ=32)

                       gAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAggg                                                                                                                      >  2:1026477/1‑141 (MQ=32)

                            aGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTa                                                                                                                 <  1:654802/141‑1 (MQ=32)

                               aGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGt                                                                                                              >  2:472235/1‑141 (MQ=32)

                               aGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGt                                                                                                              >  1:190360/1‑141 (MQ=32)

                               aGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGt                                                                                                              >  2:1384262/1‑141 (MQ=32)

                                gTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGCACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTc                                                                                                             <  2:542811/141‑1 (MQ=17)

                                  gCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTAGTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTccc                                                                                                           <  2:1193576/141‑1 (MQ=17)

                                        cGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCctc                                                                                                     <  2:897955/141‑1 (MQ=32)

                                        cGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCctc                                                                                                     >  1:1439096/1‑141 (MQ=32)

                                         gAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCctct                                                                                                    <  2:849448/141‑1 (MQ=32)

                                         gAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCctct                                                                                                    <  1:616440/141‑1 (MQ=32)

                                          aaGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTa                                                                                                   >  1:1528408/1‑141 (MQ=33)

                                              ccAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAgg                                                                                                                       >  2:445890/1‑117 (MQ=32)

                                              ccAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAgg                                                                                                                       <  1:445890/117‑1 (MQ=32)

                                                  ggCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat                                                                                           <  2:1044777/141‑1 (MQ=37)

                                                  ggCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat                                                                                           >  1:836615/1‑141 (MQ=37)

                                                  ggCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGttat                                                                                           >  1:833291/1‑141 (MQ=37)

                                                     tGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAg                                                                                        >  1:345248/1‑141 (MQ=39)

                                                     tGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAAGT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAg                                                                                        <  1:1161530/141‑1 (MQ=255)

                                                        gaCGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTa                                                                                     >  2:870077/1‑141 (MQ=255)

                                                          cGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtt                                                                                   <  2:1143700/141‑1 (MQ=255)

                                                           ggCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAttt                                                                                  >  2:1544187/1‑141 (MQ=255)

                                                            gCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtttt                                                                                 <  2:254283/141‑1 (MQ=255)

                                                            gCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTAACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAtttt                                                                                 >  2:1317109/1‑141 (MQ=255)

                                                             cACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGTTTTATGTCCCAGCCTCTAAATGTTATTAGTTAttttt                                                                                <  2:391557/141‑1 (MQ=255)

                                                             cACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAttttt                                                                                <  2:373125/141‑1 (MQ=255)

                                                             cACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTAttttt                                                                                <  1:1125338/141‑1 (MQ=255)

                                                               cccTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTAc                                                                              >  2:538375/1‑141 (MQ=255)

                                                                 cTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACtt                                                                            <  2:1642774/141‑1 (MQ=255)

                                                                    ggAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCa                                                                         >  1:627302/1‑141 (MQ=255)

                                                                              gCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGaaa                                                               <  1:590228/141‑1 (MQ=255)

                                                                               ccATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAaat                                                              <  2:833291/141‑1 (MQ=255)

                                                                                    caatacttcgTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAAttt                                                         >  1:1566906/10‑141 (MQ=255)

                                                                                    cAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAAttt                                                         >  1:165845/1‑141 (MQ=255)

                                                                                     aatacttcgTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAACAATTTc                                                        >  1:708507/9‑141 (MQ=255)

                                                                                       tACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCtt                                                      <  2:31812/141‑1 (MQ=255)

                                                                                         cGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTa                                                    <  1:1138113/141‑1 (MQ=255)

                                                                                           aaGTATTGTATAAATAGAGAACAGCAGCAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGt                                                  <  2:261452/141‑1 (MQ=255)

                                                                                            aGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGtt                                                 >  2:687914/1‑141 (MQ=255)

                                                                                             gTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGttt                                                <  1:816384/141‑1 (MQ=255)

                                                                                             gTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGttt                                                <  2:883916/141‑1 (MQ=255)

                                                                                               aTTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGttttt                                              >  1:162821/1‑141 (MQ=255)

                                                                                                   tataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTtttcttt                                          >  2:650881/1‑141 (MQ=255)

                                                                                                   tataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTtttcttt                                          >  2:289413/1‑141 (MQ=255)

                                                                                                    ataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTtttctttc                                         <  2:1128738/141‑1 (MQ=255)

                                                                                                    ataAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTtttctttc                                         >  2:1000997/1‑141 (MQ=255)

                                                                                                      ccctaaaaaaCAGCAGTAAGATAATTTTC‑AATTAGAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCaa                                       >  2:1190835/9‑141 (MQ=25)

                                                                                                       aaTAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAat                                      <  2:162821/141‑1 (MQ=255)

                                                                                                         tAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAatat                                    >  1:1160881/1‑141 (MQ=255)

                                                                                                          agagAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTCGTTTTTCTTTCAATATc                                   >  1:1335344/1‑141 (MQ=255)

                                                                                                           gagaACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCt                                  >  1:1175874/1‑141 (MQ=255)

                                                                                                              aaCAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCa                               >  2:922449/1‑141 (MQ=255)

                                                                                                                cagcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATCAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAatagaaacgtt                             >  1:359205/1‑133 (MQ=255)

                                                                                                                  gcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTCTTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAAt                           >  2:730137/1‑141 (MQ=255)

                                                                                                                  gcagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAAt                           >  2:1061823/1‑141 (MQ=255)

                                                                                                                   cagTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATc                          >  1:491239/1‑141 (MQ=255)

                                                                                                                     gTAAGAT‑AATTTTCAATT‑AGAAAATATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTAAAAATCtt                        >  2:525653/1‑141 (MQ=37)

                                                                                                                      tAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCCGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCttt                       >  1:435264/1‑141 (MQ=255)

                                                                                                                        agataatTTTCAATTA‑GAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAGATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCttttt                     >  1:1304732/6‑141 (MQ=37)

                                                                                                                          at‑atTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTa                   <  2:1377418/141‑1 (MQ=255)

                                                                                                                             tttctgataTT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGaaa                                                               >  2:550544/9‑95 (MQ=255)

                                                                                                                             tttctgataTT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGaaa                                                               <  1:550544/87‑1 (MQ=255)

                                                                                                                             atTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGTTTTTTTTGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTaaa                 >  2:563811/1‑141 (MQ=255)

                                                                                                                             atTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCtttaaaaaa                 <  2:434517/141‑7 (MQ=255)

                                                                                                                             aattttcaattagAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTaaa                 >  2:527971/13‑141 (MQ=38)

                                                                                                                               tttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTTaaaaaagattttgaaga               <  1:870077/141‑14 (MQ=37)

                                                                                                                                ttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGat              >  1:184135/1‑141 (MQ=255)

                                                                                                                                ttCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGat              >  1:276324/1‑141 (MQ=255)

                                                                                                                                  tcaattagAAAA‑TATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGatat            >  1:541276/8‑141 (MQ=255)

                                                                                                                                   tAATT‑GAAAATTATCTTAGTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGAtata           <  2:1288026/141‑1 (MQ=255)

                                                                                                                                   tAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGttatt           >  2:1112552/1‑140 (MQ=255)

                                                                                                                                      tt‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGc        <  1:1357481/141‑1 (MQ=255)

                                                                                                                                         gAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCaa      <  1:705268/141‑1 (MQ=255)

                                                                                                                                          aaaaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCaaa     <  2:1340126/141‑1 (MQ=255)

                                                                                                                                          aaaaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAACAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCaaa     <  2:708507/141‑1 (MQ=255)

                                                                                                                                            aaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCAAATg   >  1:143647/1‑141 (MQ=255)

                                                                                                                                            aaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCAAATg   <  1:279420/141‑1 (MQ=255)

                                                                                                                                            aaTTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCAAATg   <  1:1484233/141‑1 (MQ=255)

GTAGTATTAATGGCGAGACTCCTGAGGGAGCAGTGCCAGTCGAAGACCAAGGCTGAGACGGCACCCTAGGAAAGCGAAGCCATTCAATACGAAGTATTGTATAAATAGAGAACAGCAGTAAGAT‑ATTTTCTAATT‑GAAAATTATCTTACTGCTGTTTTTTAGGGATTTATGTCCCAGCCTCTAAATGTTATTAGTTATTTTTACTTTCAATATCAGAAATAATTTCTTTAGTTTTTCTTTCAATATCTTCAAAATCTTTTTTAAAGATATAAGCAAATGC  >  CP000730/2071888‑2072167

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