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%	111.7 / inf	342	R36R (AGG→AGA)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base C (84/92); new base T (84/79); total (171/171)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 347 total aligned reads displayed.

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTA  >  CP000730/36304‑36583

   tttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                           <  2:407037/140‑1 (MQ=39)

   tttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                           <  2:416068/140‑1 (MQ=11)

   aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                           <  2:475625/141‑1 (MQ=255)

    ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa                                                                                                                                          <  1:951470/141‑1 (MQ=18)

    ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTa                                                                                                                                          <  2:750102/141‑1 (MQ=11)

    ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTa                                                                                                                                          <  1:126474/141‑1 (MQ=255)

     tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCGAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa                                                                                                                                         <  2:699041/141‑1 (MQ=255)

     tttGATGTACCTCTCATCAATACGCCACTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGCGCATATTCTTGAACCCA‑ACCGTAGACCGTTGAATTATGAACGTTTACACCACGTTCCCTTaa                                                                                                                                         <  2:808891/141‑1 (MQ=14)

     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                                                                                                         <  2:21935/141‑1 (MQ=255)

      ttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAat                                                                                                                                        <  2:784533/141‑1 (MQ=11)

       tGATTTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAata                                                                                                                                       <  2:127791/141‑1 (MQ=18)

        ggcaaggatctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCCTAatat                                                                                                                                      >  2:275609/9‑141 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGTTCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                                      >  1:619626/5‑141 (MQ=11)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                                      >  2:334003/5‑141 (MQ=255)

        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTAGGGCATATTCTTGAACCCA‑ACGTTAGATCGTTGAATGATGAACGTTTACAACACGTTCTCTTAatat                                                                                                                                      >  1:156234/5‑141 (MQ=11)

         atatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACTACGTTCTCTTAATAtt                                                                                                                                     >  1:879967/4‑141 (MQ=11)

         aTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAtt                                                                                                                                     >  1:559107/1‑141 (MQ=21)

            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCCCGTTCACTTAATATTTCa                                                                                                                                  >  2:622053/1‑141 (MQ=11)

            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAGAAATACAAT‑CATCTTGGAATGATGAACGATTGCACCACGTTCCCTTAATATTTCa                                                                                                                                  >  1:1037395/1‑141 (MQ=11)

            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                                                  >  1:976893/1‑141 (MQ=40)

            aaCGTCTCATCAATACGCCATCTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATACTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                                                                                                  <  1:567834/140‑1 (MQ=11)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTTTTTTTTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTCCACCTCTTTCTCTTAATATTTCAGa                                                                                                                                >  2:31529/1‑141 (MQ=11)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                                                >  1:232992/1‑141 (MQ=255)

              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                                >  2:972092/1‑141 (MQ=21)

                gctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGata                                                                                                                              <  2:1060084/140‑1 (MQ=40)

                 ttcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                             <  2:194425/140‑1 (MQ=255)

                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                             <  2:548382/141‑1 (MQ=255)

                    ttcaATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  1:795783/140‑1 (MQ=255)

                    cacaATACGCCATGTTTAATAAGCTTTCTTTTGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCATATTCTTGAACCCA‑CCCGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  2:919179/139‑1 (MQ=11)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATc                                                                                                                          <  1:382823/141‑1 (MQ=255)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  2:442966/141‑1 (MQ=18)

                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATc                                                                                                                          >  2:843246/1‑141 (MQ=255)

                    atcaATACGCCAATTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAGTGGAGATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGATCCCCTAATATTTCAGatatcta                                                                                                                          >  2:358341/1‑138 (MQ=11)

                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCCTATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCGGATATATCa                                                                                                                         >  2:186123/1‑141 (MQ=11)

                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCa                                                                                                                         >  1:984379/1‑141 (MQ=21)

                     acaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                         <  2:869601/140‑1 (MQ=255)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAACTTGGGGCATCTTCTTGAACCCA‑AAGGTAGCTCGTTGAATGTTGGAGCTTTACACTGCGTTGTTTTAATATTTCTGatttataag                                                                                                                        >  2:466540/1‑138 (MQ=2)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                                                        >  2:144167/1‑141 (MQ=255)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATAAc                                                                                                                        >  2:943744/1‑141 (MQ=255)

                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAc                                                                                                                        <  2:247809/141‑1 (MQ=39)

                       aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGGTCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACg                                                                                                                       >  2:937515/1‑141 (MQ=11)

                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                                                      >  1:420394/1‑141 (MQ=255)

                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACATTTACACCACGTTCTCTTAATGTTTCAGATATATCACGa                                                                                                                      >  1:112365/1‑141 (MQ=11)

                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGa                                                                                                                      <  2:849154/141‑1 (MQ=18)

                         tACGCCATTTGTAATAAGCTTTTTTTTTTTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTTTTGAACGCA‑ACGGTAGATCGTTGAATGCTGAACGTGTACACCACGTTCTCTTAATATTTAAGATATGTTACGAt                                                                                                                     >  2:800466/1‑141 (MQ=11)

                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                                                     <  1:446464/141‑1 (MQ=255)

                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                                                     >  1:259792/1‑141 (MQ=255)

                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGAt                                                                                                                     >  1:652192/1‑141 (MQ=21)

                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGATTACACCACGTTCCCTTAATATTTCAGATATATCACGAt                                                                                                                     >  1:474862/1‑141 (MQ=18)

                          cccactcttgaactaaGCTTTTTCATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATa                                                                                                                    <  2:727628/128‑1 (MQ=11)

                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTACTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTGCTCTTAATATTTCAGATGTATCAGGATa                                                                                                                    >  1:434471/1‑141 (MQ=11)

                          aCGCAATTTGTACTAAGCTCTTTTATGTTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTAATATTTCAGATATATCACGATa                                                                                                                    <  1:677009/141‑1 (MQ=11)

                           cGCCATTTGTAATAAGCTTTTTTTTTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGATATATCCCGATaa                                                                                                                   >  2:654033/1‑141 (MQ=11)

                           cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                   <  1:854560/141‑1 (MQ=18)

                            gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                                                  <  1:751357/141‑3 (MQ=255)

                            gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                  <  2:278061/141‑1 (MQ=39)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATGTTCTTGAACCCA‑ACGGTAGACCGTTGAATGGTGAACGTTTACACCCCGTTCCCTTACTATTTCAGATATATGACGATAACt                                                                                                                 >  2:582954/1‑141 (MQ=11)

                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                 >  1:908149/1‑141 (MQ=255)

                              cATTTGTCATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCAGGTCCCCTTACTGTTTCTGATATATCACGATCACTc                                                                                                                >  2:596130/1‑141 (MQ=11)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                >  2:563539/1‑141 (MQ=255)

                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                <  1:54976/141‑1 (MQ=25)

                               aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                               <  1:846625/141‑1 (MQ=255)

                                tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                              <  2:738573/141‑1 (MQ=11)

                                tctGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                              <  1:281844/139‑1 (MQ=11)

                                 ttttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             <  1:247912/138‑1 (MQ=11)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             <  1:413385/141‑1 (MQ=255)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             <  2:451898/141‑1 (MQ=11)

                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGACCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                             <  1:703768/141‑1 (MQ=18)

                                  tGTAATAAGGTTTATTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  1:49450/141‑1 (MQ=11)

                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  1:1008278/141‑1 (MQ=255)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                           <  1:1001540/141‑2 (MQ=11)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                           <  1:105569/141‑2 (MQ=11)

                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATAAATCACGATAACtaccaga                                                                                                           >  1:392533/1‑135 (MQ=255)

                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCATTgca                                                                                                          >  1:303607/1‑139 (MQ=11)

                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAAAGTTTACACAACGATCTCTTAATATTTCAGATATATCACGATAGCTCaaaaca                                                                                                          >  2:897913/1‑137 (MQ=11)

                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTAGGTCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                          <  1:344916/141‑3 (MQ=11)

                                    taataaGCTTTTTTATGCTTTTTCCTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTATTATTTTACATATATCACGATAACTCAATgca                                                                                                          <  1:514735/141‑3 (MQ=11)

                                     aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTGTTAATATTACGGATATATCCAGATAGCTCGATgcaa                                                                                                         >  2:1048947/1‑138 (MQ=11)

                                     aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATGCCTCAATgcaa                                                                                                         >  1:20341/1‑138 (MQ=255)

                                            ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAg                                                                                                 >  1:661431/1‑141 (MQ=35)

                                            ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATCACGATAACTCAATG‑TATATCTTaa                                                                                                 >  1:321067/1‑140 (MQ=25)

                                                  aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAACTCAATG‑TATATCTTAGAtagta                                                                                           >  1:191467/1‑141 (MQ=25)

                                                       ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCaa                                                                                      <  2:1004698/141‑1 (MQ=35)

                                                        tttCCTCCAAATTTGATACAAAATTGGTGCTTATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGccgtt                                                                                     <  1:35865/141‑4 (MQ=255)

                                                             tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTTGCCAACGGCTa                                                                                >  1:586296/1‑141 (MQ=25)

                                                             tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTa                                                                                >  1:1025172/1‑141 (MQ=35)

                                                                  aTTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGTTAGACGGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                           <  2:515535/141‑1 (MQ=21)

                                                                    ttGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                         <  2:506800/141‑1 (MQ=35)

                                                                       aTACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                      >  1:716115/1‑141 (MQ=35)

                                                                              aTTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAAGG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                               <  1:728073/141‑1 (MQ=17)

                                                                               ttcgggaagagtCTTTAACCCA‑ACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                              <  1:87006/130‑1 (MQ=11)

                                                                                tGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                             <  1:147010/141‑1 (MQ=32)

                                                                                      atatTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGtacaaacg                                                       >  2:630935/1‑134 (MQ=255)

                                                                                        atTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                     <  1:773464/141‑1 (MQ=32)

                                                                                                              ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                <  1:66892/141‑1 (MQ=32)

                                                                                                                    aaTGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                          <  2:692814/141‑1 (MQ=17)

                                                                                                                          gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                    <  2:686530/141‑1 (MQ=32)

                                                                                                                            aCGTTCACACCACGTTCGCTTAATATTACTGATATATCA‑GCCCAGTTAACGTTAGAGGTTAGTTAATATCGATCTGCTACACTGAGTACCTCCTTGTTAAATTGTTTATATCTGAAATATTTCATACAGAAAACTCCtttt                  >  2:382093/1‑141 (MQ=2)

                                                                                                                                     ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaag         >  2:645991/1‑140 (MQ=17)

                                                                                                                                          ttgccTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc    <  2:296580/138‑1 (MQ=17)

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTA  >  CP000730/36304‑36583

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