BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,865	A→C	37.6%	intergenic (‑948/‑220)	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,865	0	A	C	37.6%	93.9 / 163.9	160	intergenic (‑948/‑220)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (44/56); new base C (36/24); total (80/80)
Fisher's exact test for biased strand distribution p-value = 7.20e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

Only 100 of 163 total aligned reads displayed.

GTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929726‑1930004

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                                                                                                                                   <  1:455596/140‑1 (MQ=255)

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

            cAATCTATCAGCTAGCATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTAt                                                                                                                                 >  1:495638/1‑140 (MQ=255)

            cAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTAt                                                                                                                                 <  1:416060/140‑1 (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:398272/140‑1 (MQ=255)

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

                       cTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTcaca                                                                                                                      <  1:195487/140‑1 (MQ=255)

                           aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCAcatcat                                                                                                                  <  1:523940/140‑1 (MQ=255)

                               atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCa                                                                                                              <  1:599815/140‑1 (MQ=255)

                               atatAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCATATTTTCGTTCACATCATAGcc                                                                                                              >  2:601680/1‑139 (MQ=255)

                                   aaCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAg                                                                                                          <  2:232125/140‑1 (MQ=255)

                                          cAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAg                                                                                                   >  1:642826/1‑140 (MQ=255)

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

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

                                            agcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                                                                                                <  2:186851/126‑1 (MQ=255)

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

                                               ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                             <  2:564867/129‑1 (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)

                                                 aagtagatatggttGTATATTCACAGATGGTATAATAACAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                                           <  2:392132/128‑1 (MQ=255)

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

                                                  agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                          <  1:512835/132‑1 (MQ=255)

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

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

                                                   gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                         >  1:97691/8‑139 (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)

                                                      gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtgatgtt                                                                                      >  2:484475/5‑135 (MQ=255)

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

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

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

                                                       atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                                                                                     >  1:236127/4‑135 (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)

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

                                                               gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat                                                                             <  1:652921/140‑3 (MQ=255)

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

                                                                  tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                          >  2:424596/9‑135 (MQ=255)

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

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

                                                                     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                       >  2:262819/6‑140 (MQ=255)

                                                                     tcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAt                                                                       >  2:223806/6‑140 (MQ=255)

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

                                                                      cacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTAtt                                                                      >  2:250494/5‑140 (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)

                                                                        cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAAAGGAGGTGCA‑TTA‑C‑TTTGCT‑AATCATTTTCGGTAACAACATAGGACCAGTCATCAGTGGCTGTGCTGTTTCGTATTATACTTATTgg                                                                    >  1:337951/3‑140 (MQ=255)

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

                                                                          gatgatATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCTCCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                   <  2:113652/140‑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‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCttag                                                               <  1:426444/138‑4 (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)

                                                                               tataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                             >  2:558378/1‑140 (MQ=255)

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

                                                                                ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCTGTGCTGTTGCGTATTATACTTCTTGGcggagaat                                                            >  2:135965/1‑133 (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)

                                                                                          agTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTGATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCATTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                   >  2:371126/1‑140 (MQ=255)

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

                                                                                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                 >  1:631400/1‑135 (MQ=255)

                                                                                           gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaat                                                 >  1:438656/1‑135 (MQ=255)

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

                                                                                            tCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaa                                                 >  1:284975/1‑140 (MQ=255)

                                                                                            tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑TCTCAATTGCCGCCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                <  1:484475/140‑7 (MQ=255)

                                                                                            tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                                                >  2:540239/1‑134 (MQ=255)

                                                                                              gCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACaaaa                                               <  2:177297/140‑1 (MQ=255)

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

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

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

                                                                                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                           >  1:418785/1‑129 (MQ=255)

                                                                                                 tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑TGCT‑AATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                           >  1:307337/1‑129 (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)

                                                                                                         aGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                    >  2:53056/1‑140 (MQ=255)

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

                                                                                                                aaTTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGc                             >  2:298503/1‑140 (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)

                                                                                                                   ttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTTACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                          >  2:355254/1‑140 (MQ=255)

                                                                                                                           agagagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaaa                  <  1:546191/140‑1 (MQ=255)

                                                                                                                            gagagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGAACCAATAAAAAt                 >  1:373733/1‑140 (MQ=255)

                                                                                                                             agagGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATc                <  1:581282/140‑1 (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‑GCTTATTCTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGAAATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTa       >  1:354392/1‑140 (MQ=255)

GTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑TGAT‑GCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCA  >  CP000730/1929726‑1930004

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