BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,865	A→C	34.8%	intergenic (‑948/‑220)	USA300HOU_1794 ← / → USA300HOU_1796	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,929,865	0	A	C	34.8%	170.8 / 255.4	201	intergenic (‑948/‑220)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (78/53); new base C (32/38); total (110/91)
Fisher's exact test for biased strand distribution p-value = 7.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 204 total aligned reads displayed.

AGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929725‑1930004

aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑Taa                                                                                                                                                >  1:82037/1‑141 (MQ=255)

 gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑g                                                                                                                                              <  2:188592/141‑1 (MQ=255)

    tAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑Tg                                                                                                                                         >  1:69780/1‑141 (MQ=255)

     aGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGa                                                                                                                                        <  1:687652/141‑1 (MQ=255)

       aCGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATg                                                                                                                                      <  2:625618/141‑1 (MQ=255)

       aCGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATg                                                                                                                                      <  1:1148875/141‑1 (MQ=255)

         gATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCt                                                                                                                                    <  1:1583504/141‑1 (MQ=255)

          aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑T‑G‑TGATGCtt                                                                                                                                   <  1:1458031/141‑1 (MQ=255)

          aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCtt                                                                                                                                   <  1:832439/141‑1 (MQ=255)

          aTTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCtt                                                                                                                                   <  2:263904/141‑1 (MQ=255)

               aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTa                                                                                                                                  >  2:596623/1‑137 (MQ=255)

               aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTAtttt                                                                                                                              >  2:110052/1‑141 (MQ=255)

                tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTc                                                                                                                             >  1:1503349/1‑141 (MQ=255)

                tCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGGATATTTGCCGATGATATAATAGTAGAGGCGCCTAGCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTc                                                                                                                             >  1:163598/1‑141 (MQ=255)

                  tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGt                                                                                                                           >  2:793150/1‑141 (MQ=255)

                     cAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTca                                                                                                                        <  2:244915/141‑1 (MQ=255)

                      aGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTAt                                                                                                                                 <  1:107119/131‑1 (MQ=255)

                      aGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTAt                                                                                                                                 >  2:107119/1‑131 (MQ=255)

                       gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGGTATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTcaca                                                                                                                      >  2:32926/1‑141 (MQ=255)

                         tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAcatc                                                                                                                    <  2:770331/141‑1 (MQ=255)

                         tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAcatc                                                                                                                    <  1:356793/141‑1 (MQ=255)

                               cATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCa                                                                                                              >  2:616746/1‑141 (MQ=255)

                                atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCGGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAc                                                                                                             <  2:759330/141‑1 (MQ=255)

                                  ataACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                           >  1:642179/1‑141 (MQ=255)

                                  ataACTTATCAGAAGCAAGTGACTGAGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                                                                                           >  1:1374335/1‑141 (MQ=255)

                                    aaCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGt                                                                                                         <  2:122840/141‑1 (MQ=255)

                                    aaCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAcctac                                                                                                         >  2:153051/1‑138 (MQ=255)

                                         aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                                                                                    <  2:884032/141‑1 (MQ=255)

                                            aGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAgccactgatg                                                                                                 <  1:110052/141‑10 (MQ=255)

                                            aGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGccactgatg                                                                                                 <  1:1121115/141‑7 (MQ=255)

                                             gAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                                                                                                >  2:753997/1‑141 (MQ=255)

                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                               <  2:1472902/127‑1 (MQ=255)

                                             atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                               <  2:548217/127‑1 (MQ=255)

                                              tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                              <  1:380260/128‑1 (MQ=255)

                                              tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                              >  1:1564238/14‑141 (MQ=255)

                                               cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                             >  1:949039/13‑141 (MQ=255)

                                               aGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                                                                                              >  2:1002215/1‑141 (MQ=255)

                                                ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                            >  2:1356337/12‑141 (MQ=255)

                                                   aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                          >  2:781169/1‑141 (MQ=255)

                                                    gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                        >  1:834694/8‑139 (MQ=255)

                                                     tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                       <  1:317160/135‑4 (MQ=255)

                                                     tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                       <  1:1525531/135‑4 (MQ=255)

                                                      gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       <  2:199577/141‑1 (MQ=255)

                                                      gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       >  1:273451/1‑141 (MQ=255)

                                                      gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       <  2:862931/141‑1 (MQ=255)

                                                      gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                       >  2:1601841/1‑141 (MQ=255)

                                                      agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                                      >  2:185500/6‑137 (MQ=255)

                                                      agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                                      <  2:715002/136‑5 (MQ=255)

                                                       gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                     <  2:1656859/137‑6 (MQ=255)

                                                       gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                     >  1:1127805/5‑136 (MQ=255)

                                                        cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                                                                                     <  2:1188427/141‑1 (MQ=255)

                                                        atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                    >  2:1558170/4‑135 (MQ=255)

                                                            gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCgat                                                                                 >  2:227279/1‑139 (MQ=255)

                                                            gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                                                                                >  2:417076/1‑139 (MQ=255)

                                                             cTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGa                                                                                                                                        >  1:1609507/1‑85 (MQ=255)

                                                                 tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                            <  1:611918/141‑1 (MQ=255)

                                                                    ttttGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCCATTGCGTTTTTTTcgt                                                                         <  2:895126/140‑3 (MQ=255)

                                                                       cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                               >  2:953128/5‑128 (MQ=255)

                                                                       cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                               <  1:953128/127‑4 (MQ=255)

                                                                       cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGCAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                     <  1:782490/137‑1 (MQ=255)

                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                   <  1:1543040/139‑1 (MQ=255)

                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                   >  1:1257520/3‑141 (MQ=255)

                                                                           gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt                                                                  <  1:172675/141‑1 (MQ=255)

                                                                              gatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                               <  1:753997/141‑1 (MQ=255)

                                                                               gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                             >  2:863701/2‑141 (MQ=255)

                                                                               atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                              >  1:562697/1‑141 (MQ=255)

                                                                               atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                              >  1:427389/1‑141 (MQ=255)

                                                                                tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGGAGAGAGGAGGT‑GTA‑TGA‑G‑G‑TGATGCTTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGCTAAGTa                                                             <  2:319493/141‑1 (MQ=255)

                                                                                tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                            >  1:224511/1‑140 (MQ=255)

                                                                                    ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                        <  2:1369602/137‑1 (MQ=255)

                                                                                    ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                        <  2:939624/137‑1 (MQ=255)

                                                                                     tagtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGc                                                        >  2:650653/1‑141 (MQ=255)

                                                                                      agtagAGTCGCCTATCTCGCAGGCGTCAAATTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                       >  2:473945/1‑141 (MQ=255)

                                                                                      aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                    <  2:586757/109‑7 (MQ=255)

                                                                                      aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                    >  1:586757/3‑105 (MQ=255)

                                                                                       atagAGTCGCCTATCTCTCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                     <  2:1651970/140‑2 (MQ=255)

                                                                                        tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                     >  1:393789/1‑141 (MQ=255)

                                                                                         agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                   >  1:483866/1‑138 (MQ=255)

                                                                                           agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                 <  1:1173820/141‑6 (MQ=255)

                                                                                            gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGATGGGG‑TAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                >  1:798539/1‑135 (MQ=255)

                                                                                               gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGCGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                              >  2:1102600/1‑141 (MQ=255)

                                                                                               gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                              >  2:1452066/1‑141 (MQ=255)

                                                                                               gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                             <  1:1499478/141‑10 (MQ=255)

                                                                                                 cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                           <  2:434342/141‑12 (MQ=255)

                                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGT‑GCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGGAAACGCaacaaataattaata                                         <  2:389714/141‑14 (MQ=255)

                                                                                                    tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGACGCAATACAAAATAGGTg                                         >  1:627953/1‑141 (MQ=255)

                                                                                                          aGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACatat                                   <  1:1058421/141‑1 (MQ=255)

                                                                                                                 aaTTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGCCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                            <  2:1303011/141‑1 (MQ=255)

                                                                                                                   tttAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                          <  2:82037/141‑1 (MQ=255)

                                                                                                                      aGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                       <  1:1087759/141‑1 (MQ=255)

                                                                                                                      aGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                       >  2:1034142/1‑141 (MQ=255)

                                                                                                                      aGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                       <  2:1503160/141‑1 (MQ=255)

                                                                                                                       gACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATa                      >  2:398111/1‑141 (MQ=255)

                                                                                                                          gCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                   <  1:1004352/141‑1 (MQ=255)

                                                                                                                          gCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                   >  1:71364/1‑141 (MQ=255)

                                                                                                                          gCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTAATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                   >  2:1318084/1‑141 (MQ=255)

                                                                                                                                   aggT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa          >  2:767599/1‑141 (MQ=255)

                                                                                                                                   aggT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCa          >  1:1449911/1‑141 (MQ=255)

                                                                                                                                             aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCa  >  1:482454/1‑141 (MQ=255)

                                                                                                                                             aa‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCATCACTACCGCa  >  2:1164721/1‑141 (MQ=255)

AGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGT‑GTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929725‑1930004

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