BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,929,797	G→A	49.0%	intergenic (‑880/‑288)	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,797	0	G	A	49.0%	15.8 / 138.7	113	intergenic (‑880/‑288)	USA300HOU_1794/USA300HOU_1796	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (27/29); new base A (31/24); total (58/55)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CTTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT  >  CP000730/1929657‑1929936

ctTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTg                                                                                                                                                   >  2:748055/1‑141 (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)

             aaTCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATataat                                                                                                                                      >  1:932464/1‑141 (MQ=255)

                cGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAAtagt                                                                                                                                   <  1:194269/141‑1 (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)

                          cTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGGAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑Ct                                                                                                                        <  1:754644/141‑1 (MQ=255)

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

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

                                       tatCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGt                                                                                                           >  2:1035769/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                                                                                       >  2:625737/1‑140 (MQ=255)

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

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

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

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

                                                                     gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑g                                                                          >  2:285229/1‑141 (MQ=255)

                                                                     gtgtAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCACAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGGGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑g                                                                          >  2:603150/1‑141 (MQ=255)

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

                                                                                tCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTAt                                                             <  2:138621/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                                             >  2:713111/1‑141 (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)

                                                                                                     tataACTTATCACAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCAcc                                        <  1:477129/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                            <  1:513320/141‑1 (MQ=255)

                                                                                                                 gAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTgg                            <  2:642341/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:690224/128‑1 (MQ=255)

                                                                                                                   tcggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCt                          >  1:779770/14‑141 (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                        <  2:289819/130‑1 (MQ=255)

                                                                                                                     ggaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgt                        <  2:692563/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                       >  2:116323/11‑141 (MQ=255)

                                                                                                                      gaagtagatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCtgtg                       <  1:660359/131‑1 (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:328843/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                 >  1:649318/5‑136 (MQ=255)

                                                                                                                            gatatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttg                 >  2:713399/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)

                                                                                                                             atatgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTA‑C‑TTTGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGctgttgc                >  2:889125/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)

                                                                                                                               atgCTTGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCCCAGGCGGCAATT‑‑GACGAAGAGAGGAGGTG‑CA‑TTACT‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGt              >  2:397605/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)

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

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

                                                                                                                                    tGTATATTCACAGATGGTATAATAATAGAGTCGC‑CTATCTCTCAGGCGTCAATT‑‑GACGAAGAGAGGAGGTG‑CATTAC‑T‑T‑TGCTAATCATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCTGTTGCGtattat         <  2:446697/141‑6 (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‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTc       >  2:1101970/1‑141 (MQ=255)

CTTGTCATAATAAAATCGCCTACTTTCTTAAATTAACAATATCTATTCTCATAGAATTTGTCCAATTAAGTGTAGACGATTCAATCTATCAGCTAGAATCATATAACTTATCAGAAGC‑AAGTGACTGTGCTTGTATATTTGCCGATGATATAATAGTAGAGTCGC‑CTATCTCTCAGGCGTCAATTTAGACGCAGAGAGGAGGTG‑TA‑TAA‑G‑G‑TGATGCTTATTTTCGTTCACATCATAGCACCAGTCATCAGTGGCTGTGCCATTGCGTTTTTTTCTTATT  >  CP000730/1929657‑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: `‑`