BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,997	T→C	29.1%	T185T (ACA→ACG)	USA300HOU_0026 ←	IS431mec transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	35,997	0	T	C	29.1%	867.9 / inf	738	T185T (ACA→ACG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (256/266); new base C (116/98); total (373/365)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

Only 100 of 738 total aligned reads displayed.

TCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35861‑36136

tcttgtgAGTGTTCGTTTAATTTCCTAGCATGATGCTAATTTCGTGGCCTGGCGAAAATCCGTAGATCTGAAGAGACTTGCGGTTCTTTTTATACAGCGCGTGAAGACACTCAATACCTTTCAAAGTATTCTTTGCTGTAt                                                                                                                                          <  2:364006/138‑1 (MQ=1)

       aGTGTTCGTTTTACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt                                                                                                                                   <  2:148783/141‑1 (MQ=9)

       aGTGTTCGTTTAACTAGCTGGCATGATGCTAATTTCGTGGCACGGCGAAAATCCGTAGATCTGATGAGACCTGCCGTTCTTTTTATATAAAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACt                                                                                                                                   >  1:577575/1‑141 (MQ=11)

       aGTGGTCGTTTAGCTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATAGATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                                                                                                                   >  1:546951/1‑141 (MQ=11)

        gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                                                  >  1:271978/1‑141 (MQ=17)

         ggtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                                                 <  1:34528/140‑1 (MQ=12)

          gtTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTg                                                                                                                                <  1:635484/141‑1 (MQ=17)

           ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                                                                               >  1:144662/1‑141 (MQ=17)

             cGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                                                                             >  1:39537/1‑141 (MQ=32)

              gtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAc                                                                                                                            <  2:434122/139‑1 (MQ=17)

               tttAACTTGCGAGCATGACGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATatc                                                                                                                           <  2:555603/140‑3 (MQ=11)

               cTTAACTTGCTAGCATGTTGGTAAGTTCTGTGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAcc                                                                                                                           <  2:540537/141‑1 (MQ=11)

                ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACTGTTCTTTTTATATAGAGCGTAAAGACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCt                                                                                                                          <  2:340342/141‑1 (MQ=11)

                    cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTCTTGATACTTTGATACCTTGtc                                                                                                                      >  1:53299/1‑141 (MQ=17)

                    cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAGAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                                                                                      >  2:223914/1‑141 (MQ=17)

                    cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                                                                                      <  2:251446/141‑1 (MQ=17)

                      tGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                    >  2:330751/1‑141 (MQ=17)

                       gCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtcttt                                                                                                                   <  2:664539/141‑1 (MQ=17)

                        cTAGCATGATGCTAATTTCGTGGCATGGGGAAAAGCCGTAGAACTCAAGAGACGTACGGTTCTTTTGATATAGAGCGTAAATGCATGCAATACCTTTTAAAGCATTCTTTGCCGCAGTGATACTGTGATATCTTgcgttgt                                                                                                                  >  2:113780/1‑135 (MQ=11)

                        cTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttt                                                                                                                  >  1:724083/1‑140 (MQ=17)

                         tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCCATACCTTTTAAAGTATTCTTTGTTGTATTGATACTTTGATACCTTGtttttct                                                                                                                 >  1:91447/1‑141 (MQ=255)

                          aGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttctt                                                                                                                <  1:68374/141‑1 (MQ=255)

                          aGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTCAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttctt                                                                                                                >  1:597208/1‑141 (MQ=17)

                           gCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATGCTTTGATATCTTGTCTTTCTTa                                                                                                               >  1:168222/1‑141 (MQ=11)

                           gCATGATGCTAATTTCATGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGCTACTTTGATACCTTGtctttgttt                                                                                                               >  1:382374/1‑140 (MQ=255)

                              tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACtt                                                                                                            >  2:313299/1‑141 (MQ=255)

                              tgatgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACtt                                                                                                            <  1:427471/141‑1 (MQ=25)

                                     aaTTTCGTGGCATGGCGAAAATCCGTATATCTCAAGGTACCTGCGGTCCTTTTTATATTGAGCGTAAATATATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                                                                                     <  2:463112/141‑1 (MQ=11)

                                        ttCGTGGCATGGCGAAAATCCGTATATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTTATACCTTGTCTCTCTTACTTTAATATGCCg                                                                                                  >  1:60573/1‑141 (MQ=11)

                                          tgTGGCATGGCGAAAATCCGGAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTTTCTTTCTTACTTTAATATGACGGt                                                                                                <  1:652504/140‑1 (MQ=255)

                                          cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt                                                                                                <  2:651099/141‑1 (MQ=37)

                                            gggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGa                                                                                              <  1:100253/140‑1 (MQ=255)

                                             ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt                                                                                             <  1:776865/141‑1 (MQ=255)

                                             ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt                                                                                             >  1:258789/1‑141 (MQ=255)

                                             ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAt                                                                                             >  2:726129/1‑141 (MQ=38)

                                                     gAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCa                                                                                    >  2:275121/1‑141 (MQ=255)

                                                      aaaaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCaa                                                                                   <  1:130817/141‑1 (MQ=255)

                                                        aaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATg                                                                                 >  1:282232/1‑141 (MQ=255)

                                                           ccGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAgg                                                                              <  1:238922/141‑1 (MQ=255)

                                                           ccGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAgg                                                                              <  2:39630/141‑1 (MQ=38)

                                                           ccGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAgg                                                                              <  2:763292/141‑1 (MQ=255)

                                                           ccGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAgg                                                                              <  1:404685/141‑1 (MQ=38)

                                                           ccGTAGATCTGAAGAGACCTACGGATCTTTTTATATAGAGCGTACATACATTCCCTACCTTTTACAGTACTCTACGCCGTATTGATACTTAGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAgg                                                                              >  2:181554/1‑141 (MQ=11)

                                                            cGTAGATCTGAAGAGTCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCACTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATAACTTGTCTTTCTTACTTTAATATGACGTTGATCTTG‑TTAAATAAGGt                                                                             >  1:49288/1‑141 (MQ=11)

                                                            cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGGGGt                                                                             >  1:316125/1‑141 (MQ=38)

                                                            cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGt                                                                             >  1:669348/1‑141 (MQ=37)

                                                            cGTAGATCTGAAGAAAGCTGCGGTTCATAACATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGt                                                                             >  2:585556/1‑141 (MQ=11)

                                                             gTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATTTGACGGTGATCTTG‑CTCAATGAGGtt                                                                            >  1:416909/1‑141 (MQ=255)

                                                               tgATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGGCGTATTGATCCGGGGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTAt                                                                          <  1:180602/140‑1 (MQ=14)

                                                                         agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCCATGAGGTTATTCAGATAttt                                                                >  1:479788/1‑141 (MQ=255)

                                                                           aCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCg                                                              <  2:380022/141‑1 (MQ=38)

                                                                            cctacGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGa                                                             >  1:358838/5‑141 (MQ=255)

                                                                            ccTGCGGTTCTTTTTAGTTATAGCGTAAAGACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTACTATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGa                                                             <  1:336836/141‑1 (MQ=11)

                                                                                cGGTTCTTTTTATATAGAGCGTAAATACACTGAATACCTTTTCAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTG‑CTCAATGACGTTATTCCGATACATCGCTGTa                                                         >  2:624885/1‑141 (MQ=11)

                                                                                    tCTGTTTATATAGAGCGTAAATACATTCAATCCCTTTTAGAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGAGCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                     <  2:21021/141‑1 (MQ=255)

                                                                                        tttATATTGAGCTTAAATCCCCTCAATACCTTTTAAAGTATGCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATCTGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCa                                                 <  2:484769/141‑1 (MQ=11)

                                                                                         ttATATAGAGCGTAAAAACCAACGATGCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAg                                                >  2:355909/1‑141 (MQ=14)

                                                                                            tataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTCCTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCa                                             >  1:295709/1‑141 (MQ=255)

                                                                                             ataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAg                                            >  1:509129/1‑141 (MQ=255)

                                                                                              taGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAgg                                           <  1:765570/141‑1 (MQ=255)

                                                                                               agagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgt                                          >  1:218452/1‑141 (MQ=255)

                                                                                                gagCGTAGATACAGTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtt                                         >  2:669933/1‑141 (MQ=255)

                                                                                                gagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATGTCGATGTACAATGGCAGTCAGgtt                                         >  2:504491/1‑141 (MQ=255)

                                                                                                 agCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttt                                        >  1:64331/1‑141 (MQ=255)

                                                                                                     tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATCTGAC‑ACCACGTTGTCGCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagt                                    >  1:761018/1‑141 (MQ=11)

                                                                                                     tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAAACCATGCCGTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTGAGgtttaagt                                    <  1:562072/141‑1 (MQ=14)

                                                                                                      aaaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATCTCGATGTACAATGACAGTCAGgtttaagtt                                   >  2:490899/1‑141 (MQ=255)

                                                                                                      aaaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGCTCTTG‑ATAACTGAAGACCCTCGGATATTTCGATGTACAATGGCAGGCAGgtttaagtt                                   >  2:234963/1‑141 (MQ=14)

                                                                                                       aaTACACTCAATACCTATTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                  >  1:459616/1‑141 (MQ=255)

                                                                                                        aTACATTCAATATCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACTGTGACCTTG‑CGCATTGAGGTTATTCATATATTTCTATGTACAATGACAGTCAGgtttaagttta                                 <  2:744821/141‑1 (MQ=14)

                                                                                                        aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttta                                 >  2:681575/1‑141 (MQ=255)

                                                                                                         tacacacAATACCTTTTAAGGTATTCTTTGCTGGATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATTAGGTTATTCAGATATGTCGGTGTACGTTGGCAGTCAGgtttaggtttaa                                <  2:103926/135‑1 (MQ=14)

                                                                                                         tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                >  1:707180/1‑141 (MQ=255)

                                                                                                         tACACTCAATACCTTTTAAAGTAGTCTTTGCTGTATTGATACTTTGATACGTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCCATGTACAATGGCAGTCAGgtttaagtttaa                                >  2:569235/1‑141 (MQ=14)

                                                                                                          aCATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTaaa                               >  2:13421/1‑141 (MQ=255)

                                                                                                           cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCATTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                              >  2:424562/4‑141 (MQ=255)

                                                                                                           cATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTaaaa                              >  1:735132/1‑141 (MQ=255)

                                                                                                             ctCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaacc                            >  1:251407/2‑139 (MQ=255)

                                                                                                              ttaaTTCCTTTTTAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCt                           <  2:289711/139‑1 (MQ=255)

                                                                                                              tCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCt                           <  1:329756/141‑1 (MQ=255)

                                                                                                                  tACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTaa                       <  1:113780/141‑1 (MQ=255)

                                                                                                                   aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      <  2:735140/141‑1 (MQ=255)

                                                                                                                   aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      >  1:623361/1‑141 (MQ=255)

                                                                                                                   aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGATTTact                      >  2:773278/1‑139 (MQ=255)

                                                                                                                   aCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTGAAt                      >  2:520271/1‑141 (MQ=255)

                                                                                                                   aCCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAACTTTAAt                      >  2:196072/1‑141 (MQ=255)

                                                                                                                    ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     >  2:34520/1‑141 (MQ=255)

                                                                                                                    ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAACTTTAAAAGATTTAAtt                     >  2:118841/1‑141 (MQ=255)

                                                                                                                     cTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    <  1:369547/141‑1 (MQ=255)

                                                                                                                       tttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGAAGTAGGGTTAAAGGTTAAAAGCTTTAATTACt                  >  2:157153/1‑141 (MQ=14)

                                                                                                                         tAACGTATTATATGATGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGTACTTG‑CTAAATGAGGTTATTCCGATATTTCGATGTACAATGGCAGTCAGGTTTAAATTTAAAAGCTTTAAtttgttt                >  2:61084/1‑136 (MQ=14)

                                                                                                                         tAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTACGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACttt                >  1:441157/1‑141 (MQ=255)

                                                                                                                           aaGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTaa              >  2:601/1‑140 (MQ=255)

                                                                                                                            agcaTTCTATGCTGTATTAAGACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTaac             <  1:208000/138‑3 (MQ=14)

                                                                                                                            aGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTaac             <  1:188120/141‑3 (MQ=255)

                                                                                                                               aTTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCAt          >  1:617176/1‑141 (MQ=255)

                                                                                                                               aTTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCAt          >  1:151281/1‑141 (MQ=255)

                                                                                                                                   tttGCTGTATTGATACTTTGAGACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCCTTAATTACTTTAACCATTGCt      >  2:70107/1‑141 (MQ=255)

                                                                                                                                    ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTa     <  2:451101/141‑1 (MQ=255)

                                                                                                                                      gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   >  2:411865/1‑141 (MQ=255)

TCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTG‑CTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35861‑36136

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