BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,913	C→G	29.3%	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	29.3%	639.3 / inf	570	S213S (TCG→TCC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (176/222); new base G (81/85); total (260/310)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 577 total aligned reads displayed.

GAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35776‑36053

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

    aaGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATc                                                                                                                                        >  1:395350/1‑141 (MQ=255)

     agagCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATcc                                                                                                                                       >  2:474122/1‑141 (MQ=255)

              acaaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGaa                                                                                                                           >  1:814929/15‑141 (MQ=11)

               caaaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAag                                                                                                                          <  2:612587/128‑1 (MQ=11)

                 aaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAagag                                                                                                                        >  1:886754/12‑141 (MQ=11)

                 aaagagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGAAC‑TGAagag                                                                                                                        >  1:241131/12‑141 (MQ=11)

                    gagtattttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGAcc                                                                                                                     >  1:881506/9‑141 (MQ=11)

                      gtgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGAcc                                                                                                                     >  2:722151/1‑141 (MQ=255)

                       tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCt                                                                                                                    <  1:298545/141‑1 (MQ=255)

                        attttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCg                                                                                                                 <  2:926098/137‑1 (MQ=11)

                           ttaataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGCAGATA‑TGAAGAGACCTGCGGtt                                                                                                              >  2:45472/2‑141 (MQ=11)

                              ataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCtt                                                                                                           >  1:471879/1‑141 (MQ=11)

                                aaaaaTAGAAGGTTCTGTTGCAAAGTCTAAAAATATAGCTAACCACTAATTAATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCtt                                                                                                           >  2:668240/1‑141 (MQ=255)

                                      aGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTTCCTGAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTtata                                                                                                     <  1:218511/141‑1 (MQ=12)

                                               tgctGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTa                                                                                            <  2:559605/138‑1 (MQ=1)

                                                                            acgaaTTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCt                                                               <  1:544719/138‑1 (MQ=1)

                                                                                          tctgtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATa                                                 <  2:593001/138‑1 (MQ=1)

                                                                                           cAGTGTTCGTTTAACGTGCTAGCATGGTGCTAATCTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTCTTATATAGAGCGTCAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAc                                                <  1:572894/141‑1 (MQ=11)

                                                                                            tgtgtatgtttAATTTGCTTGCATTATTCTTGTTTCGTGGCATGGGGAAAGTCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                               >  2:693551/10‑141 (MQ=11)

                                                                                            aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGGCGTATTGATACt                                               >  1:73670/1‑141 (MQ=11)

                                                                                            aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                               >  1:848576/1‑141 (MQ=255)

                                                                                            aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                               >  2:299885/1‑141 (MQ=255)

                                                                                             gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCAGGGGGACAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTGAAGTATTCTTTGCCGTATTGATACtt                                              >  1:366618/1‑141 (MQ=11)

                                                                                             gtgtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                              >  1:833100/1‑141 (MQ=17)

                                                                                               gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCGTGTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTATTTGCCGTATTGATACTTTg                                            >  2:544628/1‑141 (MQ=11)

                                                                                               gtTCGTTTAACTTGCTAGCATGATGCGAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TCAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAATTATTCTTTGTCGTATTGATACTTTg                                            <  2:218164/141‑1 (MQ=11)

                                                                                               gtTCGTCTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTCGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTCTTACAGCCTTCTTTGCTGTATTGATACttag                                            >  2:48978/1‑139 (MQ=11)

                                                                                               gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCG‑AAATCCGTAGGTC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCACTACCTTTTAAAGTATTCTCTGCTGTATTGATACTTTGa                                           <  1:526230/141‑1 (MQ=11)

                                                                                                 tCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                          <  2:776268/141‑1 (MQ=17)

                                                                                                  cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                         <  2:878693/138‑1 (MQ=17)

                                                                                                  cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                         <  2:515221/138‑1 (MQ=17)

                                                                                                  cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                         >  1:50235/1‑141 (MQ=17)

                                                                                                   gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc                                        <  2:764376/139‑1 (MQ=17)

                                                                                                   gCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTCTTAAAGTATTCTTTGCTGTATTGATGCTTTGATAc                                        >  1:509658/1‑141 (MQ=17)

                                                                                                   gCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc                                        <  2:638881/141‑1 (MQ=17)

                                                                                                    cgtAACGTGTTAGAATGTTGCTAATTTCGGGGCATGGCGAAAAGCCGGAGATC‑TGAAGAGACCTGCTGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                       <  1:194356/139‑1 (MQ=11)

                                                                                                    cTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATGCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                       <  2:106535/141‑1 (MQ=17)

                                                                                                     ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCt                                      <  2:84883/141‑1 (MQ=255)

                                                                                                     ttAACTTGCTAGCATGATGCTAATTGCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                      <  1:616335/141‑1 (MQ=17)

                                                                                                       accTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTg                                    <  1:168545/139‑1 (MQ=17)

                                                                                                        aCTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGt                                   <  2:895764/141‑1 (MQ=255)

                                                                                                        aCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                   <  2:537812/141‑1 (MQ=17)

                                                                                                         cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtt                                  >  1:206403/1‑140 (MQ=17)

                                                                                                          ttGCTAGCTGGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGCGCGTGAATACATTCAATACCTTTTTAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtct                                 <  2:286759/141‑1 (MQ=11)

                                                                                                          ttGCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                 <  1:181533/141‑1 (MQ=17)

                                                                                                           tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                <  1:539078/141‑1 (MQ=17)

                                                                                                            gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAGTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                               >  2:832834/1‑141 (MQ=11)

                                                                                                            gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                               <  2:801757/141‑1 (MQ=255)

                                                                                                            gCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtcttt                               <  1:279144/141‑1 (MQ=11)

                                                                                                             caaGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttc                              <  2:351128/139‑1 (MQ=11)

                                                                                                             cTAGCATGCTGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttt                              >  1:64409/1‑140 (MQ=11)

                                                                                                             cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                              >  1:629621/1‑141 (MQ=17)

                                                                                                              tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGGTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                             <  1:446522/141‑1 (MQ=11)

                                                                                                              tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                             >  1:133755/1‑141 (MQ=17)

                                                                                                                gCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTTTTCTTTGCCGTATTGATACTTTGATATCTTGtctttcttt                           >  2:628132/1‑140 (MQ=11)

                                                                                                                  atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                         >  1:60167/1‑141 (MQ=18)

                                                                                                                  atgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                         >  1:593270/1‑141 (MQ=21)

                                                                                                                  atgatgCTAATTTCGGGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACt                         <  2:231281/141‑1 (MQ=21)

                                                                                                                   tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCCATACCTTTTAAAGTCTTCTTTGCCGTATTGATACTTTGATGTATTGtctttctttctt                        >  1:900038/1‑141 (MQ=11)

                                                                                                                   tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                        <  1:68056/141‑1 (MQ=11)

                                                                                                                   tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACtt                        >  1:144926/1‑141 (MQ=21)

                                                                                                                   tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACtt                        >  1:342763/1‑141 (MQ=21)

                                                                                                                    gatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACttt                       <  1:856022/141‑1 (MQ=11)

                                                                                                                    gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                       <  2:50235/141‑1 (MQ=25)

                                                                                                                     atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATCCTTTGATACCTTGTCTTTCTTACTTTa                      >  2:395568/1‑141 (MQ=21)

                                                                                                                     atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCCTTGCTGTTTTGATGATTTGGTACCTTGTCTTTGTTACTTTa                      >  1:660844/1‑141 (MQ=11)

                                                                                                                      tgCTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGTTCCTTTGATATCTTGTCTTTCTTACTTTaa                     >  2:747522/1‑141 (MQ=11)

                                                                                                                      tgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACAATGAAAGCTCGTCTTTCTTACTTTaa                     <  1:226650/141‑1 (MQ=11)

                                                                                                                      tgCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTGTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTaa                     >  2:419364/1‑141 (MQ=21)

                                                                                                                       gCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAat                    >  1:695996/1‑141 (MQ=25)

                                                                                                                       gCTAATTTCGTGGCATAGGGAAAATCCGTAGATC‑TGAAGAGACCTCCGGTTCTTTTTATCTAGAGCGTCAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAat                    >  2:317813/1‑141 (MQ=11)

                                                                                                                        cTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAata                   <  2:860439/141‑1 (MQ=255)

                                                                                                                        cTAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAata                   >  2:372596/1‑141 (MQ=255)

                                                                                                                         tAATTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTTTCTTTCTTACTTCAatat                  >  2:253909/1‑141 (MQ=11)

                                                                                                                         tAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                  >  1:481905/1‑141 (MQ=25)

                                                                                                                          actttCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATg                 <  1:921000/139‑1 (MQ=11)

                                                                                                                          aaTTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGATACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATg                 <  2:20146/141‑1 (MQ=11)

                                                                                                                          aaTTTCGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATg                 >  2:582336/1‑141 (MQ=255)

                                                                                                                          aaTTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatctt                 >  2:134995/1‑138 (MQ=21)

                                                                                                                          aaTTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                 <  2:817999/141‑1 (MQ=25)

                                                                                                                           aTTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                <  2:438503/141‑1 (MQ=25)

                                                                                                                            tttCGTGGCGTGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGAc               <  1:87173/141‑1 (MQ=25)

                                                                                                                             ttccttgCATGGCTAAAATCCGTAGATC‑TGAAGAGACGTACCGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAACTATTCTTTGCCGTATAGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACg              <  1:35361/135‑1 (MQ=11)

                                                                                                                               cGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGt            <  2:151787/141‑1 (MQ=255)

                                                                                                                               cGTGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGCTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGt            <  2:480021/141‑1 (MQ=14)

                                                                                                                               cGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            <  1:431017/141‑1 (MQ=255)

                                                                                                                               cGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            <  1:407136/141‑1 (MQ=255)

                                                                                                                               cGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTAAATAGAGCGTAAATACACTCAATACCTTTTAAAGGATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTACCTTTAAAATGACGGt            >  1:118373/1‑141 (MQ=11)

                                                                                                                               cGTGGCATCGCGAAACTCCCTAGATC‑TGAAGAGAGCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTAGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt            <  1:784582/141‑1 (MQ=11)

                                                                                                                                 tGGCATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGa          >  1:800633/1‑141 (MQ=255)

                                                                                                                                  ggCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTCTTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt         >  2:331609/1‑141 (MQ=255)

                                                                                                                                    cATGGGGAAAATCCGTAGATC‑TGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTAAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTTATATCTTGTCTTTCTTACTTTAATATGACGGTGATCt       >  1:771370/1‑141 (MQ=11)

                                                                                                                                    cATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCt       >  2:184505/1‑141 (MQ=38)

                                                                                                                                      tGGCGAAAATCCGTAGATCTTGTTGTG‑CCTTCTGTTCGTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  2:658517/141‑1 (MQ=14)

                                                                                                                                      tGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTTTATTATACTTTAATATGACGGTGATCTTg     <  2:910169/141‑1 (MQ=14)

                                                                                                                                      tGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     >  1:706273/1‑141 (MQ=38)

                                                                                                                                      tGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg     <  1:104525/141‑1 (MQ=255)

                                                                                                                                      tGGCGAAAATCCGCAGATC‑TGACCACACCTGCGGTTCCTCTTATATAGAGCGTACATACACTCAATCCCTTTTAAAGTATTCTTTGCTGTTTTGATAATTTGATACCTTGTCTTTCTTACTTTATTATGACGGTGATCttt     >  2:270956/1‑140 (MQ=11)

                                                                                                                                         cGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTACTCTGACGGTGATCTTGcct  >  2:840880/1‑139 (MQ=255)

GAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATC‑TGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTC  >  CP000730/35776‑36053

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