BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,928	T→A	54.5%	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	54.5%	24.0 / 174.4	126	intergenic (‑1011/‑157)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (25/31); new base A (35/33); total (60/66)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

Only 100 of 130 total aligned reads displayed.

TGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAA  >  CP000730/1929788‑1930066

tGTTTATTTGCCGAGGATGTAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttt                                                                                                                                            <  2:836436/141‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

           aggtggtataATAAGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                                <  2:514985/135‑1 (MQ=255)

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

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

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

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

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

              tgttataAGAATAGAGTCGCCTAACTCTCCGGCGTCAATT‑‑GACGAAGAGAGGAGGTTCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                             <  2:303127/138‑1 (MQ=255)

              tggtataATACTAGATGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                                                                                             <  2:756122/138‑4 (MQ=255)

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

              gggtataAGAGTAGAGTCGCCTATCTTTCGTGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGTACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                                                                             <  2:498392/138‑1 (MQ=255)

               ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                            <  1:361531/139‑1 (MQ=255)

               ggtataATAATAGAGTCGCCTATCCCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCTCATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                            <  1:241552/139‑1 (MQ=255)

               ggtacaataatAGATTCGCCTATCTCTCAGGCGTCCATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                                            <  2:555807/136‑1 (MQ=255)

                gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCATTTTCGTTCACATGATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTGGCttggta                                                                                                                           >  2:197150/2‑136 (MQ=255)

                atatAATAGTAGAGTCGACTCTCTAAAGGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                                                                            >  2:385095/1‑141 (MQ=255)

                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                          <  2:233635/141‑2 (MQ=255)

                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                          >  2:365784/1‑140 (MQ=255)

                 tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                                          >  2:115692/1‑140 (MQ=255)

                  ataataataGAGTCGCCTATCTCTCCGGTGTAAA‑‑TGGACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTCTACTTATGGGCtttgtaaa                                                                                                                         <  2:800948/141‑8 (MQ=255)

                  ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                                         >  1:156826/1‑139 (MQ=255)

                  ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATGTTCGTTCGCAGCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTCAGTAg                                                                                                                          >  2:505018/1‑141 (MQ=255)

                   taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                                                                                        >  2:191068/1‑141 (MQ=255)

                    aataatAGAGTTGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                       >  1:307755/1‑141 (MQ=255)

                    aataatAGAGTGGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                                                                       <  2:239390/141‑1 (MQ=255)

                     ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                                      <  1:810716/141‑1 (MQ=255)

                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                    <  1:111478/139‑1 (MQ=255)

                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                                    >  2:650488/3‑141 (MQ=255)

                       aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCCTAGTAAACGCaa                                                                                                                    >  2:39152/3‑141 (MQ=255)

                        atatagTCCCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                   <  1:107199/137‑2 (MQ=255)

                        atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                                                                                   >  2:396951/2‑140 (MQ=255)

                          agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                                  >  2:836573/1‑141 (MQ=255)

                           gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                >  2:642553/1‑137 (MQ=255)

                           gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                                                                                <  1:376056/141‑5 (MQ=255)

                           gagTCGCCGATGTGTCAGGCTGCAATCGAGGCGCAGAGGGGAGGTGTA‑TAA‑G‑G‑GGA‑TGCTTATTTTCGTTCACATCATAGCTCCAGTCATTAGTGGCTGTGCCAGTGCGTTTTTGTCTGATGGGCTAAGTAGACGCAATAc                                                                                                                 <  2:171833/141‑1 (MQ=255)

                            agTCGCCTATCTGTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                               >  1:592453/1‑136 (MQ=255)

                            agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCGGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaact                                                                                                               >  2:90573/1‑136 (MQ=255)

                            agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                                                                               >  2:198806/1‑136 (MQ=255)

                             gTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGG‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                                                                               >  2:325926/1‑141 (MQ=255)

                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCAACAaaaa                                                                                                              >  1:13125/1‑141 (MQ=255)

                             gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                              <  2:347035/141‑7 (MQ=255)

                              tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAGTACaaa                                                                                                              >  1:26459/1‑141 (MQ=255)

                                gCCTGTCTCTCAGGCGTCAATT‑‑GGCGAGGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                           <  2:133413/141‑10 (MQ=255)

                                gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                           <  2:568671/141‑10 (MQ=255)

                                gCCTATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                                           <  2:641506/141‑10 (MQ=255)

                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                          >  1:311927/1‑131 (MQ=255)

                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                                                                                          <  1:814002/141‑11 (MQ=255)

                                 ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCCTAGTAAACGCaccaaataatta                                                                                                          >  1:191735/1‑130 (MQ=255)

                                  cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAGTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                         <  2:571334/141‑12 (MQ=255)

                                  cTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                                                                         <  1:775842/141‑12 (MQ=255)

                                  atATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTACTAAACGCaccaaataattaa                                                                                                         >  2:514129/2‑129 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                        <  2:59972/141‑13 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaaa                                                                                                        >  2:284760/1‑129 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCTTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                        >  1:781601/1‑129 (MQ=255)

                                   tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                        <  2:42152/141‑13 (MQ=255)

                                   tATCTCTCAGGCGTCAA‑‑TAGACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                                        <  1:735220/141‑13 (MQ=255)

                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                       >  2:198995/1‑128 (MQ=255)

                                    aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGC‑TAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTAGTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                                       >  2:751092/1‑128 (MQ=255)

                                      ctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑A‑G‑TGA‑TGCTTATCTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATCCAAAATACGTgg                                                                                                      >  2:139871/1‑140 (MQ=255)

                                          cAGGCTTGAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                                                                                  <  2:663794/141‑1 (MQ=255)

                                            ggCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGAGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACAtata                                                                                                >  1:788125/1‑141 (MQ=255)

                                                tCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCg                                                                                            <  2:122296/141‑1 (MQ=255)

                                                     ttAGCCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                                                                       <  1:331499/141‑1 (MQ=255)

                                                     ttAGCCGCAGAGAGGAGGTGTA‑CAA‑G‑G‑TGC‑TGCTTATTTTCGAAAGAGTCATCGTTTATGGCCCCGATCGCTGTGCCATTGCGTTTTGTGCTAGTTTTCCAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                                                                       <  2:804414/141‑1 (MQ=255)

                                                      tggACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCaa                                                                                      <  1:423127/139‑1 (MQ=255)

                                                           gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                                 <  2:515932/141‑1 (MQ=255)

                                                           gCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                                                                                 <  2:576682/141‑1 (MQ=255)

                                                                  ggaggTGTA‑TAA‑G‑G‑TGA‑TGCTGATTTTCGTTCGCATCATAGAACCAGTAATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                                                                          <  2:610102/141‑1 (MQ=255)

                                                                    aggTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTTTTTTTTGTTTAGTAGACGCAATACAAAATGGGTGACATATACCCGCGCAAGTACAAATCCCCTcc                                                                        >  2:787829/1‑140 (MQ=255)

                                                                     tgtgtA‑TAA‑C‑G‑TTA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                                                                       <  1:52103/140‑1 (MQ=255)

                                                                             aa‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTCCCGCa                                                                >  2:22551/1‑141 (MQ=255)

                                                                                      a‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGt                                                          >  2:522931/1‑141 (MQ=255)

                                                                                           ttATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgg                                                      <  2:686996/141‑1 (MQ=255)

                                                                                            tagtttCGTTCACATCATAGCACCAGTCATCAGTGGTTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAggg                                                     <  2:136150/138‑1 (MQ=255)

                                                                                              ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATGGTGAggggg                                                   >  1:669930/1‑140 (MQ=255)

                                                                                              ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                   <  2:292295/141‑1 (MQ=255)

                                                                                              ttttCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                                   <  1:314793/141‑1 (MQ=255)

                                                                                                         atcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATaa                                        <  2:713905/141‑1 (MQ=255)

                                                                                                          tcatAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATTTTTCTTATTGGCTTAGTAGACGCAATACAAAATAGGTGACATATAGCCCCACCAATAAAAATCCCCTCACTACCCACAATAGCGAGGGGATTGGTTTACaaa                                       >  2:146387/1‑140 (MQ=255)

                                                                                                               gCACCAGTCGGCACTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAAt                                  <  1:603952/141‑1 (MQ=255)

                                                                                                                cACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCGCCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATa                                 >  1:653156/1‑141 (MQ=255)

                                                                                                                cACCAGTCATCAGTGCCTGTGCCATTGCGTTTTTTTCATATTGCCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGAGATTGGTGTATAAGTAAATa                                 <  1:589259/141‑1 (MQ=255)

                                                                                                                cACCAATCATCACGGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAAAGCAATCGATTATGTTATTGGTGTATAAGTAAATa                                 <  1:405245/141‑1 (MQ=255)

                                                                                                                    aGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGGGAGGGGAGTGGTGTATAAGTAAATACttt                             >  1:548900/1‑140 (MQ=255)

                                                                                                                       catcaGTTGGTGTTCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAttt                          <  2:564504/141‑1 (MQ=255)

                                                                                                                        atcaGTGGCTGTGCCAGTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATGGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTTTAAGTAAATACTTAtttt                         >  1:393119/1‑141 (MQ=255)

                                                                                                                             tGGCTGTGCCATTGCGTTTTTTTTTTTTTTTTTAAGTAGACGCAATACAAAGTAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTTCCGCAAATAGTGAGGGGGTTGGTGTATAAGTACATACTTTTTTTCGttt                    >  2:572003/1‑140 (MQ=255)

                                                                                                                                 tgtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTg                    <  1:260277/137‑1 (MQ=255)

                                                                                                                                  gtgCCATTGCGTTTCTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTa               >  1:546464/1‑141 (MQ=255)

                                                                                                                                               tttttCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa  <  1:229317/141‑1 (MQ=255)

                                                                                                                                               tttttCTTACTGGCTAAGTACACGCAATACAACATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGTAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTaa  <  1:139871/141‑1 (MQ=255)

TGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGA‑TGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTAATTATACTGCTAA  >  CP000730/1929788‑1930066

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