BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,174	T→C	34.0%	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.0%	418.9 / inf	595	K126K (AAA→AAG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (190/199); new base C (100/104); total (292/303)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 607 total aligned reads displayed.

ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑36314

    gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGAGAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                                           >  1:512959/1‑141 (MQ=11)

    gACGGTCATCCTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCGGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                                           <  1:753633/141‑1 (MQ=11)

    gACGGGGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCACGTTTAAGTTTTACAGCTTTAATTACTGTTTCCATTTCTACCTTCGTTGAAGGTGCCTTTTCTGTAATTACCTTTTGAGGTTTAc                                                                                                                                           <  1:145932/141‑1 (MQ=11)

     aCGGTGATCTTACTCAATGAGGTTATTCAGATATATCGATGTACAATGACAGGCAGGATTAAGTTTAACAGTTTTAATTATTTTAGCCATAGCTGACTTCGTTGCAGGTGACTGAACTGTAACTACCTATTGAGGCTGAcc                                                                                                                                          >  2:210443/1‑141 (MQ=11)

        gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                                       <  2:100378/141‑1 (MQ=17)

         tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                                      >  1:297758/1‑141 (MQ=32)

          gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAtt                                                                                                                                     <  2:282232/141‑1 (MQ=255)

           aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                                                                    <  2:678660/141‑1 (MQ=17)

           aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCCTACCAAATTg                                                                                                                                    >  1:713851/1‑141 (MQ=11)

            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                   >  2:369861/1‑141 (MQ=12)

            tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGt                                                                                                                                   >  2:667415/1‑141 (MQ=255)

             cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:392852/1‑141 (MQ=255)

             cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:634321/1‑141 (MQ=11)

                gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                               >  2:298137/1‑141 (MQ=32)

                 gtCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAAt                                                                                                                              <  1:549/140‑1 (MQ=11)

                 cTCAATGAGGTTATTCAGATATTTAGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGATGGTGCCTTATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                              <  2:484216/141‑1 (MQ=11)

                  ttaaTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                                             <  1:234963/139‑1 (MQ=12)

                   cAATGAGGTTATTCCGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaaata                                                                                                                            <  2:710864/141‑4 (MQ=11)

                   cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                                            <  1:134906/141‑1 (MQ=255)

                        gggTTATTTAGATTTTTTGATGTACATTGGCTGGCAGTTTTATGTTTAAAAGGTTTAATTGCTTTATCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                       <  2:408221/140‑1 (MQ=11)

                         ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                      <  1:661625/141‑1 (MQ=17)

                         ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                                      >  2:244132/1‑141 (MQ=11)

                         ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGtt                                                                                                                      >  1:274442/1‑141 (MQ=11)

                          gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCCTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                     <  1:70107/141‑1 (MQ=12)

                          gTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGttt                                                                                                                     >  1:484295/1‑141 (MQ=255)

                           ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                    >  2:533378/1‑141 (MQ=11)

                            tATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                   <  1:345814/141‑1 (MQ=11)

                             aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAGCAAATTGTTTAATGAGACGTTTGAt                                                                                                                  >  1:617508/1‑141 (MQ=12)

                             aTTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGGTTTAATTGTTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                                  >  1:602676/1‑141 (MQ=11)

                                cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑aaa                                                                                                              >  2:134363/1‑141 (MQ=12)

                                cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑aaa                                                                                                              >  2:328882/1‑141 (MQ=12)

                                 aGATATTTCCATGGACAATGGCAGTTAGATTATAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTCTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAAc                                                                                                             <  1:692294/141‑1 (MQ=11)

                                  gttatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACg                                                                                                            <  2:569599/139‑1 (MQ=12)

                                    tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCa                                                                                                          <  1:770659/141‑1 (MQ=11)

                                       ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCctat                                                                                                       >  1:89789/1‑141 (MQ=12)

                                         cGATGTACAATGGCAGTCAGGATTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGc                                                                                                     >  1:401648/1‑141 (MQ=255)

                                          gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCt                                                                                                    >  2:620335/1‑141 (MQ=12)

                                          gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCCTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCt                                                                                                    >  2:272011/1‑141 (MQ=12)

                                             gTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGaa                                                                                                 <  1:415137/141‑1 (MQ=17)

                                               aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTCCCCTTTGGGGTGTACTCTATATTGTAATGAGACGTTTGAT‑AACCGCACTTGCTGAATg                                                                                               <  1:566500/141‑1 (MQ=11)

                                               aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATg                                                                                               <  1:249819/141‑1 (MQ=17)

                                                cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGa                                                                                              >  2:294514/1‑141 (MQ=32)

                                                 aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGAt                                                                                             <  1:104700/141‑1 (MQ=14)

                                                  aTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTGGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGCCTTACCAACTTGTTTAATGAGACGTTTGAT‑AAACGCGTATGCTGAATGAtt                                                                                            >  2:451907/1‑141 (MQ=11)

                                                    gACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTAt                                                                                          <  1:20975/141‑1 (MQ=11)

                                                                                      tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        >  1:18233/1‑141 (MQ=2)

                                                                                        tAGCCCTTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACAACATTTTGAGGCTTACCAAATTGTTTAATGAGACGTCTGAT‑AAACGAATATGCTGAATGAGTATCTCGTTGCTTACGCACCCAAATGTCTAATGTATGTc                                                      <  2:642025/141‑1 (MQ=11)

                                                                                         aaccATTGCGACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAAGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AGACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  1:473820/139‑1 (MQ=255)

                                                                                         aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  2:42139/141‑1 (MQ=12)

                                                                                              ttGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATATATGTCCctctt                                                >  1:645667/1‑140 (MQ=17)

                                                                                                   gccTTCTTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  2:655159/140‑1 (MQ=11)

                                                                                                   gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  2:634321/140‑1 (MQ=255)

                                                                                                   aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑CATCACTGAAACGCAATGATTATCTCGTTGCTTACGCAACCAAATATATAATGTATGTCCCTCTGCATCa                                           >  1:761582/1‑141 (MQ=11)

                                                                                                   aCCTTCGTGGAAGGTGCCTGTTGTGTAATGACCTTTTGAGGTTTACCAAATGGTTTAAGGAGCCGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                                           <  1:682323/141‑1 (MQ=11)

                                                                                                    cctgctctaaatgTGCCTGTTCTTAAATTACCTTTTGAGTTTTACCAAATTGTTTAATGAGACGTTTGAT‑AACCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  1:193515/129‑1 (MQ=11)

                                                                                                    ccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCaa                                          <  2:461510/141‑1 (MQ=32)

                                                                                                     cTTCGTTGAAGGTGCCTGTTCTGTAACTACCTTTTGAGGCTTCCCAGATTGTCTATTGAGTCGTTTGGT‑GAACGCGTATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         <  1:221671/141‑1 (MQ=11)

                                                                                                     cTTCGTTGAAGGTGCCTGATCTGTAACTCCCTTTTGAGGCGTACTAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCATCTGTATCAAt                                         >  1:730934/1‑141 (MQ=11)

                                                                                                      ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGTT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAATTATCTAATGTATGTCCCTCTGCATCAATg                                        >  1:340933/1‑141 (MQ=17)

                                                                                                      ttCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  2:223041/141‑1 (MQ=255)

                                                                                                         gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCACTGGCa                                     >  2:400619/1‑141 (MQ=17)

                                                                                                          tttaaGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTTTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                                    <  1:157212/138‑1 (MQ=11)

                                                                                                           tGAAGGTGTCGGATCTGTAATTTCCTTTTGAGGTTTACCGAATTGTTTCATGAGACGTTTGCT‑AAGCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACg                                   <  1:113942/141‑1 (MQ=11)

                                                                                                           tGAAGGGGCCTGACCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACg                                   <  2:148591/141‑1 (MQ=11)

                                                                                                            gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                                  <  2:353193/141‑1 (MQ=32)

                                                                                                             aaGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCACCCACATATCTCATGTATGTACCTCTGCAACAATAGCagcac                                 >  1:540632/1‑137 (MQ=11)

                                                                                                               ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatat                               >  2:561109/1‑141 (MQ=32)

                                                                                                                gggCCTGAGCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAtata                              <  1:361354/139‑1 (MQ=11)

                                                                                                                gTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATCGCACGAtata                              <  2:168833/141‑1 (MQ=17)

                                                                                                                  gCCTGATCTGTAATAACCTATTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATaaa                            >  1:425486/1‑141 (MQ=17)

                                                                                                                   ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCCAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAAt                           >  2:118240/1‑141 (MQ=17)

                                                                                                                   ccTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATAGCTAATGTATGTCCCTCGGTATCAATGGCACGATATaaac                           >  2:181663/1‑140 (MQ=12)

                                                                                                                   ccTGATCTGTAACTACATTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                           >  1:385217/1‑141 (MQ=255)

                                                                                                                    cataTCTGTAACTATCTTTTGAGGCTTACCAAATTTTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATa                          <  2:473544/138‑1 (MQ=11)

                                                                                                                     tGATCTGTAATTACCTTTTGAGGTTTACCACATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAg                         >  1:497936/1‑141 (MQ=17)

                                                                                                                     tGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATACAATGTATGTCCCTCTGTATCAATGGCACGAAATAAATAg                         >  2:393317/1‑141 (MQ=11)

                                                                                                                     tGATCTGTAACTACATTTTGAGGCTTACCAAATTGTTTAATGAGACGTATGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCACATAACTAATGTATGTACCACTGTATCAATGGAACGATGTaaaatg                         >  2:589721/1‑138 (MQ=11)

                                                                                                                      gATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGACACGATATAAATaac                        >  1:317097/1‑139 (MQ=12)

                                                                                                                      gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  1:675970/141‑1 (MQ=12)

                                                                                                                       aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGCCGTTTGCTCTTACGC‑TCAACTTAATCATTCTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt                       <  2:545221/141‑1 (MQ=11)

                                                                                                                       aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  2:758951/141‑1 (MQ=17)

                                                                                                                        tCTGTAATTAGCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTc                      <  2:68863/141‑1 (MQ=17)

                                                                                                                         cTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTCACGCAACCAAAAACCTAACGTATGTCCCTCTGTATCAATGGCACGACATAAATAGCAcc                     >  2:80251/1‑141 (MQ=11)

                                                                                                                           gTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAt                   <  2:344243/141‑1 (MQ=17)

                                                                                                                            tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtt                  <  1:287666/141‑1 (MQ=17)

                                                                                                                             aactACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAttt                 <  2:746277/138‑1 (MQ=12)

                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTGCCTCTGCATCAATGGCACGATATAAATCGCTCCAtttt                >  2:608841/1‑141 (MQ=17)

                                                                                                                               ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTc               >  1:261466/1‑141 (MQ=32)

                                                                                                                               ctACCGTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAC‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc               <  2:501487/140‑1 (MQ=255)

                                                                                                                                tACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              <  2:726396/141‑1 (MQ=17)

                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAGCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATAGCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:669933/141‑1 (MQ=12)

                                                                                                                                   cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑ACACGCATATGCTTGATGTTTATCTCGTTGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGCTATAAATAGCTCCATTTTCCttt           <  2:382814/141‑1 (MQ=11)

                                                                                                                                   cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           >  2:430388/1‑141 (MQ=17)

                                                                                                                                   cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  1:445327/141‑1 (MQ=17)

                                                                                                                                    ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCtttt          >  2:52010/1‑141 (MQ=32)

                                                                                                                                     tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATTTAAATACCAACATTTTTCttttt         >  1:522290/1‑140 (MQ=11)

                                                                                                                                     tttGAGGCTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACACATATGCTGAATGGTTTTCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttttt         >  1:369911/1‑140 (MQ=11)

                                                                                                                                       tGAGGTTTACCAAATTGTTTAATGAGACGGTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATAACTAATGTATGTCACTCTGCATCAATGGCACGAGATAAAAAGCTCCACTTTCCTTTTAtt       >  2:308156/1‑141 (MQ=11)

                                                                                                                                        gAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGACACGATATAAATAGCTCCATTTTCCTTTTAttt      >  2:249638/1‑141 (MQ=17)

                                                                                                                                            tttACCAAATTGTTTAATGAGGCGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAt  <  1:585556/141‑1 (MQ=17)

ATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAT‑AAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36034‑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: `‑`