BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,917	C→T	55.6%	intergenic (‑1000/‑168)	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,917	0	C	T	55.6%	37.8 / 199.3	163	intergenic (‑1000/‑168)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (40/32); new base T (41/47); total (84/79)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

Only 100 of 164 total aligned reads displayed.

GACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTA  >  CP000730/1929779‑1930053

gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATATTACGAGGCATCAGTGGCTGTGCCAt                                                                                                                                        <  2:646762/141‑1 (MQ=255)

gACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                                                                        <  1:691294/141‑1 (MQ=255)

   tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                                                                     <  2:779727/141‑1 (MQ=255)

   tgtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                                                                                     >  1:122779/1‑141 (MQ=255)

   tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                                                                    >  1:208309/3‑141 (MQ=255)

   gatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTGTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTCAGCATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                                                                                    <  2:56404/139‑1 (MQ=255)

    atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                   >  2:810351/2‑141 (MQ=255)

    atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                                                                                                   <  1:637536/140‑1 (MQ=255)

       cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                                                                <  1:110514/141‑4 (MQ=255)

       cTTGTATATTCACAGATGGTATAATAATAGAGGCGCCTATCTTTGAGGCGTCAAGTT‑‑ACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCTGTCATCAATGGCTGTGCTGTTGCGtatt                                                                                                                                >  1:447839/1‑141 (MQ=255)

        ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                                                               >  1:821819/1‑140 (MQ=255)

        ttGTATATTCACAGAAGGTATAATATGAGAGTCGCCTAGCTCTCAGGCTTCAAGTT‑‑ACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGGCAGCAGTGGCTGTGCTGTTGCGtagta                                                                                                                               <  2:166106/141‑5 (MQ=255)

           tataTTTGCCGAGGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt                                                                                                                             <  2:718397/141‑1 (MQ=255)

           tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                            <  1:194973/141‑8 (MQ=255)

           tataTTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                                                            <  2:307152/141‑5 (MQ=255)

            atattgacagatggtATAATAATAGAGTCGCCTAGCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                                                                           <  2:26993/132‑6 (MQ=255)

            atatTTGCCGATGATATAATAGTAGAGTGGCATATCTCTCAGGCGTCAAGTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGTTTTTTTc                                                                                                                            >  1:264965/1‑141 (MQ=255)

               ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCAATCACATCATAGCACCAGTCATCACTGTCAGTGCTGTTGCGTATTATACTTAt                                                                                                                        >  2:145154/7‑141 (MQ=255)

               ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGAAGGTGCATTAC‑T‑T‑TGCTAATCATTTTAGTTCACATCATAGCACTCGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                                                                        >  2:290243/7‑141 (MQ=255)

                 cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                                                                      >  2:115255/5‑141 (MQ=255)

                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                    <  2:914863/139‑1 (MQ=255)

                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                    <  2:337430/139‑1 (MQ=255)

                   cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                                                                    <  2:104301/139‑1 (MQ=255)

                    cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                                                                    <  1:723492/141‑1 (MQ=255)

                    agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCCATT‑‑GACGAAGGGAGGAGGTGCA‑TTT‑CTT‑TGCTAGTCATTTTCGTTGACATCGTAGGACCGGTCGTCGGTGGCTGTGCTGTTGCGTATTATACTTATTGGCt                                                                                                                   <  2:667089/140‑1 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                               >  1:104876/3‑141 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                                                               >  2:404668/3‑141 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTACTTGCTTAGt                                                                                                               >  2:491458/3‑141 (MQ=255)

                        ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATGATACTTATTGGCTTAGt                                                                                                               <  1:219423/139‑1 (MQ=255)

                        ggtacaataatAGAGTCGCCTATCTCTCAGCCGTCAATT‑‑GACGAAGAGAGGAGATGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAga                                                                                                               >  2:554136/6‑140 (MQ=255)

                          tataATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                                                                              >  1:325535/1‑141 (MQ=255)

                          tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                             >  2:133234/1‑140 (MQ=255)

                          tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTaa                                                                                                             >  2:698714/1‑140 (MQ=255)

                          tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                                                                             >  1:255172/1‑140 (MQ=255)

                           ataataataGAGTCGCCTATCTCTCAGGCGTCTATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATATCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                            <  1:609189/141‑3 (MQ=255)

                           ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                            <  1:517830/141‑3 (MQ=255)

                           ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                                                            >  1:777417/1‑139 (MQ=255)

                           ataataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAg                                                                                                             >  2:202846/1‑141 (MQ=255)

                            taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGa                                                                                                            >  2:65728/1‑141 (MQ=255)

                              ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                                                         >  1:291904/5‑141 (MQ=255)

                               taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAAGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                                                                        <  2:852428/138‑1 (MQ=255)

                                aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                       <  1:870882/139‑1 (MQ=255)

                                aatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑GTTGCTAATCATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                                                                       <  2:671873/139‑1 (MQ=255)

                                  aagagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                                                                     <  2:342829/140‑3 (MQ=255)

                                   tgagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                     <  1:846571/140‑1 (MQ=255)

                                   agagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATa                                                                                                     <  1:846460/141‑1 (MQ=255)

                                   agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                    <  1:790465/141‑4 (MQ=255)

                                   agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                                                                    >  2:630967/1‑138 (MQ=255)

                                    gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGAGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                                                                    >  1:510759/1‑141 (MQ=255)

                                      gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTTAACGCaacaaatt                                                                                                 >  1:676366/1‑135 (MQ=255)

                                      gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGATTAGTAAACGGaacaaata                                                                                                 >  2:590923/1‑135 (MQ=255)

                                      gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGTTTGCGTTTTATACTTATTGGCTTAGTAAACGCaacaaata                                                                                                 >  1:189030/1‑135 (MQ=255)

                                      gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGAATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACACaacaaata                                                                                                 >  2:597656/1‑135 (MQ=255)

                                       tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                <  2:73604/141‑8 (MQ=255)

                                       ccGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                                                                <  1:737207/140‑8 (MQ=255)

                                         gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAAt                                                                                               >  2:432488/1‑141 (MQ=255)

                                         gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                              >  1:889378/1‑132 (MQ=255)

                                         gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                                                              <  2:623208/141‑10 (MQ=255)

                                            tATCTCTCAGGCTTCAATT‑‑GGCGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTTGTTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                           <  2:365251/141‑13 (MQ=255)

                                            tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                           >  1:165115/1‑129 (MQ=255)

                                            tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                           <  2:59483/141‑13 (MQ=255)

                                            tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCGGTTGCGTATTATACTTATTGGCTGAGTAGACACaacaaataattaaa                                                                                           >  1:141676/1‑129 (MQ=255)

                                            tATCTCCCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                                                           <  1:23842/141‑13 (MQ=255)

                                             ttctctcAGGCGTCAATT‑‑GACGGAGAGAGGAGTTGCA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                          <  1:98123/140‑14 (MQ=255)

                                             aTCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGt                                                                                           >  1:5845/1‑141 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaagtaattaaaa                                                                                          >  2:535413/1‑128 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                          >  2:807944/1‑128 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                                                          >  1:606306/1‑128 (MQ=255)

                                             aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatcattaatg                                                                                          >  1:59342/1‑128 (MQ=255)

                                              tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCTTTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTg                                                                                          >  1:802774/1‑141 (MQ=255)

                                              tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattattag                                                                                         <  2:790341/141‑15 (MQ=255)

                                              tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                                                                         <  2:633684/141‑15 (MQ=255)

                                              tccctcggacGTCAGTTTAAGCGAAGAGAGGAGGCGTA‑TAA‑G‑G‑TGATGCGGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTGAGGGGTAAGGAGACGCAAGACAAAATAGGTg                                                                                          >  2:873909/10‑141 (MQ=255)

                                               ctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGATCACGCGGGTTCACGTAATAGTCCCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTGCTGATTTGCTAAGTAGACGCAATACAAAATAGGTGa                                                                                         >  1:901647/1‑141 (MQ=255)

                                                       cGTAAATTTAGACGCAGAGCGGAGGTGTA‑TAA‑G‑G‑GGATGCTTATTTTGGTTCACATCATAGCACCAGTCATCAGTGGCTGTGACCTTGCGTTTTGTTCTTAGTAGTTAAGTAGAGGCAATACAAAGACGGTGAGATATAGc                                                                                 >  2:467110/1‑141 (MQ=255)

                                                              ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCa                                                                          >  1:610001/1‑141 (MQ=255)

                                                                tgACTCAGAGAGGAGGTGTA‑TAA‑G‑A‑TGATGCTGATCTTCGTGCACACGATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAAt                                                                        <  1:802507/140‑1 (MQ=255)

                                                                     cAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                   <  2:714853/141‑1 (MQ=255)

                                                                     cAGAGAGGAGGGGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                                                                   <  1:368091/141‑1 (MQ=255)

                                                                          aggaggAGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAAATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc                                                              >  2:597906/1‑141 (MQ=255)

                                                                               gtgtA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt                                                         <  1:753867/141‑1 (MQ=255)

                                                                                tgtA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAATTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa                                                        >  2:736959/1‑141 (MQ=255)

                                                                                       a‑G‑G‑TGATGCTTATTTACGTTCACATCTTAGCACCTGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTCGTCGCAATACAAAATGGGTGACATATAGCCGCACCAATATAAATCCCCTCACTACCGCaa                                                  <  1:551202/141‑1 (MQ=255)

                                                                                               aTGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAgg                                             >  1:229113/1‑140 (MQ=255)

                                                                                                 gCTTATTTTCGTTCACATCATAGCACCAGTCATAAGTGGCTGGGGCATTGCGTGTTTTTCTTATTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCACATCACTACCGCAAATAGTGa                                           >  1:352959/1‑141 (MQ=255)

                                                                                                    tATTTTCGTTTACATCATAGCACCAGTCATCAGTGGCTTTGCCATTGCGTTTTTTCCTTATTGTCTAAGTAGACGCAATACAAAATAAGTCACATATAGCCGCACCAATAAAAAACCCCTCACTACCGCAAATAGTGAggg                                        <  1:777971/141‑1 (MQ=255)

                                                                                                    tATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTTTTGTCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATACAAATCCCCTCACTCCCCCCAACAGCGAggg                                        >  2:106029/1‑141 (MQ=255)

                                                                                                     tttttCGTTCCCATCATAGCACCGGTCGTCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAACTCCCCTCACTACCGCAAATAGTGAgggg                                       <  2:229402/140‑1 (MQ=255)

                                                                                                     aTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAgggg                                       >  2:649720/1‑141 (MQ=255)

                                                                                                      ttttCGTTCACATCAAAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                      <  2:140705/141‑1 (MQ=255)

                                                                                                      ttttCGGTCACATCATAGGACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGa                                      <  2:208213/141‑1 (MQ=255)

                                                                                                            ttCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGTTTGgtg                                >  1:248434/1‑141 (MQ=255)

                                                                                                              cacaTCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGta                              >  1:401170/1‑141 (MQ=255)

                                                                                                                    atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATGAGTa                        >  2:115433/1‑141 (MQ=255)

                                                                                                                    atAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAAAAAATAGGTGACATATAGCCGCACCAATAAAAATCACCTCACTCCCGCAAAAAGTGAGGGGATTGGACTATAAGTa                        >  2:254711/1‑141 (MQ=255)

                                                                                                                             gTCATCAGTGGCTGTGGCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAt               <  1:688181/141‑1 (MQ=255)

                                                                                                                             gTCAGCAGTGGCGGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTACACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCCCCATACCGGAACTAGTGAGGGGATTGGTGTATAAGTAAATCCTTAt               >  1:642473/1‑141 (MQ=255)

                                                                                                                              tcatcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCCCTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTAtt              <  2:610001/141‑1 (MQ=255)

                                                                                                                               catcaGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTTAATACTTGttt             >  1:242302/1‑141 (MQ=255)

                                                                                                                                          gtgCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAAGACAAAATAGGGGACATATAGCCGCACCAATAAACATCCCCTCACTACCGAAAGTAGAGAGGGGGTTGGTGTATAAGTAAATACTTATTTTACTTGGCTTa  >  1:423806/1‑141 (MQ=255)

GACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGAGGGGATTGGTGTATAAGTAAATACTTATTTTCGTTGTCTTA  >  CP000730/1929779‑1930053

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