BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,161	T→C	23.2%	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	23.2%	630.4 / 264.1	487	I131V (ATT→GTT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (181/189); new base C (45/67); total (228/259)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 494 total aligned reads displayed.

CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301

 tttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAtt                                                                                                                                             <  1:87309/141‑1 (MQ=255)

  ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                                                            >  1:404976/1‑141 (MQ=255)

  ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                                                            >  2:280861/1‑141 (MQ=255)

  ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                                                            <  2:355392/141‑1 (MQ=255)

  ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                                                            >  2:577497/1‑141 (MQ=255)

  ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                                                                                            <  1:560951/141‑1 (MQ=37)

   tcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGCTATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAc                                                                                                                                           >  1:64111/1‑141 (MQ=255)

   tcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAc                                                                                                                                           >  1:442320/1‑141 (MQ=255)

     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCt                                                                                                                                         <  1:358975/141‑1 (MQ=255)

       aCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCttt                                                                                                                                       >  2:191742/1‑141 (MQ=255)

       aCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCttt                                                                                                                                       <  2:396600/141‑1 (MQ=255)

        cTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAAGTTTAAGTTTAAAAGCTTTACTTACTTTAACCATTGCTACCTTCGTTGCAGGTGCCTGATCTGTAATTACCtttt                                                                                                                                      >  1:421251/1‑141 (MQ=11)

        cTTTAATATGACGGTGATCTTGCCCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCtttt                                                                                                                                      <  1:286649/141‑1 (MQ=11)

         tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                                                     <  2:615833/141‑1 (MQ=255)

         tttAATAAGACGGTGATCTAGCTTAATGAAGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                                                                                                     <  2:247987/141‑1 (MQ=11)

          ttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGa                                                                                                                                    <  2:43239/141‑1 (MQ=35)

            tatatGACGGTGATCTTGCTCAATGAGGTTATTCAGGTATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAgg                                                                                                                                  <  1:234775/140‑1 (MQ=255)

             atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTTAAGGTGCATCATCTTTAATTACTTTTTGAGGt                                                                                                                                 >  2:417812/1‑141 (MQ=11)

             atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTCCGTCGTAGGAGCCTGATCTGTAACTACCTTTTGAggc                                                                                                                                 <  1:660815/141‑2 (MQ=11)

              tatGGCGGTGATCTTGCTCAATGAGGTTATTCAGATCTTTCGATGTACAATGACAGTCAGGTTTAAGTTTAACAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGGTCTGTAACTACCTTTTGAggct                                                                                                                                <  2:338743/141‑3 (MQ=11)

                tGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTa                                                                                                                              <  1:507716/141‑1 (MQ=35)

                 gccGGTGATCTTTCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  2:536313/139‑1 (MQ=255)

                 gACGGTGATCTTGGTCAAGGAGTTTATTCAGATATTTCGATGGACAAGGACAGTTAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCATTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTAc                                                                                                                             <  2:275585/141‑1 (MQ=11)

                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTATCGTCGGCCACGTCAGCGGCAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  1:157430/141‑1 (MQ=11)

                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  1:568524/141‑1 (MQ=35)

                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  1:69154/141‑1 (MQ=34)

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

                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                            <  1:627492/141‑1 (MQ=25)

                        aTCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTg                                                                                                                      <  1:516434/141‑1 (MQ=12)

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

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

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGGTTACCAAATTGt                                                                                                                     >  2:120357/1‑141 (MQ=12)

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                     >  2:569943/1‑141 (MQ=32)

                         tCTTACTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTACTTACCTCTGCAGATTTACCAAATTGt                                                                                                                     >  2:112103/1‑141 (MQ=11)

                           ttGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAAATACCTTTTTTGGCTTACGAAAGTGttt                                                                                                                   >  2:669417/1‑141 (MQ=11)

                            tGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTa                                                                                                                  >  2:516704/1‑141 (MQ=12)

                            tGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTACAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTa                                                                                                                  <  1:384544/141‑1 (MQ=11)

                             gCTCAATGAGGTTATTCAAATATTTGCATGGACAAGGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                 <  1:619131/141‑1 (MQ=11)

                                cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                              <  2:626279/141‑1 (MQ=11)

                                cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATga                                                                                                              <  1:654365/141‑1 (MQ=255)

                                 aaTGAGGTTATTCAAATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgag                                                                                                             <  2:29749/141‑1 (MQ=12)

                                  aTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTCACCAGTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATgaga                                                                                                            <  1:256542/141‑1 (MQ=11)

                                  aTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTAGCCTAAGGGGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATgaga                                                                                                            >  1:455198/1‑141 (MQ=11)

                                    gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACg                                                                                                          <  2:542697/141‑1 (MQ=17)

                                      ggTTATTCAGATATTTCGATGTACAATGAAAATCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                        >  1:323821/1‑141 (MQ=17)

                                        ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                      >  1:175445/1‑141 (MQ=17)

                                        ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTCTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGCAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGTGCCGTTTg                                                                                                      >  2:325624/1‑141 (MQ=11)

                                            tCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaa                                                                                                  <  2:411092/141‑1 (MQ=32)

                                             cAGATATTTCGATGTACAAAGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCGGTTGAGGTTTACCAAATTGTTTAATGAGACGCTTGATaaa                                                                                                 >  2:476432/1‑141 (MQ=255)

                                             cAGATATTTCGATATACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                 >  2:490478/1‑141 (MQ=11)

                                                atatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                              >  1:540961/1‑141 (MQ=17)

                                                atatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCGGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                              <  2:1299/141‑1 (MQ=12)

                                                atatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                              >  1:574863/1‑141 (MQ=11)

                                                atatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGc                                                                                              >  2:155730/1‑141 (MQ=255)

                                                 tatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                             <  2:89689/141‑1 (MQ=17)

                                                 tatTTCGAGGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                             <  2:179861/141‑1 (MQ=12)

                                                  atTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATATGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCat                                                                                            <  2:240526/141‑1 (MQ=12)

                                                    tttgATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                          <  1:530245/138‑1 (MQ=12)

                                                     tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTTCCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                         <  2:152169/141‑1 (MQ=11)

                                                        aTGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                      >  2:351301/1‑141 (MQ=11)

                                                         tGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGTCTTTGCTGTCTTCGTTGCAGGTGCCTGACCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTCATGAGACGTTTGATAAACGCCTATGCTGa                                                                                     <  2:210930/141‑1 (MQ=11)

                                                         tGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                     <  1:102195/141‑1 (MQ=32)

                                                         tGTACAATGACAGGCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGa                                                                                     <  1:37548/141‑1 (MQ=17)

                                                          gTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGaa                                                                                    <  2:457603/141‑1 (MQ=17)

                                                            gcactGACAGTCAGGTTCAAGCTTAAAAGCTTTACTTATTTCAGCTATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTCATAAACGCATATGCTGAATg                                                                                  <  1:468713/137‑1 (MQ=11)

                                                                tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                              <  1:547096/141‑1 (MQ=31)

                                                                                     aatgCTTTAAGTATTTTAACCATTTCTTTCTTCTTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCaa                                                         <  2:317231/138‑1 (MQ=1)

                                                                                                   tttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatg                                           >  2:592702/1‑141 (MQ=2)

                                                                                                    ttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                          <  1:296118/141‑1 (MQ=9)

                                                                                                      aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAACCGCATATGCTGAATGATTATCTCGTTGCTTACGCACCCAAATATCTAATGCATGTcc                                        <  2:661282/139‑1 (MQ=12)

                                                                                                      aGCAATTGCTACATTCTTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        <  2:377344/141‑1 (MQ=12)

                                                                                                       gCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTtcc                                       >  1:130167/1‑141 (MQ=255)

                                                                                                         cATTGCTACCTTCGTTGAAGGTGTCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  1:244566/141‑1 (MQ=17)

                                                                                                         cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  1:226690/141‑1 (MQ=32)

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

                                                                                                             tctgccTTCTTTGAATGTTCCTGATCTTTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgta                                 <  2:198381/137‑3 (MQ=11)

                                                                                                             tccaCCTTCGTTGAAGGTGCCTGATATGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGTGACGTTTGATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                 <  2:522798/138‑1 (MQ=255)

                                                                                                             gCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                 <  1:601273/141‑1 (MQ=32)

                                                                                                                gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGCCGGTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTACctctttatac                              >  2:440432/2‑135 (MQ=11)

                                                                                                                gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                              >  1:478424/2‑141 (MQ=255)

                                                                                                                aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                              <  1:588281/141‑1 (MQ=17)

                                                                                                                aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                              <  2:80064/141‑1 (MQ=32)

                                                                                                                aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                              >  1:192524/1‑141 (MQ=32)

                                                                                                                  cttcaacgAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTTATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCCAt                            >  1:308408/8‑141 (MQ=255)

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

                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTTCGCAACCAAATATCTAATGTATGTCCCTCTTCAGCAATGGCa                        >  2:268400/1‑141 (MQ=12)

                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTCTGTCGCTCTGCATCAATGGCa                        <  1:342348/141‑1 (MQ=17)

                                                                                                                      gTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTCTCAATGCCa                        >  2:101789/1‑141 (MQ=11)

                                                                                                                         gAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGa                     >  1:112321/1‑141 (MQ=32)

                                                                                                                          aaGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCCTCAATGGCACGat                    >  2:44918/1‑141 (MQ=17)

                                                                                                                           aGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATCAACGCATACGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCTTCTGAATCAATGGCACGatt                   >  2:598800/1‑140 (MQ=255)

                                                                                                                           aGGTGCCTGATCAGTAATAACCTTTTGAGGTTTACAAAATGGTTTAATGAGACGTTGGATAAACGCATATGCTGAATGATTATCTCGTTGCTGACGCAACCAAAGATCCAATGTATGACCCTCTGCATCAATGACAGGatt                   >  2:648985/1‑140 (MQ=11)

                                                                                                                            ggTGCCTGATCTGTAATTATCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGATTATGCTGAATGATTATCTCGTTTCTTACGCAACCAAATATCTAATTTATGTCCATGTTTATTAATGGCAAGttat                  >  1:427378/1‑141 (MQ=11)

                                                                                                                             gTGCCTGATCTGTAACTACCTCTTGAGGCTTACCAAATTGTTTAATGAGACGTCTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGAtata                 <  2:218699/141‑1 (MQ=11)

                                                                                                                                ccTGAGCTGTAATTACCTTTTGCGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCGTCTGCGCAATGATTATCTCGTTGCTTAGGCAACCAAACATCTAATGTATGTCCCTCTTCATCAATGGCACGATATaaaa              >  2:202912/1‑140 (MQ=11)

                                                                                                                                  tGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGGTTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCGCTCTGCATCGATGGCACGaaattacaag            >  2:548927/1‑135 (MQ=11)

                                                                                                                                   taTCTGTAACTACCTTTTGAGGCTTACCCAATTGTTTAATGAGATGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACTAAATTTCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc           <  2:141570/140‑1 (MQ=11)

                                                                                                                                    aTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCt          <  2:503960/141‑1 (MQ=32)

                                                                                                                                       tGTACTTATCTTTTGAGGTTTTCAAAATTGTTTAATGAGACGTTTGATACACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAATTATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCa       <  1:201574/141‑1 (MQ=11)

                                                                                                                                         tAATTACCTTTTGAGGTTTATCAAATTGTTTAATGAGACGTTTGGTAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAAGCGAATATCTAATGTATGTCCCACTGCATCAATGGCACGATATAAATAGCTCCAtt     >  2:366743/1‑141 (MQ=11)

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