BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,941	A→T	38.9%	intergenic (‑1024/‑144)	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,941	0	A	T	38.9%	103.6 / 135.9	123	intergenic (‑1024/‑144)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (39/36); new base T (26/22); total (65/58)
Fisher's exact test for biased strand distribution p-value = 8.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 125 total aligned reads displayed.

TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAA  >  CP000730/1929802‑1930081

tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGGGGCGGTCCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                                                   >  2:61134/4‑141 (MQ=255)

   tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGGGGAGGTGCATTCC‑T‑TTGCTAA‑TCATTTTCACATCGCCGGA‑AGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                                >  1:177596/1‑140 (MQ=255)

   tataATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                                                <  2:747866/141‑2 (MQ=255)

     taataaTGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                                              <  2:251347/141‑4 (MQ=255)

     taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                                              <  2:657295/141‑4 (MQ=255)

     taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                              >  2:300211/1‑141 (MQ=255)

       ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                            <  1:608995/137‑1 (MQ=255)

       ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                            <  2:740733/141‑1 (MQ=255)

       ataCTAGAGTCGCCCACTTCTCACGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                            <  2:597991/141‑1 (MQ=255)

        taatagAGTCGCCTATCTCTCGGGCGGCAATT‑‑GACGAAGCGAGAAGGAATAGTTC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                                           <  1:362513/138‑1 (MQ=255)

         aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑a                                                                                                                                        <  1:323271/139‑1 (MQ=255)

          gtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa‑‑a                                                                                                                                        <  1:103884/141‑1 (MQ=255)

          atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAA‑TAATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                       >  2:195353/2‑140 (MQ=255)

            agagTTCCCTATCGCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGCCGCA‑‑ATa                                                                                                                                      <  2:515252/141‑1 (MQ=255)

            agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATGCTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                     <  1:98392/141‑4 (MQ=255)

            agagTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                     >  1:466623/1‑138 (MQ=255)

             gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATAc                                                                                                                                     <  2:708963/141‑1 (MQ=255)

             gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                    <  2:521906/141‑5 (MQ=255)

             gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                    >  1:719709/1‑137 (MQ=255)

              agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACa                                                                                                                                    >  1:132538/1‑141 (MQ=255)

              agTCGCCTATCTCTCAAGCGTCAATTTAGACGCAGAGAGGAGGTGAA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAGTGCGTTTTTTTCTTAGTGGCCAAGTAGACGCA‑‑ATACa                                                                                                                                    >  1:193572/1‑141 (MQ=255)

               gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAC‑G‑CT‑CTGCGTTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAACTAGACGCA‑‑ATACaa                                                                                                                                   <  1:402258/141‑1 (MQ=255)

               gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa                                                                                                                                   <  1:193340/141‑1 (MQ=255)

               gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaata                                                                                                                                  >  1:742323/1‑135 (MQ=255)

                tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataa                                                                                                                                 >  1:653965/1‑134 (MQ=255)

                gcGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACaa‑‑a                                                                                                                                <  1:79128/140‑1 (MQ=255)

                 cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                >  2:149529/1‑133 (MQ=255)

                 cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                >  1:645293/1‑133 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                               >  2:412835/1‑132 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                               <  2:725208/141‑10 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                               <  1:317320/141‑10 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                               >  2:499373/1‑132 (MQ=255)

                  gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                               <  2:409843/141‑10 (MQ=255)

                   ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                                             >  2:263200/1‑141 (MQ=255)

                   ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTGTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                                             >  1:639477/1‑141 (MQ=255)

                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                              >  2:266377/1‑131 (MQ=255)

                   ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                              >  2:184971/1‑131 (MQ=255)

                    cTATCTCTCAGGCGTCAATTTAGACGCCGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAg                                                                                                                            >  2:529813/1‑141 (MQ=255)

                    cTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTTAGTAGATGCA‑‑ATACAA‑‑AATAg                                                                                                                            >  2:712681/1‑141 (MQ=255)

                    cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                             >  1:140733/1‑130 (MQ=255)

                    cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                             <  1:627713/141‑12 (MQ=255)

                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTCGCGTATTATACTTATTGGCTTAGCAAACGCa‑‑acaaataattaat                                                                                                                            <  1:467103/141‑13 (MQ=255)

                     tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTTTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaagtaattaat                                                                                                                            >  1:286993/1‑129 (MQ=255)

                      aTCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                           <  1:416478/141‑14 (MQ=255)

                        ctatctggggTCCATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGa                                                                                                                        <  2:269466/132‑1 (MQ=255)

                           tcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑at                                                                                                                    <  2:36661/141‑1 (MQ=255)

                            cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAccga                                                                                                                    >  1:322778/1‑138 (MQ=255)

                            cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ata                                                                                                                   >  1:259031/1‑141 (MQ=255)

                            cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ata                                                                                                                   <  1:46988/141‑1 (MQ=255)

                            cAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ata                                                                                                                   >  2:495019/1‑141 (MQ=255)

                            cAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAC‑TCATTTTA‑GTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCTGTTCCGTATTATACTTATTGGCTTAGTAAACGCA‑‑ACACAT‑‑AAttaataggcgac                                                                                                                   >  2:111110/1‑130 (MQ=255)

                               gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑Ag                                                                                                               >  1:689926/1‑141 (MQ=255)

                                   cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gc                                                                                                          <  1:668149/141‑1 (MQ=255)

                                   cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATACGTGAC‑ATAT‑AGCC‑Gc                                                                                                          >  1:517510/1‑141 (MQ=255)

                                    aaTTTAGACGCAGATAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGTC‑GCa                                                                                                         <  1:547667/141‑1 (MQ=255)

                                       ttAGACGCAGAGGGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                                                      <  1:146770/141‑1 (MQ=255)

                                       ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCa                                                                                                      <  1:123006/141‑1 (MQ=255)

                                        tAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑Gccacaa                                                                                                     >  1:3115/1‑136 (MQ=255)

                                         aGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑GTAT‑AGCC‑GCACCAAt                                                                                                    <  1:568180/141‑1 (MQ=255)

                                         aGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAAt                                                                                                    >  1:292570/1‑141 (MQ=255)

                                               agagagGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAAt                                                                                              >  1:503593/1‑141 (MQ=255)

                                               agagagGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACAGCATAGCAACAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAAt                                                                                              >  1:369635/1‑141 (MQ=255)

                                                gagagGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AACAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                                             <  2:582040/141‑1 (MQ=255)

                                                    ggaggTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                                         >  1:653597/1‑141 (MQ=255)

                                                       ggTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGGTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                                      >  1:511171/1‑141 (MQ=255)

                                                       ggTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAc                                                                                      >  2:321042/1‑141 (MQ=255)

                                                              taA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGc                                                                                <  1:576517/141‑1 (MQ=255)

                                                               aa‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                               <  1:723494/141‑1 (MQ=255)

                                                               aa‑G‑GGGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                               <  1:233930/141‑1 (MQ=255)

                                                                             tATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                    <  1:756786/141‑1 (MQ=255)

                                                                              aTTTTC‑GTTCACATGATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGGGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                   <  1:247488/141‑1 (MQ=255)

                                                                              aTTTTC‑GATCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTATTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                   >  1:484195/1‑141 (MQ=255)

                                                                               ttttC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                  <  1:273750/141‑1 (MQ=255)

                                                                               ttttC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                  >  1:281437/1‑141 (MQ=255)

                                                                                     gTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgt                                                             >  1:90260/1‑141 (MQ=255)

                                                                                      ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTCCATTGCCAAGTCGCCGGG‑‑CCTCAA‑‑AAGAGGGGAC‑ATAT‑AGCA‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                                            >  1:754183/1‑141 (MQ=255)

                                                                                            tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTATATAAg                                                      >  1:693435/1‑141 (MQ=255)

                                                                                                aGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaaa                                                  >  2:692661/1‑141 (MQ=255)

                                                                                                 gCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcattaa                                                      >  2:23121/1‑138 (MQ=255)

                                                                                                  cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                                                >  2:432630/1‑141 (MQ=255)

                                                                                                   aCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                               >  1:699499/1‑141 (MQ=255)

                                                                                                    ccAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACt                                              >  1:77343/1‑141 (MQ=255)

                                                                                                         catcaGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                         <  1:727494/141‑1 (MQ=255)

                                                                                                            caGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAATTATTTTCg                                      <  2:260684/141‑1 (MQ=255)

                                                                                                               tGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCACATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg                                   >  1:121964/1‑141 (MQ=255)

                                                                                                                ggCTGTGCCATTGCGTTTTTTTCTTATTGGCGAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGt                                  <  2:578995/141‑1 (MQ=255)

                                                                                                                      gCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAAt                            <  2:653597/141‑1 (MQ=255)

                                                                                                                       ccATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AGACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTAGGGGGATTGGTGTATAAGTAAAGACTTATTATAGTTGTCTTAAtt                           >  2:498752/1‑141 (MQ=255)

                                                                                                                        cATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATta                          <  2:389378/141‑1 (MQ=255)

                                                                                                                          ttGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCGCTCACTACCGCAAATAGTGAGGGGATTGGTGTGTAAGTAAATACTTATTTTCGTTGTCTTAATtata                        >  2:339613/1‑141 (MQ=255)

                                                                                                                                   tttCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAtt               <  2:699499/141‑1 (MQ=255)

                                                                                                                                           tGGCTAAGTAAACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtttt       <  1:242692/141‑1 (MQ=255)

                                                                                                                                           tGGCAAAGTAGAGGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAATTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCtttt       >  1:542113/1‑141 (MQ=255)

                                                                                                                                            ggCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCttttt      <  1:687506/141‑1 (MQ=255)

                                                                                                                                            ggCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCttttt      >  1:51563/1‑141 (MQ=255)

                                                                                                                                            ggCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCttttt      >  1:150188/1‑141 (MQ=255)

                                                                                                                                            ggCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGAGTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCttttt      >  1:592960/1‑141 (MQ=255)

                                                                                                                                              cTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATGGTGAGGGGATTGGGGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGt    >  1:161978/1‑141 (MQ=255)

                                                                                                                                                aaGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaa  <  2:689926/141‑1 (MQ=255)

                                                                                                                                                aaGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTaa  >  1:226725/1‑141 (MQ=255)

TGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTC‑GTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAA  >  CP000730/1929802‑1930081

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