BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,161	T→C	30.8%	I131V (ATT→GTT)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,161	0	T	C	30.8%	618.3 / inf	675	I131V (ATT→GTT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (229/232); new base C (85/123); total (316/359)
Fisher's exact test for biased strand distribution p-value = 3.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 687 total aligned reads displayed.

CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301

 ttttttatcttAATATGACGGTGATTTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                                                             <  2:310162/132‑1 (MQ=11)

    cttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTTAAGGTGCCTGATCTGTAATTAcc                                                                                                                                          >  1:723570/1‑141 (MQ=255)

    cttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAcc                                                                                                                                          >  2:442216/1‑141 (MQ=255)

     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCt                                                                                                                                         <  2:794747/141‑1 (MQ=255)

       aCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGCCAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTCACTACCttt                                                                                                                                       >  2:357045/1‑141 (MQ=11)

         tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                                                     <  2:251842/141‑1 (MQ=25)

             atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTGTAGGt                                                                                                                                 >  2:663500/1‑141 (MQ=255)

                tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTCATTACCTTTTGAGGTTTa                                                                                                                              >  2:76227/1‑141 (MQ=255)

                tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACATTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                              >  1:186758/1‑141 (MQ=255)

                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCGAGTTCAAGTTTAAAAGCCATAATTACTTTAACCATTGCTACATTCGTAGGAGGTGCCTGATCTGTAATTACCTTTTGAGGTGTAcc                                                                                                                            >  2:177996/1‑141 (MQ=11)

                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                            >  1:604454/1‑141 (MQ=21)

                  aCGGTGATCTTGCGCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCATGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                            <  2:320354/141‑1 (MQ=255)

                   cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCa                                                                                                                           <  2:150384/141‑1 (MQ=11)

                    ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaa                                                                                                                          <  1:912971/141‑1 (MQ=16)

                    ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGTTTACCaa                                                                                                                          >  2:190530/1‑141 (MQ=255)

                     gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTTAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                         <  2:723570/141‑1 (MQ=12)

                     gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                         <  2:605808/141‑1 (MQ=12)

                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCCGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                        <  1:169407/141‑1 (MQ=12)

                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                        <  2:648897/141‑1 (MQ=17)

                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                        <  1:769822/141‑1 (MQ=17)

                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTGTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                        >  2:826752/1‑141 (MQ=12)

                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGGGCCTGATCTGTAAGTACCTTTTGAGGTTTACCAAAt                                                                                                                        >  1:12223/1‑141 (MQ=17)

                       gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAAtt                                                                                                                       <  2:552643/141‑1 (MQ=11)

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTCCCAAATTGt                                                                                                                     >  1:283105/1‑141 (MQ=12)

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                     >  1:157282/1‑141 (MQ=17)

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGt                                                                                                                     <  1:594682/141‑1 (MQ=11)

                          cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                    >  2:526952/1‑141 (MQ=11)

                          cTTACTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGtt                                                                                                                    >  1:655990/1‑141 (MQ=11)

                              gtCAATGAGGTTATTCGGATATTTCGATGTACAATGACAGTCAGGTTTAAGGTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAAt                                                                                                                <  1:316647/140‑1 (MQ=11)

                              cTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                <  2:126667/141‑1 (MQ=12)

                                   tGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                                                                                                           >  1:737577/1‑141 (MQ=17)

                                    gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTATTGAGACg                                                                                                          >  1:234659/1‑141 (MQ=12)

                                     aGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                                                                                                         <  1:443446/141‑1 (MQ=17)

                                      ggTTATTCATATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                        <  1:367494/141‑1 (MQ=12)

                                      ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTGACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                        >  1:572785/1‑141 (MQ=12)

                                      ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                        <  1:326041/141‑1 (MQ=17)

                                       gTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGttt                                                                                                       <  2:573024/141‑1 (MQ=11)

                                        ttgttCAGATATTTCGAGGTAAAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTAGCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                      <  2:207009/138‑1 (MQ=11)

                                         tATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                     <  1:783056/141‑1 (MQ=17)

                                          aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACTATTGCTACCTTCGTTGAAGGTGACTGATCTGTAATTAAGTGATTAGATTTACCAAATTGTTTAATGAGAGGTTTGAt                                                                                                    >  1:471423/1‑141 (MQ=11)

                                             cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTCACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGGTTAATGAGACGTTTGATaaa                                                                                                 >  2:862118/1‑141 (MQ=255)

                                             cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                 <  2:267958/141‑1 (MQ=17)

                                             cAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                 >  2:370758/1‑141 (MQ=255)

                                               gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATACACg                                                                                               >  2:913621/1‑141 (MQ=255)

                                               gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGACGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                               >  2:334710/1‑141 (MQ=12)

                                                  atTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                            >  2:673532/1‑141 (MQ=32)

                                                   tttCGATGTACCATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGTGGATTACCAAATTGTTAAATTAGACGTTTGATAAACGCata                                                                                           <  1:583817/141‑1 (MQ=11)

                                                   tgtCGAGGTAGAATGACAGTCAGGTTTACGTGTAAAAGCTTTAATGATTTTAGCCATTGCTGCCCTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCata                                                                                           <  1:518513/139‑1 (MQ=11)

                                                    ttCGATGTACAATGGCAGTCCGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCACTGCCACCTTCGTTGAAGGAGCCCGATCTGTAATTACCTTGTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                          <  2:28125/141‑1 (MQ=11)

                                                     tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                         <  2:739587/141‑1 (MQ=17)

                                                       gATTTAAAATGGCAGTCATGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCCTCGTTGAAGGTGCCGGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                       <  1:894953/141‑1 (MQ=11)

                                                       gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt                                                                                       >  2:595191/1‑141 (MQ=12)

                                                        aTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                      >  2:475995/1‑141 (MQ=12)

                                                           tACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAAt                                                                                   >  2:349605/1‑141 (MQ=12)

                                                            aCAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                  <  2:468502/141‑1 (MQ=255)

                                                              aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                <  1:645956/141‑1 (MQ=11)

                                                                 gACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATACACGCATATGCTGAATGATTAt                                                                             >  1:45048/1‑141 (MQ=11)

                                                                                                   tttAGCCATTGCTATCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           <  2:891007/141‑1 (MQ=11)

                                                                                                      tgacaatgctgtcTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        <  1:266465/129‑1 (MQ=11)

                                                                                                      aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        >  2:847193/1‑141 (MQ=255)

                                                                                                        ccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAGTGAGACGTTTGATAAACGCATATGCTGCATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCct                                      >  1:298984/1‑141 (MQ=17)

                                                                                                         cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     >  2:679852/1‑141 (MQ=32)

                                                                                                          aTTGCTTCCTTCGTTGAATGTGCCTGCTTTGTACTCATCTTTTGAGGTTTATCAAATTGTTTAATGTGATGTTTGATAAACGTTGGCTTTTAATTATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                    <  1:213924/141‑1 (MQ=2)

                                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                    <  2:634596/141‑1 (MQ=32)

                                                                                                            tgctgcCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                  <  2:38699/141‑2 (MQ=12)

                                                                                                            tgctgcATTCTTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                  <  2:220941/141‑2 (MQ=11)

                                                                                                            tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTCTGTCCCTCTGc                                  >  1:864922/1‑141 (MQ=17)

                                                                                                             gCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCTACCATAGATCTAATGTATGTCCCTCTgta                                 <  2:416902/141‑3 (MQ=11)

                                                                                                              cTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCAt                                >  1:824613/1‑141 (MQ=32)

                                                                                                               tgccTTCGTTGAAGGTGCCTGATGTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                               <  2:43264/139‑1 (MQ=11)

                                                                                                               tgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                               <  2:111847/139‑1 (MQ=255)

                                                                                                                 ccTTCATTGAAGGTGCCTGTTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCaa                             <  1:872324/141‑1 (MQ=17)

                                                                                                                  cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATTAGACGTGTGATAAACGCAGATGCTCAATGATTATCGCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                            >  2:620/1‑141 (MQ=11)

                                                                                                                  cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCagt                            >  2:40838/1‑139 (MQ=17)

                                                                                                                  cTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAAt                            >  1:810606/1‑141 (MQ=32)

                                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                           <  2:582527/141‑1 (MQ=32)

                                                                                                                    tCGTTGAATGTGCCTGATTTGTAATTACCTGTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATgg                          <  2:904747/141‑1 (MQ=12)

                                                                                                                     cGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGc                         <  1:838009/141‑1 (MQ=17)

                                                                                                                       ttGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCAc                       >  1:69394/1‑141 (MQ=12)

                                                                                                                         gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                     <  1:348506/141‑1 (MQ=32)

                                                                                                                         gAAGGTGCCTGATCTGGAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                     <  1:511544/141‑1 (MQ=17)

                                                                                                                          aaGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGat                    >  1:355484/1‑141 (MQ=32)

                                                                                                                            ggTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                  >  1:10389/1‑141 (MQ=17)

                                                                                                                            ggTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACTCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATAAATGGCTAGatat                  >  1:627834/1‑141 (MQ=11)

                                                                                                                            ggTGCCAGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACTTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGGAGGTGCATCTGTATCAATGGCACGagat                  >  2:502861/1‑141 (MQ=11)

                                                                                                                              tGCCTGATCTGTAACGACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATaa                >  1:675334/1‑141 (MQ=11)

                                                                                                                               gacTGATCGATAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATaaa               <  1:386134/139‑1 (MQ=12)

                                                                                                                                ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt              >  2:707751/1‑141 (MQ=255)

                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCGTCTGTATCAATGACACGATATAAATAg            >  2:166334/1‑141 (MQ=12)

                                                                                                                                   gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTCTGTATCAATGGCACGTTATAAATAGc           >  2:493250/1‑141 (MQ=11)

                                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCATTGGCACGTTATAAATAGCt          >  1:413817/1‑141 (MQ=12)

                                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCACTCTGTATCAATGGCACGATGTCAATAGCt          >  2:695715/1‑141 (MQ=255)

                                                                                                                                     tCTGTAATTACCTTTTGACGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTATGCAACAAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc         <  1:1335/141‑1 (MQ=255)

                                                                                                                                     tCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTTAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTc         <  2:779600/141‑1 (MQ=255)

                                                                                                                                      cTGTAACTACCTTTTGAGGCTTCCCAAATTGTTTAATGATACGTTTGATACACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTcc        <  1:489382/141‑1 (MQ=11)

                                                                                                                                       tGTAATTACTGTTTGAGGTTTACCAATTTGTTTAATGTGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa       <  1:747370/141‑1 (MQ=255)

                                                                                                                                       tGTAACTACCTTTTGAGGCTTATCAAATTGTTTAATGAGACGTTAGCTTTATTGGTAGGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa       <  1:253229/141‑1 (MQ=11)

                                                                                                                                        gTAATTATCTTTTGTGGTTTGCTTAGTTGTTTAATGAGACGTTTGATGAAGGCTTATGCTGAATGATTGTCTCGTTGCTTACGCAACCAAATTTCTAATGTATGTCCCTCTGCATCAGTGGCACGATATAAATAGCTCCAt      <  1:54379/141‑1 (MQ=11)

                                                                                                                                         tAACTACCTTTTGAGGCTTACCAAATAGTTTAATGAGACGTTTGATAAACGCATATGCTGAATCATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGAACGCTATAAATAGATCCact     >  2:545969/1‑139 (MQ=11)

                                                                                                                                           aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt   >  2:80426/1‑141 (MQ=17)

CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301

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