BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,378	C→T	92.6%	L58L (TTG→TTA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,378	0	C	T	92.6%	1656.0 / 104.4	811	L58L (TTG→TTA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (34/28); new base T (376/366); total (415/396)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 814 total aligned reads displayed.

CAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG  >  CP000730/36240‑36510

   ccAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATaaatt                                                                                                                                 >  2:238445/1‑139 (MQ=255)

    cAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTTTACGTCTCATCAATACTCAATTTGTAATAAGATTTTTTATGCTTTCTCTTCCAAGTTTGATACAACAtt                                                                                                                                >  1:57802/1‑141 (MQ=14)

      aaTATCTAGTGGATGTGGCTCTGGATCAATGGCGCGATATAAATAGCTCCATTTTGCTTTTATTTGGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                                                                                                                              <  2:31392/141‑1 (MQ=11)

      aaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCCATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTTATATAAATTTgg                                                                                                                              >  2:405244/1‑141 (MQ=11)

      aaTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCAATTTGTAATAAGCTTTTTTATACTTTTTATTCCAAATTTGATATAACATTgg                                                                                                                              >  1:79684/1‑141 (MQ=11)

      aaTATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTgg                                                                                                                              >  2:689836/1‑141 (MQ=38)

       atatCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTggg                                                                                                                             >  2:469341/1‑141 (MQ=255)

       atatCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGTTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTCTGCTTTTGCTTCCAAATTTGATATAAAATTggg                                                                                                                             >  1:154955/1‑141 (MQ=255)

            tactgtatctCCCCTTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                        <  1:710601/138‑1 (MQ=11)

            tAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGGTATACGTCTCTTCAATACGCCATTTGTGATAAGCTTTTTTATGTTTTTTCTTCCAAATATGGTATAAAATTGGGGcttt                                                                                                                        >  1:233610/1‑138 (MQ=11)

            tAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                        <  1:684672/141‑1 (MQ=37)

            tAATGCATGTCCCTCTGTATTAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCata                                                                                                                        <  1:390602/141‑1 (MQ=11)

             aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTTATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatat                                                                                                                       >  1:659089/1‑141 (MQ=255)

             aaTGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCatat                                                                                                                       <  1:305171/141‑1 (MQ=25)

              atgtatgtCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGTTATAAAATTGGGGCATAtt                                                                                                                      >  2:482884/1‑141 (MQ=255)

                 tatgtCCCTCTGTATCAATGGCACGATATAAATTGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCtt                                                                                                                   <  2:424146/141‑1 (MQ=255)

                  ctgtCCCTCTGTATCAATTGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATATTCTTg                                                                                                                  <  1:119731/140‑1 (MQ=11)

                   tgtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGa                                                                                                                 <  1:730842/141‑1 (MQ=255)

                    gtCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGAATATTATTGaa                                                                                                                >  2:567838/1‑141 (MQ=11)

                      cccTCTGTATCAATGGCACGATATAAATAGCTCCATTTTGCTTTGATTGTGGTATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGTTATAGAATTGGGGCATATTCTTGAAcc                                                                                                              >  2:518306/1‑141 (MQ=11)

                       ccTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAATTTGGGGCATATTCTTGAAccc                                                                                                             >  2:113676/1‑141 (MQ=21)

                         tctGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACTCTATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCCAATTTGATATAACGTTGGGTCTTATTCtttattcgac                                                                                                           >  1:274273/1‑133 (MQ=11)

                         tctGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCaa                                                                                                           >  2:54544/1‑141 (MQ=255)

                         tctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATAAAATTGGGGCATATTCTTGAAACCaa                                                                                                           >  2:178596/1‑141 (MQ=255)

                         tctGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTCGGGCATATTCTTGAACCCaa                                                                                                           <  1:479470/141‑1 (MQ=21)

                         tctGCATCAATGACACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCaa                                                                                                           >  2:539143/1‑141 (MQ=25)

                            gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTGCGTATTCTTGAACCCAACgg                                                                                                        >  1:167272/3‑141 (MQ=11)

                            gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                                                                        >  2:226975/3‑141 (MQ=255)

                            gtaTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACgg                                                                                                        <  2:746145/139‑1 (MQ=255)

                             taTCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGt                                                                                                       <  2:785867/140‑1 (MQ=19)

                                cactGGCACGACATAAACAGCTTCTTTTTCCTTTTATTTTGATGTACTTCTCATCAATACGCCATTTGTAATCAGCTTTTCTCAACTCTTTCTGAACACTTAAAGAAAATATTGGGGCATATTCTTGAACCCAACGGTAGa                                                                                                    <  1:778463/138‑1 (MQ=2)

                                cAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGa                                                                                                    <  2:620068/141‑1 (MQ=18)

                                cAATGGCACGATATAAATAGATCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCTTATTCTTGAACTCAACGGTAGa                                                                                                    >  2:156821/1‑141 (MQ=11)

                                cAATGCCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCACCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGa                                                                                                    <  1:703076/141‑1 (MQ=255)

                                 aagggCACTATATAAATTGCTCCTTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTTTAATAAGTTTTTTTATGCTTTTTCTTCTAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGat                                                                                                   <  1:355525/138‑2 (MQ=11)

                                 aaTGGCACGATATAAATAGCTCCTTTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGat                                                                                                   <  2:49238/141‑2 (MQ=255)

                                 aaTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGat                                                                                                   >  1:285221/1‑140 (MQ=21)

                                  aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCACATTTGATACACAATTGGGGCATATTCTTGAACCCAACGGTAGAcc                                                                                                  >  2:66445/1‑141 (MQ=255)

                                  aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatc                                                                                                  >  2:396105/1‑139 (MQ=255)

                                  aTGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGatc                                                                                                  >  1:459713/1‑139 (MQ=255)

                                   tGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAGTCCATATAAAACTGGCGCATATCCTTGGACGCAACGGTAGACCg                                                                                                 >  1:550147/1‑141 (MQ=11)

                                    ggCACGCTATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGt                                                                                                <  1:521904/141‑1 (MQ=21)

                                    ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGt                                                                                                >  2:638343/1‑141 (MQ=39)

                                    ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTGGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAAACCAACGGTAGACCGt                                                                                                >  2:144508/1‑141 (MQ=255)

                                    ggCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCGCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGt                                                                                                >  1:542316/1‑141 (MQ=255)

                                       ccGATATAAATAGCTCCATTTTCCGTTTCTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGa                                                                                             <  2:575408/140‑1 (MQ=255)

                                       aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCCTATTCTTGAACCCAACGGTAGATCGttca                                                                                             >  2:721756/1‑139 (MQ=11)

                                       aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGa                                                                                             <  1:367592/141‑1 (MQ=255)

                                        cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCTAAAATTGGGGCATATTATTGAACCCAACGGTAGATCGTTGaa                                                                                            >  2:766887/1‑141 (MQ=11)

                                         gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAat                                                                                           <  1:447733/141‑1 (MQ=255)

                                           tataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTGATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAatga                                                                                         >  2:512389/1‑141 (MQ=25)

                                            aaaaaTAGCCCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgat                                                                                        <  2:279769/139‑1 (MQ=11)

                                             taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgatg                                                                                       >  1:86554/1‑141 (MQ=255)

                                             taAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAatgatg                                                                                       <  2:376858/141‑1 (MQ=255)

                                                   gCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                 >  1:295180/1‑141 (MQ=255)

                                                     tCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                                                                               <  1:469341/141‑1 (MQ=255)

                                                      ccATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTac                                                                              >  1:588448/1‑141 (MQ=255)

                                                       cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTCGATCGTTGAATGATGAACggtaccc                                                                             >  1:764061/1‑135 (MQ=11)

                                                         ttttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacc                                                                           <  1:564467/141‑1 (MQ=255)

                                                         ttttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGCGGCATATTCTTGAACCCAACGGTAGATCGTTGAAGGATGACCATTTAcacc                                                                           >  2:445235/1‑141 (MQ=11)

                                                             ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGCTGAACGTTTACACCACGt                                                                       >  2:713354/1‑141 (MQ=255)

                                                              cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt                                                                      >  1:127458/1‑141 (MQ=25)

                                                               ttttATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTTCACCAGGTTc                                                                     >  1:183825/1‑141 (MQ=255)

                                                                 ttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctc                                                                   <  1:426604/141‑3 (MQ=255)

                                                                   aTTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                 <  2:599218/141‑1 (MQ=18)

                                                                     ttttATGTACATCTCATCAATACGCCATTTGTAATAAGCTCTTTTATGCTTTTTCTTCCAAATCTTATATAAAATTGGCCCATATTCTTGAACCCAACGGTAGACCGTTGAACGATTAACGTTTACACCACGTTCCCTTaa                                                               <  2:303389/141‑1 (MQ=11)

                                                                     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                               <  1:622828/141‑1 (MQ=255)

                                                                        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                            >  2:125817/5‑141 (MQ=255)

                                                                        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                            >  2:293694/5‑141 (MQ=255)

                                                                        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGGTTACACCCCGTTCTCTTAatat                                                            >  1:101883/5‑141 (MQ=11)

                                                                            tACGTCTCATCAATTCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCATGCAGTAATTTTAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                        <  1:375128/141‑1 (MQ=11)

                                                                            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                        <  1:2265/141‑1 (MQ=18)

                                                                            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                        <  1:603880/141‑1 (MQ=21)

                                                                                tctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCTACTTTCATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGata                                                    <  1:537262/141‑1 (MQ=11)

                                                                                  tcatcaATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGGATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatata                                                  <  1:748994/141‑1 (MQ=255)

                                                                                  tcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTCCTCTTAATATTTCAGatata                                                  >  1:638205/1‑141 (MQ=255)

                                                                                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                <  1:719125/141‑1 (MQ=21)

                                                                                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCa                                               >  1:644462/1‑141 (MQ=25)

                                                                                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATAAc                                              >  1:349139/1‑141 (MQ=255)

                                                                                      caATACGCCATTTGTAATAAGATTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACATTTACAATAGGTTCCCTTAATATTTCAGatatatacc                                              >  2:356891/1‑138 (MQ=14)

                                                                                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGCTATATCATGa                                            >  2:610358/1‑141 (MQ=11)

                                                                                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGGATATAGATTAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                            <  2:166956/141‑1 (MQ=11)

                                                                                        aTACGCCATTTGTAATAAGCTTTTGTATGCTTTTTCTTTCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                            >  1:506144/1‑141 (MQ=11)

                                                                                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                           >  1:213866/1‑141 (MQ=255)

                                                                                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTAGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTGCAGATATATCACGAt                                           <  1:325324/141‑1 (MQ=11)

                                                                                          aCGCCATTTGTAATACGCTTTTGTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATACCACGATa                                          <  1:491481/141‑1 (MQ=11)

                                                                                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATa                                          <  1:670016/141‑1 (MQ=18)

                                                                                           cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTAAGATATATCACGATag                                         >  2:631795/1‑140 (MQ=255)

                                                                                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAAGTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                       >  1:416904/1‑141 (MQ=11)

                                                                                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCGGATATTTTACGGTaatt                                       >  1:264366/1‑139 (MQ=11)

                                                                                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                      <  1:436345/141‑1 (MQ=11)

                                                                                               ctttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTGGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                     <  2:384586/140‑1 (MQ=11)

                                                                                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                  <  2:253603/141‑1 (MQ=40)

                                                                                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                >  2:103200/1‑139 (MQ=255)

                                                                                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                <  1:242959/141‑3 (MQ=11)

                                                                                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                >  1:198885/1‑139 (MQ=11)

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

                                                                                                      ataaGCTTTTTTATGCTTTATATTGCACATTTGATATAAACTTGGGGGAGATTAGAGAACCCGAAGATGGATCGGTGAGTGATGAACGTTTAAGAGACGTTATCTTAATATTTCAGATACATCACGATAGCGCaagatatt                              >  2:204848/1‑135 (MQ=255)

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

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

CAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTG  >  CP000730/36240‑36510

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