BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	36,161	T→C	32.0%	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	32.0%	374.9 / inf	431	I131V (ATT→GTT)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (140/150); new base C (75/62); total (216/215)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 436 total aligned reads displayed.

CTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTC  >  CP000730/36021‑36301

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

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

     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTCCCt                                                                                                                                         >  1:143219/1‑141 (MQ=255)

     ttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCt                                                                                                                                         >  1:23956/1‑141 (MQ=255)

       acttcgatatGAAGGTGATCTTGCTCAATGAGGTTACTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTACAAGCTTTAATTACTTTACCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCttt                                                                                                                                       <  1:95860/135‑1 (MQ=11)

           tAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTGACTACCTTTTGAg                                                                                                                                   >  1:10135/1‑141 (MQ=255)

            aaTATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAgg                                                                                                                                  <  1:295419/141‑1 (MQ=255)

            aaTATGACGGTGATCTTGCTCAATGAGGGTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTATAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAgg                                                                                                                                  >  1:246777/1‑141 (MQ=11)

              tatGACGGTTACCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAggct                                                                                                                                <  1:43884/141‑3 (MQ=11)

               atGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGttt                                                                                                                               <  1:261537/141‑1 (MQ=255)

                 gacggtttctttccTCAATGAGGTTATTAAGATATTTCGATGCACAATGGCAGTCATGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  2:347503/128‑1 (MQ=11)

                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  1:366709/141‑1 (MQ=25)

                 gACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAc                                                                                                                             <  1:3377/141‑1 (MQ=21)

                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGTAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                            >  2:361634/1‑141 (MQ=255)

                  aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGCTGCCTGAACTGTAATTACGTTTCGCGGATTAcc                                                                                                                            >  2:46234/1‑141 (MQ=11)

                   cGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCa                                                                                                                           <  1:85201/141‑1 (MQ=18)

                    ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaa                                                                                                                          <  1:205183/141‑1 (MQ=33)

                     gttaTCTTGCTTAATGAGGTTATTCAGATATTTCGATGTACAATGTCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                         <  2:254915/138‑1 (MQ=255)

                     gTGATCTTGCCCAAGTAGGTTATTCAGATATTTCTCTGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                                                                                                         <  2:302710/141‑1 (MQ=11)

                      tGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAt                                                                                                                        >  2:250325/1‑141 (MQ=255)

                       gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAtt                                                                                                                       >  2:27439/1‑141 (MQ=12)

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                     <  2:330339/141‑1 (MQ=17)

                         tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGt                                                                                                                     >  1:250388/1‑141 (MQ=255)

                          cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                    >  2:197913/1‑141 (MQ=12)

                            tGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTa                                                                                                                  <  2:170267/141‑1 (MQ=11)

                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                 <  1:53574/141‑1 (MQ=11)

                             gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTaa                                                                                                                 >  2:340000/1‑141 (MQ=255)

                               tCAATGGGGTTATTCAGATATTTCGATGGACATTGACAGTCAGTTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCTTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                               <  2:79936/141‑1 (MQ=11)

                               tCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                                                                                               <  2:637/141‑1 (MQ=17)

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

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

                                    gAGGTTATTCAGATATTTCGATGGACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGTTTGAGCAAATTGTTTAATGAGACg                                                                                                          <  2:176341/141‑1 (MQ=11)

                                      ggTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                        >  1:235151/1‑141 (MQ=17)

                                       gTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGttt                                                                                                       <  1:19621/141‑1 (MQ=17)

                                       gTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTTCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTACTGATACGttt                                                                                                       >  1:173427/1‑141 (MQ=11)

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

                                         tATTTAGATTTTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                     <  1:7485/141‑1 (MQ=255)

                                         tATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                     <  1:264742/141‑1 (MQ=11)

                                              aGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                                                                <  2:107856/141‑1 (MQ=17)

                                              aGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAAc                                                                                                <  1:342130/141‑1 (MQ=255)

                                               gATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                               >  2:272516/1‑141 (MQ=12)

                                                 tatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                             >  1:154000/1‑141 (MQ=11)

                                                  atTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGACGTTTACCAAATTGTTTAATGAGCCGTTTGATAAACGCat                                                                                            >  2:134979/1‑141 (MQ=255)

                                                   tttCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCata                                                                                           >  1:156172/1‑141 (MQ=17)

                                                    ttCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCatat                                                                                          <  2:199249/141‑1 (MQ=32)

                                                     tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACGACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                         <  2:297202/141‑1 (MQ=11)

                                                     tCGATGTACAAGGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAACTTTTGAATGAGACGTTTGATAAACGCATATg                                                                                         <  1:91244/141‑1 (MQ=11)

                                                      ctctGTAAAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTTATTTCTTCAGCCATTGCTACGTTCGTTGAAGGTGCCTTAACTTTAATTACCTTTCTAGGTTTACCAAATTTTTTAATGAGACGTTTGATAAACGCGTATGc                                                                                        <  2:75182/138‑1 (MQ=11)

                                                      cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  1:199667/141‑1 (MQ=17)

                                                      cGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  1:346050/141‑1 (MQ=17)

                                                      cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc                                                                                        <  2:320738/141‑1 (MQ=32)

                                                       gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATTAACGCATATGCt                                                                                       >  1:180921/1‑141 (MQ=12)

                                                       gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTGGATAAACGCATATGcg                                                                                       >  1:47152/1‑140 (MQ=255)

                                                        gttttcAATGACAGTCGGGTTTACGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACGAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                      <  1:79680/136‑1 (MQ=11)

                                                         tGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTGGCCATTGCTACCTTCGTTGAAGGTGCCCGATCTGTAATTCCCTTTTCAGGTTTACCATATTGTTTACTGAGGCGTGTGGTAAACGAATAGGCTgc                                                                                     >  2:30958/1‑140 (MQ=11)

                                                            aCAATGGCAGTCAGGATTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGATGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg                                                                                  >  1:282832/1‑141 (MQ=255)

                                                             cAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGa                                                                                 >  2:68547/1‑141 (MQ=9)

                                                              aaTGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                <  1:225685/141‑1 (MQ=14)

                                                              aaTGGCAGGCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt                                                                                <  2:17885/141‑1 (MQ=9)

                                                               aTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAtt                                                                               <  2:355797/141‑1 (MQ=255)

                                                                           ttgaagtttaaAATCTTTAATTATTTTAGCCATTGCTGCCGTCGCTGAATGTGCCTGAACTGTAACTACCTTTTTAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCGGAATGATTATCTCGTTGCtt                                                                   <  2:294118/138‑1 (MQ=1)

                                                                                              aatattttACCCATTGCTGCCTTCGTGGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAat                                                <  1:232251/136‑1 (MQ=1)

                                                                                               ttgcTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGTTCTGTGATTACCGGTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatg                                               <  2:268577/138‑1 (MQ=1)

                                                                                                tattttAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgt                                              <  1:255940/138‑1 (MQ=0)

                                                                                                     taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  2:284427/138‑1 (MQ=17)

                                                                                                     taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTc                                         <  1:288775/138‑1 (MQ=17)

                                                                                                      aaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                        >  2:235418/3‑141 (MQ=17)

                                                                                                      aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTATAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTTATGTATGTcc                                        >  1:378585/1‑141 (MQ=11)

                                                                                                       accATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTACTGTATGTctc                                       >  2:117078/2‑139 (MQ=17)

                                                                                                         cATTGCTAGCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTTAGGTTTACCAAATTGTTTTATGAGTCGTTTGATAAACGTATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  2:50678/141‑1 (MQ=11)

                                                                                                         cATTGCTACCTGCGTTGAAGGTGCTTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                     <  2:140395/141‑1 (MQ=17)

                                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                    >  1:308666/1‑141 (MQ=32)

                                                                                                          aTTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctct                                    >  1:304769/1‑141 (MQ=32)

                                                                                                            tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCACTCTgt                                  >  2:227129/1‑140 (MQ=17)

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

                                                                                                               tgccTTGGTTGAACGTGCCTGATCTGTAATCAATTTTTTAGGCTTACCAAATTGGTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCGCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATc                               <  1:32921/139‑1 (MQ=11)

                                                                                                               tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATc                               <  2:188896/141‑1 (MQ=32)

                                                                                                                gccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                              <  2:154000/140‑1 (MQ=255)

                                                                                                                aCCTTCGTTGAAGGTGCCTGTTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCa                              <  2:86753/141‑1 (MQ=17)

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

                                                                                                                 ccTTCGTTGAAGGTGCCTGATCTGTACCTACCTTTTGTGGCTTGCCATATTGTGTACTGAGACGGTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                             <  1:270457/141‑1 (MQ=11)

                                                                                                                 ccTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCaa                             <  2:215240/141‑1 (MQ=32)

                                                                                                                   ttCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCCCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                           <  2:150437/141‑1 (MQ=17)

                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCa                        <  2:335405/141‑1 (MQ=32)

                                                                                                                      gTTGAAGGTGCCTGATCTGTAATTACCTATTGAGGTTTACCAAATTGTTTAATGAGGGGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCAAAGGTATGTGCCTCTACAtgagaggca                        >  2:160783/1‑133 (MQ=11)

                                                                                                                        tGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACg                      >  1:165714/1‑141 (MQ=17)

                                                                                                                           aGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAACTTGTTTACTAACGAGATTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAGCCAATTATCTAATGTATGTCCCTCTGTATCAATGGCACGtta                   >  1:318037/1‑141 (MQ=11)

                                                                                                                           aGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGata                   >  2:182769/1‑141 (MQ=32)

                                                                                                                            ggTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGatat                  >  1:117906/1‑141 (MQ=12)

                                                                                                                               gCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATaaa               >  1:109129/1‑141 (MQ=17)

                                                                                                                               gCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATcaa               >  1:216029/1‑141 (MQ=255)

                                                                                                                                ccTGTTCTGTTATTTCCTTTTGAGGTTTACCAAATTGTTTAATGCGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCTCTCTGCATCAATGGCACGATATAAAt              <  2:9119/141‑1 (MQ=11)

                                                                                                                                  tGATCTGTAATTACCTTTTGAGGAATACCAAATTGTTTAATGAGACGTTTGATAAACACATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAAAGTAAGACCCCCTGTATCAATGGCCCGATATaaaaaa            >  2:242616/1‑138 (MQ=11)

                                                                                                                                   gATTTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc           <  2:381568/141‑1 (MQ=17)

                                                                                                                                    aTCTGTGATTACCGTTTGAGTTTTACCATATTTTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt          <  1:131051/141‑1 (MQ=11)

                                                                                                                                      cTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATCTCTCACGTATGTCCCCCCCAAACACTGCCCCGCTATAAATCGCTcc        <  2:223562/141‑1 (MQ=11)

                                                                                                                                       tGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCa       <  1:61370/141‑1 (MQ=32)

                                                                                                                                          aaTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACACAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAttt    <  1:293323/141‑1 (MQ=17)

                                                                                                                                           aTTACCTTTTGAGGTTTACCAAATAGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCAtttt   >  1:278941/1‑141 (MQ=17)

                                                                                                                                            ttACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTc  <  1:187844/141‑1 (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: `‑`