BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,501	G→A	4.4%	A17A (GCC→GCT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,501	0	G	A	4.4%	367.6 / 3.6	148	A17A (GCC→GCT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base G (68/71); new base A (2/5); total (71/77)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 148 total aligned reads displayed.

TCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTAT  >  CP000730/36365‑36640

tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTTGCCAACGGCTa                                                                                                                                         >  1:586296/1‑141 (MQ=25)

     aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAAGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                                                    <  1:580186/141‑1 (MQ=25)

       ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                                                  <  2:506800/141‑1 (MQ=35)

        tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGCTGCCCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGAGAACTCAATG‑TATAGCTTAGATAATAGCCAACGGCTACCGTGATa                                                                                                                                 >  2:692055/1‑141 (MQ=11)

        tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCACCGGCTCCAGTGATa                                                                                                                                 >  1:475285/1‑141 (MQ=21)

        tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                                                 >  1:228455/1‑141 (MQ=35)

          aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                                               >  1:716115/1‑141 (MQ=35)

             cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACTTc                                                                                                                            >  1:656641/1‑141 (MQ=17)

                 aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                                                        <  1:432282/141‑1 (MQ=32)

                 aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAAGG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                                                        <  1:728073/141‑1 (MQ=17)

                  ttcgggaagagtCTTTAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                                       <  1:87006/130‑1 (MQ=11)

                  ttGGGGCATATTCTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                                       <  1:871889/141‑1 (MQ=17)

                  ttGGGACATATTCTTTAACCCAACGCTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGATATATCACGTTGTCCCATTG‑TATGTCTGACATAGTATCCAACGGCTAGAGTGAGAACATCCTTGt                                                                                                                       <  2:40639/141‑1 (MQ=11)

                   tGGGGCATATTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                                                      <  2:777507/141‑1 (MQ=17)

                     gggCATATTCTTGAACCCAACGGCAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                                                    <  1:1035448/141‑1 (MQ=17)

                       gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                                  >  1:548503/1‑141 (MQ=32)

                         atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTg                                                                                                                <  1:188848/141‑1 (MQ=32)

                         atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGtacaaacg                                                                                                                >  2:630935/1‑134 (MQ=255)

                           atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                                              <  1:709239/141‑1 (MQ=32)

                           atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                                              <  1:773464/141‑1 (MQ=32)

                                    cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                                                     >  2:906165/1‑141 (MQ=32)

                                        aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                                                 >  1:414258/1‑141 (MQ=32)

                                           ggtgACCGTTAAATGATTAACGGTTACACCACGTTCCCTTCATATCTCAGATATATCACGACACCTCAATG‑TATATCCTAGATAGTAGCCAACGGCTAAAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                                              <  1:489085/138‑1 (MQ=11)

                                              gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                                           >  1:1031138/1‑141 (MQ=32)

                                                ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                                         <  1:66892/141‑1 (MQ=32)

                                                 cGTTGAATTATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACa                                                                                        <  1:387279/141‑1 (MQ=17)

                                                      aaggatgtACGGTTACACCACGTTCTCTCAGTACTTCAGTGATATCTCGCTAGCTATATG‑TATCGCTTGGATAGTAGCCAACGTCTACAGTGATGTCAGCTTTGTTAAATTGTTTATATCTGAAATAGTTCATGCagaagt                                                                                   <  2:467281/138‑2 (MQ=2)

                                                      aaTGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                                                   <  2:692814/141‑1 (MQ=17)

                                                      aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                                                   <  1:312597/141‑1 (MQ=32)

                                                        tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                                                 >  2:742127/1‑141 (MQ=32)

                                                              aCGTTCACACCACGTTCGCTTAATATTACTGATATATCA‑GCCCAGTTAACGTTAGAGGTTAGTTAATATCGATCTGCTACACTGAGTACCTCCTTGTTAAATTGTTTATATCTGAAATATTTCATACAGAAAACTCCtttt                                                                           >  2:382093/1‑141 (MQ=2)

                                                               cGTTTACACCACGTACCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGGAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGATGAATCCttttt                                                                          >  2:461282/1‑141 (MQ=255)

                                                                  tcacatcgcGTTCCCTTTATATTTCAGATATATCACGATAACTCTATG‑GATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                                       <  2:901211/133‑1 (MQ=1)

                                                                   tACAGCACGTTCCCTTAATATTTCAGATATGTCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGttc                                                                      >  2:646080/1‑140 (MQ=17)

                                                                      accccGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                                                   <  2:581429/141‑1 (MQ=17)

                                                                       ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaag                                                                  >  2:645991/1‑140 (MQ=17)

                                                                       ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGATACAGTGATAACTTCATTGTTACATTGTTTATATATGGGTTAGTTCAGGCGGAATAGACGTTATTGTTAAAAt                                                                  >  2:728509/1‑141 (MQ=11)

                                                                        cacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAtt                                                                 <  1:1003147/141‑1 (MQ=32)

                                                                          cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttttttttttttttt                                                               >  2:496266/1‑133 (MQ=255)

                                                                            ttgccTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                                                             <  2:296580/138‑1 (MQ=17)

                                                                              cccTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACta                                                           <  2:985869/141‑1 (MQ=32)

                                                                                 ttAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                                                        >  2:955833/1‑141 (MQ=33)

                                                                                          cAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                                               >  1:344295/1‑141 (MQ=35)

                                                                                              tatatCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                                           >  2:880522/1‑141 (MQ=35)

                                                                                               atatCACGATAGCTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCaa                                          >  1:156417/1‑141 (MQ=25)

                                                                                                   cACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa                                      >  2:289140/1‑141 (MQ=34)

                                                                                                      gATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc                                   <  1:418981/141‑1 (MQ=34)

                                                                                                      gATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc                                   <  1:449144/141‑1 (MQ=35)

                                                                                                        tagcTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATCGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTCTAACTTCAACTTTGCAACGGAACCGt                                 <  2:619899/138‑1 (MQ=11)

                                                                                                        tagcTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                 <  1:520365/138‑1 (MQ=9)

                                                                                                        tAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt                                 <  2:441751/141‑1 (MQ=37)

                                                                                                          gcTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGCAATAGGTCATACCGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtct                               >  2:519792/2‑139 (MQ=11)

                                                                                                          gcTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAACCGtat                               >  2:345531/2‑141 (MQ=12)

                                                                                                          gcTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                               <  1:1021151/140‑1 (MQ=16)

                                                                                                          gcTCAATG‑CATATCTTAGATAGTAACCAAAGGCTACAGTGATAACATCCTTGTTAAATTGTATATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATATTGTAACTTCAACTTTGCGACAGAACCGtat                               >  2:569962/2‑141 (MQ=11)

                                                                                                           cTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                              >  2:111099/1‑141 (MQ=16)

                                                                                                           cTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                              >  2:956179/1‑141 (MQ=16)

                                                                                                            tCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta                             >  1:1011714/1‑141 (MQ=40)

                                                                                                            tCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAACCGtattt                             <  2:1052761/141‑2 (MQ=255)

                                                                                                            tCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                             <  1:106745/141‑2 (MQ=16)

                                                                                                            tCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                             <  1:223217/141‑2 (MQ=12)

                                                                                                            tCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                             <  1:107569/141‑2 (MQ=19)

                                                                                                            tCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatgt                             >  1:47710/1‑139 (MQ=12)

                                                                                                            tCAATG‑CATATCTTAAATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                             >  1:787303/1‑140 (MQ=255)

                                                                                                            cacagtgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctccta                              >  2:535934/9‑136 (MQ=11)

                                                                                                             cAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt                            <  1:52945/141‑3 (MQ=255)

                                                                                                             cAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt                            <  2:192865/141‑3 (MQ=255)

                                                                                                              aaTG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGTCTTTTg                           >  1:1033368/1‑141 (MQ=255)

                                                                                                                agacatatCTTAGATCGAAGCCAACACCTAAAGTGATAACATCCAACTTCAATAGCTTCTCTCTAACCTACTTCCTACAGAAGACTCATTTCTGTTACCATTATACTATAAATTCAACTTTGCAACAGAAcctcatttact                          >  2:129945/5‑132 (MQ=2)

                                                                                                                 tgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                         >  1:269761/4‑136 (MQ=255)

                                                                                                                  gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                        >  2:456595/3‑135 (MQ=255)

                                                                                                                   catatCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                       <  1:204202/140‑13 (MQ=11)

                                                                                                                   catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                       <  1:576011/140‑13 (MQ=11)

                                                                                                                    atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtttt                      >  2:989078/1‑128 (MQ=11)

                                                                                                                    atatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTGTTGTTAACATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatat                      <  2:436535/141‑14 (MQ=11)

                                                                                                                     tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAACTTGTATATATCTGAAATCGTCCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCgccatctatgagta                     >  2:874529/1‑128 (MQ=11)

                                                                                                                     tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGCAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCATATTATGGAATag                     >  2:598639/1‑141 (MQ=255)

                                                                                                                     tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTTATTTTATTATAAATTCAACTTTCCAACAGAAcctcctaactgtatta                     >  2:214881/1‑127 (MQ=11)

                                                                                                                     tatCGTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                     <  1:862672/141‑15 (MQ=11)

                                                                                                                     tatCGTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGCTATTATTTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                     >  1:723699/1‑127 (MQ=11)

                                                                                                                      atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtac                    >  2:229069/1‑131 (MQ=11)

                                                                                                                      atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaa                    >  2:661945/1‑131 (MQ=255)

                                                                                                                      atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTAATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCACCTTTGAACCAGAACAGcattttataagtac                    >  2:825977/1‑131 (MQ=11)

                                                                                                                        cTTAGTTAGTAGCAAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacac                  <  1:919351/141‑13 (MQ=11)

                                                                                                                        cTTAGATAGTAGCCGACGGGTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtaatgtatgggtacac                  >  1:583934/1‑127 (MQ=11)

                                                                                                                        cTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAgaga                  <  2:238954/141‑1 (MQ=255)

                                                                                                                        cTTAGATAGTAGCCAACGGCTACAGTGATAACATACTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTTTAAATTCAACTTTGCAACAGAACCGTATTATGTAGTagata                  >  2:828521/1‑139 (MQ=255)

                                                                                                                         ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                 <  1:27082/141‑14 (MQ=11)

                                                                                                                         ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                 >  1:851587/1‑128 (MQ=11)

                                                                                                                         ttAGATAGTAGCCAACGGCTACAGTGACAACAGCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACCGCACCGTATTATGGAATAGAGAt                 <  2:48523/141‑1 (MQ=255)

                                                                                                                          tAGATTGTATGCAACGGCTACAGTGATAACATCCCTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                <  2:750601/141‑15 (MQ=11)

                                                                                                                          tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                <  1:278204/141‑15 (MQ=11)

                                                                                                                          tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAATTTTTTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                <  1:205135/141‑15 (MQ=11)

                                                                                                                            gATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGtt              <  2:548503/141‑1 (MQ=255)

                                                                                                                               agtagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGCCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGt           >  2:377102/1‑141 (MQ=255)

                                                                                                                                gtagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTa          >  1:955216/1‑141 (MQ=255)

                                                                                                                                       aCGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTTTTATGGAATAGAGAGGTTGGTAACATTta   <  1:566526/141‑1 (MQ=255)

                                                                                                                                       aCGGCTACAGTGATAACATCCTTGTTAAACTGTTTATATCTGAAATAGTTCCTACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTta   <  2:562022/141‑1 (MQ=255)

                                                                                                                                        cGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtat  <  2:663039/141‑1 (MQ=255)

                                                                                                                                        cGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTtat  <  1:1029555/141‑1 (MQ=255)

TCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTAT  >  CP000730/36365‑36640

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