BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,804	A→G	42.7%	intergenic (‑887/‑281)	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,804	0	A	G	42.7%	35.3 / 96.3	89	intergenic (‑887/‑281)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (22/29); new base G (14/24); total (36/53)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

ATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGT  >  CP000730/1929666‑1929944

    aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                              >  1:136150/1‑141 (MQ=255)

      tCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtag                                                                                                                                            >  1:83127/1‑141 (MQ=255)

       tggcTACTTTCTTAAATTAACAATATTTATTCTCATACAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                                                           <  2:162808/138‑1 (MQ=255)

          cTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTaga                                                                                                                                        >  2:357100/1‑141 (MQ=255)

                tCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGcc                                                                                                                                  >  1:686996/1‑141 (MQ=255)

                     aaTTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTATACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtct                                                                                                                             <  2:533084/141‑1 (MQ=255)

                     aaTTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGGGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTAtct                                                                                                                             <  2:362503/141‑1 (MQ=255)

                            aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                      <  2:517654/141‑1 (MQ=255)

                            aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                      >  1:564504/1‑141 (MQ=255)

                            aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                      <  2:517621/141‑1 (MQ=255)

                            aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                      <  1:535998/141‑1 (MQ=255)

                            aaTATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGc                                                                                                                      <  1:348710/141‑1 (MQ=255)

                                tCTATTCTCATAGAATTTGTCCAATTGAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGGATTTTTGCCGATGATATAATAGTAGAGTCGCCTATCTTCCAGGCGTCa                                                                                                                  <  1:613299/141‑1 (MQ=39)

                                  tATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGACATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAt                                                                                                                >  2:589259/1‑141 (MQ=255)

                                  tATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAAATTATCAGAAGAAAGTTACTGTGCTTGTATATTTGACGATGATATAATAGTAGAGTCGCATATCTCTCAGGGGTCAAt                                                                                                                >  2:405245/1‑141 (MQ=38)

                                  tATGCTCAAAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAt                                                                                                                <  2:811567/141‑1 (MQ=255)

                                            gcaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                                      <  1:677700/139‑1 (MQ=255)

                                            aaaTTTGTCCAATTAAGTGTAGACGATTCAATCAATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTACATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                                      <  1:160852/140‑1 (MQ=255)

                                            aaaTTTGTCCAATTAAGTGAAAACGATACAACCCATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCag                                                                                                      <  1:160462/140‑1 (MQ=255)

                                             aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                     <  1:354427/141‑1 (MQ=255)

                                             aaTTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCaga                                                                                                     >  2:630651/1‑141 (MQ=255)

                                               tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                   <  2:43878/141‑1 (MQ=255)

                                               tttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagaga                                                                                                   <  1:357100/141‑1 (MQ=255)

                                                ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                                  <  2:815550/141‑1 (MQ=255)

                                                  gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAgagga                                                                                                >  1:816617/1‑141 (MQ=255)

                                                     cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGaggaggg                                                                                             >  1:122296/1‑140 (MQ=255)

                                                          aaGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAGTTTTGAAGCAGAGATGAGGTGtc‑ta                                                                                       >  2:336648/1‑141 (MQ=255)

                                                           aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGGGGAGGAGTA‑Taa                                                                                      >  1:334782/1‑141 (MQ=255)

                                                                          aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGGGGTGTA‑TAA‑G‑G‑T‑GA‑TTGTTGtttt                                                                  >  1:720683/1‑141 (MQ=255)

                                                                          aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTCtttt                                                                  >  2:819802/1‑141 (MQ=255)

                                                                                cAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCGGAGAGGGGGTGTG‑TAA‑G‑G‑T‑GT‑TGCTTTTTTTCGTTca                                                            >  1:594507/1‑141 (MQ=255)

                                                                                 aGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTcac                                                           <  2:816617/141‑1 (MQ=255)

                                                                                  gCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTcaca                                                          >  2:35556/1‑141 (MQ=255)

                                                                                    tAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCAcatc                                                        <  1:808235/141‑1 (MQ=255)

                                                                                      gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCAcatcat                                                      >  1:698946/1‑141 (MQ=255)

                                                                                        aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTAGTTCACATCATAg                                                    >  2:603952/1‑141 (MQ=255)

                                                                                           atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCAc                                                 >  2:327764/1‑141 (MQ=255)

                                                                                            tataACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGAAcc                                                >  2:791623/1‑141 (MQ=255)

                                                                                                    aTCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGtcatca                                        <  1:138099/141‑1 (MQ=255)

                                                                                                        atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                   >  2:192086/15‑141 (MQ=255)

                                                                                                        atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                   >  2:247314/15‑141 (MQ=255)

                                                                                                         tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCCCCAGCCATCAGTGCCt                                  >  1:42930/14‑141 (MQ=255)

                                                                                                          cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                 >  2:136253/13‑141 (MQ=255)

                                                                                                           ggaagtagatatgCTTTTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                <  2:433576/130‑1 (MQ=255)

                                                                                                           ggaagtagatatgCTTGTATTTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCAGTTGCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                <  2:158240/130‑1 (MQ=255)

                                                                                                           ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                <  1:488992/130‑1 (MQ=255)

                                                                                                           ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                <  1:803573/130‑1 (MQ=255)

                                                                                                           gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                 <  1:167235/141‑1 (MQ=255)

                                                                                                            gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTTTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                               <  1:649193/131‑1 (MQ=255)

                                                                                                            gaagtagataggCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                               <  1:817627/130‑1 (MQ=255)

                                                                                                             gagttgatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                              <  2:321725/132‑1 (MQ=255)

                                                                                                             aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                              <  2:558803/132‑1 (MQ=255)

                                                                                                             aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TGGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                              >  1:337450/10‑141 (MQ=255)

                                                                                                             aaGTGACTGTGCTTGGATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                               <  1:59458/141‑1 (MQ=255)

                                                                                                              agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                             <  2:673894/133‑2 (MQ=255)

                                                                                                              agtagatatgCTTGTATATTCACAGATGGGATAATAATAGAGTCGCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                             <  2:26005/133‑2 (MQ=255)

                                                                                                               gtagatttgCTTGTATATTCACAGATGGTATAATAATGGAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                            <  2:11155/134‑3 (MQ=255)

                                                                                                               gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                            >  1:349963/8‑139 (MQ=255)

                                                                                                               gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                            >  1:760732/8‑139 (MQ=255)

                                                                                                               gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTACGGCGACTAATT‑‑ACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                            <  1:823123/134‑3 (MQ=255)

                                                                                                                tagatatgCTTGTATATTCAAAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                           <  2:774787/135‑4 (MQ=255)

                                                                                                                 agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAGTT‑‑GCCGACGATAGACGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTAGTTCGCATCATAGCCCCCGTCAGCAGTGGCTgcgctgct                          >  2:219394/6‑137 (MQ=255)

                                                                                                                 agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATAAGTGGCTGTGctgtt                          >  2:352146/6‑137 (MQ=255)

                                                                                                                  gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                         >  1:208871/5‑136 (MQ=255)

                                                                                                                  gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                         <  1:297781/137‑6 (MQ=255)

                                                                                                                   cttttcTTGTATATTTGGCGATGATATAATAGTCGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTg                         <  2:784971/136‑1 (MQ=255)

                                                                                                                    tatgCTTGTATATTCACAGAAGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                       <  1:467640/139‑1 (MQ=255)

                                                                                                                     gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                       <  1:632163/141‑1 (MQ=255)

                                                                                                                     gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                       <  1:35556/141‑1 (MQ=255)

                                                                                                                     gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACAAGTCATCAGTGGCTGTGCCATTACg                       >  1:220442/1‑141 (MQ=255)

                                                                                                                     atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                      >  2:601971/2‑141 (MQ=255)

                                                                                                                      tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                     <  1:479830/141‑2 (MQ=255)

                                                                                                                      tgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCACTTC‑A‑T‑T‑GC‑AAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGta                     >  1:822430/1‑140 (MQ=255)

                                                                                                                       gCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat                    <  1:550332/141‑3 (MQ=255)

                                                                                                                       gCTTGGATATTGGCCGATGATAGAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt                     <  2:199097/141‑1 (MQ=255)

                                                                                                                        cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                   <  1:120611/141‑4 (MQ=255)

                                                                                                                         ttGTATATTCACAGATGGTGTAATAATAGAGTCGCCTATCTCTCAGGCGTCAACT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                  >  2:511428/1‑140 (MQ=255)

                                                                                                                           gTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt                 <  1:704293/141‑1 (MQ=255)

                                                                                                                            tataTTCACAGATGGTATAATAATAGATTCGCCTATCTTTAAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac               <  2:278859/141‑5 (MQ=255)

                                                                                                                            tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac               <  2:356748/141‑8 (MQ=255)

                                                                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt             >  2:55213/9‑141 (MQ=255)

                                                                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt             >  1:618788/9‑141 (MQ=255)

                                                                                                                              tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑TACTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACtt             >  2:605023/9‑141 (MQ=255)

                                                                                                                              tatTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt              <  2:255223/141‑1 (MQ=255)

                                                                                                                                ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt           >  2:810716/7‑141 (MQ=255)

                                                                                                                                ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGGGGCTGTGCTGTTGCGTATTATACTTAt           >  1:143126/7‑141 (MQ=255)

                                                                                                                                 tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt          <  1:361057/136‑1 (MQ=255)

                                                                                                                                  cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg         >  1:790995/5‑141 (MQ=255)

                                                                                                                                  cacagatggtAGAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg         >  2:616332/5‑141 (MQ=255)

                                                                                                                                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       <  1:413275/139‑1 (MQ=255)

                                                                                                                                    cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc       <  1:563988/139‑1 (MQ=255)

                                                                                                                                     agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      <  1:786131/140‑1 (MQ=255)

                                                                                                                                     agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt      <  1:786025/140‑1 (MQ=255)

                                                                                                                                      gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATCATTCTTGGCtt     >  2:752379/1‑140 (MQ=255)

                                                                                                                                      gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     >  1:48936/1‑140 (MQ=255)

                                                                                                                                      gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCt      <  1:619346/141‑1 (MQ=255)

                                                                                                                                      aatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GC‑TAATCATTTTCGTTCACGTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt     <  2:50584/140‑2 (MQ=255)

                                                                                                                                       atgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTa     <  1:490140/141‑1 (MQ=255)

                                                                                                                                        tgatATAATAGTAGAGTCGCCTATATCTCAAGCTTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTaa    >  2:347823/1‑141 (MQ=255)

                                                                                                                                          atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt  >  2:41776/1‑141 (MQ=255)

ATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGT  >  CP000730/1929666‑1929944

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