BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,413	C→T	51.1%	V47I (GTC→ATC)	tpn ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,413	0	C	T	51.1%	16.5 / inf	514	V47I (GTC→ATC)	tpn	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (127/122); new base T (126/135); total (257/257)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-03

Only 100 of 517 total aligned reads displayed.

AATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAC  >  CP000730/36273‑36553

    gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTAGGGCATTATCTTGAACCCAACTGTAGATCGtt                                                                                                                                          >  2:569222/1‑141 (MQ=11)

     cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTg                                                                                                                                         <  1:619896/141‑1 (MQ=18)

     cACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTg                                                                                                                                         <  1:336088/141‑1 (MQ=18)

      aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGa                                                                                                                                        >  2:327384/1‑141 (MQ=255)

       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGaa                                                                                                                                       >  1:361299/1‑141 (MQ=21)

        gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGt                                                                                                                                           <  2:647551/136‑1 (MQ=255)

         atatatAGAGCCCCATTTTCCTTTTATTTTGATATACGCCTCATCATAACTCCATTTGTAATAATCTTT‑TTTATGCTTTTTCTTCCAAATTTGCTATAAAATTGGGGCAAACTGGTTAACCCAACTCTAGATCGTTGAatg                                                                                                                                     <  2:141797/141‑1 (MQ=2)

         atatAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                                                     >  1:442796/1‑141 (MQ=255)

          tataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatga                                                                                                                                    >  1:760606/1‑141 (MQ=255)

           ataAATAGCTCTCTTTTTCTTGTATTTGGAAATAAGTCTCATCAATACGCCATCGTTCATAACTTCT‑TTTTTGCTTTCTCTTCCAAATTTCATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgat                                                                                                                                   <  2:716521/141‑1 (MQ=2)

            taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgatg                                                                                                                                  >  2:491375/1‑141 (MQ=255)

              aaTAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGaa                                                                                                                                >  2:783862/1‑141 (MQ=255)

                tAGCTCCATTTTCCTTTTTTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACg                                                                                                                              <  1:527070/141‑1 (MQ=255)

                tAGCTCCATTTTCCTTTTATTTTGATGTACGTGTCATCAATACGCCATTTGTAATAAGGTTT‑TTTATGCTTTTTCTTTCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCACCGGTAGACCGTTGAATGATGAACg                                                                                                                              >  2:48185/1‑141 (MQ=11)

                 agttcctttttccttttATTTTGATGTCCGACTCATCAATACTCCATTTGTTATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGt                                                                                                                             <  1:185327/138‑1 (MQ=14)

                    tCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTa                                                                                                                          >  2:122849/1‑141 (MQ=25)

                      cATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca                                                                                                                        >  1:834306/1‑141 (MQ=21)

                       aTTTTCCTTTTATTTTGATGTACGTCTCATCAAGGCGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGACCCCAACGGTAGACCGTTGAAGGATGAACGTCTacgc                                                                                                                       <  2:173378/141‑3 (MQ=11)

                       aTTTTCCTTTTATTTTCATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTacac                                                                                                                       <  1:194859/141‑1 (MQ=11)

                        ttttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacc                                                                                                                      <  1:648720/141‑1 (MQ=255)

                         tttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                                                     <  1:55169/141‑1 (MQ=40)

                         tttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTCCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                                                     <  2:606315/141‑1 (MQ=18)

                         gttCCTTTTATTTTGATGTACGTCTCATCACTACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                                                     <  2:508659/140‑1 (MQ=16)

                          ttccagttattaGGATGGACGTCGTATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGATCCCAACGGTAGACCGTTGAATGATGAACGTTTAcaccac                                                                                                                    <  1:165554/135‑1 (MQ=11)

                          ttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacgac                                                                                                                    >  2:536356/1‑141 (MQ=255)

                          ttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcaccac                                                                                                                    >  1:803444/1‑141 (MQ=255)

                          ttCCTTTTATTTTGATATACGTCTCATCAATACGCCATGTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcaccac                                                                                                                    <  2:513502/141‑1 (MQ=255)

                            ccTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCACCTTTGAGAGACAATTGGATCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGt                                                                                                                  <  1:401318/141‑1 (MQ=11)

                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTCATCAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATGCGGACCCATTATTGAACCCAACGGTCGACCGATAAATAATGAACgatcagccaacact                                                                                                                 >  2:155086/1‑128 (MQ=11)

                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  2:42192/1‑141 (MQ=40)

                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCGTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtc                                                                                                                 >  1:576039/1‑140 (MQ=9)

                             cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  1:466579/1‑141 (MQ=255)

                             cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAAGAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 <  1:102679/141‑1 (MQ=255)

                               tttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTcc                                                                                                               >  1:677631/1‑141 (MQ=21)

                                 tATTTTGATATACGTCTTATCAATACGCCATCTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATGTGATCTAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctct                                                                                                             <  2:664153/141‑1 (MQ=11)

                                  tttttGAGATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  2:454776/140‑1 (MQ=11)

                                  aTTTTGTTATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATCCTTTTTCTTCCAAAGTTGATATAAAATTGAGGCATGTACTTGATCCAACCGGTGGGTGGTTGAATGATGAAGATTTACACCACGATctttt                                                                                                            >  2:341807/1‑137 (MQ=2)

                                  aTTTTGTTATACGTCTCATCAAGACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  2:348734/141‑1 (MQ=11)

                                  aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            <  1:752655/141‑1 (MQ=18)

                                  aTTTTGATATACGTCTCATCAATACGTCACTTTGAATAAGGTTTAAGTAT‑CTTAGCCATCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  1:604789/141‑1 (MQ=2)

                                  aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  1:477312/141‑1 (MQ=255)

                                  aTTGTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            <  1:26777/141‑1 (MQ=21)

                                  aTTGTGATGTACGTCGCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            <  2:631671/141‑1 (MQ=255)

                                   ttttGATTGACGTCGCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACCAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa                                                                                                           <  1:207852/141‑1 (MQ=255)

                                    tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTGCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGATCCCTTCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa                                                                                                          <  1:827755/141‑1 (MQ=11)

                                     ttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAat                                                                                                         <  2:619666/141‑1 (MQ=18)

                                     ttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAat                                                                                                         <  1:246872/141‑1 (MQ=255)

                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCCCGTTCTATTTatat                                                                                                       >  2:584256/5‑141 (MQ=11)

                                        aTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAtt                                                                                                      <  2:708341/141‑1 (MQ=39)

                                           tACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                                                                   <  2:286093/141‑1 (MQ=255)

                                           tACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGATCCCTTAATATTTCa                                                                                                   >  1:51341/1‑141 (MQ=21)

                                            aCGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAg                                                                                                  >  1:504738/1‑141 (MQ=255)

                                             cgtctacccaATACGCTATTTGTAATAAGCTTT‑GTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGTATATTCTTGAACCCAACGGTAGATCGCTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCCGa                                                                                                 <  1:636782/133‑1 (MQ=11)

                                             cGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTACTATTTCAGa                                                                                                 >  2:747121/1‑141 (MQ=255)

                                              ctctcATCAATACGCCATTTGTAATCAGCTTT‑TTTCTGCTCTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGat                                                                                                <  1:154985/140‑1 (MQ=255)

                                                ctcATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAAGACACCTTCTCTTAATATTTCAGatat                                                                                              <  2:594257/141‑1 (MQ=11)

                                                 tcatcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAACTTGAGGCATCTTCTTGAAGCCCACCGTAGAACCTTGAATGATGAACCTTTACACCACGACCCCTTAATCTTTCAGttata                                                                                             >  2:90053/1‑141 (MQ=11)

                                                 tcatcaATACGCCATTTGTAATAAGCTTT‑GTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatata                                                                                             <  1:87242/141‑1 (MQ=17)

                                                 tcatcaATACGCCATTTGTAATAAGCTTA‑TTTATGCTTTTTCTTCCAAATTTGACACAAAATTGGGGCATATTCTTGAAGCCAACGGTAGACCGTTGAATGATGAACGTTTACNCCANGNTCCCTTAATATTTCAAatata                                                                                             >  1:541288/1‑141 (MQ=14)

                                                   atcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                           <  2:135360/141‑1 (MQ=21)

                                                   atcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGATGAATGAGGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                           >  1:567782/1‑141 (MQ=255)

                                                    tcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          >  2:768561/1‑141 (MQ=255)

                                                    tcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          >  2:736945/1‑141 (MQ=255)

                                                    tcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          >  1:666001/1‑141 (MQ=255)

                                                    tcaATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCa                                                                                          >  1:38510/1‑141 (MQ=21)

                                                     caATACGCCATTTGTAATAAGCTTT‑TTTATGGTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                         >  2:366153/1‑141 (MQ=255)

                                                      aaTACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCg                                                                                        >  1:175172/1‑141 (MQ=21)

                                                      aaTACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACg                                                                                        >  2:195651/1‑141 (MQ=21)

                                                       aTACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                       >  2:102738/1‑141 (MQ=255)

                                                       aTACGCCATTTGTAACAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGa                                                                                       <  1:226734/141‑1 (MQ=21)

                                                        tACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                      <  1:13402/141‑1 (MQ=255)

                                                         aCGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                                                                     <  2:308163/141‑1 (MQ=18)

                                                         aCGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                                                                     >  2:685115/1‑141 (MQ=40)

                                                         aCGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                                                                     >  1:795335/1‑141 (MQ=40)

                                                          cGCCATTTGTATTAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATag                                                                                    <  2:526084/141‑2 (MQ=11)

                                                          cGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATag                                                                                    >  1:713007/1‑140 (MQ=255)

                                                          cGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                    <  1:391525/141‑1 (MQ=40)

                                                           gCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                   <  1:43232/141‑3 (MQ=255)

                                                           cccAGTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                   <  1:306443/140‑1 (MQ=21)

                                                            ccATTTGTAATAAGCTTT‑TTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACgtttgct                                                                                  >  2:10875/1‑135 (MQ=11)

                                                            ccATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                  >  2:227140/1‑141 (MQ=255)

                                                            ccATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                  <  1:327425/141‑1 (MQ=18)

                                                            ccATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                  >  1:258046/1‑141 (MQ=21)

                                                              aTTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATTTTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGCTAACTCa                                                                                >  2:540692/1‑141 (MQ=255)

                                                               tttGTAATCAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTATCACGATAACTCaa                                                                               <  1:416217/141‑1 (MQ=12)

                                                                ttGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                              <  2:795412/141‑1 (MQ=21)

                                                                 tGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCACCGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                             >  2:258692/1‑141 (MQ=11)

                                                                  gTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTTCACAGCGGCCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                            >  1:781269/1‑140 (MQ=11)

                                                                   taataaGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                           >  1:127246/1‑139 (MQ=11)

                                                                   taataaGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                           >  1:119254/1‑139 (MQ=11)

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

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

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

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

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

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

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

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

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

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

AATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTT‑TTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAC  >  CP000730/36273‑36553

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