BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,174	T→C	34.5%	K126K (AAA→AAG)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,174	0	T	C	34.5%	594.5 / inf	785	K126K (AAA→AAG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (255/254); new base C (130/141); total (387/398)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 807 total aligned reads displayed.

TATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36035‑36314

tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGtt                                                                                                                                             <  1:529228/141‑1 (MQ=25)

tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGtt                                                                                                                                             >  1:568346/1‑141 (MQ=25)

  agACTGTGATATTGCTCAATTAAGTTATTTAGATATTTCGATGCACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTa                                                                                                                                           <  2:360919/140‑1 (MQ=11)

   gACGGTGACCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTAc                                                                                                                                          <  1:100669/141‑1 (MQ=255)

    aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTAcc                                                                                                                                         >  1:311126/1‑141 (MQ=255)

     cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCa                                                                                                                                        <  1:260114/141‑1 (MQ=255)

      ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCaa                                                                                                                                       <  2:131338/141‑1 (MQ=18)

       gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCaaa                                                                                                                                      >  2:11101/1‑141 (MQ=255)

       gTGATCTTGCTCAATGAGGTTATTCAGATACTTCGATTTACAATGGCAGTCACGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCGACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCaaa                                                                                                                                      <  2:417854/141‑1 (MQ=11)

          atgtggCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTGTACCAAATTg                                                                                                                                   <  2:241024/136‑1 (MQ=255)

          aTCTTGCTCAATGAGGCTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTg                                                                                                                                   <  1:11876/141‑1 (MQ=11)

            cTTGCTCAATGAGGTTCTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  1:395579/1‑141 (MQ=12)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGCTCTGTAATG‑ACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 <  1:158948/141‑1 (MQ=255)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 <  2:76823/141‑1 (MQ=17)

            cTTGCTCAATGAGGNTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  2:379754/1‑141 (MQ=12)

               gcccAATGTGGTTATTCATATATTTCGCTGTACAATGACAGTCAGGTTTTAGTTTAAAAGCTTTAATTCTTTTAGCGATTGATGCCTTCGTGGAAGGTGCCTGTTCAGTAGCT‑ACCTTTTGTGGCTATCCACATTGttgac                                                                                                                              <  2:210154/138‑4 (MQ=255)

                 tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTACTTCCTTTAACCATTGCCACCTTCGTTGCAGGTGCCTGATCCGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTattc                                                                                                                            <  2:588886/141‑4 (MQ=11)

                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                                           >  2:357499/1‑141 (MQ=255)

                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           <  1:55809/141‑1 (MQ=32)

                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           <  1:598979/141‑1 (MQ=32)

                   aaTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATgag                                                                                                                          <  2:568885/141‑1 (MQ=11)

                    aTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATgaga                                                                                                                         <  1:478991/141‑1 (MQ=11)

                      gAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAACTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGCAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACg                                                                                                                       <  2:446821/141‑1 (MQ=11)

                       gggTTATGCAGGTATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      <  1:83461/140‑1 (MQ=255)

                         gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  1:245270/141‑1 (MQ=17)

                         gTTATTCAGATATTTCGATGTACAATGACAGGCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  1:507389/141‑1 (MQ=11)

                          ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtttt                                                                                                                   >  2:39758/1‑140 (MQ=12)

                          ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  2:642057/141‑1 (MQ=17)

                            aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGTTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATTAGACGTTTGAt                                                                                                                 >  1:279871/1‑141 (MQ=255)

                            aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAGTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGCCGTTTGAt                                                                                                                 >  2:272665/1‑141 (MQ=255)

                             ttCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTTAGGTTTACCAAATTGGTTAATGAGAAGTTTGATa                                                                                                                >  1:45936/1‑141 (MQ=11)

                               cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              >  1:96091/1‑141 (MQ=255)

                                 gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:555837/141‑1 (MQ=17)

                                       tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      >  1:228382/1‑141 (MQ=17)

                                       tCGATGTACAATGACAGTTAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGAGCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      <  1:576847/141‑1 (MQ=11)

                                       tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTATCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      <  1:382940/141‑1 (MQ=11)

                                        tgATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAGTATGTTAGCCAGTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  1:635599/140‑1 (MQ=11)

                                        cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTAGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATCTGc                                                                                                     >  2:461435/1‑141 (MQ=255)

                                        cGATGTACAATGGCAGTCAGGTCTAAGTTTAAAAGCTTTAATTACTTCAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAAGTGTACCAAATTGTGTAATGACACGTTGGATAAACGCatatac                                                                                                     >  2:486807/1‑139 (MQ=11)

                                        cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGTTGTTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAACTTCTTTAACTAACCGTTTGATaatagcggatgc                                                                                                     >  1:301171/1‑131 (MQ=11)

                                        cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  1:471334/141‑1 (MQ=11)

                                         gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                                    >  1:414477/1‑141 (MQ=255)

                                         gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                                    >  2:349560/1‑141 (MQ=32)

                                          aTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                                   >  2:542001/1‑141 (MQ=12)

                                          aTGTACAATGGCAGCCAGGTTAAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGAGCTGTAATT‑ACCTGTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCta                                                                                                   >  2:455874/1‑140 (MQ=11)

                                           tGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGACGATGCCAGATCTGTAATT‑ACCTTTTGAGGTTCACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                                  <  1:497305/141‑1 (MQ=11)

                                             tACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAACTTGTTTAATGAGACGTTTGGTAAACGCATATGCTGaag                                                                                                >  1:190638/1‑140 (MQ=17)

                                               cAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                              <  1:542160/141‑1 (MQ=9)

                                               cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                              >  2:494316/1‑141 (MQ=11)

                                                aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             >  1:373387/1‑141 (MQ=2)

                                                  tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           <  1:148379/141‑1 (MQ=255)

                                                                                        aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTACTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTAGTTGCTTATGATACAATATATATAatgtatgttc                                                     >  1:13325/1‑139 (MQ=11)

                                                                                        aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     >  1:443732/1‑141 (MQ=255)

                                                                                          ccATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCct                                                   >  1:550243/1‑141 (MQ=12)

                                                                                           ctttGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  2:390455/139‑1 (MQ=255)

                                                                                           cATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  2:578123/141‑1 (MQ=12)

                                                                                              tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGCGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTga                                               >  2:86094/1‑140 (MQ=17)

                                                                                               gctgccTTCGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCATATATCTAATGTATGTCCCTCTgta                                              <  1:426647/137‑3 (MQ=11)

                                                                                                 tACCTCCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAAGCACAACCCCTTGCTAACGCAAACGAATTTCTAatgtatgtgcttttgtaac                                            >  2:596854/1‑130 (MQ=11)

                                                                                                  tccTTCGTTGAAGGTGCCTTATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAATCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                           <  2:437205/140‑1 (MQ=17)

                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                           <  2:161893/141‑1 (MQ=32)

                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTAGGTCCCTCTGTATCa                                           >  1:588515/1‑141 (MQ=17)

                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTAGGTCCCTCTGTATCa                                           <  1:368251/141‑1 (MQ=17)

                                                                                                  aCCTTAGTTGAAGGTGCCTGATCTGTAATT‑ACATTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  1:213909/141‑1 (MQ=12)

                                                                                                   ccTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          >  2:190722/1‑141 (MQ=17)

                                                                                                       cGTTGAAGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGc                                      <  1:389058/141‑1 (MQ=255)

                                                                                                         ttGAAGGTGCCTGATCTGTCATTGATCTTGT‑AAGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                                    >  1:489911/1‑141 (MQ=11)

                                                                                                         ttGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:200321/141‑1 (MQ=32)

                                                                                                            aaGGTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTTTGTATCAATGGAACGat                                 >  2:371623/1‑141 (MQ=11)

                                                                                                            aaGGTGCCTGATCCGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGACACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGat                                 <  1:287840/141‑1 (MQ=12)

                                                                                                              ggTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCAAATGTGTGTACCTATGTATCAATGGCACGctat                               >  2:291453/1‑141 (MQ=11)

                                                                                                               gTGCCTGATCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAtcta                              >  2:333906/1‑141 (MQ=255)

                                                                                                               gTGACTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAGTGAGACGTTTGATAAACGCATTTTCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTACTTCAGTggacccattta                              >  1:329301/1‑132 (MQ=11)

                                                                                                                tGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATaa                             <  2:514312/141‑1 (MQ=17)

                                                                                                                  ccTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCCAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                           >  2:532959/1‑141 (MQ=17)

                                                                                                                  ccTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTCTGTCACTGTGCATCAGTgggccgtaattaat                           >  2:667855/1‑129 (MQ=11)

                                                                                                                     gATCTGTAATT‑AACTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGGAATTGGTGTTTGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCATGGCTCGATATAAATAGc                        >  2:303438/1‑141 (MQ=11)

                                                                                                                      aTCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATCCGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACAAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  2:170124/141‑1 (MQ=11)

                                                                                                                       tCTTTAATT‑ACCTTTTGAGGTTTACTAAATTGTTTAATTAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc                      <  1:68981/141‑1 (MQ=255)

                                                                                                                       tCTGTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATAGAACGAGATAAGTACCTc                      <  2:475670/141‑1 (MQ=11)

                                                                                                                        cTGTAATT‑CCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                     <  2:627492/141‑1 (MQ=17)

                                                                                                                         tGTAACT‑ACCTTTTCAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa                    <  2:123105/141‑1 (MQ=255)

                                                                                                                          gTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAt                   >  1:342769/1‑141 (MQ=32)

                                                                                                                          gTAACT‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACCATATAAATAGACCCAt                   >  2:655494/1‑141 (MQ=11)

                                                                                                                            aact‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGAATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAttt                 <  1:252878/138‑1 (MQ=255)

                                                                                                                             act‑ACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                >  1:360781/3‑141 (MQ=12)

                                                                                                                             aTT‑ACCTTTTGAGGTTTACCAAATTGTTTATTGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAtttt                <  1:170078/141‑1 (MQ=17)

                                                                                                                              tt‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc               <  2:154557/141‑1 (MQ=17)

                                                                                                                              tt‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTc               >  1:106249/1‑141 (MQ=32)

                                                                                                                               t‑ACCGTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTTCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              <  2:221003/141‑1 (MQ=12)

                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:517968/141‑1 (MQ=17)

                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCCTTTTCCtt            >  2:316868/1‑141 (MQ=17)

                                                                                                                                    ttttGAGGCTTACCAAGTTGTTTAATGAGACGTTTTATAAAGGCATAGGCTGAATGATTCTCTCCTTGCTTACGCAACCAAATATCCAATGTATGTCCCTCTGTATCAATGGCACGATATAAACAGCTCCATTTTCCtttt          <  2:503570/141‑1 (MQ=11)

                                                                                                                                     tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTACTTACGCAACCAAATAGCTAATGTATGTACCTCTGTATCAATGGCACGATGTAAATAGCCACATTTTCCATTTa         >  2:263494/1‑141 (MQ=11)

                                                                                                                                     tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAACGATTATCTCGTTGCTTACGCAACCAAATATCTATTTTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTa         <  2:366326/141‑1 (MQ=12)

                                                                                                                                      ttGAGGTTTACCAAATTGTTTAATGAGACGTTTTATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAt        <  1:652326/141‑1 (MQ=17)

                                                                                                                                       tGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtt       <  2:615541/141‑1 (MQ=17)

                                                                                                                                           gTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCTACCAAATATCTAATGTATGTCCCTCTGCATCATTGGCACGATATAAATAGCTCGATTTTCCTTTTAtttttt   >  1:18013/1‑139 (MQ=12)

                                                                                                                                           gTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAAGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGCTATAAATAGATCCATTTTCCTTTTATTTTGa   >  2:431344/1‑141 (MQ=255)

                                                                                                                                            tttACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt  >  1:616633/1‑141 (MQ=32)

TATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATT‑ACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36035‑36314

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