BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,940	G→A	32.4%	R204R (CGC→CGT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,940	0	G	A	32.4%	215.7 / 141.6	196	R204R (CGC→CGT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (75/57); new base A (27/37); total (102/94)
Fisher's exact test for biased strand distribution p-value = 6.74e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 198 total aligned reads displayed.

TAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGT  >  CP000730/35804‑36080

 aataaAAAT‑GATGGTTCTGTTGCAAAGT‑AAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTGGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt                                                                                                                                        <  2:53817/141‑1 (MQ=11)

       aaTAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttt                                                                                                                                    <  1:133753/141‑1 (MQ=34)

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

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

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

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

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

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

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

                                                                     tcgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                                                      >  2:289801/5‑141 (MQ=12)

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

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

                                                                      cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGCTACTTTGATa                                                                     >  2:304484/4‑141 (MQ=12)

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

                                                                       ttttAACTTGCCAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATat                                                                    <  2:310489/139‑2 (MQ=11)

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

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

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

                                                                         ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatct                                                                  >  2:183310/1‑138 (MQ=255)

                                                                         ttAACTTGCTAGCATGATGCTAAGTTCTTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCt                                                                  <  2:260539/141‑1 (MQ=11)

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

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

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

                                                                             cTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATATATGTAACTCCGTGGAAATTATTCTTTGCCGTATTGATACTTTGATATCTTGtc                                                              >  1:25654/1‑141 (MQ=11)

                                                                             cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                              >  1:278805/1‑141 (MQ=17)

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

                                                                                gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt‑t                                                          <  1:312004/141‑1 (MQ=32)

                                                                                gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt‑t                                                          >  1:199532/1‑141 (MQ=17)

                                                                                 cTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctt‑tc                                                         <  1:90009/141‑1 (MQ=255)

                                                                                 cTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctt‑tc                                                         >  1:132981/1‑141 (MQ=255)

                                                                                 cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt‑tc                                                         <  1:72235/141‑1 (MQ=17)

                                                                                   aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGGTGTATTGATACTTTGATACCTTGtctt‑tctt                                                       >  1:45867/1‑141 (MQ=17)

                                                                                    gCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTa                                                      >  2:190284/1‑141 (MQ=17)

                                                                                       tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACtt                                                   >  2:238728/1‑141 (MQ=25)

                                                                                        gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACttt                                                  >  2:166491/1‑141 (MQ=25)

                                                                                        gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACttt                                                  <  2:265513/141‑1 (MQ=25)

                                                                                        gatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACttt                                                  <  2:274575/141‑1 (MQ=25)

                                                                                        gatgCTAATTTCGGGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACttt                                                  <  2:120776/141‑1 (MQ=11)

                                                                                          tgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTaa                                                >  1:247614/1‑141 (MQ=21)

                                                                                           gCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAat                                               >  1:299381/1‑141 (MQ=21)

                                                                                             tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAatat                                             >  2:171515/1‑141 (MQ=255)

                                                                                             tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAatat                                             >  2:73882/1‑141 (MQ=25)

                                                                                             tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAatat                                             <  1:117968/141‑1 (MQ=25)

                                                                                              aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGCGACCTGCAGTTCTTTTTATATAGAGCGTTAATACCGCCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATg                                            <  1:159455/141‑1 (MQ=11)

                                                                                               aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGa                                           >  1:274263/1‑141 (MQ=21)

                                                                                                tttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGATCTTTTTATATAGAGCGTAAATACACTCAATAACTTTTAAAGTATTCTTTGCTGTATTGGTACTTTGATACCTTTTCTT‑TCTTACTTTAATATgcc                                          >  2:123672/1‑139 (MQ=11)

                                                                                                tttCGTGACATGGCGAAAATCCGGAGATCCCAAGAGACCTGCGGATATTTTTATATAGAGCGTAAATACACCCAATACCTTTTACAGTAGTCTGAACTGAATGGATACCTTACTACCTCGTCTT‑TCTCACttaaccatgcc                                          >  2:47485/1‑132 (MQ=2)

                                                                                                 ttCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACg                                         >  1:107022/1‑141 (MQ=255)

                                                                                                  tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACgg                                        <  1:193910/141‑1 (MQ=38)

                                                                                                   cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTTAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGt                                       <  1:171511/141‑1 (MQ=255)

                                                                                                   cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATCCCTTGTCTT‑TCTTACTTGAATATGACGGt                                       <  2:138357/141‑1 (MQ=255)

                                                                                                   cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGt                                       <  1:183660/141‑1 (MQ=38)

                                                                                                     tGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCAGTCTTGTTACTTTGTTATCTTGTCTT‑TCTTACTTTAATATGACGGTgt                                     >  2:234114/1‑140 (MQ=11)

                                                                                                     tGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATCACGGTGa                                     <  1:209436/141‑1 (MQ=255)

                                                                                                        cATGGCGAAAATCCGTAGATCTGACGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGATGATCt                                  >  2:44394/1‑141 (MQ=255)

                                                                                                        cATGGCGAAAATCCGGAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTACAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCt                                  <  1:32406/141‑1 (MQ=255)

                                                                                                         aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCtt                                 >  2:281442/1‑141 (MQ=11)

                                                                                                         aTGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCtt                                 >  2:164842/1‑141 (MQ=255)

                                                                                                         aTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAAGACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTAATAATTTGATAACTTCTCTA‑TCTTACCTTAATAGGACGGTGATCtt                                 >  2:214219/1‑141 (MQ=11)

                                                                                                          tGGGGAAAATCCGTAAATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTATTATGACGGTGATCTTg                                <  1:165294/141‑1 (MQ=14)

                                                                                                           ggggAAAGTGCGTAGATCGGAAGGGGCCTGCGGTTGTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATAGGACGGTGATCTTGc                               <  2:198663/138‑1 (MQ=11)

                                                                                                           ggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACGTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGc                               <  1:105510/138‑1 (MQ=14)

                                                                                                           ggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAACGTAAATACATTCAATCCCTTTTTAAGTATTCTTTGCCGTATTGAGACTTTGATATCTTGTCTT‑TATTACTTTAATATGACggggtacttac                               >  2:9294/4‑132 (MQ=11)

                                                                                                            gCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTACAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCt                              >  1:63394/1‑141 (MQ=255)

                                                                                                              gAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCa                            >  1:237640/1‑141 (MQ=38)

                                                                                                                aaaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTATGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTCGTCAAt                          >  1:224549/1‑141 (MQ=255)

                                                                                                                aaaTCCGCAGATCTAAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAAt                          <  2:151065/141‑1 (MQ=11)

                                                                                                                 agtCCGTAGATATGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACTTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTATAATTTGACGGTGATCTTGCTCAATg                         <  2:223075/139‑1 (MQ=11)

                                                                                                                 aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATg                         <  2:65858/141‑1 (MQ=38)

                                                                                                                 aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATg                         >  2:222587/1‑141 (MQ=38)

                                                                                                                  aTCCGTAGATCTGAATAGACCTGCCGTTCTTTTTATATAGTGCGTAAATACACTCAATACCTTTTACAGTATTCTTTGCTGTATTTATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGa                        <  2:169436/141‑1 (MQ=11)

                                                                                                                  aTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGa                        <  1:137689/141‑1 (MQ=255)

                                                                                                                   tCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAg                       <  1:5660/141‑1 (MQ=38)

                                                                                                                     cGTAGATCTGAAGGGGACTACGCGTCTTTTTATATAGAGCGTAAATACATTTAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTAGCTCGATGAGGt                     <  2:39438/141‑1 (MQ=11)

                                                                                                                     cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                     >  2:205719/1‑141 (MQ=38)

                                                                                                                     cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                     <  1:289801/141‑1 (MQ=255)

                                                                                                                     cGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAGTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGGTCTTGCTCAATGGGGt                     >  1:91761/1‑141 (MQ=11)

                                                                                                                     cGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                     >  2:296928/1‑141 (MQ=255)

                                                                                                                       tAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTa                   <  2:144501/141‑1 (MQ=255)

                                                                                                                          aTCTGTAGAGATCTGCTGTTCTTTTTATATATGGGGTATATACACCCAATACCTTTTAAAGTATTCTTTGCTGTACTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                <  2:319474/141‑1 (MQ=11)

                                                                                                                          aTCTGACGCGACCTGCGCTCCTTTTTATATAGACCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                <  1:189798/141‑1 (MQ=14)

                                                                                                                          aTCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                >  1:43951/1‑141 (MQ=38)

                                                                                                                          aTCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                <  2:248584/141‑1 (MQ=255)

                                                                                                                           gcTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCa               >  2:161626/2‑141 (MQ=37)

                                                                                                                             tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGTTCAATGAGGTTATTCAGa             <  2:131659/141‑1 (MQ=255)

                                                                                                                             tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa             <  1:139769/141‑1 (MQ=255)

                                                                                                                             tGAACAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa             <  2:10169/141‑1 (MQ=255)

                                                                                                                              gAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGat            <  1:296632/141‑1 (MQ=255)

                                                                                                                              gAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGat            <  1:109372/141‑1 (MQ=255)

                                                                                                                               aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGata           <  1:1121/141‑1 (MQ=38)

                                                                                                                               aaGAGACCTGCGGGTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATTTGACGGGTATCTTGCTCAATGAGGtactccagact           >  2:50374/1‑131 (MQ=11)

                                                                                                                                 gagaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt         <  2:103487/141‑1 (MQ=255)

                                                                                                                                  agaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAttt        <  2:150899/141‑1 (MQ=255)

                                                                                                                                    aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATGCGGATATTTCg      >  1:12270/1‑141 (MQ=255)

                                                                                                                                    aCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGAATTATGTTGC‑TGAAAATGACGGTGATCTTGCTCAATGCGGTTATTCAGATATTTCg      >  1:326883/1‑141 (MQ=11)

                                                                                                                                     ccTGCGGTTCTTTTTAGATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGa     <  2:262062/141‑1 (MQ=255)

                                                                                                                                      cTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGAt    <  2:329679/141‑1 (MQ=38)

                                                                                                                                      cTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGAt    >  2:104684/1‑141 (MQ=38)

                                                                                                                                        gCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTTTTGTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGt  <  1:75815/141‑1 (MQ=255)

                                                                                                                                        gCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGt  >  1:234454/1‑141 (MQ=38)

TAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTT‑TCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGT  >  CP000730/35804‑36080

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