BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,873	C→A	55.5%	intergenic (‑956/‑212)	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,873	0	C	A	55.5%	33.4 / 155.0	134	intergenic (‑956/‑212)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base C (29/30); new base A (45/28); total (75/59)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 134 total aligned reads displayed.

AATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGA  >  CP000730/1929739‑1930012

 aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGGGGAGGTGTA‑TAA‑G‑G‑G‑GGTGC‑TTGtttt                                                                                                                                      >  1:67697/1‑141 (MQ=255)

 aTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTAtttt                                                                                                                                      >  1:617662/1‑141 (MQ=255)

             gAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCAcatcat                                                                                                                          >  2:589440/1‑141 (MQ=255)

              aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGCTAAAAGATTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATa                                                                                                                         <  1:625249/141‑1 (MQ=255)

              aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATa                                                                                                                         <  1:40282/141‑1 (MQ=255)

               aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAg                                                                                                                        >  1:104477/1‑141 (MQ=255)

                    ataACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCAATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCa                                                                                                                   >  1:456894/1‑141 (MQ=255)

                       aCTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCTGtc                                                                                                                >  2:546668/1‑141 (MQ=255)

                             cAGAAGCAAGTGACTGTGCTTGTTCTTGAACCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGt                                                                                                          >  1:661663/1‑141 (MQ=255)

                             cAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGt                                                                                                          <  2:447099/141‑1 (MQ=255)

                             cAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGt                                                                                                          >  2:655242/1‑141 (MQ=255)

                               atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                                                                                                       >  1:621917/15‑141 (MQ=255)

                                  tgaagtcgatatccTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                    <  1:509801/128‑1 (MQ=255)

                                  ggaagtaggtatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                    <  1:82708/130‑1 (MQ=255)

                                  ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                                                                                                    >  1:168329/12‑141 (MQ=255)

                                  ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTAATAAGTGGCtgt                                                                                                    >  1:368623/12‑141 (MQ=255)

                                  gCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACGTAATAGCAGCAGTCATCAGTGGCtg                                                                                                     >  2:189474/1‑141 (MQ=255)

                                  gCAAGTGACTGTGCTTGTATATTTACCGATGATATAATAGTAGAGTCGCCTATCTCGCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                                                                                                     >  2:490596/1‑141 (MQ=255)

                                    aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                                                                                                  <  1:442291/132‑1 (MQ=255)

                                     agtagatattcTTCTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                 <  2:500018/131‑2 (MQ=255)

                                     agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTGTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                                                 <  1:430107/133‑2 (MQ=255)

                                      gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                <  1:543834/134‑3 (MQ=255)

                                      gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                <  1:170476/134‑3 (MQ=255)

                                      gcagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                                                                                                >  2:578301/8‑139 (MQ=255)

                                        agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgct                                                                                              >  1:476726/6‑137 (MQ=255)

                                          atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                                                            <  1:603429/138‑7 (MQ=255)

                                           tatgattGTATATTCACAGATGGTATAATAATAGAATCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑T‑GGAGCATCATTGTCGGTCACAGCATAGCCCCAGTCATCAGTGGCCGTGCCGGTGCg                                                                                           >  2:41985/6‑141 (MQ=255)

                                           tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCGGTGGCTGTGCTGTTGCg                                                                                           >  1:117086/3‑141 (MQ=255)

                                            gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                                                           >  1:416956/1‑141 (MQ=255)

                                               cTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt                                                                                       <  1:557506/141‑4 (MQ=255)

                                                ttGTATATTCACAGATGGGATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCCGTCATCAGTGGCTGTGCTGTTGCGtatta                                                                                      <  1:594485/141‑2 (MQ=255)

                                                  gTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATAATAGCACCAGTCATCAGTGGCCGTGCTGTTGCGtattata                                                                                    >  2:374439/1‑138 (MQ=255)

                                                   tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtgttatac                                                                                   >  1:187644/11‑137 (MQ=255)

                                                   tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatat                                                                                   >  2:605751/11‑137 (MQ=255)

                                                   tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑T‑TGGCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                   >  1:451585/11‑137 (MQ=255)

                                                   tatattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGACGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TAATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCTGTTGCGtattatac                                                                                   >  2:605744/11‑137 (MQ=255)

                                                   tataTTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATTC‑TTATCTTCGTTCGCATCATAGCACCAGTCATCAGTGCCTGTGCCATTGCGttttttt                                                                                    <  1:407313/141‑1 (MQ=255)

                                                    atattcacagatggtATAATAATAGAGTCGCGTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGAGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTACACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattatact                                                                                  >  1:491760/10‑136 (MQ=255)

                                                     tattcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGTATTATACtt                                                                                 >  1:303874/9‑141 (MQ=255)

                                                     aatTTGCCGATGATATAATAGGAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCt                                                                                  <  2:262162/140‑1 (MQ=255)

                                                      attcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                                                >  2:121783/8‑141 (MQ=255)

                                                      atTTGCCGATGATATAATAGTAGAGTCGACTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCtt                                                                                 <  1:463886/141‑1 (MQ=255)

                                                          acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGGTGCATTAC‑T‑T‑T‑GCTAG‑TCCTTGTCGCTCACATCAGAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                            >  2:640275/4‑141 (MQ=255)

                                                          acagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTgg                                                                            <  1:252100/138‑1 (MQ=255)

                                                           cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                                                           >  1:9688/3‑141 (MQ=255)

                                                            cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTgtt                                                                           >  2:606001/1‑139 (MQ=255)

                                                            cGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGc                                                                           <  2:40546/141‑1 (MQ=255)

                                                               tggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAg                                                                       >  1:183012/4‑141 (MQ=255)

                                                                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTCCGTATTGTACTTCTTGGCTTAgg                                                                      >  2:2184/3‑140 (MQ=255)

                                                                ggtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                      >  2:621136/3‑141 (MQ=255)

                                                                ggtacaatgttagAGTCGCCTATCTCCTAGGCGGCACTT‑‑GCCGAAGAGAGGAGGTGCATTCC‑T‑T‑T‑GCTAG‑TCATTTCCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGt                                                                      <  2:112381/131‑1 (MQ=255)

                                                                 gtataATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGCTCACCTCAAAACACCTGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTa                                                                     >  1:126114/2‑141 (MQ=255)

                                                                 atatAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGt                                                                      >  2:647888/1‑141 (MQ=255)

                                                                   ctaataaGAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                   <  1:533308/140‑3 (MQ=255)

                                                                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                   >  1:328466/1‑139 (MQ=255)

                                                                   ataataataGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaa                                                                   >  1:145633/1‑139 (MQ=255)

                                                                    taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑GCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATAATTATTGGCttggtaaac                                                                  >  1:91844/1‑133 (MQ=255)

                                                                     aataatACAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACg                                                                 <  2:430551/141‑1 (MQ=255)

                                                                      ataatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGc                                                                >  2:336579/5‑141 (MQ=255)

                                                                      ataGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTacctcgacc                                                                 >  2:429641/1‑133 (MQ=255)

                                                                       taatagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TAATTTTCGTTCACATCATAGCAGCAGTCATCAAAGGATGTGCTGTTGCGTCTTATACTTATTGGCTTAGTAACCGCa                                                               >  2:281208/4‑141 (MQ=255)

                                                                       aaatagAGGCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCa                                                               <  2:581824/138‑1 (MQ=255)

                                                                          tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                            >  1:293777/1‑139 (MQ=255)

                                                                          tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                            <  2:299261/141‑3 (MQ=255)

                                                                          tagAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑CTT‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaaca                                                            >  1:2819/1‑139 (MQ=255)

                                                                           agagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaa                                                           >  1:40321/1‑138 (MQ=255)

                                                                            gagTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATAc                                                           <  1:459516/141‑1 (MQ=255)

                                                                            gagTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaa                                                          >  2:282973/1‑137 (MQ=255)

                                                                              gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAACCGCaacaaata                                                        >  1:213076/1‑135 (MQ=255)

                                                                               tCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataa                                                       >  2:497279/1‑134 (MQ=255)

                                                                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                     <  1:464365/141‑10 (MQ=255)

                                                                                 gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTA‑C‑TTT‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGGTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                                                     <  1:456075/141‑10 (MQ=255)

                                                                                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCGGTCGTCATCGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                   >  1:154206/1‑130 (MQ=255)

                                                                                   cTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGGGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaa                                                   <  2:45839/141‑12 (MQ=255)

                                                                                    tacccctcAGGCGTAAATTTAGACGCAGAGAGGAGGTGGA‑TAA‑G‑G‑T‑GAGGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAgg                                                   <  1:190978/136‑1 (MQ=255)

                                                                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGCGGAGTGGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATTAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                  <  2:372869/141‑13 (MQ=255)

                                                                                    tATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaat                                                  >  2:17659/1‑129 (MQ=255)

                                                                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAC‑TCATTTTCTTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                 <  1:349945/141‑14 (MQ=255)

                                                                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaatccatactc                                                 >  2:328789/1‑128 (MQ=255)

                                                                                     aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                                                 >  2:101483/1‑128 (MQ=255)

                                                                                      tctctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCCCATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGCCTAAGTAGCAGCCATACAACATAGGTg                                                 >  2:528254/1‑141 (MQ=255)

                                                                                      tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAA‑TCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattcatag                                                >  2:195512/1‑127 (MQ=255)

                                                                                         cccacgCGTCAATTTATCCGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACa                                              <  1:401522/136‑1 (MQ=255)

                                                                                           cAGGCGTCAATTTAGACGAAGAGAGGAGGTGTA‑GAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACata                                            <  2:33137/141‑1 (MQ=255)

                                                                                                     tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                                  <  1:425465/141‑1 (MQ=255)

                                                                                                     tttAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTATGTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCAcc                                  >  2:368015/1‑141 (MQ=255)

                                                                                                            gCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATaaaa                           >  2:659479/1‑141 (MQ=255)

                                                                                                                  aggaggTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGACATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATcccc                     >  1:628334/1‑141 (MQ=255)

                                                                                                                   ggaggTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCt                    <  2:543100/141‑1 (MQ=255)

                                                                                                                      ggTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCAc                 <  1:589440/141‑1 (MQ=255)

                                                                                                                       gtgtA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACt                >  1:231736/1‑141 (MQ=255)

                                                                                                                          ta‑taA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTAcc             >  1:281535/1‑141 (MQ=255)

                                                                                                                             taA‑G‑G‑T‑GATTC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc           <  1:24990/141‑1 (MQ=255)

                                                                                                                             taA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGc           <  1:494761/141‑1 (MQ=255)

                                                                                                                                 g‑g‑T‑GAGTC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCaaa        <  2:232352/141‑1 (MQ=255)

                                                                                                                                       gATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAACATAGGAGTCATAGAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg     <  1:293784/141‑1 (MQ=255)

                                                                                                                                       gATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGACATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg     <  1:437436/141‑1 (MQ=255)

                                                                                                                                       gATGC‑TCATTTTCGTTCACATCAGAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAg     <  1:2760/141‑1 (MQ=255)

                                                                                                                                          gC‑TTATTTTCGTTGCAATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa  <  1:2023/141‑1 (MQ=255)

                                                                                                                                          gC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGa  >  1:65744/1‑141 (MQ=255)

AATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGC‑TTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACCAATAAAAATCCCCTCACTACCGCAAATAGTGA  >  CP000730/1929739‑1930012

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