BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,479	T→C	22.9%	T25A (ACA→GCA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,479	0	T	C	22.9%	176.5 / 58.7	110	T25A (ACA→GCA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (46/38); new base C (10/15); total (56/54)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 112 total aligned reads displayed.

AATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36341‑36616

aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTgaatg‑aat                                                                                                                                         >  1:693710/1‑141 (MQ=11)

       ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAg                                                                                                                                  >  1:661431/1‑141 (MQ=35)

       ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATCACGATAACTCAATG‑TATATCTTaa                                                                                                                                  >  1:321067/1‑140 (MQ=25)

             aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAACTCAATG‑TATATCTTAGAtagta                                                                                                                            >  1:191467/1‑141 (MQ=25)

               gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGc                                                                                                                          >  2:1053063/1‑141 (MQ=35)

                  ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCaa                                                                                                                       <  2:1004698/141‑1 (MQ=35)

                   tttCCTCCAAATTTGATACAAAATTGGTGCTTATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGccgtt                                                                                                                      <  1:35865/141‑4 (MQ=255)

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

                        tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTa                                                                                                                 >  1:1025172/1‑141 (MQ=35)

                             aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGTTAGACGGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                            <  2:515535/141‑1 (MQ=21)

                             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)

                                         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)

                                           tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                              <  1:147010/141‑1 (MQ=32)

                                             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)

                                                         gagccgaaCGGTAGACCGTTGAATGATGAACGTTTACACGACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                <  2:894835/135‑1 (MQ=255)

                                                            cccAACTGTAGACCGTTGAATGATTAACGTTTACACCACGTTCCCTTAATATTTCAGATATAGCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                             <  2:923171/141‑1 (MQ=12)

                                                            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)

                                                                        ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAAGGCAATC‑TCCATTTTATATAGTAGCCAACGGCTACAGTGATAACATCCGTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                 >  1:508680/1‑141 (MQ=11)

                                                                        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)

                                                                                         tttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                <  2:847029/141‑1 (MQ=32)

                                                                                          ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAACCGTTCAGACGGAAGACTCCTTTTTGtt                                               <  1:1053063/141‑1 (MQ=255)

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

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

                                                                                           tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                              >  2:449144/1‑141 (MQ=32)

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

                                                                                              accacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                           >  2:284929/1‑141 (MQ=32)

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

                                                                                               ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                          >  2:569891/1‑141 (MQ=32)

                                                                                               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)

                                                                                                      cccTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTGTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAATTTTTTACtt                                   >  2:703258/1‑140 (MQ=11)

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

                                                                                                                 tCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                        >  2:379504/1‑141 (MQ=35)

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

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

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

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

                                                                                                                      aatatTACGATAACTCAATG‑TATATCTTAGATAGTAGACAACGGCTAGAGTAATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                   <  1:248578/140‑1 (MQ=255)

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

                                                                                                                          tCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCACCAg               >  2:644793/1‑141 (MQ=25)

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

                                                                                                                              gATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATAAGACTCTTTTTTGTTAAAATTATACTATAATTTCAACTTTGCAACAGAAcc           <  1:646080/141‑1 (MQ=255)

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

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

                                                                                                                               aTAGCTCAATG‑CATATCTTAGATAGTAGCCAACGGCTACAGTGAGAACAGCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGACGCCTCCGTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg          <  1:205571/141‑1 (MQ=11)

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

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

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

                                                                                                                                  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‑CATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAACCGtatg      <  2:650714/141‑2 (MQ=255)

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

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

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

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

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

                                                                                                                                    tCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta     <  1:906165/141‑1 (MQ=40)

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

                                                                                                                                    tCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGAAACAGAACCGtatta     >  1:810811/1‑141 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                                                                                       aTG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAAAGAACCGTATTATgg  >  2:463064/1‑141 (MQ=255)

AATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36341‑36616

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