BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,479	0	T	C	26.2%	107.9 / 46.3	86	T25A (ACA→GCA)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (34/27); new base C (9/13); total (43/43)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618

taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtt                                                                                                                                             >  2:535678/1‑140 (MQ=255)

 aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCATTGtaa                                                                                                                                            >  2:268896/1‑140 (MQ=255)

 aataaGCTTCTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtat                                                                                                                                            >  1:577815/1‑141 (MQ=255)

     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAGCGTTTACACCACGTTCCCTTAATATTTCAGATATATCATGATAACTCAATGTATATCt                                                                                                                                        <  1:473852/141‑1 (MQ=255)

     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                        <  2:533888/141‑1 (MQ=38)

            ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACAGCACGTTCACTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                                 >  1:582588/1‑141 (MQ=25)

            ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                                 >  2:337639/1‑141 (MQ=35)

             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:467483/1‑141 (MQ=35)

             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:96053/1‑141 (MQ=35)

             tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACAGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:512139/1‑141 (MQ=25)

              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                               >  1:94844/1‑141 (MQ=35)

              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                               >  2:332236/1‑141 (MQ=35)

                gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                             >  1:564426/1‑141 (MQ=35)

                  tcttgcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCa                                                                                                                           <  1:241332/136‑1 (MQ=25)

                          ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                   <  2:49111/141‑1 (MQ=21)

                          ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTTTATCTTAGATAGTAGCCAACGGCTAc                                                                                                                   >  1:599114/1‑141 (MQ=21)

                          ccAAAGTTGATACAGAGTGTGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                   <  2:123005/141‑1 (MQ=11)

                           ttttttttATACAAAATTGGGGCATAATCTTGTACCCAACGGTAGACCGTTGAATGGTGAACGTTTACACCACTTTTTTTTATTATTTCAGATATATGATGATAACTCAATTTATATTTTAGATAGTAGCCAACgggtgaa                                                                                                                  >  2:495916/5‑136 (MQ=11)

                                  gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGAAGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAGGGCTACAGTgctca                                                                                                           >  1:521911/1‑137 (MQ=255)

                                         aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATGCtt                                                                                                   >  2:313297/1‑141 (MQ=17)

                                           ttGGGGCATATTCTTGAACCCAACGGTATACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGt                                                                                                 <  1:447623/141‑1 (MQ=17)

                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGt                                                                                                 <  1:85203/141‑1 (MQ=32)

                                           ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCCACAGTG‑ATAACATCCTTGt                                                                                                 <  1:55735/141‑1 (MQ=17)

                                           tgggggCATATTCTTGAACCCACCCGTAGACCGTTGAATGATGAACCTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGt                                                                                                 <  1:162038/139‑1 (MQ=12)

                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTaa                                                                                              >  1:320265/1‑141 (MQ=32)

                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATTTCAGATATATCACGATAACTCAATTTAGATCTTAGATAGTAGCCAACGGCTACAGTG‑CTACGATCCTTGTTaa                                                                                              >  1:22184/1‑141 (MQ=11)

                                               ggCCTATTCTTGAACCCAACGGTAGCCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATACCACGATAACTCAATGTATGTCTTAGATAGTAGCCACCGGCCACAGTG‑ATAACCTCCCTGTTaca                                                                                             >  2:92202/1‑139 (MQ=11)

                                               ggCAAATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGGTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGCTaaa                                                                                             >  2:545530/1‑141 (MQ=12)

                                               agCATATTCTTGAACCCAACGGTAGACAGTTGAATGATGAACGTTTACACCACGATCCCTTCATATTTCAGAAATATCACGATAACTCACTGTAAATCTTAGATCATGGACAACGCCT‑CAATGAATAACATCCTTGTTcaa                                                                                             >  2:92214/2‑141 (MQ=11)

                                                 cATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACACAATATATGTCTTAGATATTAGCCAACGGCTACAGTG‑ATAACATCCTAGTTAAAtt                                                                                           >  2:754851/1‑141 (MQ=11)

                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTCTATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGGTtat                                                                                     >  2:2533/1‑141 (MQ=17)

                                                       cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATACCATCGTTGTTAAATTGTTtat                                                                                     <  1:405336/141‑1 (MQ=17)

                                                           aaCCCATCCGTAGATTGTTTAATGATGAACGTTTACACGACGTCCTCCTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGTCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCt                                                                                 <  2:537633/141‑1 (MQ=11)

                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGa                                                                               >  2:394832/1‑141 (MQ=32)

                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAGATATTGCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGAG‑ATAACATCCTTTTTAACTTGTTTATATCTGa                                                                               >  1:327676/1‑141 (MQ=11)

                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACAACACGTTCCCTTAATATTTCAGATATATCAAGATAACTCAATGTATACTTTAGATAGTAGGCAACGGCTACAGTG‑ATAACATCCTTGTTAACTTGTTTAGATCTGa                                                                               >  2:329954/1‑141 (MQ=11)

                                                               cAACCGTAGACCCTTGAATAATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                             <  2:468732/141‑1 (MQ=12)

                                                                      aGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                                                                      <  2:96053/141‑1 (MQ=32)

                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                     >  1:303520/1‑141 (MQ=32)

                                                                            ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                                >  2:773191/1‑141 (MQ=32)

                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttt                                                     >  2:581495/1‑141 (MQ=17)

                                                                                           ttACACCACGTTCCTTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  2:610387/141‑1 (MQ=17)

                                                                                           ttACACCACCTTCCCTTAACATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  1:28769/141‑1 (MQ=17)

                                                                                           tgacacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  1:14448/139‑1 (MQ=17)

                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  1:124833/141‑1 (MQ=32)

                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  1:129676/141‑1 (MQ=32)

                                                                                               cccccGTTCTCTTAGTATTTCAAATATACCACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                             <  1:266322/137‑1 (MQ=1)

                                                                                                ccacGTTCCCTTAATATTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                            <  1:23867/141‑1 (MQ=17)

                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                         >  2:18857/1‑141 (MQ=32)

                                                                                                      tCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                                      >  2:497075/1‑141 (MQ=32)

                                                                                                      tCCCTTAACATATCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACt                                      <  2:507431/141‑1 (MQ=17)

                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                <  1:632673/141‑1 (MQ=35)

                                                                                                            aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                <  2:157457/141‑1 (MQ=35)

                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                          >  1:440501/1‑141 (MQ=35)

                                                                                                                   cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                         >  2:773679/1‑141 (MQ=35)

                                                                                                                     aatatatCACGATAACTCAATGGATATCTTAGATAGTAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                       <  2:144553/140‑1 (MQ=17)

                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAATTTATAATATAAATTTAACTTTGCa                     >  2:471689/1‑141 (MQ=255)

                                                                                                                       tatatCACGATAACTCAACGGATATCGTAGATAGAAGCCAACAGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATGCAGAAGACTCCTTTTTGTTAAAATTATACTCTAAATTCAACTTTGCa                     <  1:202110/141‑1 (MQ=11)

                                                                                                                       taaatCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  1:353346/138‑1 (MQ=21)

                                                                                                                         tatCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc                   <  1:416965/141‑1 (MQ=255)

                                                                                                                          atCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACa                  <  2:577815/141‑1 (MQ=35)

                                                                                                                              cGATAACTCAATGTATATCTTAGATAGTAGCCATCGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc              >  2:520093/1‑141 (MQ=25)

                                                                                                                               gatagcgccaTGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcc             <  2:456042/132‑1 (MQ=1)

                                                                                                                               gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAccc             >  2:193791/1‑141 (MQ=25)

                                                                                                                                 tagcTCAGTGCATATCTGATATAGGAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGGTTGGATCTTAAATAGTTCAGACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  1:550511/138‑1 (MQ=11)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  2:242193/138‑1 (MQ=12)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           <  2:76936/138‑1 (MQ=9)

                                                                                                                                 tagcTCAATGCACATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACCCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAccacacctt           <  2:683328/138‑9 (MQ=11)

                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta          >  1:756098/3‑141 (MQ=12)

                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta          >  1:640460/3‑141 (MQ=9)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  2:231193/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  2:788733/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCAAACGGCTACAGTG‑ATAACATCCTGGTTAAATTGTTTATATCTGAAATAGATCATACACAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACTGAACCGtat         >  2:190875/2‑141 (MQ=11)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:772242/1‑141 (MQ=17)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        <  2:702947/141‑1 (MQ=18)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:171321/1‑141 (MQ=16)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:710573/1‑141 (MQ=12)

                                                                                                                                     ccAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta       <  1:66445/140‑1 (MQ=40)

                                                                                                                                      cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATAATTGCTGGACGTTGGATGTCTGAATCCTTGCATACAGAAGACTCCTTTTGGTTAGAAGGATACTATGAATTCAACTTTGCAACAGAACCGtattat      >  1:777847/1‑141 (MQ=9)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:60161/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  2:101883/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:518306/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:129214/141‑3 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTACTTTTATTTGAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt      <  1:346955/141‑3 (MQ=11)

                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     >  2:583819/1‑138 (MQ=255)

                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    >  1:451434/1‑141 (MQ=255)

                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    <  1:713354/141‑1 (MQ=255)

                                                                                                                                        aTGCATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACCGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat    <  1:65944/141‑5 (MQ=255)

                                                                                                                                          gTACATCTTAGATAGTAGCCAACCGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATACTTCATACAGAAGACTCCTTCTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgta  <  2:410062/141‑12 (MQ=11)

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG‑ATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618

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