BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,444	C→T	47.9%	R36R (AGG→AGA)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,444	0	C	T	47.9%	98.3 / 248.6	236	R36R (AGG→AGA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (64/58); new base T (52/60); total (118/118)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.48e-01

Only 100 of 241 total aligned reads displayed.

TTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTA  >  CP000730/36305‑36578

 tATTTTGATGTACGTCTTCTCAATACGCCATTTGTAATAAGCTTTTGTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTccc                                                                                                                                        <  1:578253/140‑1 (MQ=255)

 tATTTGGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCATATTCTTGAACCCAACGGTAGACCGGTGAATGATGAACGTTTACACCACGTTccc                                                                                                                                        <  1:336495/140‑1 (MQ=255)

   ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGTTGAACGTTTACACCACGTTCCCtt                                                                                                                                      >  1:415688/1‑140 (MQ=25)

   ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                      >  2:531781/1‑140 (MQ=21)

   ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                      <  1:31153/140‑1 (MQ=21)

   ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                      <  2:234596/140‑1 (MQ=255)

   ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                      <  2:654941/140‑1 (MQ=11)

   ttttGATATACGTCTCATCAAGACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCTAAGTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                      <  1:281771/140‑1 (MQ=11)

   ctctGATATACGTCTCATCAATACGCCATTCGTAAAAAGCTTTTTTAAGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAAATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                      <  1:2508/137‑1 (MQ=11)

     ttGATATACGTCTCATTAACAGGTCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                                                                                                    <  2:268186/140‑1 (MQ=11)

     ttGATATACGTCTCAGCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAGTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                                                                                                    <  1:120608/140‑1 (MQ=11)

      tGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGACGAACGTTTACACCACGTTCTCTTAat                                                                                                                                   <  2:412886/140‑1 (MQ=11)

       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTCata                                                                                                                                  >  2:568054/5‑140 (MQ=11)

       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAata                                                                                                                                  >  2:478917/5‑140 (MQ=11)

       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTAATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATTGTTGAATTATGAAAGTTTACACCACGTTGTCTTAata                                                                                                                                  >  1:385211/5‑140 (MQ=11)

       gATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGCCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAata                                                                                                                                  >  2:578867/1‑140 (MQ=25)

       gATATACGTCTCACTAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACAACACGTTCTCTTAata                                                                                                                                  <  1:494185/140‑1 (MQ=11)

       caTGTACGTCTCTTCAAACCGCCATGTGTAATAAGCTTCTTTATGCTTTTTCTTTCAAGATTGAGATAAAATAGAGCCATATTTTTGAACCGCACTTTACACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAata                                                                                                                                  >  2:219227/2‑140 (MQ=2)

        aTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                                                 >  2:194368/1‑140 (MQ=25)

         tatACGTCTTATCAATACGACATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATAtt                                                                                                                                <  1:300599/138‑1 (MQ=11)

         tatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTTATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATAtt                                                                                                                                <  1:533600/138‑1 (MQ=11)

           tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTACAAATTTGATATAAAATTGGGGCATATTCTTGAACACAACGGTAGACCGTTGAATGATGCACGTTTACACCACGTTCCCTTAATAttca                                                                                                                              >  2:491171/1‑138 (MQ=11)

            aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                                             <  1:436918/140‑1 (MQ=18)

            aCGTCTCATCAATACGCCATTTGTAATAAGCTCTTTTATGCTTTTTCTTCCAAATTTGATACAAATTTGGGGCATATTCTAGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                                             <  2:371815/140‑1 (MQ=255)

             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTGCACCACGTTCTCTTAATATTTCAg                                                                                                                            >  2:611476/1‑140 (MQ=255)

             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAg                                                                                                                            <  2:48495/140‑1 (MQ=18)

             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAACGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTTAg                                                                                                                            >  1:439514/1‑140 (MQ=255)

               tctcATCAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGAGAGTAGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGat                                                                                                                          <  1:129740/140‑1 (MQ=11)

                ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACTGTTTCTTAATATTTCAGata                                                                                                                         >  2:217646/1‑140 (MQ=11)

                  catcaATACGCCATTTGTAATACGCCTTTTCACGCTTCTTCTACCAAAAAAGATAAAAAAAAGGGACATAATCTTGAACCCAACGGCAAATCGCTGAATGATGAACGTTTACACCACGTCCTCTAACTATTTCAGatata                                                                                                                       >  1:172068/1‑140 (MQ=2)

                  catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatata                                                                                                                       <  1:383038/140‑1 (MQ=255)

                   atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAAATAGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatat                                                                                                                      >  1:217568/1‑140 (MQ=11)

                    tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCCCATACTCCTAATATTTCAGATATATc                                                                                                                     <  2:643241/140‑1 (MQ=11)

                    tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATc                                                                                                                     >  1:147980/1‑140 (MQ=255)

                    tcaATACGCCATTCGTAATAAGCTTTTTTATGCTTTTTCTTGCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGACCGTAGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                     <  2:430956/140‑1 (MQ=11)

                     caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATTTCTATTACACAAAGTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                    <  1:316844/140‑1 (MQ=11)

                     caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCa                                                                                                                    <  1:192247/140‑1 (MQ=21)

                     caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGTCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCa                                                                                                                    <  1:450087/140‑1 (MQ=21)

                     caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTAATATTTCAGATATATCa                                                                                                                    >  2:601144/1‑140 (MQ=19)

                      aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                                                   <  1:564137/140‑1 (MQ=255)

                      aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGCTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                                                   <  1:662027/140‑1 (MQ=255)

                       aTACGTCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACg                                                                                                                  <  1:257318/140‑1 (MQ=21)

                       aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTCGACCGTTCAATGATGAACGTTTTCACCGCGTTCCATTAATATTTCAGatatacaacg                                                                                                                  >  2:18589/1‑135 (MQ=14)

                        tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                                                 >  1:402853/1‑140 (MQ=255)

                        tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATATGGGCATATTCTTGAACCCAACGGTAGATAGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATAAATCACGa                                                                                                                 >  1:422022/1‑140 (MQ=11)

                        tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGAGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGa                                                                                                                 >  1:45790/1‑140 (MQ=255)

                         aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTTTTAATATTTCAGATATATCACGAt                                                                                                                <  1:485848/140‑1 (MQ=11)

                           gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATag                                                                                                              >  1:512503/1‑139 (MQ=255)

                            ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTATTATTTCATATATATCACGATagc                                                                                                             >  2:520696/1‑138 (MQ=11)

                            ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                                             <  2:357473/140‑3 (MQ=255)

                             cATTTGTAATAAGCTTTTTTTTTTTTTTTCTTCCAAAAAAGATATAAAATTGGGTCATATTCTTGAACCCAACGGAAGATCGTTGAATGACGAACGTTTACACCACGTTCTCTTAAAATTTCAGATATATCACTATAGCt                                                                                                            >  2:32446/1‑140 (MQ=11)

                              aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                           <  1:630100/140‑1 (MQ=11)

                              aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                           <  2:528791/140‑1 (MQ=11)

                                 tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCCGATATATCACGATAGCTCAAt                                                                                                        >  1:461802/1‑140 (MQ=11)

                                 tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                        >  2:131162/1‑140 (MQ=21)

                                 tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTCCCCTTAATAGTTCATAGATATAACGATAACTCAAt                                                                                                        >  2:293139/1‑140 (MQ=11)

                                  gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCCATg                                                                                                       >  2:1552/1‑140 (MQ=19)

                                   taataaGCTTTTTTATGCTTTTTTTTTTTTTTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCCATgc                                                                                                      >  2:539164/1‑139 (MQ=11)

                                   taataaGCTTTTTTATGCTTTTTCTTCTAAATTTGATATAAAAGTGGGGCATATTCTTGAACCTAACGGTAGACCGTTGAATGATGAACGTCTACACCACGTTCCCTTAATATTTCCGATATATCACGATAACTCAATgc                                                                                                      <  2:165868/140‑2 (MQ=11)

                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                      >  2:570834/1‑139 (MQ=11)

                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                      >  1:222695/1‑139 (MQ=255)

                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                      <  2:157739/140‑2 (MQ=18)

                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgc                                                                                                      <  2:451913/140‑2 (MQ=18)

                                    aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCTTATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                     <  1:169307/140‑3 (MQ=11)

                                    aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                     <  1:53867/140‑3 (MQ=11)

                                    aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                     <  1:628070/140‑3 (MQ=11)

                                    aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTGAACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                     >  2:409238/1‑138 (MQ=255)

                                    aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                     <  2:619212/140‑3 (MQ=17)

                                     ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTAATCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCCTTGt‑at                                                                                                   >  1:303678/1‑140 (MQ=255)

                                       aaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATCGCTCAATgactat                                                                                                  >  2:664917/1‑135 (MQ=11)

                                          cTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCtt                                                                                              <  2:476748/140‑1 (MQ=37)

                                            tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAg                                                                                            >  2:51894/1‑140 (MQ=34)

                                               ttATGCTTTTTCTTCCAAATTTGATACATAATTGGGGCATATTCTTGTTCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAta                                                                                         <  1:184666/140‑1 (MQ=255)

                                                 aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtagt                                                                                       <  1:68217/140‑1 (MQ=35)

                                                 aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtagt                                                                                       <  1:163321/140‑1 (MQ=35)

                                                 aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGAtagt                                                                                       <  2:105841/140‑1 (MQ=35)

                                                    cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATATTAGc                                                                                    >  2:28334/1‑140 (MQ=21)

                                                    cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGc                                                                                    >  2:28385/1‑140 (MQ=34)

                                                           ttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGc                                                                             <  2:609109/140‑1 (MQ=37)

                                                                         cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAt                                                               >  2:118336/1‑140 (MQ=33)

                                                                           aaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAAACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                             >  1:14716/1‑140 (MQ=17)

                                                                           aaaTTGGGGCAAATTCTTGAACTCAACGTGAAACCGATGAATGATGAACGTTTACACCACGTTCGCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGACTACAGTGATCACCTcc                                                             >  2:659706/1‑140 (MQ=11)

                                                                              ttGGGGCATATTCTTGAACCCATCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                          <  2:630548/140‑1 (MQ=18)

                                                                              ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                          <  1:551904/140‑1 (MQ=33)

                                                                               tGGGGCATATTCTTGAACCCAACGGTAGACGGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTTAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                         <  1:40507/140‑1 (MQ=16)

                                                                                  ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                      >  1:109408/1‑140 (MQ=32)

                                                                                  ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACACAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                      >  1:207991/1‑140 (MQ=17)

                                                                                  ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCTTGTTaa                                                      >  1:37007/1‑140 (MQ=17)

                                                                                   gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑TAATCTTGAATAGTAGCCAACGGCTACAGTGGTAAAATCCTTTTTaaa                                                     >  2:86028/1‑140 (MQ=11)

                                                                                     atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAtt                                                   >  2:404552/1‑140 (MQ=33)

                                                                                               aCCCAACGGTAGACCGTTGAATGACGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                         <  2:204354/140‑1 (MQ=17)

                                                                                                         ctaCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                               >  1:413571/3‑140 (MQ=17)

                                                                                                          gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCa                              >  1:530069/1‑140 (MQ=32)

                                                                                                              gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAAATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATATGAAATAGTACATAcc                          >  1:84679/1‑139 (MQ=12)

                                                                                                                 gAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                       <  1:608582/140‑1 (MQ=33)

                                                                                                                  aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                      >  1:535945/1‑140 (MQ=32)

                                                                                                                    tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                    <  1:275542/140‑1 (MQ=33)

                                                                                                                               tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt         <  2:246545/140‑1 (MQ=33)

                                                                                                                                acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACTAGAGACGCCTTTTTGTTa        >  1:582957/1‑140 (MQ=255)

                                                                                                                                      ctttCTCTTAATTATTCAGTGATAGCACGATAGCTC‑ATGTCATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta  <  2:582835/138‑1 (MQ=11)

TTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGT‑ATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTA  >  CP000730/36305‑36578

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 13 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`