BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	36,174	0	T	C	33.3%	760.4 / inf	984	K126K (AAA→AAG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (323/323); new base C (116/208); total (445/539)
Fisher's exact test for biased strand distribution p-value = 2.95e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.19e-01

Only 100 of 1005 total aligned reads displayed.

AATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36033‑36314

 atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTACTTACCTTTTGAGGt                                                                                                                                              >  2:438598/1‑141 (MQ=255)

 atatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAACTTTAATTACTTTAACCATTGCTACCTTCGTTGAGGGTGCCAGATCTGTAATTACATTTTGAGGt                                                                                                                                              >  2:321736/1‑141 (MQ=11)

  tatGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGtt                                                                                                                                             >  2:625794/1‑141 (MQ=21)

   atGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGttt                                                                                                                                            <  2:403809/141‑1 (MQ=21)

      tcGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  1:377645/140‑1 (MQ=18)

      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCCGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAAGCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         >  2:695391/1‑141 (MQ=255)

      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  1:106163/141‑1 (MQ=19)

      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCCTTTGAGGTTTAcc                                                                                                                                         >  2:53840/1‑141 (MQ=255)

      aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                                                                                                         <  2:167480/141‑1 (MQ=34)

        ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTATGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTCATTACCTTTTGAGGTTTACCaa                                                                                                                                       >  1:708609/1‑141 (MQ=255)

        ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaa                                                                                                                                       <  1:541971/141‑1 (MQ=11)

         gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGACCTGTAATCACCTTTTGAGGTTTACCaaa                                                                                                                                      <  1:505449/141‑1 (MQ=255)

         gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTGCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      >  1:682082/1‑141 (MQ=11)

         gTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGGACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCaaa                                                                                                                                      <  1:582302/141‑1 (MQ=11)

           gATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAAtt                                                                                                                                    >  1:406723/1‑141 (MQ=17)

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

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

             tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                                                                                                                  >  2:740833/1‑141 (MQ=12)

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

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

              gttGGTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGCTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAGCTTCGTTGAAGGTGCCTGATCTCTAATTACCTTTTGAGGTTTACCAAAGTGtt                                                                                                                                 <  2:574275/140‑1 (MQ=11)

              cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  1:744842/1‑141 (MQ=17)

              cTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGtt                                                                                                                                 >  2:481569/1‑141 (MQ=32)

              cTTGCACCATGAGGCTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTGTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTCAGGATAACCAAATtttt                                                                                                                                 <  2:111009/141‑1 (MQ=11)

                 gCTCAATGTGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                                                                                                              <  1:477998/141‑1 (MQ=255)

                  cTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAAt                                                                                                                             <  2:883022/141‑1 (MQ=17)

                    cAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                                                                                                           >  1:813856/1‑141 (MQ=17)

                       tGAGGGTATTCAGATATGTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGAACAATGACTTGCTCGTCCCTTTTGAGGCTTACCAAATTGTTTAATGAGAc                                                                                                                        >  2:195502/1‑141 (MQ=11)

                          ggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                                                                                                                     <  1:733975/141‑1 (MQ=32)

                            ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  1:943341/1‑141 (MQ=12)

                            ttATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   <  1:793025/141‑1 (MQ=12)

                            ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGTAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAACTTGTTTAATGAGACGTTTg                                                                                                                   <  1:567065/141‑1 (MQ=11)

                            ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  2:639514/1‑141 (MQ=255)

                            ttATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTg                                                                                                                   >  1:30909/1‑141 (MQ=255)

                             tATTCAGATATTTCGATGTACAATGGCAGTCATGTTTAAGTTTAAAAGCTTTAATTCCATTATCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGa                                                                                                                  <  1:448796/141‑1 (MQ=11)

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

                              aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGaa                                                                                                                 >  1:832948/1‑140 (MQ=12)

                              aTTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGAt                                                                                                                 >  1:458669/1‑141 (MQ=17)

                               ttCAGATATTTCGATGGACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATa                                                                                                                >  1:468857/1‑141 (MQ=12)

                                 cAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTATAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATaaa                                                                                                              <  1:668093/141‑1 (MQ=12)

                                   gATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:936112/141‑1 (MQ=255)

                                   aatatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACg                                                                                                            <  2:280023/140‑1 (MQ=17)

                                     tatTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                                                                                                          <  2:758916/141‑1 (MQ=12)

                                         tCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAAGTTTGATAAACGCATATg                                                                                                      >  2:800686/1‑141 (MQ=12)

                                         tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      >  2:361359/1‑141 (MQ=255)

                                         tCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATg                                                                                                      <  1:359604/141‑1 (MQ=255)

                                           gATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGTGACGTTTGATAAAAGCATATGCt                                                                                                    >  2:507597/1‑141 (MQ=255)

                                            gtGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg                                                                                                   <  2:296153/140‑1 (MQ=11)

                                                aCAATGTCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGATGAAtc                                                                                               <  2:244951/141‑2 (MQ=12)

                                                    tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa                                                                                           <  2:411506/141‑1 (MQ=255)

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

                                                                                    tactgtaaacATTGCTACCTTCGTTGAAGGTCCCTGATCTGTAACTACCTTTTGAGGTTTACCCAATTGTTTAATGCGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgt                                                           <  1:393291/132‑1 (MQ=1)

                                                                                       tttAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCCTTGCTTATGCAACCAAATATCTAatgtatg                                                        >  1:875549/1‑141 (MQ=2)

                                                                                        ttAACCATTGCTACCTTCGTTGAATGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAatgtatgt                                                       <  1:215119/141‑1 (MQ=14)

                                                                                         taaccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAACATCTAATGTATGTc                                                      <  1:493428/138‑1 (MQ=17)

                                                                                          aGCCATTGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTcc                                                     <  2:869232/141‑1 (MQ=12)

                                                                                          aGCCATTGCTGCCTTCGTTAGAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACACAACCAAATATCTAATGTATGTcc                                                     <  1:600493/141‑1 (MQ=11)

                                                                                            ccATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGATTGATTATTTCGTTGCTTATGCAACGAAGTATCTAATGTATGTctct                                                   >  1:267343/1‑141 (MQ=11)

                                                                                             cattgtgatcTTGGTTGAAGGTGCCTGATCTTTAATTACCTTTTTAGGTTTACCAAATTGTTTAATGAGACGTTTGATAACCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCctc                                                  <  1:34912/132‑1 (MQ=11)

                                                                                             cATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTCTGTccggc                                                  <  2:816210/141‑4 (MQ=17)

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

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

                                                                                               ttGCTGCCTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTg                                                >  2:228065/1‑141 (MQ=11)

                                                                                                tGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTgt                                               >  1:864622/1‑140 (MQ=17)

                                                                                                 tcTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCa                                              <  2:389014/140‑1 (MQ=17)

                                                                                                  ctgccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGTGGCTGCAGGTCGTGGCTTCCGCAACCAAATATCTAATGTATGTCCCTCTgtat                                             >  1:923128/4‑138 (MQ=11)

                                                                                                   tACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCAGATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAAGGTATGTCCCTCTGTATc                                            >  1:608365/1‑141 (MQ=12)

                                                                                                    aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAACTATCTAATGTATGTCCCTCTGTATCa                                           <  2:523849/141‑1 (MQ=17)

                                                                                                    aCCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCa                                           >  1:624478/1‑141 (MQ=17)

                                                                                                     ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:178206/141‑1 (MQ=255)

                                                                                                     ccTTCGTTGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTATTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          >  2:668002/1‑141 (MQ=255)

                                                                                                     ccTTCGTTGAAGGTGCCTGATCTGGAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCaa                                          <  2:208183/141‑1 (MQ=255)

                                                                                                      ctgcggtGAAGGTGCCTTATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         <  2:460441/135‑1 (MQ=255)

                                                                                                      cTTCGGTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTATGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAAt                                         <  2:683227/141‑1 (MQ=255)

                                                                                                       ttCGTGGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATg                                        <  1:281026/141‑1 (MQ=17)

                                                                                                          gTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGcc                                     >  2:838263/1‑140 (MQ=17)

                                                                                                           ttGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    >  1:244375/1‑141 (MQ=32)

                                                                                                           tcgAAGTTGCATGATTTGTAATTACCTTTTGACGTTTACTACATTGTTTAATGCGACGTTTGATACACGCATATGCGTAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCAc                                    <  1:321736/139‑1 (MQ=11)

                                                                                                            tGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATTAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAAGATCTAATGTATGTCCCTCTGCATCCATGGCACg                                   >  1:129055/1‑141 (MQ=12)

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

                                                                                                            tGAAGGTGCCTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCGCTCTGTATCAATGGCACg                                   >  1:189830/1‑141 (MQ=255)

                                                                                                                   gCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATaaa                            >  1:413758/1‑141 (MQ=32)

                                                                                                                    ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATCTGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                           >  1:436979/1‑141 (MQ=255)

                                                                                                                    ccTGATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAAt                           >  1:853246/1‑141 (MQ=12)

                                                                                                                       gATCTGTAACTACCTTTTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGc                        <  2:30909/141‑1 (MQ=12)

                                                                                                                        acttGTAACTACCTTTTGAGGCTTACCAAATTGTTGAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCt                       <  1:230531/138‑1 (MQ=11)

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

                                                                                                                          ttGTAATTACCTTTTTATGTTTACCAAATTGTTTAATGAGACGTTTGATAACCGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTcc                     <  1:780956/140‑1 (MQ=12)

                                                                                                                              aaTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCAttt                 <  2:787996/141‑1 (MQ=32)

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

                                                                                                                               aTTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGATTACGCAACCAAATATCTAATGTATGTCCCTCTGGATCAATGGGACGATATAAATAGCGCCActtt                >  2:913769/1‑141 (MQ=11)

                                                                                                                                  tccTTGTGAGGCTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  1:150237/140‑1 (MQ=255)

                                                                                                                                  aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  2:143432/141‑1 (MQ=17)

                                                                                                                                  aCCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAACGGCACCACATAAATAGcccaatctccct             >  2:56186/1‑130 (MQ=11)

                                                                                                                                  aCCTTTTGAGGCTGACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACAAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCt             <  2:848576/141‑1 (MQ=255)

                                                                                                                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCCTATGCTGAATGATTATCTCGTTGCTTACGCAACCACATATCTAATGTATGTCCCTCTGTATCCATGCCCCGATATAAAACGCCCCCTTGTCCttt           >  2:128111/1‑141 (MQ=11)

                                                                                                                                    cTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttt           >  1:893515/1‑141 (MQ=32)

                                                                                                                                      tttGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCttttt         >  2:350328/1‑140 (MQ=17)

                                                                                                                                         gAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAttt      >  1:198646/1‑141 (MQ=17)

                                                                                                                                          aGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGTATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTAtttt     >  1:75830/1‑141 (MQ=17)

AATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTATCTCGTTGCTTACGCAACCAAATATCTAATGTATGTCCCTCTGCATCAATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGAT  >  CP000730/36033‑36314

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 22 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`