BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,479	T→C	15.9%	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	15.9%	164.8 / 24.6	89	T25A (ACA→GCA)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (37/36); new base C (9/5); total (46/43)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618

taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                            >  2:83144/1‑141 (MQ=255)

taataaGCTTTTTGATGCTTTTTCTTCCAAATTTCATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGTACTTTTACACCACGATCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                            <  1:325140/141‑1 (MQ=14)

 aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCAATGaat                                                                                                                                           >  2:150208/1‑141 (MQ=12)

 aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtaa                                                                                                                                           >  1:455946/1‑140 (MQ=255)

  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGaata                                                                                                                                          >  2:5688/1‑141 (MQ=17)

     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                       <  1:727624/141‑1 (MQ=39)

     aGCTTTTTTAGTCTTTTTCTTCCATATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                       <  2:694149/141‑1 (MQ=255)

         tttttATGCTTTTTCTTCCAAATTTGATTCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                                   <  2:690943/141‑1 (MQ=25)

          ttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAt                                                                                                                                  <  1:42192/141‑1 (MQ=35)

           tttATGCTTTTTCTTCCAAATTTGATACCAAATTGTTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                                 >  1:404323/1‑141 (MQ=255)

           tttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                                 >  1:188222/1‑141 (MQ=35)

              ttttttttttttACAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGAACGTTGAATGATGAACGTATCCACTACGTTCCCTTAATTTTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                              >  2:61471/5‑141 (MQ=11)

               tGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagtag                                                                                                                             >  1:699634/1‑141 (MQ=34)

                 cTTTTTCTTCCAAATTTGATACAAAATTTGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGcc                                                                                                                           <  1:837981/141‑1 (MQ=25)

                         tcCAAATTTGATACAAACTTGGAAATGCTATTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                   >  1:830863/1‑141 (MQ=11)

                         tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                   >  1:762679/1‑141 (MQ=35)

                             aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                               <  2:556495/141‑1 (MQ=35)

                                ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                            <  1:295031/141‑1 (MQ=35)

                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                           <  1:271504/141‑1 (MQ=35)

                                  gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATaa                                                                                                          >  2:112495/1‑141 (MQ=35)

                                     aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGGATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAt                                                                                                       >  2:92717/1‑141 (MQ=21)

                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTTCAGTGATAACATCCTTGTTa                                                                                               <  2:595522/141‑1 (MQ=17)

                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCGTAGATAGTAGCCAACGGCTACAGGGATAACATCCTTGTTa                                                                                               >  1:71894/1‑141 (MQ=17)

                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                              >  1:18646/1‑141 (MQ=32)

                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                             >  1:2874/1‑141 (MQ=32)

                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACATTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                            >  2:834053/1‑141 (MQ=17)

                                                    atTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTTCTTCCCGTTGTCTTAATATTTCTGATCTCTTACGATAACAAAAACTACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                        <  1:401720/141‑1 (MQ=11)

                                                      tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                      <  1:136086/141‑1 (MQ=32)

                                                      tCTGGAACCCAACGTGAGGCCCTTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAGTATTTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                      <  1:650565/141‑1 (MQ=11)

                                                       cTTGAACCCAACGGAAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtat                                                                                     >  1:437071/1‑141 (MQ=17)

                                                          gAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATGTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATc                                                                                  >  1:581478/1‑141 (MQ=255)

                                                          gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                  >  1:255631/1‑141 (MQ=32)

                                                          gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                  >  2:686882/1‑141 (MQ=32)

                                                          gAACCCAACGGGAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACATTGATAACATCCTTGTTAAATGGTTTATATc                                                                                  >  2:284525/1‑141 (MQ=12)

                                                            aCCCAAACGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                                <  1:179927/141‑1 (MQ=17)

                                                             cccAACGGTAGACCGTTGAATTCTGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTAGATCTGa                                                                               <  2:464436/141‑1 (MQ=12)

                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                               >  1:837695/1‑141 (MQ=32)

                                                              tcAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATTTTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAGAACATCCTTGTTAAATTGTTTATATTTGaa                                                                              <  2:508739/140‑1 (MQ=12)

                                                               cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATGGTTTATATCTGaaa                                                                             <  1:547741/141‑1 (MQ=17)

                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        <  1:786618/141‑1 (MQ=32)

                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        >  1:470960/1‑141 (MQ=32)

                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        >  1:470750/1‑141 (MQ=32)

                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        >  2:645567/1‑141 (MQ=32)

                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        <  2:762679/141‑1 (MQ=32)

                                                                         ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                   >  2:171911/1‑141 (MQ=32)

                                                                            gtGAATGATGAACGGTTAAACCCAGTTCTCTTAAAATTACCGGTATATCACGATCGTTCAATGTCGATCTGGGATATTAGACGTGGGCCACAGGTAGAATAGCCTTGGTATATTGGTTATATCTGAAATAGTTCATACaga                                                                <  2:198023/140‑1 (MQ=2)

                                                                             tGAATGATTAACGTTTACACCACGTGCCCTGAAGATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                               <  2:54295/141‑1 (MQ=255)

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

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

                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACGTAGTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                                                     >  1:478394/1‑141 (MQ=255)

                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                                                     >  2:691288/1‑141 (MQ=32)

                                                                                         gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                   >  2:425974/1‑141 (MQ=32)

                                                                                          tttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                  >  2:56874/1‑141 (MQ=32)

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

                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 >  2:149682/1‑141 (MQ=32)

                                                                                             tcaccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  2:349532/140‑1 (MQ=17)

                                                                                              caccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaa                                              >  1:211831/1‑141 (MQ=32)

                                                                                                  acGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                          >  2:704761/1‑141 (MQ=32)

                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                         >  2:765449/1‑141 (MQ=17)

                                                                                                        ccTTAATAGTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                    <  1:266046/141‑1 (MQ=17)

                                                                                                          ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                                  >  1:30189/1‑141 (MQ=33)

                                                                                                           tAATATTTTATATAAATCACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATaaa                                 <  1:412/141‑1 (MQ=11)

                                                                                                           tAATATTTCAGATATAGCACGATAACTCAATGTATATCTTAGATAATAGCCAACTGCTACAGTTATAATATCCTTGTTAAATTGTTTATATCTCAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                                 <  2:247909/141‑1 (MQ=11)

                                                                                                               atgtCAGCTATATCACGACAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                             <  1:108250/138‑1 (MQ=255)

                                                                                                               atTTCCGATATATCACGATAACTCCATGTACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                             <  1:358286/141‑1 (MQ=255)

                                                                                                               atTTCAGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                             <  1:17475/141‑1 (MQ=25)

                                                                                                                  ttaGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                          <  1:27951/139‑1 (MQ=25)

                                                                                                                    aGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACttt                        >  1:102941/1‑141 (MQ=35)

                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     >  2:333053/1‑141 (MQ=34)

                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  1:83144/141‑1 (MQ=35)

                                                                                                                       taaatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  2:699634/138‑1 (MQ=21)

                                                                                                                             aCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGaa               >  1:96305/1‑141 (MQ=35)

                                                                                                                             aCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACACAAGACTCCTTTTTATTAAAATTGTACTATAAATTAACCTTGGGACCAGaa               >  1:195328/1‑141 (MQ=11)

                                                                                                                              cGATAACTCAATGTATATCTTAGATCGTACCTATCGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCACCAGAAc              >  1:404794/1‑141 (MQ=11)

                                                                                                                               gagaaCTCAATGCATATATTAGATAGTAGCCAACAGCGAAAGTGATAACATCCTTGTTAAATTGTTAATAACTGAAATAGTTCATACAGAAGACCCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc             <  1:535071/138‑1 (MQ=11)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCATCGGCTACAGTGATATCGACCTTGTTAAATTGTTTATATCTCAAATAGTTTATACATAATACTCCTTTTTTTTAAAATTATATTATTAATTCAACTTTGCAACAGAACCGt           <  2:31402/138‑1 (MQ=11)

                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           >  2:252535/4‑141 (MQ=2)

                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt           <  1:143588/141‑1 (MQ=37)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATATTATAAATACAACTTTGCAACAGAACCGtat         >  2:824909/2‑141 (MQ=255)

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

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTACTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  1:662117/2‑141 (MQ=12)

                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAGCCTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  1:795425/2‑141 (MQ=255)

                                                                                                                                    cTCAATGTATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatt        <  1:168676/141‑1 (MQ=38)

                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:86602/1‑141 (MQ=12)

                                                                                                                                    atCAAAGCATATCGTAGATGGTCGACAACGGCTACAGTGATAACATCCTAGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTTGTTTTTGTTAGAATTATATTATAAATTCAACTTTCCAACAGAACGGtatt        >  2:391030/2‑141 (MQ=11)

                                                                                                                                     tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta       >  1:332464/1‑141 (MQ=40)

                                                                                                                                     tCAATGCATATCTTAGATAGTAGCCAACCGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt       <  1:10620/141‑2 (MQ=255)

                                                                                                                                     tCAATGCATATCTTAGATAATAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt       >  1:469163/1‑140 (MQ=255)

                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTAATTTAAAGTATATTATAAATTCACCTTTGCAACAGAACCGtatttt      >  1:567633/1‑139 (MQ=11)

                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     >  1:762678/1‑138 (MQ=255)

                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     <  1:344343/141‑4 (MQ=255)

                                                                                                                                       aaTGAATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaact     <  2:345749/141‑9 (MQ=11)

                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    <  1:766548/141‑1 (MQ=255)

                                                                                                                                          gTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGaa  <  1:92717/141‑1 (MQ=255)

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