BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,174	T→C	30.3%	K126K (AAA→AAG)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,174	0	T	C	30.3%	569.1 / inf	625	K126K (AAA→AAG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (233/195); new base C (85/102); total (322/303)
Fisher's exact test for biased strand distribution p-value = 4.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.63e-02

Only 100 of 642 total aligned reads displayed.

TATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36035‑36314

tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGtt                                                                                                                                              >  2:119578/1‑141 (MQ=35)

 atGACGGTGATCTTGCTCAAGGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGttt                                                                                                                                             <  2:813996/141‑1 (MQ=21)

  tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTa                                                                                                                                            <  1:117907/141‑1 (MQ=25)

  tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTa                                                                                                                                            >  2:835126/1‑141 (MQ=255)

    gcGGTGATCTTGCTCAATGAGTTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTAcc                                                                                                                                          <  1:189159/140‑1 (MQ=11)

    aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTAcc                                                                                                                                          <  2:695532/141‑1 (MQ=21)

    aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTAcc                                                                                                                                          >  1:265534/1‑141 (MQ=255)

    aCGGTGATCTTGCTAAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGCGGTTTAcc                                                                                                                                          >  1:457263/1‑141 (MQ=255)

     cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCa                                                                                                                                         <  2:109052/141‑1 (MQ=11)

          aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAATTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAAttt                                                                                                                                    >  2:504416/1‑140 (MQ=12)

           tcacgCTCAATGAGGTAATAACAATGTTGCGATGTACAACGGCAGTCAGGTTTAAGTTTCAAAGGTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGAGCCTGATCTGTA‑ATTACCTTTGGAGGATTACCAAATTGt                                                                                                                                   >  2:385761/5‑141 (MQ=11)

           tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                   >  1:514919/1‑141 (MQ=17)

           tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                   >  2:742396/1‑141 (MQ=12)

           tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCGGGTTTAAGTTTAAGAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAATGTGCCTGATCTGTG‑ACTACCTTTTGAGGCTTAGCAAATTGt                                                                                                                                   >  2:408084/1‑141 (MQ=11)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCATGGTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                  >  1:788821/1‑141 (MQ=255)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:305584/1‑141 (MQ=255)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:675770/1‑141 (MQ=11)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAACCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                                  >  1:618625/1‑141 (MQ=11)

            cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAgttgtt                                                                                                                                  >  2:609937/1‑141 (MQ=17)

               gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGGTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                               >  1:73176/1‑141 (MQ=255)

               gCTCAATGAGGTTATTCAGATAGTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                               <  1:9171/141‑1 (MQ=12)

                  cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                            <  1:520892/141‑1 (MQ=12)

                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                                            >  1:293221/1‑141 (MQ=11)

                  cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                                            <  2:792591/141‑1 (MQ=11)

                   aaTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                                                                           <  2:324384/141‑1 (MQ=32)

                      gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGGG‑ACg                                                                                                                       >  2:5758/1‑141 (MQ=12)

                       atgTTATGCATATATTTCGATGTACAATGACAGTCATGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTTAAGTTGCCTATTCTGTA‑ACTTCCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGt                                                                                                                      <  2:602823/139‑1 (MQ=11)

                        ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAGGGTGGCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGtt                                                                                                                     <  2:341094/141‑1 (MQ=11)

                         gTTATTCAGATATTTCGATGTACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGttt                                                                                                                    <  2:457499/141‑1 (MQ=12)

                         gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGttt                                                                                                                    >  1:348580/1‑141 (MQ=17)

                          ttagtgagagATTTCGATGTACAATGACAGTGAGATTTAAGTTTAAAAGCTTTATTTATTTTCGCCATTGCTGCCTTTGTTGAAGGTGCCTGATCTGTA‑ACTAGCTTTGGAGGCTTACGAAATTGTAAAATCAG‑ACGTTTg                                                                                                                   <  2:183703/135‑1 (MQ=2)

                          ttATTCGAATATTTCGATTTACAAGGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTg                                                                                                                   <  2:444276/141‑1 (MQ=11)

                          ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTg                                                                                                                   <  2:362691/141‑1 (MQ=17)

                               cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTTACCATTTCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTGGATatg                                                                                                              >  1:81684/1‑139 (MQ=11)

                                    atTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCat                                                                                                         >  2:160577/1‑141 (MQ=12)

                                     ttgcGATGTAAAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAACTACTTTAACCATTGCTACCTTCGGTGAAGGTGCCGGATCTGGA‑ACTACCTTTTGAGGGTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACTCata                                                                                                        <  2:215974/138‑1 (MQ=11)

                                     tcacGACGTACCATGACAGTCAAGTTTAAGTTTAAAAGCTTTAAATACTTTAGCCATTGCTAACTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAC‑ACGTTTGATAAACGCata                                                                                                        <  2:752406/138‑1 (MQ=11)

                                      ttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTGTGACCATTGCTACCTTCGTTGAAGGTGACTGATCTGTA‑ATTACCTTTTGATGTTTCCCTAATTGTTTAAAGAGTACGttattaaatggatat                                                                                                        >  2:23086/1‑128 (MQ=11)

                                        cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACGTTTTGAGGTTTACCAAATTGTTTAATGAA‑GCGTTTGATAAACGCATAAGc                                                                                                     >  2:548021/1‑141 (MQ=11)

                                        cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGc                                                                                                     <  2:778947/141‑1 (MQ=17)

                                         gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCt                                                                                                    >  1:481828/1‑141 (MQ=255)

                                          aTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTg                                                                                                   >  1:2328/1‑141 (MQ=12)

                                            gTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGaa                                                                                                 <  1:419923/141‑1 (MQ=32)

                                                aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGAt                                                                                             <  1:717763/141‑1 (MQ=14)

                                                  tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           >  1:443034/1‑141 (MQ=255)

                                                   gACAGGCCGGGTTAAGTTTAAAAGCTGTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACTTTTTGAGGGTTACCAAATTGTTTAATGAG‑AGGTTTGATAAACGCATCTGCtaccggattat                                                                                          >  2:260001/1‑131 (MQ=11)

                                                                                     tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        <  2:91474/141‑1 (MQ=2)

                                                                                     tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                                        <  1:383292/141‑1 (MQ=2)

                                                                                       taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                                      <  1:804006/138‑1 (MQ=17)

                                                                                        aaccATTGCTATCTTCGTTGAAGGCGCCTGATCTGTA‑ATTACCTTTTAAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGAGCAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  2:300272/139‑1 (MQ=11)

                                                                                        aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTCCCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTACGTcc                                                     <  1:591834/139‑1 (MQ=9)

                                                                                        aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     >  2:332934/3‑141 (MQ=17)

                                                                                        aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTACGCAACCAAATATTTAaagtatgtac                                                     >  2:43232/1‑139 (MQ=11)

                                                                                           cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAACTTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  >  1:741243/1‑141 (MQ=17)

                                                                                            aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTTTCTCGTTGCTTACGCAAACAAATATCTAATGTATGTCCctct                                                 >  1:171141/1‑141 (MQ=17)

                                                                                            aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAAGCAAATATCTAATGTATGTCCctct                                                 >  2:523133/1‑141 (MQ=17)

                                                                                             ttGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTg                                                <  1:579765/141‑1 (MQ=32)

                                                                                             ttGCTACCTTCGTTGAAAGTGCCTCACCTCAG‑GTGACGGTTTGAGGGTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTCACTTACGCAACAAAATATCTAATCTAAGTCCctcta                                                <  2:173588/141‑2 (MQ=11)

                                                                                              tGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               <  2:136123/141‑1 (MQ=32)

                                                                                              tGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               >  2:277428/1‑141 (MQ=32)

                                                                                              tGCTACCTTCGTTGAAGGTGCCTG‑GCGATAGAGGACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGc                                               >  1:830313/1‑141 (MQ=11)

                                                                                               tgtggcTTCGTTGAAGGTGCCTGATCCGTA‑ACTACCTTTCGAGGCTTACCAAATTGTTTATTGAG‑TAGTTTGATAAACGCATATGCTGAATTATTATCTCGTTGCTTACGCAGCCAAATATCTAATGTATGACCCTCTgta                                              <  2:67873/136‑3 (MQ=11)

                                                                                               gcgtccTTCGGTGAGGGTGCTTGATGTGCA‑TTTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTGGATAAACGCATATGCGGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                              <  1:101720/137‑1 (MQ=11)

                                                                                                ttACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             <  2:582052/140‑4 (MQ=17)

                                                                                                cTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                             <  1:252591/141‑1 (MQ=32)

                                                                                                 tACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                                            >  1:690943/1‑141 (MQ=32)

                                                                                                  gccTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTCTCTCGTTGCTTACGCAACCAGATATCTAATGTATGTCCCTCTGTATCa                                           >  2:668396/2‑141 (MQ=11)

                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTCTGCATCa                                           >  1:135663/1‑141 (MQ=17)

                                                                                                  aCCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           <  2:115843/141‑1 (MQ=17)

                                                                                                   ccTTCGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  1:631987/141‑1 (MQ=12)

                                                                                                     ttCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCACATTTTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATg                                        <  1:373940/141‑1 (MQ=12)

                                                                                                     ttCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                                        <  1:551829/141‑1 (MQ=32)

                                                                                                       cGTTGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGc                                      >  2:124834/1‑141 (MQ=12)

                                                                                                         ttGAAGGTGGCTGCTCTGTA‑ATTAGCTTTTGAGTTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:316545/141‑1 (MQ=255)

                                                                                                         ttGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                                    >  1:754477/1‑141 (MQ=17)

                                                                                                          tGAAGGTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACg                                   <  2:558472/141‑1 (MQ=12)

                                                                                                             aGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGata                                >  2:648720/1‑141 (MQ=17)

                                                                                                              ggTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAGTGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  2:623405/1‑141 (MQ=17)

                                                                                                              ggTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  1:665852/1‑141 (MQ=17)

                                                                                                              ggTGCCTGATCTGTA‑ACTACCTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                               >  1:772190/1‑141 (MQ=12)

                                                                                                                  ccTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGAGATAAAt                           >  2:310603/1‑141 (MQ=17)

                                                                                                                  ccTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCACTCTGCATCAATGGCACGATATAAAt                           >  2:686619/1‑141 (MQ=17)

                                                                                                                    tGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCCACTATCTAATGAATGGCCCTCTACATCAATGGCACGCTATAAATac                         >  2:418195/1‑140 (MQ=11)

                                                                                                                      aTCTGTA‑AATACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  1:746443/141‑1 (MQ=17)

                                                                                                                              actACCTGTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTATTTTATGCTCCTCGGTATCAATGGCACGATATAAATAGCTCCAtttt                <  1:538199/139‑1 (MQ=11)

                                                                                                                              aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt                >  1:331475/1‑141 (MQ=17)

                                                                                                                                tACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTcc              >  1:651498/1‑141 (MQ=17)

                                                                                                                                 aCCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTAACGCTGCATAAATGTCACGATATAAATAGATCCATTTTGCt             >  2:759318/1‑141 (MQ=11)

                                                                                                                                  ccTTTTGAGGTTTACCAAATTGTTGAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  2:730044/141‑1 (MQ=17)

                                                                                                                                  ccTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCtt            <  1:127401/141‑1 (MQ=17)

                                                                                                                                   ctgttcaGGTCTTCCAAATTGTTTAATGAG‑ACGTTTGATAAACGAATATGCTGAATTATTATCTCGTTGCTTACGCAACCAAATCTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  2:474933/135‑1 (MQ=11)

                                                                                                                                   cTTTTTAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt           <  1:177128/141‑1 (MQ=17)

                                                                                                                                   cTTTTGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCttt           >  1:220048/1‑141 (MQ=12)

                                                                                                                                     tttGAGGCTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTa         <  2:481828/141‑1 (MQ=12)

                                                                                                                                      ttGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAt        >  1:464637/1‑141 (MQ=17)

                                                                                                                                       tGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTCCGCAACCAAATTTCTGTTGGAGGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtt       <  1:401333/141‑1 (MQ=11)

                                                                                                                                        gAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAATTAAATAGCTCCATTTTCCTTTTAttt      >  1:141752/1‑141 (MQ=12)

                                                                                                                                          ggTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTg    <  2:246837/141‑1 (MQ=32)

                                                                                                                                          ggTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCCCCATTTTCCTTTTATTTTg    >  2:201900/1‑141 (MQ=17)

                                                                                                                                          ggTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTACCTCTGCATCAATGGCACGATATAAGTAGGTCCATTTTCCTTTTATTTTg    >  1:708341/1‑141 (MQ=12)

TATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTA‑ATTACCTTTTGAGGTTTACCAAATTGTTTAATGAG‑ACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36035‑36314

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