BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,913	C→G	25.9%	S213S (TCG→TCC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,913	0	C	G	25.9%	193.8 / 101.5	143	S213S (TCG→TCC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (53/52); new base G (20/17); total (74/69)
Fisher's exact test for biased strand distribution p-value = 8.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 145 total aligned reads displayed.

TTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35774‑36053

ttGAAAAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGTAATGGct                                                                                                                                              <  2:289270/141‑2 (MQ=255)

     cAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAt                                                                                                                                         <  1:58720/141‑1 (MQ=255)

     cAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAt                                                                                                                                         <  1:323583/141‑1 (MQ=255)

         agCTAAAATAAAATTGTGATCTAATAGAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGACACTCCGt                                                                                                                                     >  2:306491/1‑141 (MQ=255)

              aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                                                                <  1:18041/141‑1 (MQ=255)

              aaataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATc                                                                                                                                >  2:1121/1‑141 (MQ=255)

               aataaaatTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACAACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGAGGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCt                                                                                                                               >  2:76765/1‑141 (MQ=255)

                acaaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                                                                            >  1:133051/15‑141 (MQ=11)

                acaaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                                                                                                                            >  1:138666/15‑141 (MQ=11)

                   aaatTGTGATCTAATAAGAATAGAAGGTGCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCCTGTCAGTGTTCCTTTAACTTGCTAGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTGAag                                                                                                                           <  2:113847/141‑1 (MQ=14)

                     atagtattttaataaAAAT‑GATGGTTCCGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                       <  1:247329/132‑1 (MQ=11)

                     agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                       <  1:56848/132‑1 (MQ=11)

                     agagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc                                                                                                                       <  1:247202/132‑1 (MQ=11)

                      gagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                                                                                                                      <  1:276739/133‑1 (MQ=11)

                           ttttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATGATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCGTGGCGAAAATCAGTAGATCTGAAGAGGCCTGCgg                                                                                                                 >  1:133254/4‑141 (MQ=11)

                              taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACGTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt                                                                                                                >  2:167419/1‑141 (MQ=255)

                                                                              accaatatctCATGGTAGTGTTCGTTGAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCTAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAGAGTATTCt                                                                <  2:333292/132‑1 (MQ=1)

                                                                                caagataccATGTTAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCttt                                                              <  1:249583/133‑1 (MQ=1)

                                                                                              aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACt                                               >  2:58900/1‑141 (MQ=11)

                                                                                              aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACt                                               >  2:2305/1‑141 (MQ=14)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTTGTGGTATGGTGAAAATCCGTATATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACtt                                              >  2:279205/1‑141 (MQ=11)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACtt                                              <  2:242828/141‑1 (MQ=255)

                                                                                               gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACtt                                              <  2:165168/141‑1 (MQ=12)

                                                                                               gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGCACAGACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACtt                                              <  1:175654/141‑1 (MQ=255)

                                                                                               gggtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACtt                                              <  1:15005/139‑1 (MQ=255)

                                                                                                tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACttt                                             >  1:115136/1‑141 (MQ=11)

                                                                                                tgtTCGTCTAACTTGCTAGCATGATGCTAATTTCGGGGCATGGGTAAACTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTGTTTGCCGTATTG‑ATACttt                                             >  2:235536/1‑141 (MQ=11)

                                                                                                 gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTg                                            >  1:218223/1‑141 (MQ=17)

                                                                                                 gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACATTTTAAAGTATTCTTTGCTGTATTG‑ATACtttt                                            >  1:277969/1‑140 (MQ=17)

                                                                                                   tCGTTTAACGTGCTAGCATGATGCTAATTGCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGAt                                          <  1:254666/141‑1 (MQ=255)

                                                                                                   tCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGAt                                          >  1:34949/1‑141 (MQ=32)

                                                                                                    cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGTTa                                         >  2:62272/4‑141 (MQ=12)

                                                                                                     gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATat                                        <  1:279851/139‑2 (MQ=255)

                                                                                                     gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATat                                        <  2:59971/139‑2 (MQ=11)

                                                                                                     gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTGACGTAGTCTTTGCTGTATTG‑GTACTTGGATAc                                        >  1:21620/3‑141 (MQ=11)

                                                                                                     gCTTGACTGGCTAACATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCGGGCGAGACGAGGGGTTCTTTTTAGATATGGCGTAAATACACTCCATACCTTATAAAGTATTCCTTGCTATCTTG‑ATAcgtggagac                                        >  2:18553/1‑133 (MQ=2)

                                                                                                      tttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATatc                                       >  1:108439/2‑139 (MQ=11)

                                                                                                        tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTTAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCtt                                     <  1:124226/141‑1 (MQ=11)

                                                                                                         aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTg                                    >  2:123543/1‑141 (MQ=32)

                                                                                                         aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACtt                                              >  2:168075/1‑131 (MQ=17)

                                                                                                         aaCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACtt                                              <  1:168075/131‑1 (MQ=17)

                                                                                                          acgtGCTAGCATGATGCTAATGTCGTGGCATGGCGAAAATCCGTAGATCTGAAGTGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGt                                   <  2:296264/138‑1 (MQ=255)

                                                                                                          aCTTGCAAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGt                                   >  2:184078/1‑141 (MQ=11)

                                                                                                           cTTGCTAGCATGATGCTAATTTCTTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTACTG‑ATACTTTGATACCTTGtc                                  <  1:210202/141‑1 (MQ=12)

                                                                                                            ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtct                                 <  1:196818/141‑1 (MQ=17)

                                                                                                            ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtct                                 <  2:326566/141‑1 (MQ=17)

                                                                                                             tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtctt                                <  2:257462/141‑1 (MQ=17)

                                                                                                             tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACATTTtctt                                >  2:264823/1‑141 (MQ=12)

                                                                                                              gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtcttt                               <  1:312004/141‑1 (MQ=32)

                                                                                                              gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtcttt                               >  1:199532/1‑141 (MQ=17)

                                                                                                               cTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGtctttc                              <  1:90009/141‑1 (MQ=255)

                                                                                                               cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtctttc                              <  1:7345/141‑1 (MQ=17)

                                                                                                               cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtctttc                              >  2:39557/1‑141 (MQ=17)

                                                                                                                tcgCATGGTGCTAATTTCGTGGGATGGGGAAAATCCGTAGATCTTAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGtctttct                             <  2:148056/139‑1 (MQ=11)

                                                                                                                tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGtctttct                             <  1:294507/141‑1 (MQ=255)

                                                                                                                tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGtctttct                             <  1:50387/141‑1 (MQ=17)

                                                                                                                  gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTa                           <  1:136055/141‑1 (MQ=17)

                                                                                                                  ccATGATGCTAATTTCGGGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTa                           <  1:18461/140‑1 (MQ=17)

                                                                                                                   cATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTAc                          >  2:55812/1‑141 (MQ=33)

                                                                                                                      ggtgCTAATGTCCTCGCATGGGGCACATCCGTAGATCTGAAGAGACCTACGGTTCCTTTTATATAGAGCGTAAATACATTCAATACCTTTTACAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTTTCTTTCTTACttt                       <  2:120777/139‑1 (MQ=11)

                                                                                                                      gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACttt                       <  2:265513/141‑1 (MQ=25)

                                                                                                                         gCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTATTTCTTACTTTAat                    >  2:160235/1‑141 (MQ=11)

                                                                                                                         gCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAat                    >  1:299381/1‑141 (MQ=21)

                                                                                                                           taagtgcGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAGTGTCTTTTAAAGTCTTCTTTGCTGTCTTT‑ATACTTAGATACTTTGTCTTTCTTACTGTCatat                  >  2:168161/7‑141 (MQ=11)

                                                                                                                           tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTAatat                  >  2:171515/1‑141 (MQ=255)

                                                                                                                           tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAatat                  >  2:52750/1‑141 (MQ=25)

                                                                                                                           tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAatat                  >  2:100829/1‑141 (MQ=21)

                                                                                                                            aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGCGACCTGCAGTTCTTTTTATATAGAGCGTTAATACCGCCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 <  1:159455/141‑1 (MQ=11)

                                                                                                                             aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                >  1:274263/1‑141 (MQ=21)

                                                                                                                              tttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGAc               <  2:38184/141‑1 (MQ=25)

                                                                                                                              tttCGTGACATGGCGAAAATCCGGAGATCCCAAGAGACCTGCGGATATTTTTATATAGAGCGTAAATACACCCAATACCTTTTACAGTAGTCTGAACTGAATGG‑ATACCTTACTACCTCGTCTTTCTCACttaaccatgcc               >  2:47485/1‑132 (MQ=2)

                                                                                                                               ttCGTGGCATGGGGAAAATCCGTAGATCTGAATAGACCTACGGTTCTTTTTATCTATAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTT‑ATAATATCAAATCTTGTCTTACTTCCTTTacactgccc              >  2:91279/1‑133 (MQ=11)

                                                                                                                               ttCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTAATATGACg              >  1:107022/1‑141 (MQ=255)

                                                                                                                               ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACg              >  2:240515/1‑141 (MQ=25)

                                                                                                                                tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg             >  1:36322/1‑141 (MQ=25)

                                                                                                                                tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg             >  1:302755/1‑141 (MQ=25)

                                                                                                                                tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg             <  2:299381/141‑1 (MQ=37)

                                                                                                                                 cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTTAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            <  1:171511/141‑1 (MQ=255)

                                                                                                                                 cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATCCCTTGTCTTTCTTACTTGAATATGACGGt            <  2:138357/141‑1 (MQ=255)

                                                                                                                                  gTGGCATTGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTg           <  2:260163/141‑1 (MQ=255)

                                                                                                                                  gTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTg           >  2:310215/1‑141 (MQ=14)

                                                                                                                                  gTGGCATGGCGAAAATCCGTAGATCTTAATAGACCTGCGGGTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATTCCTTGTCTTTCTTACTTTAATATGACGGTg           <  1:70347/141‑1 (MQ=11)

                                                                                                                                  gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATCGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgtgg           >  1:192783/1‑138 (MQ=25)

                                                                                                                                   tGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCAGTCTTG‑TTACTTTGTTATCTTGTCTTTCTTACTTTAATATGACGGTgt          >  2:234114/1‑140 (MQ=11)

                                                                                                                                   tGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATCACGGTGa          <  1:209436/141‑1 (MQ=255)

                                                                                                                                    ggCATGGGGAAAATCCGTAGATCTGAAGAGCCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGttag         >  2:203015/1‑138 (MQ=11)

                                                                                                                                    ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCCATCCATTTTAAAGTTGTCGTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTggt         >  2:118692/1‑139 (MQ=11)

                                                                                                                                      cATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt       >  2:224842/1‑141 (MQ=255)

                                                                                                                                      cATGGCGAAAACTCGTAGATCTGAAGAGACCCGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt       <  1:119867/141‑1 (MQ=255)

                                                                                                                                       aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCtt      >  2:281442/1‑141 (MQ=11)

                                                                                                                                       aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCtt      >  2:164842/1‑141 (MQ=255)

                                                                                                                                       aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGGTGTTTTGATATAGAGCATAAATACATTCAATACCTTTTAAAGTAGTCTTTGCCGTATTGCATA‑TTCGATGTCTTGTCTTTCTTACTTTAATATGTCGGGGCTCtt      >  2:103706/1‑141 (MQ=2)

                                                                                                                                       aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCtt      <  2:312732/141‑1 (MQ=255)

                                                                                                                                       aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCtt      <  2:173174/141‑1 (MQ=255)

                                                                                                                                       aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAAGACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTA‑ATAATTTGATAACTTCTCTATCTTACCTTAATAGGACGGTGATCtt      >  2:214219/1‑141 (MQ=11)

                                                                                                                                        tGGGGAAAATCCGTAAATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTG‑ATACTTTGATATCTTGTCTTTCTTACTTTATTATGACGGTGATCTTg     <  1:165294/141‑1 (MQ=14)

                                                                                                                                        tGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATCCCTTGTCTTTCTTACTAAAATACGaacatgctgttg     <  2:195056/141‑12 (MQ=11)

                                                                                                                                         ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGACCTTGc    >  2:134818/1‑141 (MQ=255)

                                                                                                                                          gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTACAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt   >  1:63394/1‑141 (MQ=255)

                                                                                                                                          gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTAATCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCt   >  1:224828/1‑141 (MQ=255)

TTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTG‑ATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35774‑36053

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