BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,174	T→C	31.4%	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	31.4%	191.8 / 165.9	213	K126K (AAA→AAG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (54/91); new base C (32/35); total (87/126)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

Only 100 of 217 total aligned reads displayed.

TGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36037‑36314

tGACGGTGATCTTGCTCAATGAGGTTATACAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTa                                                                                                                                           >  1:211029/1‑141 (MQ=11)

 gatggTGATCTTGCTAAAGGTGGTTTTTCTGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAATGTGCCTGATCCGTAACTACCTTTTGAGGCTTAc                                                                                                                                          <  1:161808/138‑1 (MQ=11)

 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                                          <  2:185460/141‑1 (MQ=21)

 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCAATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTCAc                                                                                                                                          <  2:255820/141‑1 (MQ=255)

  gcGGTGATCTTGCTCAATGAGGTTATTTAGTTATTTCGATGTACAATGGCAGTCACGTGTAAGTTTAATAGCTTTAATTTTTTGCATCATTGCTACCTTCGTTGAAGGTGCTTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  2:201432/140‑1 (MQ=11)

  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  1:285519/1‑141 (MQ=25)

  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  2:42179/141‑1 (MQ=21)

  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  2:126194/1‑141 (MQ=255)

    ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaa                                                                                                                                       >  1:166910/1‑141 (MQ=18)

     gagATCGTTCTCAATGATGTTATTCAGCTATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGTTTTACCaaa                                                                                                                                      <  1:41541/139‑1 (MQ=11)

     gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                                      >  2:173201/1‑141 (MQ=17)

         tCTTGCTCAATGAGGTTACTCCTGTAGTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  2:322887/1‑141 (MQ=11)

           ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGttt                                                                                                                                >  1:72429/1‑141 (MQ=255)

             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                              >  1:143608/1‑141 (MQ=255)

              cTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAAt                                                                                                                             >  2:283287/1‑141 (MQ=11)

               tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGGGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                                            >  1:168292/1‑141 (MQ=12)

                 aaTGAGGTTATTCAGCTATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                                                                          >  1:247914/1‑141 (MQ=255)

                  aTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgaga                                                                                                                         >  1:109228/1‑141 (MQ=32)

                   tGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACTAACTTGTTTAATGAGAc                                                                                                                        >  2:301521/1‑141 (MQ=255)

                     aGGTTATTCAGATATTTCGATGTACAATGCCGGTGAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                                      >  1:297244/1‑141 (MQ=255)

                       gTTATTCATAAATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  1:18553/141‑1 (MQ=255)

                       gTTATTCAGATATTTGGATGTACTATGGCAGTCAGTTTTAAGTTTAAAAGCTTTAATTACTTTATCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  2:254785/141‑1 (MQ=11)

                       gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                                    >  1:107183/1‑141 (MQ=17)

                       gTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTCCGTTGAAGGTGCCTTATCTTTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGttt                                                                                                                    <  2:152037/141‑1 (MQ=11)

                         tactCATATAATTAGATGTACGATGACAGTTAGGTTTAAGGTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  2:197584/138‑1 (MQ=11)

                         tATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAACAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  1:190284/141‑1 (MQ=12)

                         tATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  1:52750/141‑1 (MQ=12)

                         tATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  >  1:280509/1‑141 (MQ=17)

                           ttCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                                >  2:60567/1‑141 (MQ=11)

                             cAGATATTTCGATGTACCCTGGCAGTCAGGTTTAAGTTTGGCGGCTTTGATGCCATTAAGAATTGCTGCCTGGGGTGAATGTTTCTGATCTGTAAGTAGCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              <  2:277969/141‑1 (MQ=2)

                             cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              <  1:125785/141‑1 (MQ=32)

                              aGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                                                                             <  1:266861/141‑1 (MQ=17)

                                atatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                                           <  1:3961/141‑1 (MQ=255)

                                 tatTTTGATGTATAATGGCAGTCTGGTTTAAGTTTAAAAGCTTTAATTACTTTTACCATTGCTACCTTCGTTGAAGGTGCTTGATCTGTATTTACCTTTTGAGGTTTCCCAAACTTTTCGATGAGGCGTTTGATAAACGCa                                                                                                          <  2:72046/141‑1 (MQ=11)

                                 tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  1:59329/141‑1 (MQ=255)

                                  atTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         <  1:257062/141‑1 (MQ=17)

                                  atTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                                         <  1:169761/141‑1 (MQ=255)

                                    ttCTATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                                       <  2:9152/141‑1 (MQ=11)

                                      cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     >  1:182092/1‑141 (MQ=17)

                                      cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                                     <  2:247544/141‑1 (MQ=32)

                                         tGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                                  >  1:77801/1‑141 (MQ=12)

                                          gTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGccgct                                                                                                 >  2:81741/1‑137 (MQ=11)

                                           tACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAt                                                                                                <  2:237640/141‑1 (MQ=17)

                                            aCAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCATTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                               >  2:267872/1‑141 (MQ=12)

                                             cAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                              <  1:67638/141‑1 (MQ=17)

                                             cAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGGTTACCAAATTGTTTATTGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                              >  1:288173/1‑141 (MQ=255)

                                              aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                             >  1:1887/1‑141 (MQ=14)

                                               aTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAtt                                                                                            <  2:312530/141‑1 (MQ=11)

                                                tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           <  2:169799/141‑1 (MQ=31)

                                                 gACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTAt                                                                                          <  1:153061/141‑1 (MQ=11)

                                                                    aacaGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGCTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCa                                                                       <  2:214394/138‑1 (MQ=0)

                                                                                   tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        >  1:175105/1‑141 (MQ=2)

                                                                                     taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                                      <  2:205712/138‑1 (MQ=17)

                                                                                         cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  2:324879/141‑1 (MQ=32)

                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                                 <  1:10054/141‑1 (MQ=32)

                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                                 <  2:25642/141‑1 (MQ=32)

                                                                                              cTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGATGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             <  1:142533/141‑4 (MQ=17)

                                                                                               tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                                            <  2:334687/141‑1 (MQ=17)

                                                                                               tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                                            <  2:297244/141‑1 (MQ=17)

                                                                                                gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  1:321185/140‑1 (MQ=255)

                                                                                                gccTTCGGTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  1:196145/140‑1 (MQ=255)

                                                                                                aCCTTCGTTGAAGGTGCCTGATGTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  1:231297/141‑1 (MQ=12)

                                                                                                 ccgtctttcaaGGTGCCGGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:63394/132‑1 (MQ=11)

                                                                                                 ccTTCGTTGAAGGTGGCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:224828/141‑1 (MQ=17)

                                                                                                  cTTCGTTGAAGGTGCATGATCTGTAAGTACCTTTGGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGCTGCTTACGCAACCAAATATCTAACGTATGTCCCTCTGCATCAAt                                         <  1:202218/141‑1 (MQ=11)

                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  1:304484/141‑1 (MQ=17)

                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                                        <  1:315650/141‑1 (MQ=32)

                                                                                                       ttGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                                    >  1:253605/1‑141 (MQ=12)

                                                                                                         taaGGTGCCTGAGACGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                                  <  1:100829/140‑1 (MQ=12)

                                                                                                         gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGa                                  <  1:324005/141‑1 (MQ=17)

                                                                                                          aaGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGat                                 >  1:203965/1‑141 (MQ=12)

                                                                                                           agttGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGACACGata                                <  2:263008/138‑1 (MQ=17)

                                                                                                           aGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGgta                                >  2:167456/1‑141 (MQ=17)

                                                                                                            ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGatat                               >  1:11741/1‑141 (MQ=32)

                                                                                                             gTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAtata                              >  1:45367/1‑141 (MQ=32)

                                                                                                                 ccgATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATa                          <  2:199076/139‑1 (MQ=17)

                                                                                                                 cTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTCCGCAACCTCATATCGAATGTATGTCCCTCTGTATCAATGGCACGATATAAATa                          <  2:321413/141‑1 (MQ=11)

                                                                                                                  tGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAg                         >  1:173125/1‑141 (MQ=255)

                                                                                                                   gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  1:317900/141‑1 (MQ=255)

                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  2:154675/141‑1 (MQ=17)

                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCACTGGCACGATATAAATAGcc                       >  2:160589/1‑140 (MQ=17)

                                                                                                                     tCTGTAATTACCTTTTGAGGTTTATCCAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTc                      <  1:62123/141‑1 (MQ=17)

                                                                                                                       tGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa                    >  1:205685/1‑141 (MQ=32)

                                                                                                                       tGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACAAAATATGTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa                    <  2:219076/141‑1 (MQ=17)

                                                                                                                       tGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGGATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCGTCTGTATCAATGTCACGATATAAATAGCTcaa                    >  1:222863/1‑139 (MQ=11)

                                                                                                                        tgaactACCTTTTGGCGCTTACCAAATTGTTTAATGACACATTCGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATGTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAt                   <  1:235536/136‑1 (MQ=11)

                                                                                                                         tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAtt                  >  2:240746/1‑141 (MQ=32)

                                                                                                                         tAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTACCTCTGCATCAATGGCACGATCTAAATGGATCCAtt                  >  2:258564/1‑141 (MQ=255)

                                                                                                                           actACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                <  2:107022/139‑1 (MQ=12)

                                                                                                                            ctcccTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATACCTAATGTATGTCACTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc               <  2:117726/138‑1 (MQ=11)

                                                                                                                            ctACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc               <  1:63283/140‑1 (MQ=12)

                                                                                                                              aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             >  2:152545/1‑141 (MQ=17)

                                                                                                                              aCCTTTTGAGGCTTACCAAATTGTTTAATGATACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:1357/141‑1 (MQ=12)

                                                                                                                              aCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATGTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:96914/141‑1 (MQ=12)

                                                                                                                                cctttGAGGCTTACCAAATTGTTTAATGAGACTGTTGATAAACGCATATCCCGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  1:9294/139‑1 (MQ=11)

                                                                                                                                cTTTTGAGGCTTACCAAATAGTTTAATGAGACGTTTGATAAGCGCATATGCTGAATGTTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATATATAGCACCATTTTCCttt           >  1:200078/1‑141 (MQ=11)

                                                                                                                                 ttttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCAAATGTATGTCCCTCTGCATCAATGGCACGATATAAATCGCCACATTTTCCtttt          >  2:43220/1‑141 (MQ=255)

                                                                                                                                  tttGAGGTTTACCAACTTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTTCATCAATGGCACGATATAAATAGCTCCATTTTCCttttt         >  1:34856/1‑140 (MQ=17)

                                                                                                                                     gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCGGTCCTTTTCATGGGCACGATATAAATAGCTCCATTTTCCTTTTAttt      >  1:225880/1‑141 (MQ=11)

                                                                                                                                     gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttt      >  2:304574/1‑141 (MQ=17)

TGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36037‑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: `‑`