BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,946	G→A	47.5%	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	47.5%	9.6 / 149.5	125	intergenic (‑1029/‑139)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (32/33); new base A (36/23); total (68/57)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.92e-01

Only 100 of 127 total aligned reads displayed.

TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAATAT  >  CP000730/1929807‑1930086

taataagagcggcGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                                                   <  1:42065/129‑1 (MQ=255)

  ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTCC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCGGTCAGCAGCGGCTTTGCTGTTGCGTATTATTCTTATTGGCTTAGTAAACGc                                                                                                                                                 >  2:198069/5‑141 (MQ=255)

  ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                                                 <  1:185103/137‑1 (MQ=255)

  ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATAATACGGATTGGCTTAGTAAACGc                                                                                                                                                 <  1:226877/141‑1 (MQ=255)

  ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTG‑C‑TTTACTAAGCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGTCTTAGTAAAccg                                                                                                                                                 >  2:261080/1‑139 (MQ=255)

  ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                                                  >  2:70133/1‑141 (MQ=255)

     gtagAGACGACTATCTCTCAGGCGCCAAGTTAGACGCAGAGGGGATGTG‑TA‑TAA‑G‑GGGA‑TTCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa‑‑a                                                                                                                                             <  2:414113/141‑1 (MQ=255)

     atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑ac                                                                                                                                            >  2:633990/2‑140 (MQ=255)

      tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTACTTAGACGCA‑‑At                                                                                                                                            >  2:259973/1‑141 (MQ=255)

       agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATa                                                                                                                                           <  1:913395/141‑1 (MQ=255)

       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGTACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                          <  2:932192/141‑4 (MQ=255)

       agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaa                                                                                                                                          >  2:682593/1‑138 (MQ=255)

        gagTCGCCTATCTTTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                         <  2:378383/141‑5 (MQ=255)

        gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAGGCa‑‑acaaa                                                                                                                                         >  1:68318/1‑137 (MQ=255)

        gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑agaat                                                                                                                                         >  1:442004/1‑137 (MQ=255)

        gagTCGCCTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                         <  1:733450/141‑5 (MQ=255)

        aagTCGCCTATCTCTCAGGCGTCAAGT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaa                                                                                                                                         <  1:14138/140‑5 (MQ=255)

         agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACa                                                                                                                                         >  2:60435/1‑141 (MQ=255)

         agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACa                                                                                                                                         >  1:707462/1‑141 (MQ=255)

         agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                                        >  1:704799/1‑136 (MQ=255)

         agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCACTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                                        <  1:625419/141‑6 (MQ=255)

         agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaat                                                                                                                                        <  1:687557/141‑6 (MQ=255)

            cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                     >  2:272426/1‑133 (MQ=255)

            cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                     >  2:139467/1‑133 (MQ=255)

            cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataat                                                                                                                                     >  1:929096/1‑133 (MQ=255)

             gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                    >  2:743584/1‑132 (MQ=255)

             gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                    <  2:134626/141‑10 (MQ=255)

             gCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatt                                                                                                                                    <  1:49383/141‑10 (MQ=255)

             gCCGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TTA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AAt                                                                                                                                   <  2:494823/141‑1 (MQ=255)

              ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATa                                                                                                                                  >  2:891788/1‑141 (MQ=255)

              ccTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGCAGACGCA‑‑GTACaa‑‑aaga                                                                                                                                  >  2:891759/1‑139 (MQ=255)

              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                   >  2:945098/1‑131 (MQ=255)

              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataatta                                                                                                                                   >  1:548759/1‑131 (MQ=255)

              ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTC‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGCTGTGCTGTTGCGTATTATGCATATTGGGTTAGTAAACGCA‑‑ACAaaa‑‑actga                                                                                                                                 >  1:82009/1‑137 (MQ=255)

               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaccaaataattaa                                                                                                                                    >  1:354737/1‑129 (MQ=255)

               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattca                                                                                                                                  >  2:689590/1‑130 (MQ=255)

               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                                  <  2:864551/141‑12 (MQ=255)

               cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaa                                                                                                                                  >  1:322458/1‑130 (MQ=255)

                tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCAA‑TAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaaa                                                                                                                                 >  1:526362/1‑129 (MQ=255)

                 ttctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTCCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                <  2:695177/140‑14 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGTTTAGTAAACGCa‑‑acaaataattgata                                                                                                                                >  1:714068/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatg                                                                                                                                >  1:699730/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                <  2:800637/141‑14 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  1:167979/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  1:632144/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGACTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  1:283630/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGCAGGCATTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  1:448853/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTGCACATCATAGCACCAGTCATGAGTGGCTGTGCTGTTGCGTATTATACTGATTGGCTTAGTAGACGCa‑‑acaaataattaaaa                                                                                                                                >  1:95698/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattcata                                                                                                                                >  2:327059/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  1:230164/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑TTGA‑TAATCATTATCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGCTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  2:649384/1‑128 (MQ=255)

                 aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaata                                                                                                                                >  1:839580/1‑128 (MQ=255)

                  tctctcATGGGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                                               <  2:365915/141‑15 (MQ=255)

                  tctctcAGGCGTCAATTTAGACGCAGAGAGGCGGCG‑TA‑CAA‑G‑GTGA‑TGCTTATTTTAGTTCACAACATAGAACCAGTCATCAGTGGATGTGCCATTGCGTTTTTTGCTTATTGGCTAAGTTGGCGCA‑‑ATCCAA‑‑ACTAGGTg                                                                                                                              >  2:453324/1‑141 (MQ=255)

                  tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTGGACGCA‑‑ATACAG‑‑AAtttgtg                                                                                                                              >  2:657552/1‑136 (MQ=255)

                  tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑TTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa‑‑acaaataattaatag                                                                                                                               >  1:553255/1‑127 (MQ=255)

                   ctttcAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTA‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAc                                                                                                                            <  2:172643/138‑1 (MQ=255)

                       cAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑AGACAA‑‑AATAGGTGAC‑ata                                                                                                                        >  1:736943/1‑141 (MQ=255)

                        atgCGTAAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCCATTGCGGTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑atat                                                                                                                       <  1:550462/139‑1 (MQ=255)

                         ggCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑Atat‑a                                                                                                                     <  2:398467/141‑1 (MQ=255)

                          gCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTTGATGCA‑‑GTACAA‑‑AATAGGTGAC‑ATAT‑Ag                                                                                                                    >  1:679796/1‑141 (MQ=255)

                           cGTTAATTTAGACGCAGGGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGc                                                                                                                   <  1:309534/141‑1 (MQ=255)

                            gTCAATTTTGACGCAGGGAGGAAGTT‑GA‑TAT‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTATTTCGTTTTTTTCTTATTGGCTAAGTATACGCA‑‑ATACAA‑‑TATAGGTGAC‑ATAT‑AGcc                                                                                                                  <  2:90529/141‑1 (MQ=255)

                               aaTTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCa                                                                                                              <  1:101620/141‑1 (MQ=255)

                                 tttAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GGGA‑TGCTTATGTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGAA‑‑ATACAC‑‑CATAGGAGAC‑ATAT‑AGCC‑GCAcc                                                                                                            <  2:693638/141‑1 (MQ=255)

                                   tAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑CTACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCCCCaa                                                                                                          >  1:377853/1‑141 (MQ=255)

                                     gACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                        <  1:883763/141‑1 (MQ=255)

                                     gACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATa                                                                                                        <  1:781378/141‑1 (MQ=255)

                                       cGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATaaa                                                                                                      <  1:725120/141‑1 (MQ=255)

                                           gagagGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATc                                                                                                  >  1:441603/1‑141 (MQ=255)

                                               ggaggTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCt                                                                                              >  1:348371/1‑141 (MQ=255)

                                                   gtg‑tA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCAca                                                                                          >  2:185922/1‑140 (MQ=255)

                                                            a‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCaa                                                                                   >  1:757352/1‑141 (MQ=255)

                                                                  gA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATGAAA‑‑AATAGGTGAC‑AAAT‑AGCC‑GCACCAATAAAAATCCCCTCACTAACGCAAATAg                                                                               >  2:345692/1‑141 (MQ=255)

                                                                     tGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTg                                                                             <  2:500526/141‑1 (MQ=255)

                                                                       cTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAACAGTGAg                                                                           <  2:520782/141‑1 (MQ=255)

                                                                         tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         >  1:599473/1‑141 (MQ=255)

                                                                         tATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         <  1:494724/141‑1 (MQ=255)

                                                                         cagtttCGTTCACAGCATAGCACCATGCATCAGTGGCTGTGTCATTGCGTTCTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                                         <  2:134095/138‑1 (MQ=255)

                                                                                 ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                                                 >  1:623820/1‑141 (MQ=255)

                                                                                     catcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑GATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATCGTGCGGGGAATGGTGtata                                                             >  2:41198/1‑141 (MQ=255)

                                                                                         atAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAAATAATTAATAGGCGACTTTATAAGCCAGCATCGA‑AAAAATCCCCACACTACTGC‑AATAGTGTGGGGATTGGTGcat                                                              >  2:20374/1‑141 (MQ=255)

                                                                                         atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTTTCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTa                                                         >  2:78495/1‑141 (MQ=255)

                                                                                              cccAGTCATCAGTGGCTGGGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATAc                                                    <  1:495684/140‑1 (MQ=255)

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

                                                                                                    catcaGTGGCTGTGCCATTGCGTTTTTTTCTTGTTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                                              <  2:22524/141‑1 (MQ=255)

                                                                                                        aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGt                                          >  2:573681/1‑141 (MQ=255)

                                                                                                        aGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGt                                          >  1:11169/1‑141 (MQ=255)

                                                                                                            gCTGTGCCATTGCGTCTTTTTCTCATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTc                                      <  1:764350/141‑1 (MQ=255)

                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                                   >  2:286405/1‑141 (MQ=255)

                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                                   <  1:677448/141‑1 (MQ=255)

                                                                                                               gtgCCATTGCGTTTTTTTCTTATTGGCTAAGGAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa                                   <  1:829688/141‑1 (MQ=255)

                                                                                                                          tttttttCTTATTGGCTAAGTAGACGCA‑‑ATACCA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAGTAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCt                        >  2:288488/1‑141 (MQ=255)

                                                                                                                           ttttttttTATTGGATAAGTAGACGCA‑‑ATACAT‑‑AATAGGTGAC‑ACAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTa                       <  2:377853/141‑1 (MQ=255)

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

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

                                                                                                                                   taattgCTAAGTAAACGCA‑‑ATACAA‑‑AATAGGGGAC‑ATAT‑AGCC‑GCGCTGAGGAAAATCGCCATACTATCGCATATTTTCAGGGGATTGATGTATTATGAAATACTTATTTTCCTTGTCTTAATTATACTGCTAATATTTCt               <  2:394821/136‑1 (MQ=255)

                                                                                                                                        gCTAAGTAGACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTg          <  2:937040/141‑1 (MQ=255)

                                                                                                                                                gACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  >  1:336487/1‑141 (MQ=255)

                                                                                                                                                gACGCA‑‑ATACAA‑‑AATAGGTGAC‑ATAT‑AGCC‑GCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAATTTTTCTTTTTGTAAAatat  >  1:311616/1‑141 (MQ=255)

TAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑GTGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCA‑‑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` `< 13 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`