BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,843	Δ1 bp	58.8%	intergenic (‑926/‑242)	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,843	0	A	.	58.8%	70.7 / 255.7	131	intergenic (‑926/‑242)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (23/30); new base . (47/30); total (71/60)
Fisher's exact test for biased strand distribution p-value = 5.16e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 131 total aligned reads displayed.

AGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACC  >  CP000730/1929709‑1929980

aGAATTTGTCCAATTAAGTGTAGCCGATTCAATCTCTCAGCTAGAATCTTATGACTGGTCAGAAGCAATTGACTGGTCTTGTAGATTTGCGGATGGTATAATAGTAGAGGCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                                                          <  2:223840/141‑1 (MQ=255)

aGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCa                                                                                                                                          <  2:120252/141‑1 (MQ=255)

   atgtGTCCAATTAAGGGTAGAGGATTCAACCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagag                                                                                                                                       <  2:302552/138‑1 (MQ=255)

     ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                                                                     >  2:52323/1‑141 (MQ=255)

     ttGGCCAATGAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTTCTTGTATATGCGCCGAGGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                                                                     <  2:198903/141‑1 (MQ=255)

          cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgt                                                                                                                                <  2:359265/141‑1 (MQ=255)

          cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAAGTTAGACGCGGAGGGGAGgt                                                                                                                                >  1:518861/1‑141 (MQ=255)

               aaGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGta‑ta                                                                                                                          <  1:385609/141‑1 (MQ=255)

                 gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                       <  2:104716/141‑1 (MQ=255)

                 gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGGGTA‑TAA‑g                                                                                                                       >  1:33137/1‑141 (MQ=255)

                 gtggAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑g                                                                                                                       <  2:469274/141‑1 (MQ=255)

                              aaTCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTAttt                                                                                                       <  2:55256/141‑1 (MQ=255)

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

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

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

                                            aaTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATa                                                                                         <  1:40471/141‑1 (MQ=255)

                                             aTCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAg                                                                                        >  2:651616/1‑141 (MQ=255)

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

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

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

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

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

                                                            aGAAGCAAGTGACTGTGTTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTg                                                                         <  1:501325/141‑1 (MQ=255)

                                                            aGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTGGTTCACATCATAGCACCAGGAATCAGTg                                                                         >  1:255828/1‑141 (MQ=255)

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

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

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

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

                                                                  aaGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                                                                   >  1:291877/1‑141 (MQ=255)

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

                                                                   agtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGct                                                                 <  2:618990/133‑2 (MQ=255)

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

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

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

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

                                                                      gACTGTGCTTGTATATTTGCCGATCATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCCCATCATAGCACCAGTCAGCAGTGGCTGTGCCAt                                                               <  1:134979/141‑1 (MQ=255)

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

                                                                        atatgagtGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                                                            <  1:71702/134‑7 (MQ=255)

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

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

                                                                          ttgCTTGTATAGTTGCCGATGATATAATAGTAGCGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg                                                           <  2:101255/140‑1 (MQ=255)

                                                                          atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt                                                          <  2:516675/140‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                    tatcacagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCA‑TTACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTa                                                >  2:337992/8‑141 (MQ=255)

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

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

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

                                                                                         cagatggtATAATAATAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGc                                           <  2:591541/139‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                                                                  taataaTAGAGTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTaaac                                  <  1:220128/141‑4 (MQ=255)

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

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

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

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

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

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

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

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

                                                                                                         agagTCGCCTCTCTCTCCGGCGTCAATT‑‑GACGACGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCAGTGCTGTTGCGTCTTATACTTATTGGCTTAGTAAACGccacaa                           >  1:190893/1‑136 (MQ=255)

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

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

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

                                                                                                           agTCCCCTATCTCTCAGGCGTCAATTTATACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACa                          <  1:65601/141‑1 (MQ=255)

                                                                                                            gTCGCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaata                        >  2:430302/1‑135 (MQ=255)

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

                                                                                                               gCCTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataatt                     >  2:154913/1‑132 (MQ=255)

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

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

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

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

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

                                                                                                                  tATCTCTCAGGCGTAAATT‑‑GACGAAGAGAGGAGGTGCA‑TAACT‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCCCCAGTCATCAGGGGCTGTGCTGTTGCGTATTATACTTATTGGCTGAGTAAACGCaacaaataattaat                  <  1:603319/141‑13 (MQ=255)

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

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

                                                                                                                   aTCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTAGTAAACGCaacaaataattaata                 >  1:377099/1‑128 (MQ=255)

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

                                                                                                                    tctctcAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCATTAC‑T‑T‑T‑GCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGTATTATACTTATTGGCTTTGTAAACGCaacaaataattaatag                <  1:426032/141‑15 (MQ=255)

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

                                                                                                                       ctcAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACa              >  2:389607/1‑141 (MQ=255)

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

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

                                                                                                                             cggcAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGc        <  2:463808/138‑1 (MQ=255)

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

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

AGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGCAAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGCCTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTGTA‑TAA‑G‑G‑T‑GATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATTGGCTAAGTAGACGCAATACAAAATAGGTGACATATAGCCGCACC  >  CP000730/1929709‑1929980

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