BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,161	T→C	28.6%	I131V (ATT→GTT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,161	0	T	C	28.6%	322.1 / 219.9	310	I131V (ATT→GTT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (103/117); new base C (40/49); total (143/167)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 316 total aligned reads displayed.

TTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36022‑36301

  tcttACTTTAATATGTCGGTGATCTTGCTCAATGCGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTAc                                                                                                                                            <  2:29230/141‑1 (MQ=11)

   cttACTTTAATATTACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTACCCATTGCTACCTTCGTTGAAGATGCCTGATCTGTAATTAcc                                                                                                                                           >  1:34098/1‑141 (MQ=255)

    ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                                                          >  1:181127/1‑141 (MQ=255)

     ttcTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTCCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATAACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTGATTACCtt                                                                                                                                         <  2:293991/139‑1 (MQ=11)

      aCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCttt                                                                                                                                        <  1:351610/141‑1 (MQ=255)

        tttAATTTGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                                                      <  2:394024/141‑1 (MQ=255)

           aaTATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑Agg                                                                                                                                  <  2:51079/141‑1 (MQ=25)

             tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑Aggct                                                                                                                                >  2:89165/1‑139 (MQ=255)

               tGACGGTGATCTTGCTCACTTAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCGGTAACTACCTTTTG‑AGGcttc                                                                                                                              >  2:385470/1‑140 (MQ=11)

               tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtt                                                                                                                                         <  2:171302/131‑1 (MQ=35)

                tcagctGATCTTGCTCAATTATGTTAGTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTCTAATTACCTTTTG‑AGGTTTAc                                                                                                                             <  2:273905/136‑1 (MQ=11)

                gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTAc                                                                                                                             <  1:380012/141‑1 (MQ=21)

                 gcGGTGATCTTGCTTAATGATGTTATTCATATATTTCGATGTACACTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCttttcaagcgtacc                                                                                                                             <  2:47245/140‑11 (MQ=11)

                 aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTAcc                                                                                                                            >  2:263950/1‑141 (MQ=25)

                 aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTAcc                                                                                                                            <  1:20364/141‑1 (MQ=11)

                   ggTGATCTGGCTCAATGAGGCTGTTCGGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCaa                                                                                                                          <  2:399790/141‑1 (MQ=255)

                      gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTCCCTTTTG‑AGGTTTACCAAAtt                                                                                                                       >  2:10635/1‑141 (MQ=12)

                      gATCTTGCTCAATGAAGTGATACAGATATTTCGATGTACAATGACGGTCAGGTTTAAGTTTAAAAGCTTCACTTATTTTAGCCCTTGCTGCCTTCGTTCACGGTGCATGCTCTGTACCTGCCTTTGG‑Agaccccccaactt                                                                                                                       >  2:3354/1‑129 (MQ=2)

                        tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTTATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTT‑AGGCTTACACAATTGt                                                                                                                     >  1:102962/1‑141 (MQ=11)

                        tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTCACTACCTTTTG‑AGGCTGACCagcgtgt                                                                                                                     >  2:56558/1‑135 (MQ=11)

                        tCTTGCTCAATGAGGGTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGt                                                                                                                     <  2:22050/141‑1 (MQ=12)

                         cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGtt                                                                                                                    >  1:162852/1‑141 (MQ=17)

                         cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGtt                                                                                                                    >  1:199995/1‑141 (MQ=255)

                          ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGttt                                                                                                                   >  1:55641/1‑141 (MQ=17)

                             ttCAATGGAGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTGCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAAt                                                                                                                <  2:226175/140‑1 (MQ=11)

                               cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATga                                                                                                              <  2:250904/141‑1 (MQ=12)

                               aaaTGTGGTTATTCAGATATTTCGATGTACAATGACAGTCAGTTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGTTGCCTGATCTGTAACT‑CCTTTTGTATCCTTACCAAATTGTTTAATga                                                                                                              <  2:85482/140‑1 (MQ=11)

                                    aGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTCATGAGACGt                                                                                                         >  1:106989/1‑141 (MQ=11)

                                     ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                        <  1:132108/141‑1 (MQ=11)

                                      gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                       <  1:316058/141‑1 (MQ=12)

                                      gTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGttt                                                                                                       <  2:144231/141‑1 (MQ=255)

                                       ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGtttt                                                                                                      >  2:340049/1‑140 (MQ=12)

                                       ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                      >  1:392516/1‑141 (MQ=17)

                                       ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                      >  2:159497/1‑141 (MQ=255)

                                       ttATTCAGATATTTCCATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTCCTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                      <  1:348859/141‑1 (MQ=255)

                                        tATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTAACAAATTGTTTAATGAGACGTTTGa                                                                                                     >  2:68824/1‑141 (MQ=11)

                                         aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTTCTTTAACCATTGCTACCTTCGTTGAAGGATCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGaa                                                                                                    >  1:172714/1‑140 (MQ=11)

                                          ttCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                   <  2:330255/141‑1 (MQ=17)

                                          ttCAGAATTTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                   <  2:116758/141‑1 (MQ=11)

                                             aGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                                                                >  2:351429/1‑141 (MQ=17)

                                             aGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                                                                >  2:23891/1‑141 (MQ=32)

                                               atatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                              >  1:330742/1‑141 (MQ=17)

                                                 atTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                            <  1:325226/141‑1 (MQ=11)

                                                   ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTAATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                          <  1:113770/141‑1 (MQ=12)

                                                    tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGTCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                         >  1:73561/1‑141 (MQ=255)

                                                    tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                         <  1:178690/141‑1 (MQ=17)

                                                     cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  2:43526/141‑1 (MQ=17)

                                                     cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  2:38041/141‑1 (MQ=17)

                                                      gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                       >  2:395809/1‑141 (MQ=17)

                                                       aTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                      >  2:181301/1‑141 (MQ=17)

                                                       aTGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATACGCTg                                                                                      >  2:148591/1‑141 (MQ=17)

                                                          tACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAt                                                                                   <  1:233291/141‑1 (MQ=17)

                                                          tACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAt                                                                                   <  1:364936/141‑1 (MQ=17)

                                                            cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                 >  2:155486/1‑141 (MQ=32)

                                                              aTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTGGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGATTACCCAATTCTTTAATGAGAAGTTTGATAAACGCATATGCTGAAGGAtt                                                                               <  2:36211/141‑1 (MQ=11)

                                                                gACAGTCAGGTTTAAGTTTAAATATAGTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATTGGCTGAATGATTAt                                                                             >  1:95792/1‑141 (MQ=11)

                                                                 gcaatCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTAGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTATAATGAGACGTTTGATACACCCATATGCTGAATGATTAtc                                                                            <  2:285543/137‑1 (MQ=1)

                                                                                               tattttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgt                                              <  2:307588/138‑1 (MQ=0)

                                                                                                  tttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           <  2:215189/141‑1 (MQ=11)

                                                                                                    tcgctaTTGCATCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAACTTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  1:25274/136‑1 (MQ=11)

                                                                                                    taaccATTGCTACCTTCGTTGAAGGTTCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  1:26428/138‑1 (MQ=17)

                                                                                                     aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        <  2:397963/141‑1 (MQ=12)

                                                                                                     aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        >  1:46284/1‑141 (MQ=12)

                                                                                                      gCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTccc                                       >  1:139572/1‑141 (MQ=255)

                                                                                                       ccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcccc                                      >  2:63103/1‑140 (MQ=17)

                                                                                                       ccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCct                                      <  1:370201/141‑1 (MQ=32)

                                                                                                        cATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  1:246197/141‑1 (MQ=12)

                                                                                                        cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTAGCAAATTGTTTAATGAGACGTTTGATCAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  1:286310/141‑1 (MQ=17)

                                                                                                         aTTGCTACCTTCGTTGCAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                    >  1:198023/1‑141 (MQ=17)

                                                                                                          ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                   >  1:384584/1‑141 (MQ=32)

                                                                                                           tGCTACCTTCGTTTAAGGTGCCTGATCTGTAATTATCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                  <  1:102138/141‑1 (MQ=17)

                                                                                                           tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGTCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                  >  1:33480/1‑140 (MQ=17)

                                                                                                            gCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATCTATGTCCCTCTGCa                                 <  1:208430/141‑1 (MQ=17)

                                                                                                             ccaCCTTCGTCGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCACATTGTTTAATGCGACGTTTGATAAACTCAGATGCTGAATTAGTATCTCGTTGCTGACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                <  2:94132/139‑1 (MQ=11)

                                                                                                              tgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCACTCTgtgtc                               >  2:138665/3‑137 (MQ=11)

                                                                                                                 tttCGTTTAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                            <  2:77912/140‑1 (MQ=17)

                                                                                                                  ttCTTTGAAGGTGCCTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                           <  2:301824/141‑1 (MQ=255)

                                                                                                                    cGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGCCGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTCCCTCAATGGc                         >  2:227023/1‑141 (MQ=12)

                                                                                                                    cGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGc                         <  1:398934/141‑1 (MQ=32)

                                                                                                                     gTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGGTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCATTGGCa                        >  2:192682/1‑141 (MQ=17)

                                                                                                                     gTTGAAGGTGCCTGATCTGTAACCACCTTTTG‑ATGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCa                        <  1:47687/141‑1 (MQ=11)

                                                                                                                      ttGAAGGTGCCTCATCTGTAATTACCTTTTG‑AGCTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATAATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGGCCCTCTGCATCAATGGCAc                       <  2:223034/141‑1 (MQ=255)

                                                                                                                         aaGGTGCCTGATCTGTAACTGCCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTTTATCAATGGCACGat                    <  2:343074/141‑1 (MQ=11)

                                                                                                                          aGGTACCTGATCTGTAATTCGCTTTTG‑AGGTTTACCAAATTGTTTAGTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCGACCAAGTATCTAATGTATGTCCCTCTGGATCAATGGCACGata                   >  1:280822/1‑141 (MQ=11)

                                                                                                                           ggTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatat                  >  1:35830/1‑141 (MQ=32)

                                                                                                                               ccTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATTAATGTCACGATTTAAAt              >  1:148790/1‑141 (MQ=11)

                                                                                                                               ccTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt              >  1:311478/1‑141 (MQ=12)

                                                                                                                               ccTGATCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt              >  1:57520/1‑141 (MQ=12)

                                                                                                                                cTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATa             <  2:325634/141‑1 (MQ=32)

                                                                                                                                  gATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAGTGGCACGATATAAATAGc           <  1:318952/141‑1 (MQ=17)

                                                                                                                                   aTCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt          <  2:253442/141‑1 (MQ=17)

                                                                                                                                   aTCTGTAACTACCTTTTG‑AGGCTTACCAAATTGTTTATTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt          >  1:191930/1‑141 (MQ=255)

                                                                                                                                    tCTGTAATTACCTGTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc         <  1:121836/141‑1 (MQ=17)

                                                                                                                                     cTGTACTTACCTTTGT‑ATGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTACCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACCATATAAATAGCTcc        <  1:213130/141‑1 (MQ=11)

                                                                                                                                      tGTAATTACCTTTTT‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa       <  2:198305/141‑1 (MQ=17)

                                                                                                                                        tAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCACCAtt     >  2:48715/1‑141 (MQ=17)

                                                                                                                                          aTTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGAACGATATAAATAGCACCAttat   >  2:39095/1‑139 (MQ=12)

                                                                                                                                           ttATCTTTTG‑AGGCTTACCTAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc  <  1:89165/141‑1 (MQ=255)

                                                                                                                                           ttACCTTTTG‑AGGTTTACCTAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTc  <  1:155470/141‑1 (MQ=17)

                                                                                                                                           ttACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTTCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTc  >  1:353752/1‑141 (MQ=17)

TTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTG‑AGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36022‑36301

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