BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,015	C→T	30.7%	R179R (AGG→AGA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,015	0	C	T	30.7%	469.7 / 294.6	380	R179R (AGG→AGA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (120/142); new base T (50/67); total (171/209)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-02

Only 100 of 386 total aligned reads displayed.

GCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATC  >  CP000730/35875‑36155

gCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAACATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc                                                                                                                                              >  1:279781/1‑141 (MQ=17)

gCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc                                                                                                                                              >  1:182974/1‑141 (MQ=17)

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

     aCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                                                                                                                         <  1:270219/141‑1 (MQ=17)

      cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATACAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                                                                                                        >  2:259482/1‑141 (MQ=17)

       ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                                                                                                       <  1:228445/141‑1 (MQ=17)

        tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAACTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                                      <  1:90825/141‑1 (MQ=17)

           tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                                                                                                                                   <  1:238270/141‑1 (MQ=255)

           tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                                                                   <  1:212475/141‑1 (MQ=17)

           tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttct                                                                                                                                   <  2:44956/141‑1 (MQ=17)

             gCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGTTACTTTGATATCTTGtctttcttt                                                                                                                                 >  2:69645/1‑140 (MQ=11)

             gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTa                                                                                                                                 >  2:118520/1‑141 (MQ=32)

               atgatgCTAATTTCGCGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTAcc                                                                                                                               <  1:329572/141‑2 (MQ=11)

                 gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                                                                                                                             <  2:21154/141‑1 (MQ=21)

                  atgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCCACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTa                                                                                                                            >  2:303305/1‑141 (MQ=11)

                   tgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTaa                                                                                                                           >  1:316558/1‑141 (MQ=21)

                   tgCCAATTGCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTCCGTTTGATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTaa                                                                                                                           <  1:10270/141‑1 (MQ=11)

                      tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAatat                                                                                                                        >  1:208824/1‑141 (MQ=11)

                      tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTCATACTTTGCTATCTTGTCTTTCTTACTTTAatat                                                                                                                        >  2:169074/1‑141 (MQ=11)

                       aaTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATg                                                                                                                       <  1:37917/141‑1 (MQ=11)

                       aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                                                                                                       >  1:295500/1‑141 (MQ=25)

                           tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                                                                                                   <  1:52068/141‑1 (MQ=38)

                           tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                                                                                                   <  1:303012/141‑1 (MQ=38)

                            cGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGt                                                                                                                  <  1:287888/141‑1 (MQ=14)

                                   tGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATCGAGCTTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                                           <  2:109257/141‑1 (MQ=14)

                                      cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGATGGTGATCGTGCTc                                                                                                        >  1:346409/1‑141 (MQ=255)

                                       gAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCa                                                                                                       >  1:61600/1‑141 (MQ=38)

                                         aaaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTACAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTCCTTTAATATGACGGTGATCTTGCTCAAt                                                                                                     >  2:39348/1‑141 (MQ=11)

                                         aaaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTATGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAAt                                                                                                     <  2:352330/141‑1 (MQ=14)

                                          aaTCCGTAGATCTGAAGAGAGCTGCGGTTCTTTTTATATAGAGCATAAATACGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTCTCTTCCTTTACTATGACGCTGATCTTGCTCAATg                                                                                                    >  1:217718/1‑141 (MQ=11)

                                          aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                                                                    <  2:232589/141‑1 (MQ=255)

                                          aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACGCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTCTCTTTCTTACTTTAATATGACCGTGATCTTGCTCAATg                                                                                                    <  1:200051/141‑1 (MQ=255)

                                            tCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAAGTGTTCTTCCTTACTTTAATATGACGGTGATCTTGCTCAATGAg                                                                                                  >  1:333852/1‑141 (MQ=11)

                                             ccGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAgg                                                                                                 <  1:343845/141‑1 (MQ=255)

                                              cgcaGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                                                                                                <  1:98835/138‑1 (MQ=255)

                                              cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                                                                                                <  2:86950/141‑1 (MQ=255)

                                                  gATCGGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                                                                                            <  2:206686/141‑1 (MQ=255)

                                                     cTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATTAGGTTATTCAg                                                                                         >  1:267034/1‑141 (MQ=255)

                                                     cTGAAGAGACCTACGGTTCTTTTTATATAGAGCGGAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCGTGCTTACTTTAGTATGACGGTGATCTTGCTCAATGAGGTTATTCaa                                                                                         >  2:161741/1‑140 (MQ=11)

                                                      tGAAGAGACCTGCGGTTCTTTTTATATATAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTGTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa                                                                                        <  1:202905/141‑1 (MQ=255)

                                                       gAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGat                                                                                       <  1:106199/141‑1 (MQ=255)

                                                        actagaCCTACGGTTCTTTTTATATAGAGAGTGAATGCTTTCAGGACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTTATATCTTGTCTTCCTTACTTTGATATGACGGTGATCTTGCTCAATGAGGTTATTCAGata                                                                                      <  1:153910/138‑1 (MQ=11)

                                                        aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCCATGAGGTTATTCAGata                                                                                      >  2:321095/1‑141 (MQ=255)

                                                         agagACCTACGGTTCTTTTTATATAGGGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGattt                                                                                     >  2:200518/1‑139 (MQ=255)

                                                          gagaCCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTGAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                                                                    <  1:250356/141‑1 (MQ=38)

                                                               cTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGAt                                                                               >  2:138293/1‑141 (MQ=255)

                                                                  cGGTTCTTTTTATAGAGGGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTCGTAGGATATGTCTATGtc                                                                            <  2:166320/141‑2 (MQ=14)

                                                                    gTTCTTGTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCCGATATTTCGATGTACa                                                                          <  1:177437/141‑1 (MQ=255)

                                                                     ttCTTTTTATATAGGGCGTAAGTACCTTTAAGACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACaa                                                                         <  2:78814/141‑1 (MQ=11)

                                                                      tCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                                        >  2:329177/1‑141 (MQ=38)

                                                                      tCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                                        >  1:312073/1‑141 (MQ=38)

                                                                       ttttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCCTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATg                                                                       <  2:133491/140‑1 (MQ=255)

                                                                         ttttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGAc                                                                     <  1:133359/141‑1 (MQ=255)

                                                                          tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCGTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCa                                                                    >  2:269471/1‑141 (MQ=255)

                                                                             ttataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTc                                                                 <  1:7238/140‑1 (MQ=255)

                                                                              tataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCa                                                                <  1:195761/141‑1 (MQ=255)

                                                                                taGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAgg                                                              <  1:118006/141‑1 (MQ=40)

                                                                                taGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAgg                                                              <  2:295500/141‑1 (MQ=255)

                                                                                aagcgCGTACATAAACTCAATACCTTTTAAAGTATTCTTTACTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAgg                                                              <  2:182341/137‑1 (MQ=14)

                                                                                  gagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtt                                                            >  2:196466/1‑141 (MQ=255)

                                                                                  gagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtt                                                            >  2:156725/1‑141 (MQ=255)

                                                                                    gCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttta                                                          >  2:147406/1‑141 (MQ=255)

                                                                                       tAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagt                                                       >  2:31533/1‑141 (MQ=255)

                                                                                        aaaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt                                                      <  1:314456/141‑1 (MQ=255)

                                                                                         aaTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttt                                                     <  2:175811/141‑1 (MQ=40)

                                                                                         aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACTGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                                     <  1:288736/141‑1 (MQ=255)

                                                                                          aTACATCCACTACCTTTTAAAGTACTAATTGACGTATTGAGACTTTGATATCTTGTCTTTCTTACTGTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttta                                                    >  1:263322/1‑141 (MQ=14)

                                                                                             cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAGGTTTaaaa                                                 >  1:36809/4‑141 (MQ=255)

                                                                                             cactCAATACCTTTTAAAGTATTCTTTGCTGTAATGATAGTGGGATGCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                                                 <  1:8743/138‑1 (MQ=14)

                                                                                             cATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGTCAGTCAGTTTTAAATTTaaaa                                                 >  1:144785/1‑141 (MQ=14)

                                                                                               ttCAATACCTTTTAACGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTGAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGc                                               <  1:68729/141‑1 (MQ=255)

                                                                                               ttCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAgt                                               >  2:48226/1‑140 (MQ=255)

                                                                                                 cAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCtt                                             <  1:321095/141‑1 (MQ=255)

                                                                                                  aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCttt                                            >  1:298542/1‑141 (MQ=255)

                                                                                                   aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTCATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGCTGTACAATGCCAGTCTGGTTTAAGTTTAAAAGCTTTa                                           >  1:126063/1‑141 (MQ=255)

                                                                                                     aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACTGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                                         <  1:169264/141‑1 (MQ=255)

                                                                                                     aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAt                                         <  2:66459/141‑1 (MQ=40)

                                                                                                      ccTTTTAAAGTATTTTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                                        <  2:36045/141‑1 (MQ=255)

                                                                                                      ccTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGCCCTGCAGCGATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                                        <  1:139426/141‑1 (MQ=14)

                                                                                                       cTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGATTATTCAGATATTTCGATGTACCATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                                       >  2:246634/1‑141 (MQ=255)

                                                                                                       cTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACATTGTCTTTCTTACTTTAATATGACGGTGTTCTTGCTCACTGAGGTTATTCCGATATTTCGATGTACATTGGCGTTCAGGTTTAAGTTTAAAAGCTTTAATTa                                       >  2:31065/1‑141 (MQ=14)

                                                                                                       cTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                                       >  2:125065/1‑141 (MQ=255)

                                                                                                         tttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCTATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACt                                     <  1:188446/141‑1 (MQ=255)

                                                                                                          ttAAAGTATTCTTTGCTGTATGTATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACtt                                    <  2:300353/141‑1 (MQ=255)

                                                                                                             aaGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAg                                 >  2:257114/1‑141 (MQ=255)

                                                                                                                 aTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAt                             >  1:59802/1‑141 (MQ=255)

                                                                                                                 aTTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCAt                             <  1:147911/141‑1 (MQ=255)

                                                                                                                  ttCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAtt                            >  2:288291/1‑141 (MQ=255)

                                                                                                                    cTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGc                          <  2:61759/141‑1 (MQ=255)

                                                                                                                    cTTTGCCGTATTGATACTTTGATATCTTGTCTTACTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGc                          <  1:243094/141‑1 (MQ=14)

                                                                                                                     tttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAACTTGCTCAATGAGGTTATTCAGATATTTCGTTGTACAATGGTGGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt                         >  2:22007/1‑141 (MQ=14)

                                                                                                                        gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc                      <  1:29142/141‑1 (MQ=255)

                                                                                                                        gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAAAATGACGGTGATCTTTCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTACAAGCTTTCATTACTTTAACCATTGCTAcc                      <  2:342920/141‑1 (MQ=14)

                                                                                                                              ttGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGtt                <  1:350985/141‑1 (MQ=255)

                                                                                                                              ttGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGtt                >  2:227267/1‑141 (MQ=255)

                                                                                                                                gATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTACTTACTTTAGCCATTGCTACCTTCGTTGa              >  2:214587/1‑141 (MQ=255)

                                                                                                                                 aTACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGaa             >  1:19233/1‑141 (MQ=255)

                                                                                                                                    cTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAgtt          >  1:309371/1‑139 (MQ=255)

                                                                                                                                    cTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGATTAAGTTTAAAAGCTTTAATTACTTCAAACATGGCTACCTTCGTTGAAGGt          <  2:254896/141‑1 (MQ=14)

                                                                                                                                      ttGATCCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGTAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGc        <  2:170842/141‑1 (MQ=255)

GCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATC  >  CP000730/35875‑36155

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