BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,471	A→G	29.8%	S27S (AGT→AGC)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,471	0	A	G	29.8%	232.2 / 130.6	187	S27S (AGT→AGC)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base A (67/64); new base G (21/35); total (88/99)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 197 total aligned reads displayed.

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTAT  >  CP000730/36331‑36611

cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                                              <  1:365093/141‑1 (MQ=21)

cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTACTATTTCAGATATATCACGATaa                                                                                                                                              >  2:763279/1‑141 (MQ=255)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                            >  1:114409/1‑141 (MQ=255)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                                            >  1:831767/1‑141 (MQ=255)

  ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACt                                                                                                                                            >  2:659888/1‑141 (MQ=39)

   cATTTGTAATAGGCTTTTTTATGCTTCTTCTTGCCAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACCCCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                           <  2:514156/141‑1 (MQ=14)

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

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                           <  1:794981/141‑1 (MQ=21)

   cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                                           >  1:66601/1‑141 (MQ=21)

    ctttGTAATCAACTTTTTTATGCTTTGTCTTCCAACTTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTTCACCACGTTCCCTTAATATTTAAGATATATCACTATAACTCa                                                                                                                                          <  2:475585/140‑1 (MQ=11)

    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                                                          >  2:49018/1‑141 (MQ=255)

    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                          <  1:904008/141‑1 (MQ=18)

    aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                                                                          <  2:740273/141‑1 (MQ=18)

     ttttTAATAAGCTTTTTTATACTTTTTCTTCCAAATTTGATAGAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                         <  2:174476/141‑1 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGCACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGCTCCCTTAATATTTCAGATATATACCGATAACTcca                                                                                                                                         >  2:87586/1‑139 (MQ=11)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                         <  2:6070/141‑1 (MQ=255)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                         <  2:403747/141‑1 (MQ=255)

     tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                                                         <  1:781708/141‑1 (MQ=255)

     tttGTAATAAGCGTTTTTATGCTTTTTCTTCCAACTTTGATATAAAATTGGGGCATATTTTTGAACCCAACGGTAGACCGTTGAATGATGTACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCaa                                                                                                                                         <  2:941906/141‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                                                        <  2:297182/141‑1 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCCCAAt                                                                                                                                        >  1:900912/1‑141 (MQ=11)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                        <  1:246687/141‑1 (MQ=21)

      ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                                                        <  1:937925/141‑1 (MQ=21)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaagg                                                                                                                                       >  1:846773/1‑139 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       <  2:223622/141‑1 (MQ=255)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       >  1:454022/1‑141 (MQ=255)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       <  1:412182/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAAGATTTCAGATATATCACGATAGCTCAATg                                                                                                                                       <  1:298819/141‑1 (MQ=11)

       tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                                                                                       <  2:903491/141‑1 (MQ=39)

        gTAATAAGCTTTTTTATGCTTTTGCTTCCAAATTTGATATAAAATGAGGGCATGTTCTTGAACCCAACGGTTGCTCGATCAATGATGAAAGCTTACACCACGATCTCTTACTCTTTAAGATATATCACGATAGCTCactcc                                                                                                                                      >  2:185216/1‑137 (MQ=2)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAACTTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                     <  1:406528/141‑3 (MQ=255)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGGATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                     <  1:590879/141‑3 (MQ=12)

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

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGTTCCCTTAATATTTCAGATATATCAGGATAACTCAATgca                                                                                                                                     >  2:385080/1‑139 (MQ=255)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                                                     <  1:569164/141‑3 (MQ=12)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                     <  2:627971/141‑3 (MQ=17)

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                     <  2:226426/141‑1 (MQ=255)

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

         taataaGCTTTTTTATGCTTTTTCTTCCAAATTCGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATgca                                                                                                                                     <  1:377304/141‑3 (MQ=12)

         taagaaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACACAACGGTAGATCGTTGAATGATGAAAGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTAAATggt                                                                                                                                     >  1:677369/1‑139 (MQ=11)

          aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgcat                                                                                                                                    >  1:758683/1‑138 (MQ=11)

              aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTGACACAACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                >  2:519723/1‑141 (MQ=255)

                 ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAg                                                                                                                             >  1:649056/1‑141 (MQ=37)

                  tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                            >  2:484704/1‑141 (MQ=34)

                       tttttttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCAATGTATATCTTAGAtagta                                                                                                                       >  2:563336/5‑141 (MQ=21)

                       aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGCATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                       >  1:868803/1‑141 (MQ=25)

                             tttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGccgtt                                                                                                                 <  1:544293/141‑4 (MQ=25)

                                cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                              <  2:378390/141‑1 (MQ=35)

                                cttcCAAAGTTGATACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                              <  1:277769/141‑1 (MQ=25)

                                            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTCTCACGATAACTCAATGTATATCTTAGATAGTAGCCAGCGGCTACAGTGATaaa                                                                                                  >  1:26528/1‑140 (MQ=255)

                                            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAAAATGATGCAAAATGAACTTTTAGACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                  <  1:443933/141‑1 (MQ=11)

                                                    ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                          <  2:714914/141‑1 (MQ=32)

                                                      ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATGACATCCTTGTTa                                                                                        >  2:134463/1‑141 (MQ=17)

                                                      ggggCATATTCTTGAACCCAACGGTAGACCGCTGAATGATGAAAGTTTACGCAACGCAACTTTAATACTTTTGTTATATGCTGATAAATTGATCTAAAGCTTAAGGGGTAGCAAACTGCTACAGTTTTAACATCCTTGTTa                                                                                        >  2:785553/1‑141 (MQ=2)

                                                       gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                       <  2:84650/141‑1 (MQ=32)

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

                                                         gCATATTCTTGAACCCAACGGTAGACCGTTGATTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTACATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                     <  1:630090/141‑1 (MQ=12)

                                                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                     >  2:653764/1‑141 (MQ=32)

                                                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                     >  2:880897/1‑141 (MQ=32)

                                                           atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTg                                                                                   <  1:286313/141‑1 (MQ=17)

                                                               tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCGTTAATATTTCAGATATATCACGATAACTAAATTTATTTTTTAGTTAGTGGCTAACAGCTACAGTGGTAACATACTAGTTAAATTGTTta                                                                               >  2:287431/1‑141 (MQ=11)

                                                                     aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                         >  1:819449/1‑141 (MQ=32)

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

                                                                       ccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                       <  1:691818/141‑1 (MQ=32)

                                                                        cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                      <  2:367838/141‑1 (MQ=32)

                                                                          aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                    >  2:152787/1‑141 (MQ=32)

                                                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                  >  2:291693/1‑141 (MQ=32)

                                                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                  <  1:535179/141‑1 (MQ=32)

                                                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                  <  1:462828/141‑1 (MQ=32)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAACTCGTTCAt                                                              >  2:752821/1‑141 (MQ=17)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                              >  1:755353/1‑141 (MQ=32)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                              >  2:817193/1‑141 (MQ=32)

                                                                                gACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATTTAAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTGAATTGTTTATATCTGAAATAGTTCAt                                                              >  2:431876/1‑141 (MQ=12)

                                                                                  tcGTTGAAGGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                            <  2:701723/140‑1 (MQ=12)

                                                                                      tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                        >  2:100273/1‑141 (MQ=32)

                                                                                      tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                        >  1:224502/1‑141 (MQ=32)

                                                                                               taCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                               <  2:411153/140‑1 (MQ=17)

                                                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATGTTGTTCAGACAGAAGACTCCttttt                                             <  2:669330/141‑1 (MQ=255)

                                                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                             <  2:73262/141‑1 (MQ=32)

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

                                                                                                    ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCGTACAGAAGGCTCCTTTTTGtt                                          <  1:489987/141‑1 (MQ=17)

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

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

                                                                                                     tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCCTACAGAAGCCTCCTTTTTGTTa                                         >  2:303268/1‑141 (MQ=17)

                                                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                        >  1:858728/1‑141 (MQ=32)

                                                                                                            cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTTTTTATATTTGACATAGTTCATAAAGACGACTACTTTTTGTTTAAATtat                                  >  1:523321/1‑141 (MQ=11)

                                                                                                                cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtt                              >  1:810710/1‑140 (MQ=17)

                                                                                                                 ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                             >  1:137314/1‑141 (MQ=32)

                                                                                                                   ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAGGACTCCTTTTTGTTAACATTATACTCTaa                           >  1:74860/1‑141 (MQ=12)

                                                                                                                      atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                        <  1:864276/141‑1 (MQ=35)

                                                                                                                      atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                        <  2:367354/141‑1 (MQ=35)

                                                                                                                      atatTTCAGATATAACACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                        <  2:120770/141‑1 (MQ=25)

                                                                                                                       aatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATACTATAAATTc                       <  2:641276/140‑1 (MQ=25)

                                                                                                                           actgATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGTAGACTCCTGTTTGTTAAAATTATACTATAAATTCAACt                   <  1:187929/138‑1 (MQ=255)

                                                                                                                              gATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAATTTCAACTTTg                <  1:683885/141‑1 (MQ=255)

                                                                                                                               atatatCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc               <  1:587874/141‑1 (MQ=255)

                                                                                                                                     cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATACTATAAATTCAACTTTCCacccga         >  1:116292/1‑136 (MQ=255)

                                                                                                                                            aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat  >  2:828061/1‑141 (MQ=38)

                                                                                                                                            aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat  >  2:271117/1‑141 (MQ=38)

                                                                                                                                            aCTCAATGTATATCTAAAATAGTAGCCAACGGCTACAGTGATCACATCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat  <  1:880897/141‑1 (MQ=255)

CGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTAT  >  CP000730/36331‑36611

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