BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,862:1	+T	37.9%	intergenic (‑945/‑223)	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,862	1	.	T	37.9%	128.4 / 294.5	161	intergenic (‑945/‑223)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (45/55); new base T (36/25); total (81/80)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.76e-01

Only 100 of 162 total aligned reads displayed.

TAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTA  >  CP000730/1929723‑1929999

   gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAATAGAATCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGAGTA‑Taa                                                                                                                                            >  1:324369/1‑140 (MQ=255)

    tgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                                          >  1:177297/1‑140 (MQ=255)

          cGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑g                                                                                                                                 <  1:107495/140‑1 (MQ=255)

             ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCtt                                                                                                                              >  2:116665/1‑140 (MQ=255)

             ttCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCACCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCtt                                                                                                                              >  2:180334/1‑140 (MQ=255)

               cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTAt                                                                                                                            <  1:416060/140‑1 (MQ=255)

                aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCTTATCTCTCCGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTAtt                                                                                                                           <  1:28187/140‑1 (MQ=255)

                aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGCGTA‑TAC‑G‑G‑TGAT‑GCTTAtt                                                                                                                           >  2:420496/1‑140 (MQ=255)

                 aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTAttt                                                                                                                          >  2:127954/1‑140 (MQ=255)

                  tCTATCAGCTAGAATTATATAACTTATCAGAAGCAATTTACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTAtttt                                                                                                                         <  1:285190/140‑1 (MQ=255)

                    tATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCg                                                                                                                       <  1:327701/140‑1 (MQ=255)

                                              aGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTA‑TTTCGTTCACATCATAGCACCAGTCATCAGTg                                                                                            <  2:561362/140‑1 (MQ=255)

                                               atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAATCAACAGTgg                                                                                           >  2:194350/15‑140 (MQ=255)

                                               atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCAACAGTgg                                                                                           <  2:592135/126‑1 (MQ=255)

                                               atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                                                                                           >  2:553749/15‑140 (MQ=255)

                                                tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATAATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                          >  1:80227/14‑140 (MQ=255)

                                                  ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                        >  1:609342/12‑140 (MQ=255)

                                                   gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGTGTTGACTAACTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCAGTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                       <  1:314332/130‑1 (MQ=255)

                                                    aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                      <  1:506592/131‑1 (MQ=255)

                                                     agtacatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                     <  2:484059/132‑1 (MQ=255)

                                                     aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                      <  1:519469/140‑1 (MQ=255)

                                                     aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                      >  2:506101/1‑140 (MQ=255)

                                                     aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCggtg                                                                                      >  2:343325/1‑140 (MQ=255)

                                                     aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGAGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GATTGTTTTCGTTCACATCATAGCGCCAGACATCAGTGGCtgtg                                                                                      >  1:179563/1‑140 (MQ=255)

                                                     aGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGACGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                      >  2:240119/1‑140 (MQ=255)

                                                      gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCCCCCCAAACCACCAGTCATCAGTGGCTGTGct                                                                                    <  1:490424/133‑2 (MQ=255)

                                                      gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑A‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                    >  1:478801/8‑139 (MQ=255)

                                                        agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GTCGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                  <  1:230962/135‑4 (MQ=255)

                                                        agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                  <  1:118944/135‑4 (MQ=255)

                                                        agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgt                                                                                  <  2:25474/135‑4 (MQ=255)

                                                         gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                                 >  2:103259/5‑136 (MQ=255)

                                                         gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtt                                                                                 <  1:246185/136‑5 (MQ=255)

                                                         aCTGCGCTTCTATATTTGCCGATGATATAATAGTAGAGTCGCCTATTTCTCAGGCGTTAATTTAGACGTAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAt                                                                                  <  1:315532/140‑1 (MQ=255)

                                                          cTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCAtt                                                                                 >  1:407211/1‑140 (MQ=255)

                                                          atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                >  1:602345/4‑135 (MQ=255)

                                                          atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                >  1:236127/4‑135 (MQ=255)

                                                           tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAA‑TTA‑ACGAAGAAAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgtttc                                                                               >  1:554713/3‑134 (MQ=255)

                                                           tatgCTTGTATATTCACAGAGGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                               <  2:267487/138‑7 (MQ=255)

                                                            gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTTTCGCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGc                                                                               <  1:123323/140‑1 (MQ=255)

                                                            atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGAACCAGTCATCAGTGGCTGTGCTGTTGCg                                                                              >  1:336112/2‑140 (MQ=255)

                                                                ttGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                          <  1:500520/140‑4 (MQ=255)

                                                                 tGTATATTCACAAATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                         <  1:59706/140‑5 (MQ=255)

                                                                   tataTTTGCCGATGAAATAATAGGAGAGTTGCCTATCTCTCAGGCGTCAATTTAGTCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCCATTGCGtttttt                                                                        <  1:261145/140‑1 (MQ=255)

                                                                       ttcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                   >  1:378478/7‑140 (MQ=255)

                                                                        tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TGACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                  <  2:226637/135‑1 (MQ=255)

                                                                         cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                 >  1:262654/5‑140 (MQ=255)

                                                                          acagatagtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTg                                                                <  2:152124/132‑1 (MQ=255)

                                                                            agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                              <  2:36775/139‑1 (MQ=255)

                                                                            agatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                              <  2:383188/139‑1 (MQ=255)

                                                                             gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                              <  1:437332/140‑1 (MQ=255)

                                                                             gatgatATAATAGTAAAGTCGCCTATCTCTCAGGCGTCAATTTATACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                              <  1:110705/140‑1 (MQ=255)

                                                                               tgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCtg                                                            <  2:576950/140‑2 (MQ=255)

                                                                                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTACTATACTTATTGGCTTAg                                                          <  2:335683/138‑1 (MQ=255)

                                                                                 ttataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                         <  2:615958/139‑1 (MQ=255)

                                                                                 gttaaataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCCTAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                         <  2:332591/136‑1 (MQ=255)

                                                                                 atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                          <  2:658037/140‑1 (MQ=255)

                                                                                 atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAg                                                          <  2:495638/140‑1 (MQ=255)

                                                                                  tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAA‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCATTCATCAGTGGCTTTTCTGTTGCGTATTATACTTATTGGCTTAGTa                                                        >  2:157312/1‑140 (MQ=255)

                                                                                   ttaataGTAGAGCCGCTGATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTa                                                        >  2:146383/2‑140 (MQ=255)

                                                                                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaa                                                       >  1:374131/1‑139 (MQ=255)

                                                                                    taataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATAATAGCACCAGTTATCAGTGGCTGTGCCAATGCATTTTTTTCTTATTGGCTAAGTAg                                                       <  1:28062/140‑1 (MQ=255)

                                                                                     aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                                     <  1:540592/140‑4 (MQ=255)

                                                                                     aataatAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAAc                                                     >  1:667840/1‑140 (MQ=255)

                                                                                      ataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCT‑TTA‑G‑T‑TCTTAAAGCAATTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                    <  1:154182/140‑1 (MQ=255)

                                                                                      ataCTAGAGTAGCCTATCTCTCAGGCAAAAATT‑‑GAAGAAGAGAGGAGGTGAATTAA‑T‑T‑TGAA‑AAGCATTTTCAGTTGCATCAGATCACCAGTCATCAGTGGCTGTGGTGTTGCGTATGATACTTATTGGCTTAGTAAACg                                                    <  2:657606/140‑1 (MQ=255)

                                                                                         atagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGACGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaa                                                 >  2:433154/2‑140 (MQ=255)

                                                                                          tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaac                                                >  2:353430/1‑139 (MQ=255)

                                                                                             agTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                             >  1:644172/1‑136 (MQ=255)

                                                                                             agTCGCCTACCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTTATCAGTGGCTGTGTCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                              <  1:173036/140‑1 (MQ=255)

                                                                                              gTCGCCTATCTCTCAGACGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGATCACATCATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGTATTATCCTTATTGGCTTAGTAAACGCaaaaaat                                            >  2:188322/1‑135 (MQ=255)

                                                                                              gTCGCCTATATCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                            >  2:251746/1‑135 (MQ=255)

                                                                                                 gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGTCTTTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACTCAATACaaaa                                          >  1:467819/1‑140 (MQ=255)

                                                                                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATAATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataat                                         >  1:186500/1‑132 (MQ=255)

                                                                                                  cctctctttacggCGTAAATTGAGACCCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGTT‑GGTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGCAATACAAAAt                                         <  2:179563/129‑1 (MQ=255)

                                                                                                  ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                        >  1:53523/1‑131 (MQ=255)

                                                                                                  ccTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                        >  1:30862/1‑131 (MQ=255)

                                                                                                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatta                                       >  2:100951/1‑130 (MQ=255)

                                                                                                   cTAGCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCCGTCGATGTGCTGTTGCGTATTATACGTATTGGCTTAGTAAACGCaacaaataatta                                       <  2:304254/140‑11 (MQ=255)

                                                                                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTCC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGATGTGCTGTTGCGTATTATACTTATTAGCTTAGTAAACACaacaagtaattta                                      >  1:649359/1‑129 (MQ=255)

                                                                                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGATGTTGCGTATAATACTTATTGGCTTGGTAAACGCaacaaataattaa                                      >  1:451416/1‑129 (MQ=255)

                                                                                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                      >  1:307337/1‑129 (MQ=255)

                                                                                                     aTCTCTCATGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                     <  2:214678/140‑13 (MQ=255)

                                                                                                     aTCTCTCAGTCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑C‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                     <  2:532574/140‑13 (MQ=255)

                                                                                                       ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                   <  2:660952/140‑15 (MQ=255)

                                                                                                       ctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaatag                                   >  2:607649/1‑126 (MQ=255)

                                                                                                                 tccaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGcc                          <  2:654316/137‑1 (MQ=255)

                                                                                                                    aattgacgcctgagagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                       <  1:45457/129‑1 (MQ=255)

                                                                                                                    aTTTAGGCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCa                       <  2:117695/140‑1 (MQ=255)

                                                                                                                            gCAGATACGAGTTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATCGGCTAAGTCGCCGCAATACAAAATAGGTGACATATAGCCGCACCAATaaa               <  1:657444/140‑1 (MQ=255)

                                                                                                                                agagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATc           >  1:442154/1‑140 (MQ=255)

                                                                                                                                  gggaggGGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCCAAAAAAATcac         >  2:541456/2‑138 (MQ=255)

                                                                                                                                  aggaggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCCATACAAAATAGGTGACATATAGCCGCAAACATAAAAATccc         >  2:228734/1‑140 (MQ=255)

                                                                                                                                   ggaggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc        <  1:420496/140‑1 (MQ=255)

                                                                                                                                    gaggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt       <  2:375683/140‑1 (MQ=255)

                                                                                                                                      ggTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCCTTGCGTTTTTTTCTTATTGGCTAAGTAGAGGCAATACAAAATAGGTGACATATAGCCGCACCACTAAAAATCCCCCCa     >  2:293383/1‑140 (MQ=255)

                                                                                                                                       gtgtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc    >  1:556422/1‑140 (MQ=255)

                                                                                                                                        tgtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt   <  1:208517/140‑1 (MQ=255)

                                                                                                                                         gtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa  >  1:634105/1‑140 (MQ=255)

                                                                                                                                         gtA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCAATTGCGTTTTTTTCTTATTGGCTAAGTAGAAGAAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCACCTCACTa  >  1:433186/1‑140 (MQ=255)

                                                                                                                                         gtA‑TAA‑G‑G‑TGAT‑GCTTATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGAAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa  >  1:354392/1‑140 (MQ=255)

TAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTA  >  CP000730/1929723‑1929999

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