BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,928	T→A	67.8%	intergenic (‑1011/‑157)	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,928	0	T	A	67.8%	129.6 / 96.5	107	intergenic (‑1011/‑157)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (16/18); new base A (42/29); total (60/47)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.48e-01

Only 100 of 109 total aligned reads displayed.

TGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAT  >  CP000730/1929788‑1930067

tGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAT‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                                                                             <  1:424911/141‑1 (MQ=255)

      ttcacagatggtATAATGATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTGCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                      >  1:469921/7‑141 (MQ=255)

      ttcacagatggtATAATAATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑GACGAAAAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                                      >  1:429667/7‑141 (MQ=255)

       tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                                                                                     >  2:91718/6‑141 (MQ=255)

        cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTCTCCTTCTTg                                                                                                                                    >  1:458604/5‑141 (MQ=255)

        cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                                    >  2:609170/5‑141 (MQ=255)

        cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTACACTTATTg                                                                                                                                    >  2:534203/5‑141 (MQ=255)

         acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                                                                                   >  2:456364/4‑141 (MQ=255)

          ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                   <  2:161790/141‑1 (MQ=255)

          ccGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgg                                                                                                                                   <  1:448480/141‑1 (MQ=255)

          caggtggtataATAATAGAGTCGCCTATCTCTCATGCGTCAATT‑‑GACGGCGAGAGGGGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTTCGTATTATACTTATTGGc                                                                                                                                  <  2:285528/134‑1 (MQ=255)

           agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                 <  1:437605/140‑1 (MQ=255)

            tatggtATAATAATGGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCTTAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                <  2:320702/140‑2 (MQ=255)

            gattgaataacaATATAGTCGCCTATCTCTCAGGCGTCACTT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTACTC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                <  2:588350/135‑2 (MQ=255)

            gatggtATAATAATAGAGTCGCCTATCTCTCAGGCTTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                <  2:536543/141‑2 (MQ=255)

            gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                >  2:508681/1‑140 (MQ=255)

            gatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCtt                                                                                                                                >  2:272694/1‑140 (MQ=255)

              tgttataATAATAGAGTCGTCTATCTCTCAGGCGTCAA‑‑TGGACGAATAGAGGAGGTCCATTAC‑T‑T‑TGCTAATC‑CT‑‑TTTCGTTCGCATCATAGCACCAGTCAGCAGTGGCTGTGCTGTGGCGTATTATACTTATTGGCTTAg                                                                                                                              <  2:270883/138‑1 (MQ=255)

              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                              <  1:612933/138‑4 (MQ=255)

              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                              <  2:197083/138‑4 (MQ=255)

              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                              <  2:415240/138‑4 (MQ=255)

              tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                              >  1:649279/4‑141 (MQ=255)

              tggtataATAATAGAGTCGCCTATCTCTCAGGAGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATTATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                              <  1:446155/138‑4 (MQ=255)

              gagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                              <  1:638521/138‑4 (MQ=255)

               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTGTACTTATTGGCttggt                                                                                                                             >  1:573848/3‑137 (MQ=255)

               gatATAATAGTAGAGCCGCCTATCTCTCAGGCGACAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTAGTTCACACCATAGAACCAGGCAACAGTGGCTGTGGCATTGCGTGTTTTTcatcttcactaat                                                                                                                              >  2:625437/1‑129 (MQ=255)

                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                                                                            >  1:496047/2‑141 (MQ=255)

                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGCCATCAGTGGCTGTGCTGTTGCGTATTATACTGATTGGCTTAGTa                                                                                                                            >  1:74474/2‑141 (MQ=255)

                atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                             >  1:352569/1‑141 (MQ=255)

                atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                             >  1:652215/1‑141 (MQ=255)

                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCGATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                           >  2:280738/1‑140 (MQ=255)

                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                          >  2:493284/1‑139 (MQ=255)

                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                          >  1:354903/1‑139 (MQ=255)

                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑‑TGCTA‑ATCATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                         >  2:388492/1‑141 (MQ=255)

                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                         <  2:29742/141‑4 (MQ=255)

                   aaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                                                                                         <  1:99481/140‑4 (MQ=255)

                    aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaaag                                                                                                                        >  2:647690/1‑137 (MQ=255)

                    aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                        >  2:111394/1‑141 (MQ=255)

                     ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTGAGTAGACg                                                                                                                        >  2:122135/1‑141 (MQ=255)

                     ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACg                                                                                                                        <  2:290132/141‑1 (MQ=255)

                      taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                      >  1:566730/4‑141 (MQ=255)

                      taatagAGCCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                                      <  2:376500/138‑1 (MQ=255)

                      gaatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGGACATTAC‑T‑T‑TGCTAATC‑AT‑‑TTGCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTGGTAAACgga                                                                                                                      >  2:96135/4‑139 (MQ=255)

                       agtagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCa                                                                                                                      <  1:206125/141‑1 (MQ=255)

                       aatagGGTCGCCTATATCTCAGGCGTCAATT‑‑GACGAACAGAGGCGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAATGCaa                                                                                                                     <  1:469100/139‑1 (MQ=255)

                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                     <  1:613537/139‑1 (MQ=255)

                        atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                    >  2:498300/2‑140 (MQ=255)

                         tagAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                    <  1:407826/141‑1 (MQ=255)

                         tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGCGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                                   <  1:292403/141‑3 (MQ=255)

                         tagAGTCGCCTATCCCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAt                                                                                                                    <  1:53886/141‑1 (MQ=255)

                          agctgctccTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAAGC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                  <  1:77642/134‑4 (MQ=255)

                          agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                  >  1:270264/1‑138 (MQ=255)

                          agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                  >  2:297534/1‑138 (MQ=255)

                          agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                  >  2:571661/1‑138 (MQ=255)

                          agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TTGCGAATCAT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                                  >  2:629138/1‑138 (MQ=255)

                           gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                 <  1:617180/141‑5 (MQ=255)

                             gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                                                                                <  2:400338/141‑1 (MQ=255)

                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                               <  1:14081/141‑7 (MQ=255)

                              tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                              >  1:162660/1‑134 (MQ=255)

                              tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                              <  1:618995/141‑8 (MQ=255)

                              tCGCCTATCTCCCAGGCCTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                              <  1:8777/141‑8 (MQ=255)

                               cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                             >  1:261191/1‑133 (MQ=255)

                               cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                             >  2:379385/1‑133 (MQ=255)

                               cGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                                                                                             >  1:358919/1‑133 (MQ=255)

                                cccTATCTCTCGGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCAAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                            <  2:59850/140‑10 (MQ=255)

                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTTTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                           >  2:173791/1‑131 (MQ=255)

                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                           >  2:641759/1‑131 (MQ=255)

                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCACCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                           <  1:97113/141‑11 (MQ=255)

                                 ccTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                           <  1:602673/141‑11 (MQ=255)

                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                          >  1:144924/1‑130 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatttag                                                                                                         >  2:164795/1‑129 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                         >  1:508557/1‑129 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTATTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                         >  2:13627/1‑129 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                         >  1:93114/1‑129 (MQ=255)

                                    ttctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAT‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACGTCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                        <  1:450432/140‑14 (MQ=255)

                                    aTCTCTCAGGCGTCAATT‑‑GCCGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatg                                                                                                        >  1:172730/1‑128 (MQ=255)

                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                        >  1:168804/1‑128 (MQ=255)

                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                        >  2:580444/1‑128 (MQ=255)

                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATC‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                        >  2:577827/1‑128 (MQ=255)

                                             gCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAg                                                                                                >  2:526166/1‑141 (MQ=255)

                                                 cAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                                                                                            >  1:464158/1‑141 (MQ=255)

                                                  aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                                                                                           <  2:454536/141‑1 (MQ=255)

                                                    tttGGACGCAGACAGGAGGTGTA‑TCA‑G‑G‑TGACGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                                                                                         <  2:63656/141‑1 (MQ=255)

                                                                  ggaggTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                           >  2:26504/1‑141 (MQ=255)

                                                                           ttaA‑G‑G‑TGATGCGT‑CT‑‑TTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATCGGCTAAGTAGACGCAATACAAAATCGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCg                                                                   <  1:86504/140‑1 (MQ=255)

                                                                                g‑g‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa                                                               <  1:280709/141‑1 (MQ=255)

                                                                                      aTGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTGATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGt                                                           >  2:190265/1‑141 (MQ=255)

                                                                                             aT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGcccctcgagagggc                                                     >  2:624889/1‑128 (MQ=255)

                                                                                             aT‑‑TTTCGTTCACATCATAGCACCAATCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAAAAGTGAgggg                                                     >  1:248261/1‑141 (MQ=255)

                                                                                                   tctgtcatatcataGCACCAGTCATCATTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTg                                                 <  2:663060/137‑1 (MQ=255)

                                                                                                      ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGgtg                                              <  1:338440/141‑1 (MQ=255)

                                                                                                           atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTCAGTAGACGCAATACAAAATAGGTGACATTTAGCCGCACCAATAAGAATCCCCTCACTTCCGCAAATAGTGAGGGGGTTGGTGTATaa                                         >  2:135692/1‑141 (MQ=255)

                                                                                                                       gTCACCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAACATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt                             <  1:133167/141‑1 (MQ=255)

                                                                                                                              gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGtt                      <  2:352569/141‑1 (MQ=255)

                                                                                                                              gTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGtt                      <  2:652215/141‑1 (MQ=255)

                                                                                                                                  cTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCt                  >  1:226948/1‑141 (MQ=255)

                                                                                                                                   tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATCGCCGCACCAATAAAAATCCCCTCACTACCGCATATAGTGAGGGGTTTGGTGTATAAGTAAATACTTATTTTCGTTGTCtt                 >  2:658784/1‑141 (MQ=255)

                                                                                                                                   tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCGGTCCCGGCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGGCtt                 >  1:153638/1‑141 (MQ=255)

                                                                                                                                         aTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATtat           <  1:358101/141‑1 (MQ=255)

                                                                                                                                                  ttttCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAt  <  2:507327/141‑1 (MQ=255)

TGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTT‑AT‑‑TTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAAT  >  CP000730/1929788‑1930067

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