BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,970	T→C	62.5%	E194E (GAA→GAG)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,970	0	T	C	62.5%	475.1 / inf	572	E194E (GAA→GAG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (109/104); new base C (165/189); total (275/297)
Fisher's exact test for biased strand distribution p-value = 2.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

Only 100 of 578 total aligned reads displayed.

AAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGC  >  CP000730/35838‑36110

                             cAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTATAGTATTCTTTGCTGTATTGATAc                                                                                                           >  1:743823/1‑141 (MQ=2)

                               gtgtTCGTTTACCTTTCGAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                                                                                         <  1:352092/141‑1 (MQ=11)

                                tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                        <  1:170252/141‑1 (MQ=17)

                                tgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                        >  1:714598/1‑141 (MQ=17)

                                 gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCGTTTTATATAGAGCGTAAATACACTCAGTACCTGTTAACGTATTCTTGTCTGTATTGATACGTTg                                                                                                       >  1:190730/1‑141 (MQ=11)

                                 gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg                                                                                                       >  1:408400/1‑141 (MQ=17)

                                  gtCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGa                                                                                                      <  1:116320/140‑1 (MQ=11)

                                   tCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                                                                                     >  1:461646/1‑141 (MQ=17)

                                   tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATCGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGCGCGTAAATACACTCAATACCTTTTAAAGTATTCTTAGCTGTATTGATACTTTGAt                                                                                                     <  1:499328/141‑1 (MQ=255)

                                    cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATa                                                                                                    >  2:15175/4‑141 (MQ=255)

                                    cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                                                    >  1:19951/1‑141 (MQ=17)

                                       ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatct                                                                                                 <  2:426290/141‑4 (MQ=255)

                                       ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                                                                                 <  2:81161/141‑1 (MQ=17)

                                        tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGGATTCTTTGCCGTATTGATACTTTGATATCtt                                                                                                >  1:211393/1‑141 (MQ=11)

                                        tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCtt                                                                                                <  2:606278/141‑1 (MQ=32)

                                         aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTTATACCTTg                                                                                               >  1:281380/1‑141 (MQ=17)

                                           cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACAAACCATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                                                             <  1:148808/141‑1 (MQ=12)

                                              gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                                                          <  2:743629/141‑1 (MQ=17)

                                               cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                         <  2:13942/141‑1 (MQ=17)

                                               cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                         <  1:436112/141‑1 (MQ=17)

                                                tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTACATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                        >  2:579637/1‑141 (MQ=17)

                                                   cATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTCTTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑Ac                                                                                    >  2:374693/1‑141 (MQ=11)

                                                   cATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑Ac                                                                                    >  1:670484/1‑141 (MQ=255)

                                                   cATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑Ac                                                                                    >  2:30319/1‑141 (MQ=17)

                                                         gCTAATTTCGTGGCATGGGGAAAATCCGTAGATCCGAAGAGACCTACGGTTCTGTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATCTCTTGTCTTTCTT‑ACTTTAat                                                                              >  2:414963/1‑141 (MQ=11)

                                                           tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGCTATCTTGTCTTTCTT‑ACTTTAat‑at                                                                           >  2:513995/1‑141 (MQ=11)

                                                           tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAat‑at                                                                           >  1:140191/1‑141 (MQ=21)

                                                            aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTGTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATg                                                                          <  2:565983/141‑1 (MQ=255)

                                                                  gtcgCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTg                                                                    <  1:338560/138‑1 (MQ=14)

                                                                  gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTg                                                                    >  2:117079/1‑141 (MQ=11)

                                                                  gTGGCATGGCGAAAATGCGGAGATCTGAAGAGACCTGCAGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATACTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTg                                                                    >  1:429928/1‑141 (MQ=11)

                                                                  gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTg                                                                    >  2:7194/1‑141 (MQ=37)

                                                                  gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGACACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTg                                                                    <  1:598631/141‑1 (MQ=255)

                                                                  gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTGATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGaggttg                                                                    >  2:495778/1‑136 (MQ=255)

                                                                  gTGGCATGGCGAAAATCCGTAGATCCGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTg                                                                    >  1:100074/1‑141 (MQ=255)

                                                                       aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCtt                                                               <  2:442051/141‑1 (MQ=14)

                                                                        tGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTg                                                              <  1:521522/141‑1 (MQ=255)

                                                                         ggggAAAATCCGTAAATCTGAAGAGACCTACGGTTCTTTTTATATGTAGCGTCAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACAGTGATCTAGc                                                             >  2:692684/4‑141 (MQ=11)

                                                                         ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGGGATCTTGc                                                             >  1:658843/1‑141 (MQ=255)

                                                                           cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGGTCTTGCTc                                                           >  1:235813/1‑141 (MQ=255)

                                                                           cGAAAATCCGGACATCCGAAGAGACCTGCGGTTCTTTTTATATAGCGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATCCCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTc                                                           <  1:298916/141‑1 (MQ=14)

                                                                            gAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCa                                                          <  2:503782/141‑1 (MQ=38)

                                                                            gAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCa                                                          <  1:323355/141‑1 (MQ=255)

                                                                              aaaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAAt                                                        >  2:339669/1‑141 (MQ=38)

                                                                              aaaTCCGTAGATCTGAAGAGACCTACTGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑‑CATTATTCAGGACGGTGATCTTGCTCAAt                                                        >  2:535340/1‑141 (MQ=11)

                                                                                aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGa                                                      <  1:84032/141‑1 (MQ=255)

                                                                                 tCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAg                                                     <  2:288263/141‑1 (MQ=38)

                                                                                 tCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTTATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAg                                                     >  2:112696/1‑141 (MQ=255)

                                                                                   cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGt                                                   >  1:579133/1‑141 (MQ=255)

                                                                                   cGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGt                                                   >  2:646673/1‑141 (MQ=255)

                                                                                    gTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGtt                                                  <  1:761895/141‑1 (MQ=255)

                                                                                    gTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGtt                                                  >  1:222452/1‑141 (MQ=255)

                                                                                     tAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGttt                                                 >  1:164626/1‑140 (MQ=255)

                                                                                      aGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTAt                                                >  1:537073/1‑141 (MQ=37)

                                                                                       gATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAATGTATTCTTTGGCGTATTGCTACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCCCAATGAGGTTAtt                                               <  2:357552/141‑1 (MQ=14)

                                                                                       aaTCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTAtt                                               <  2:209776/140‑1 (MQ=37)

                                                                                        aTCTGAAGTGAGCTACGGTTCTTTTTATAGAAGGCGTAAATACTTTCAATACTGTTTAAGGTATTCTTTGCCGTATGGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTc                                              <  2:18908/141‑1 (MQ=11)

                                                                                          cTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                            <  1:373796/141‑1 (MQ=255)

                                                                                          cTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                            <  1:491903/141‑1 (MQ=255)

                                                                                           tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGAGCTTGCTCAATGAGGTTATTCAGa                                           >  2:469812/1‑141 (MQ=38)

                                                                                           tGAAGAGACCTACGGTTCTTTTTATAGAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGa                                           <  2:324596/141‑1 (MQ=255)

                                                                                            gAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTTCTCAATGAGGTTATTCAGat                                          >  1:688164/1‑141 (MQ=255)

                                                                                             aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGata                                         <  1:128739/141‑1 (MQ=38)

                                                                                             aaGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGata                                         >  2:708921/1‑141 (MQ=255)

                                                                                             aaGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTAT‑TTTCTTGTCTGGAAT‑ATAACGGTGATCTTGCTCAATGAGGTTATTCAGata                                         >  1:553944/1‑141 (MQ=11)

                                                                                              agagACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGatat                                        <  2:458020/141‑1 (MQ=255)

                                                                                               gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                       <  1:647034/141‑1 (MQ=255)

                                                                                               gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                       <  1:306045/141‑1 (MQ=38)

                                                                                                agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTATTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAttt                                      >  2:655729/1‑141 (MQ=38)

                                                                                                 gaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTc                                     >  1:91304/1‑141 (MQ=255)

                                                                                                 gaCCTACGGTTCTTTTTATATAAAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTTATTTTGCTCACTGCGGTTACTCAGATATTTc                                     >  2:27745/1‑141 (MQ=11)

                                                                                                  aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCg                                    <  2:677110/141‑1 (MQ=255)

                                                                                                   ccTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGa                                   <  1:474199/141‑1 (MQ=255)

                                                                                                   ccTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGa                                   <  2:626174/141‑1 (MQ=255)

                                                                                                     tGCGGTGCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATg                                 <  1:392136/141‑1 (MQ=255)

                                                                                                        ggTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTAc                              >  1:620653/1‑141 (MQ=38)

                                                                                                         gTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGCTGTACa                             >  1:252966/1‑141 (MQ=255)

                                                                                                              ttttATATAGAGCGTAAATACGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGCCGGTGCTCTTTCTCACTGCGGTTATTCAGATATTTCGGTGTACTATggc                        >  2:618764/1‑139 (MQ=11)

                                                                                                              ttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATggc                        <  1:703625/141‑3 (MQ=255)

                                                                                                               tttATATAGAGCGTAAATACACTCAATTCCTTTTAAAGTATTCTTTGCTGTCTTGATACGTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATggca                       <  1:333521/141‑4 (MQ=14)

                                                                                                                  atataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTc                    >  1:635767/1‑141 (MQ=255)

                                                                                                                    ataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAg                  <  2:263344/141‑1 (MQ=255)

                                                                                                                    ataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAg                  <  1:725606/141‑1 (MQ=255)

                                                                                                                         gCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttta             >  1:384516/1‑141 (MQ=255)

                                                                                                                            tAAATACATTCAATACCTTTTAAAGTAATCTATGCCGTATTGATACTATGACATATTGTCTTTCTT‑AATTTACT‑ATGACGGAGATATTGCTCAATGAGGGAATTCAGATATTCCGAAGTCAGATGACAATCAGttttaagt          >  2:252737/1‑141 (MQ=2)

                                                                                                                             aaaTACACTCAATACCTTTTACAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtt         <  1:77445/141‑1 (MQ=255)

                                                                                                                             aaaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtt         <  2:748206/141‑1 (MQ=255)

                                                                                                                              aaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttt        <  1:314708/141‑1 (MQ=255)

                                                                                                                              aaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttt        >  2:478947/1‑141 (MQ=255)

                                                                                                                              aaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATCTCAATTTATAATGACAGTCAGgtttaagttt        >  2:729571/1‑141 (MQ=14)

                                                                                                                              aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt        <  2:545638/141‑1 (MQ=255)

                                                                                                                              aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt        <  2:178057/141‑1 (MQ=255)

                                                                                                                               aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta       >  1:568601/1‑141 (MQ=255)

                                                                                                                                tACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATAACAGTCAGgtgtaagtttaa      >  2:235170/1‑141 (MQ=21)

                                                                                                                                tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa      >  2:396373/1‑141 (MQ=255)

                                                                                                                                tACACTCAATACCTTTTAAAGTATTCGTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa      >  1:449469/1‑141 (MQ=255)

                                                                                                                                tACACTCAATACCTTTCACAGTTTTCTTTGCTGTATTGATACTTTGATAACTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCTATTTGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa      <  1:256176/141‑1 (MQ=14)

                                                                                                                                 acacTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaa     <  2:667390/141‑1 (MQ=255)

                                                                                                                                 aCATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTtaa     >  2:221433/1‑141 (MQ=255)

                                                                                                                                 aCATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTaaa     >  1:246896/1‑141 (MQ=255)

AAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTT‑ACTTTAAT‑ATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGC  >  CP000730/35838‑36110

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