BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,946	G→A	34.6%	intergenic (‑1029/‑139)	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,946	0	G	A	34.6%	124.5 / 100.9	118	intergenic (‑1029/‑139)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (41/36); new base A (23/18); total (64/54)
Fisher's exact test for biased strand distribution p-value = 8.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

Only 100 of 121 total aligned reads displayed.

TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929807‑1930086

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

 aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                                  >  1:81441/1‑141 (MQ=255)

 aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                                                  >  2:173057/1‑141 (MQ=255)

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

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

    aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑a                                                                                                                                             <  1:702817/139‑1 (MQ=255)

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

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

     atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                            <  1:230485/140‑2 (MQ=255)

     atagAGTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTACACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                            <  1:492174/140‑2 (MQ=255)

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

       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAAGCa‑‑acaa                                                                                                                                          >  2:283390/1‑138 (MQ=255)

       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                          >  1:13729/1‑138 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

             gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                    >  2:596056/1‑132 (MQ=255)

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

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

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

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

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

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

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

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

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

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

               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑G‑TTGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                                  >  1:32603/1‑130 (MQ=255)

                tATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAgg                                                                                                                                >  1:454/1‑141 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                      aCGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaa                                                                                                       <  1:27185/141‑1 (MQ=255)

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

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

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

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

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

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

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

                                                          aa‑G‑GTGATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCa                                                                                    <  2:676258/141‑1 (MQ=255)

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

                                                                     c‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAg                                                                           >  2:255398/1‑141 (MQ=255)

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

                                                                         aTTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                                                        <  1:165511/141‑1 (MQ=255)

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

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

                                                                          ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                                       <  1:273750/141‑1 (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)

                                                                                         tAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTaa                                                        >  1:686362/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)

                                                                                             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)

                                                                                                                          tgttttCTCGTTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCGCCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTa                       <  2:711883/139‑1 (MQ=255)

                                                                                                                           tttttCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa                      >  1:536770/1‑141 (MQ=255)

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

                                                                                                                               ccTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTATtttt                  <  2:465782/140‑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:592457/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)

                                                                                                                                               gACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGGAAAatct  >  2:147902/1‑139 (MQ=255)

                                                                                                                                               gACGCA‑‑ATACAA‑‑AATAGAGGAC‑CTAT‑ACAC‑GCACAAAAAAAAAACCCCTCACTACCACAAATAGTGAGGGGATTAGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAActat  <  2:415104/141‑1 (MQ=255)

                                                                                                                                               gACGCA‑‑ATACAA‑‑AAAAGGAGAC‑ACAT‑AGCC‑GAACCAAAAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  <  1:712681/141‑1 (MQ=255)

TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑GTGATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929807‑1930086

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