BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,796	T→C	48.9%	intergenic (‑879/‑289)	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,796	0	T	C	48.9%	‑3.7 / 148.2	113	intergenic (‑879/‑289)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (27/30); new base C (31/23); total (58/55)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 116 total aligned reads displayed.

TCTTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT  >  CP000730/1929656‑1929936

tctTGTCTGAATAAAATCGCCTTCTGTCTTAAATTAACAATATCAAGTCTCACATAATTTGGCTAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTCTTCAGAAGC‑AAGTGACTGTGCTTGTATAttt                                                                                                                                                    <  2:146098/141‑1 (MQ=25)

tctTGTCATAATAAAATCGCCTGCTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATAttt                                                                                                                                                    <  1:765943/141‑1 (MQ=255)

   tGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTACTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCCAGAATCATATACCTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGcc                                                                                                                                                 <  2:123669/141‑1 (MQ=255)

   tGTCACAATAAAATCGATGACTTTCTGAAATTAACAATATCTATTCTCATAGAATTTGTCCAAGTAAGTGTAGACGATTCAATTTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGcc                                                                                                                                                 <  2:198251/141‑1 (MQ=25)

      cATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGCAGCTGAGT‑TCTTTTCTTGTATATTTGCCgat                                                                                                                                              <  2:1068457/141‑1 (MQ=21)

         aataaAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCgatgat                                                                                                                                           <  1:818537/141‑1 (MQ=255)

              aaTCGCCTACTTTCTTAAAGTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                      >  1:719774/1‑141 (MQ=255)

                 cGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                                                   <  2:1019512/141‑1 (MQ=255)

                   ccTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagtag                                                                                                                                 >  2:245631/1‑141 (MQ=255)

                      aCTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGGGCTTGTATATTTGCCGAGGATAGAATAGTAGAGt                                                                                                                              >  1:458261/1‑141 (MQ=255)

                        tttCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCg                                                                                                                            >  1:504767/1‑141 (MQ=255)

                            ttAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTa                                                                                                                       <  1:362728/141‑1 (MQ=255)

                             tAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTAt                                                                                                                      >  1:542137/1‑141 (MQ=255)

                                aTTAACAAGATCTATTCTCATAGAATGTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTAtctc                                                                                                                   <  2:353183/141‑1 (MQ=255)

                                    aCAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTGTCAg                                                                                                               >  2:739895/1‑141 (MQ=255)

                                        tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGTTTTTGTTTTTATTTTTGCCGATGATATAATAGTAGAGTCGC‑CTATGTCTCAGGCGt                                                                                                           >  1:857658/1‑141 (MQ=38)

                                        tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGt                                                                                                           >  2:1035769/1‑141 (MQ=255)

                                         atCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTc                                                                                                          >  2:874452/1‑141 (MQ=255)

                                          tCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCa                                                                                                         >  1:371047/1‑141 (MQ=255)

                                             aTTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGACAAtt                                                                                                      >  2:390148/1‑141 (MQ=255)

                                                    tAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGc                                                                                               <  2:1102013/141‑1 (MQ=255)

                                                          ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCagagag                                                                                         <  1:245631/141‑1 (MQ=255)

                                                          ttGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTAAAAGATTAGATAATCGTTTA‑ACGCGCGAGCTCTCAGGCGTCAATTTAGACGCagagag                                                                                         <  1:1093049/141‑1 (MQ=255)

                                                           tGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAgagg                                                                                        >  1:997173/1‑141 (MQ=255)

                                                            gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAgaggg                                                                                       >  1:347522/1‑140 (MQ=255)

                                                            gTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAgaggg                                                                                       >  2:625737/1‑140 (MQ=255)

                                                              ccAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGaggagg                                                                                     <  1:335820/141‑1 (MQ=255)

                                                               cAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGgt                                                                                    >  1:201777/1‑141 (MQ=255)

                                                                     aGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑Taa                                                                            <  2:719774/141‑1 (MQ=255)

                                                                                 tCAATCTATCAGCTAGAATCATATAACTTATCAGCCGC‑CAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                             <  1:748055/141‑1 (MQ=255)

                                                                                               aGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCAcatca                                               <  1:611307/141‑1 (MQ=255)

                                                                                                 aaTCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATa                                             >  1:176627/1‑141 (MQ=255)

                                                                                                     ctataACTTGTCAGGAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAc                                         <  1:649469/140‑1 (MQ=255)

                                                                                                     atatAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAc                                         <  1:616798/141‑1 (MQ=255)

                                                                                                       ataACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCTGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                       <  2:562550/141‑1 (MQ=255)

                                                                                                       ataACTTATCAGAAGC‑AAGGGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCa                                       <  2:385224/141‑1 (MQ=255)

                                                                                                                 aGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATCGTCGATGCGC‑CTATCTCTCAGGCGTCTATTTAGACGCAGAGAGGAGGGG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTGCACATCATAGCACAAGTCATGAGTg                             >  1:213763/1‑141 (MQ=255)

                                                                                                                 aGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATAAGTg                             >  2:342894/1‑141 (MQ=255)

                                                                                                                  gAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                            <  2:642341/141‑1 (MQ=255)

                                                                                                                  gAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                            <  1:513320/141‑1 (MQ=255)

                                                                                                                   atcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGc                           >  1:298789/15‑141 (MQ=255)

                                                                                                                    tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:779770/14‑141 (MQ=255)

                                                                                                                    tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          <  1:690224/128‑1 (MQ=255)

                                                                                                                    tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:694400/14‑141 (MQ=255)

                                                                                                                    tcggaagtagatatgCCTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:654921/14‑141 (MQ=255)

                                                                                                                     cggaagtagatatgCTTGTATATTCACAGATGGTATGATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                         <  2:1075602/129‑1 (MQ=255)

                                                                                                                     cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGTGCCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                         >  1:483337/13‑141 (MQ=255)

                                                                                                                     cggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGCGAGGCGCTC‑CAGTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtg                         >  1:753601/13‑141 (MQ=255)

                                                                                                                      tgaaggtgatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:936672/130‑1 (MQ=255)

                                                                                                                      ggaagtagttatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:920351/130‑1 (MQ=255)

                                                                                                                      ggaagtagatatgCTTGTATATTTGCAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:825765/130‑1 (MQ=255)

                                                                                                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCCGTGGCtgg                        >  1:755295/12‑140 (MQ=255)

                                                                                                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        >  1:49264/12‑141 (MQ=255)

                                                                                                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:692563/130‑1 (MQ=255)

                                                                                                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:289819/130‑1 (MQ=255)

                                                                                                                      ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        >  1:755035/12‑141 (MQ=255)

                                                                                                                      c‑AAGTGATTGTGCTTGTATATTTGCAGATGATAGAATAGTAGAGTCGC‑CTACCTCTCAGGCGGCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:451705/141‑1 (MQ=255)

                                                                                                                       gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       <  1:932688/131‑1 (MQ=255)

                                                                                                                       gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       >  1:341501/11‑141 (MQ=255)

                                                                                                                       gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       >  2:116323/11‑141 (MQ=255)

                                                                                                                       gaagtagatacgCTTGTAGATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTGCCG‑TAC‑G‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       <  2:721744/130‑1 (MQ=255)

                                                                                                                        aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                      <  1:1076765/132‑1 (MQ=255)

                                                                                                                        aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGc                      <  1:175702/132‑1 (MQ=255)

                                                                                                                        aagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGGGc                      >  1:561739/10‑141 (MQ=255)

                                                                                                                          gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                    <  1:990400/134‑3 (MQ=255)

                                                                                                                          gtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGGCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGGTCACATCATAGCACCAGTCATCAGTGGCTGTGctg                    >  1:173457/8‑139 (MQ=255)

                                                                                                                           tagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GGCGAAGAGAGGAGGGG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGctgt                   >  1:563036/7‑138 (MQ=255)

                                                                                                                            agatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCCCCAGTCACCAGTGGCGGTGcttgt                  >  2:168351/6‑137 (MQ=255)

                                                                                                                             ggtttgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 <  2:615632/137‑6 (MQ=255)

                                                                                                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  2:713399/5‑136 (MQ=255)

                                                                                                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  1:649318/5‑136 (MQ=255)

                                                                                                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  2:687062/5‑136 (MQ=255)

                                                                                                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGctgttg                 >  1:96340/5‑136 (MQ=255)

                                                                                                                             gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GAAGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTAGTTCACATCATATCACCAGTAATCAGTGGCCGTGctgtgg                 >  2:989420/5‑136 (MQ=255)

                                                                                                                              atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                >  2:811502/4‑135 (MQ=255)

                                                                                                                               tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGCGGCC‑AACTGC‑G‑T‑TTCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCg               >  1:255225/3‑141 (MQ=255)

                                                                                                                               tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTAGCtgtgatgttacg               >  2:436172/3‑133 (MQ=255)

                                                                                                                               tatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATAATTTTCGTCCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgcg               >  2:114621/3‑134 (MQ=255)

                                                                                                                                ttgCTTGGATATTTTCCTATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATCGCACCAGTCATCAGTGGCTGTGCCATTGCg               <  2:892052/140‑1 (MQ=255)

                                                                                                                                gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGT‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg               <  1:718776/141‑1 (MQ=255)

                                                                                                                                gtgCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCg               <  2:733135/141‑1 (MQ=255)

                                                                                                                                atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt              <  2:745872/140‑1 (MQ=255)

                                                                                                                                atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAAGCCTTTTCGTTCACATCATAGCACCAGTCATCAGTGGCGGTGCTGTTGCGt              >  1:478695/2‑141 (MQ=255)

                                                                                                                                  gcgtGTATATTAACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat            <  1:245101/138‑3 (MQ=255)

                                                                                                                                  gCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtt             >  2:951174/1‑141 (MQ=255)

                                                                                                                                  gCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat            <  1:112256/141‑3 (MQ=255)

                                                                                                                                  gCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtat            <  1:176831/141‑3 (MQ=255)

                                                                                                                                   cTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt           >  2:159904/1‑141 (MQ=255)

                                                                                                                                   cTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatt           >  2:712100/1‑138 (MQ=255)

                                                                                                                                    ttGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta          <  2:455323/141‑5 (MQ=255)

                                                                                                                                    ttGTATATTCACAGACGGTATAATAATAGAGTCGC‑CTATCTCTCTGGCGTCAATGT‑‑ACGATGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtatta          <  2:759559/141‑5 (MQ=255)

                                                                                                                                     tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         <  2:446697/141‑6 (MQ=255)

                                                                                                                                     tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         >  2:84588/1‑139 (MQ=255)

                                                                                                                                     tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCAGCAGTGGCTGTGCTGTTGCGtattat         <  2:395724/141‑6 (MQ=255)

                                                                                                                                     tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑CTT‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         <  2:761835/141‑6 (MQ=255)

                                                                                                                                      gTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttt         <  2:932464/141‑1 (MQ=255)

                                                                                                                                       tataTTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCATATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGttttttt        <  2:504767/141‑1 (MQ=255)

                                                                                                                                        atatTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGtttttttt       >  2:145239/1‑140 (MQ=255)

                                                                                                                                        atatTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc       >  2:1101970/1‑141 (MQ=255)

                                                                                                                                             tGCCGATGATATAATAGTAGAGTCGC‑ATATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGGCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTAtt  >  2:687741/1‑141 (MQ=255)

TCTTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT  >  CP000730/1929656‑1929936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`