BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,479	T→C	32.4%	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	32.4%	131.7 / 83.8	113	T25A (ACA→GCA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (35/40); new base C (18/18); total (54/59)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 115 total aligned reads displayed.

ATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36342‑36616

ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCATTGtata                                                                                                                                        >  2:54104/1‑141 (MQ=255)

ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtcta                                                                                                                                        >  2:378104/1‑141 (MQ=255)

       tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACAAAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                                 >  1:625013/1‑141 (MQ=25)

        ttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCAGGTTCCCTTAATATTTAAGATATATCACGATAACTCAATGAATATCTTAGaa                                                                                                                                >  2:130693/1‑140 (MQ=12)

          ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                              <  1:100280/141‑1 (MQ=35)

          ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAACTTAGAtag                                                                                                                              >  1:597979/1‑141 (MQ=21)

                    tcttcCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                                                                    >  2:547771/1‑141 (MQ=35)

                     cttcCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                   <  2:663807/141‑1 (MQ=35)

                     cttcCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTACATCTTAGATAGTAGCCAACGGc                                                                                                                   >  1:460684/1‑141 (MQ=25)

                        ccAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAAATAGTAGCCAACggatgc                                                                                                                >  2:478196/1‑137 (MQ=25)

                        ccAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                >  1:197381/1‑141 (MQ=25)

                        ccAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAAATATATCACGATAACTCAATGTATATCTTAGATGGTAGCCAACGGCTAc                                                                                                                >  1:286692/1‑141 (MQ=255)

                            atccgATACAAAATTGGTGCT‑TATTCTTTAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATGTCATATATATCACGATAACTCAATGTATATCTTAGTTAGTCGCCAACGGCTACAGTg                                                                                                            <  1:342458/137‑1 (MQ=11)

                              ttGATACAAAATTNGGGCA‑TATTCTTGAACNCAGCGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTCATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  1:244094/141‑1 (MQ=19)

                              ttGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTATAGTGAt                                                                                                          <  2:368021/141‑1 (MQ=25)

                              ttGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  1:265276/141‑1 (MQ=35)

                              ttGATACAAAATTGGGGAATTGTTCTT‑AACCCAACGGTAGCTAGACCGAACG‑GGAACGTTCACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  1:651257/141‑1 (MQ=11)

                               tGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                         >  1:146886/1‑141 (MQ=35)

                               tGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCATAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                         >  1:465783/1‑141 (MQ=25)

                                 aTACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                       <  2:502788/141‑1 (MQ=35)

                                  tACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                      >  2:563279/1‑141 (MQ=35)

                                  tACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATAATTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCCACGGCTACAGTGCTAACa                                                                                                      >  2:482998/1‑141 (MQ=255)

                                       aaTTGGTGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGAAGACCGTTTTCACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                                                                                 <  2:270466/141‑1 (MQ=255)

                                          tGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTGACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                              <  1:123703/141‑1 (MQ=17)

                                                a‑tatTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTg                                                                                        >  1:263384/1‑141 (MQ=32)

                                                  tatTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGCTAACATCCTTGTTCAATTGt                                                                                       >  1:599675/1‑141 (MQ=17)

                                                     tCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                    <  1:547771/141‑1 (MQ=32)

                                                          aaCCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATGTATCTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTCAATTGTTTATATCt                                                                               >  2:24713/1‑141 (MQ=12)

                                                          aaCCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                               <  1:519093/141‑1 (MQ=32)

                                                            cccAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                             <  1:37194/141‑1 (MQ=32)

                                                                  ggTAGAGCG‑TTGAATGATGAACGTTTACAGCCCGTGCCCTAAATATTTCAGATATATCACGATAACACAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATATCTGAAATAGt                                                                       >  1:489431/1‑141 (MQ=11)

                                                                   gTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                      <  1:481845/141‑1 (MQ=32)

                                                                      gACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATAAAGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAAATCCTTGTTAAATTGTTTATATCTGAACTAGttact                                                                   >  1:136728/1‑138 (MQ=11)

                                                                        ccG‑TTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                 <  1:167627/141‑1 (MQ=17)

                                                                          g‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                               >  2:591781/1‑141 (MQ=32)

                                                                          g‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                               >  1:667046/1‑141 (MQ=32)

                                                                          g‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTCTTTATATCTCAAATAGTTCATACag                                                               <  1:334069/141‑1 (MQ=17)

                                                                               aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                           <  1:355511/141‑1 (MQ=32)

                                                                                     gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                                                     >  1:29860/1‑141 (MQ=32)

                                                                                     gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCTTGTTAAATTGTTTAGATTTGAATTATTTAACAAAGATAACTCCtt                                                     >  2:520407/1‑141 (MQ=11)

                                                                                      aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                    <  2:258131/141‑1 (MQ=32)

                                                                                      aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                    <  2:275691/141‑1 (MQ=32)

                                                                                           ttACACCACGTTCCCTTAATATTTCCGCTATATTACGATCACTCAATGTAGATCTTCGATAGTAGCCAACGGCTACAGGGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACGGAAGACTCCTTTTTGtt                                               <  1:137284/141‑1 (MQ=11)

                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                               <  1:368126/141‑1 (MQ=32)

                                                                                             acacCACGTTCCCTTAATATTTCAGATATATGCCGATAACTCAGTGTATGTCTTAGATAGTAGCCAACGGCTACAGTGATGACAGCCTTGTGAAATTGTTTATATCTGCAATAGTTCATACAGAATCAGCATTTTTGTTaa                                             >  2:282174/1‑141 (MQ=11)

                                                                                               accacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                           >  2:530222/1‑141 (MQ=32)

                                                                                                          ttAATATTTCAGATATATCACGCTAACTCAATGTATAACCTAGATAGTAGCCAACGGCTACAGTGTTAACATCCTTGTCAAATTGTTTATATCTAAAATAGTTCATACAGAAGTGTCGTTTTTGTTAAAATAGTACTATaa                                >  2:48732/1‑141 (MQ=11)

                                                                                                             atctttCAGCCATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTCATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAAtt                             <  1:554685/138‑1 (MQ=1)

                                                                                                              tatTTCAGATATATCACGATAGCTCAACGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTc                            <  1:385854/141‑1 (MQ=255)

                                                                                                               atTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                           <  1:201487/141‑1 (MQ=35)

                                                                                                                tttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCaa                          <  2:496702/141‑1 (MQ=35)

                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                        >  1:174276/1‑141 (MQ=35)

                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                        >  1:461650/1‑141 (MQ=35)

                                                                                                                    ggATCTATCACGATAACACAAGTCATAACGTGGAGAGTAGCCAATAGCTACAGTGATAACAGCCTTGTTAAATTGTTTACATCTGAAATAGGTCATACAGAAGACTCCTTCTTGTTAAAATTATACTATAAATTCAACttt                      <  2:19385/140‑1 (MQ=255)

                                                                                                                     gATATATCACGATAGCTAAGTGTAGATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTg                     >  1:490584/1‑141 (MQ=11)

                                                                                                                        atacgcggctagcTCAATGAATAGCTTAGATAGTAGCCAACGGCTATAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCaa                  <  2:520474/129‑1 (MQ=11)

                                                                                                                         tagcACGATAACTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc                 <  2:183814/138‑1 (MQ=2)

                                                                                                                           tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg               >  1:288255/1‑141 (MQ=35)

                                                                                                                             aCGATAACTCAATGTATATCTTAGATAGTAGACAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTTTTTTTTCTATAAATTCGACTTTTGATCAGaa             >  2:140991/1‑141 (MQ=11)

                                                                                                                               gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc           >  2:550425/1‑141 (MQ=34)

                                                                                                                                aTAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCg          >  1:573116/1‑141 (MQ=35)

                                                                                                                                 tagcccAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:227555/136‑1 (MQ=255)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAATGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:57855/138‑1 (MQ=255)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTGAACTTTGCAACAGAACCGt         >  1:424805/4‑141 (MQ=255)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         >  1:350448/4‑141 (MQ=9)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:135747/138‑1 (MQ=9)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:262110/138‑1 (MQ=9)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         >  1:163003/4‑141 (MQ=9)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGATACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGATATAACACTATACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:126915/138‑1 (MQ=11)

                                                                                                                                 tagcTCAATGCATATCTAAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATCGTTTATATCTGAAATAGTTCCTACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:447164/138‑1 (MQ=11)

                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt         <  1:638325/141‑1 (MQ=35)

                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTCTGCAACAGAACCGt         <  1:356719/141‑1 (MQ=37)

                                                                                                                                 aagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:421961/138‑1 (MQ=12)

                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta        >  1:581353/3‑141 (MQ=12)

                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta        >  2:243286/3‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  1:35693/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  1:607485/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  2:350177/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       <  2:133010/140‑1 (MQ=16)

                                                                                                                                   gcTCAATGCATATCTCAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  1:333509/2‑141 (MQ=12)

                                                                                                                                   aCTCAATGTATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat       <  2:42401/141‑1 (MQ=38)

                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTCAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAACGtat       >  1:486273/1‑141 (MQ=255)

                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat       >  2:349187/1‑141 (MQ=38)

                                                                                                                                    cacAATGCATATCTTAAATAGTCGCCAACGGCTACAGTGATAACACCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAAAGAACCGtatt      >  2:584778/3‑141 (MQ=11)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTATAATTATTTTATAAATTCAACTTTGCAGCAGAACCGtatt      >  1:247940/1‑141 (MQ=11)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  1:509275/1‑141 (MQ=16)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  2:620215/1‑141 (MQ=16)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAACTCAACCTTGCAACAGAACCGtagt      >  1:447866/1‑139 (MQ=11)

                                                                                                                                    atCAATGTATATGTTAGATGGCAGCCAACGGCGACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatt      <  1:24713/140‑1 (MQ=255)

                                                                                                                                     tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt     >  1:77760/1‑140 (MQ=16)

                                                                                                                                      cAATGCATATCTTAGTTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:582522/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGCGAGTAGCCAACGGCTACAGTGATAACATCCGTTTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:166017/141‑3 (MQ=11)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  2:247758/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  2:171153/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:199074/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    >  1:332808/1‑139 (MQ=16)

                                                                                                                                      cAATGCATATCTTAGATAGTACAAAACGTCTGCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    >  1:659988/1‑139 (MQ=11)

                                                                                                                                       aaTGTATATCTTAGATAGTAGCCAACGTCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg   <  2:507576/141‑1 (MQ=255)

                                                                                                                                       aaTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg   <  2:613547/141‑1 (MQ=255)

                                                                                                                                        aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat  <  1:136447/141‑5 (MQ=255)

ATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCA‑TATTCTTGAACCCAACGGTAGACCG‑TTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36342‑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: `‑`