BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	35,997	T→C	31.9%	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	31.9%	386.4 / 308.6	379	T185T (ACA→ACG)	USA300HOU_0026	IS431mec transposase
Reads supporting (aligned to +/- strand): ref base T (118/138); new base C (48/72); total (167/212)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 379 total aligned reads displayed.

TGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35864‑36136

    ggttttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATACGTAGATCTGAAGAGCTCTGCGGTTCTTTTTATATAGAGCGTAAATACAATCAATACCTTTTATTGTAGCCGTTGCTGTATTGATACt                                                                                                                                  <  2:318458/137‑1 (MQ=11)

     gtgtGCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACtt                                                                                                                                 <  2:301204/141‑1 (MQ=11)

      tgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                                                >  2:71151/1‑141 (MQ=17)

        ttcttttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGa                                                                                                                              <  1:209732/136‑1 (MQ=11)

        ttCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGa                                                                                                                              <  1:213419/141‑1 (MQ=11)

         tcgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGAt                                                                                                                             >  1:147007/5‑141 (MQ=11)

         tCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATAAATAGGTCTTATGAGACCTACGGTTCTTTTTATATAGAGCGTAAATTGATTCAAGACCTTTTTAATGATTCTTTTCCGTATTGATACTTTggt                                                                                                                             >  1:207518/1‑139 (MQ=11)

         tCGTTTAACTTGCTAGCATGATGCTAATTTAGTGGCATGGCGAAACCCCGAGGATCGGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGAt                                                                                                                             >  2:252975/1‑141 (MQ=11)

          cgtttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATa                                                                                                                            >  1:140120/4‑141 (MQ=17)

             ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                                                                                                         >  2:352824/1‑141 (MQ=17)

             ttAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCt                                                                                                                         <  1:306084/141‑1 (MQ=17)

              tAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCtt                                                                                                                        <  1:181970/141‑1 (MQ=11)

                aCTTGCTAGCATGATGGTAATTTCGTGGCATGGCGAAAGTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                                                                                                      <  2:20427/141‑1 (MQ=17)

                aCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTTTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                                                                                                      <  2:291571/141‑1 (MQ=17)

                  ttGCTAGCATGATGCTAATTTCGTGGCGTGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGTGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                                                                                    <  2:120120/141‑1 (MQ=12)

                  ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtct                                                                                                                    <  1:228445/141‑1 (MQ=17)

                  atGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTCGtct                                                                                                                    <  2:218539/140‑1 (MQ=17)

                   tGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctt                                                                                                                   <  2:257423/141‑1 (MQ=255)

                     cTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTCtctttc                                                                                                                 >  1:144153/1‑141 (MQ=11)

                      tAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGtctttct                                                                                                                <  2:225577/141‑1 (MQ=255)

                              tgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTaa                                                                                                        >  1:316558/1‑141 (MQ=21)

                              tgCCAATTGCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTCCGTTTGATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTaa                                                                                                        <  1:10270/141‑1 (MQ=11)

                                 tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTCATACTTTGCTATCTTGTCTTTCTTACTTTAatat                                                                                                     >  2:169074/1‑141 (MQ=11)

                                 tAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAatat                                                                                                     >  2:145104/1‑141 (MQ=25)

                                  aaTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                                                                                    >  1:295500/1‑141 (MQ=25)

                                  aaTTTCGTGGCATGGCGAAAATCCGGAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                                                                                    <  1:279712/141‑1 (MQ=25)

                                   aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGAGGATCTTGTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                                                                                                   >  2:326502/1‑141 (MQ=11)

                                     ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATacgacg                                                                                                 <  1:243375/141‑1 (MQ=255)

                                      tCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGTTACTTTGATATTTTGTCTTTCTTACTTTAATATGACgg                                                                                                >  1:59867/1‑141 (MQ=11)

                                      tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                                                                                <  1:52068/141‑1 (MQ=38)

                                      tCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACgg                                                                                                <  2:174076/141‑1 (MQ=37)

                                       cGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGt                                                                                               <  1:336728/141‑1 (MQ=255)

                                        gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGtt                                                                                              >  1:74761/1‑140 (MQ=37)

                                              tggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                        <  1:319359/137‑1 (MQ=14)

                                              tggggAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                        <  2:45155/137‑1 (MQ=14)

                                               ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGc                                                                                       >  2:286425/1‑141 (MQ=38)

                                                 cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTGTTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc                                                                                     >  1:302332/1‑141 (MQ=255)

                                                 cGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGATGGTGATCGTGCTc                                                                                     >  1:346409/1‑141 (MQ=255)

                                                     aaTCCTGAGCTCTGAAGCGCCCTGTGGTTCTTTTTATATAACGAGTTAAAACTTGCACAACCGTTTGCAGTGCTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                                                 <  1:155770/141‑1 (MQ=2)

                                                     aaTCCGTAGATCTGAAGAGAGCTGCGGTTCTTTTTATATAGAGCATAAATACGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTCTCTTCCTTTACTATGACGCTGATCTTGCTCAATg                                                                                 >  1:217718/1‑141 (MQ=11)

                                                     aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATg                                                                                 <  2:232589/141‑1 (MQ=255)

                                                      aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTACTTTTTCGGTGATCTTGCTCAATGa                                                                                >  1:332096/1‑141 (MQ=11)

                                                        ccGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAgg                                                                              <  1:343845/141‑1 (MQ=255)

                                                         cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGt                                                                             <  2:86950/141‑1 (MQ=255)

                                                          gTAAATCTGAAGAGCCCTACGGTTCTTTGTATATAGAGCGTATATACATTCAATACCTTCTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATAGGACGGTGATCTTGCTCAATGAGGtt                                                                            <  1:17279/141‑1 (MQ=14)

                                                             gATCTGAAGAGACCTGCGGTTCTTTTTATCTAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                                                                         >  1:6207/1‑141 (MQ=255)

                                                             gATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAtt                                                                         <  1:91582/141‑1 (MQ=255)

                                                              aTCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTc                                                                        <  1:187992/141‑1 (MQ=255)

                                                                cTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                                                      <  2:332823/141‑1 (MQ=255)

                                                                cTGAAGAGACCTACGGTTCTTTTTATATAGAGCGGAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCGTGCTTACTTTAGTATGACGGTGATCTTGCTCAATGAGGTTATTCaa                                                                      >  2:161741/1‑140 (MQ=11)

                                                                 tGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa                                                                     >  1:37486/1‑141 (MQ=38)

                                                                  gAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGat                                                                    <  1:106199/141‑1 (MQ=255)

                                                                   aaGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGTTAATCTTTCTCACTGAGGTTCTTCAAatt                                                                   >  1:151172/1‑140 (MQ=11)

                                                                           tacGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATg                                                           >  2:34518/3‑141 (MQ=255)

                                                                             cGGTTCTTTTTATAGAGGGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTCGTAGGATATGTCTATGtc                                                         <  2:166320/141‑2 (MQ=14)

                                                                               gTTTTTTTTATATAGAGCGAAAATACATTCAATACCTCTTAAAGTATTCTTTGCTGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACa                                                       <  1:268054/141‑1 (MQ=255)

                                                                                 tCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAAt                                                     >  1:312073/1‑141 (MQ=38)

                                                                                  ttttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCCTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATg                                                    <  2:133491/140‑1 (MQ=255)

                                                                                  cTTTTTATATAGGGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATg                                                    <  1:118939/141‑1 (MQ=255)

                                                                                  cTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATg                                                    <  2:304893/141‑1 (MQ=255)

                                                                                   tttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATgg                                                   <  1:20224/141‑2 (MQ=255)

                                                                                   tttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTGTAATATGGCGGTGCTCTTGCTCAATGAGGTTATTCAGTTATTTCGATGTACAATgg                                                   >  1:264602/1‑140 (MQ=14)

                                                                                    ttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATggc                                                  <  2:211641/141‑3 (MQ=255)

                                                                                     tttATATAGAGCGTAAATACATTCAATACCTTTTACAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACa                                                 <  2:102572/141‑1 (MQ=255)

                                                                                     tttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCGTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCa                                                 >  2:269471/1‑141 (MQ=255)

                                                                                       tatataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGGTATTTCGATGTCCAATGGCAgc                                               >  2:207599/1‑140 (MQ=255)

                                                                                        atataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTATTTTAAAATTGCTGTTATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTc                                              <  1:70696/141‑1 (MQ=11)

                                                                                         tataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCa                                             <  1:195761/141‑1 (MQ=255)

                                                                                            agagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTGCTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAGTGGCAGTCAGgt                                          >  2:179833/1‑141 (MQ=255)

                                                                                             tagtgTAAATTCACTCAATTCCTTTTAAAGTATTCTTTGCACGTTTCATACTTTAATGCTTTGTCTTTCTTACTCTAATATGACGGTGATCTTGCTCAATGAGGTTAGTCAGATATTTCGATGTACAATGGCCGTCAGgtt                                         <  2:111247/137‑1 (MQ=14)

                                                                                             gagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtt                                         >  2:156725/1‑141 (MQ=255)

                                                                                               gCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttta                                       >  2:147406/1‑141 (MQ=255)

                                                                                                   aaaTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGCTACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt                                   <  1:132492/141‑1 (MQ=255)

                                                                                                   aaaTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATCCCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTTAGATATTTTGATGTACAATGACAGTCAGgtttaagtt                                   >  2:312158/1‑141 (MQ=255)

                                                                                                   aaaTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtt                                   <  1:314456/141‑1 (MQ=255)

                                                                                                    aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACTGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                  <  1:288736/141‑1 (MQ=255)

                                                                                                    aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                  <  1:285358/141‑1 (MQ=255)

                                                                                                     aTACATTCATTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttta                                 <  2:138321/141‑1 (MQ=255)

                                                                                                     aTACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttta                                 >  1:87768/1‑141 (MQ=255)

                                                                                                      tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTAAAATGACAGTCAGgtttaagtttaa                                >  2:55538/1‑141 (MQ=255)

                                                                                                      tACATTCAATACCTTTTAAAGTATTCTATGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATCAGGTTATTCAGATATTTCGATGTACGATGACCGTCAGgtttaatgttaa                                >  2:29711/1‑141 (MQ=14)

                                                                                                      tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                >  1:254589/1‑141 (MQ=255)

                                                                                                        cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAGTGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                              <  1:240341/138‑1 (MQ=255)

                                                                                                             aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCttt                         >  1:298542/1‑141 (MQ=255)

                                                                                                              aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTCATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGCTGTACAATGCCAGTCTGGTTTAAGTTTAAAAGCTTTa                        >  1:126063/1‑141 (MQ=255)

                                                                                                              aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTa                        >  1:356407/1‑141 (MQ=255)

                                                                                                              aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTa                        <  2:7828/141‑1 (MQ=255)

                                                                                                              aTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCATTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTa                        <  2:199498/141‑1 (MQ=255)

                                                                                                                aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      <  2:66459/141‑1 (MQ=40)

                                                                                                                 ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     <  1:321733/141‑1 (MQ=255)

                                                                                                                 ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     <  2:194648/141‑1 (MQ=255)

                                                                                                                  cTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGATTATTCAGATATTTCGATGTACCATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    >  2:246634/1‑141 (MQ=255)

                                                                                                                  cTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    >  2:125065/1‑141 (MQ=255)

                                                                                                                   ttttAAAGTATTCTGTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAc                   >  1:191322/1‑141 (MQ=255)

                                                                                                                        aaGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTaa              <  2:37486/141‑2 (MQ=255)

                                                                                                                            aTTCTTTGCTGTATTGATACTTTGATACCTTTTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAt          <  1:45739/141‑1 (MQ=255)

                                                                                                                               cTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGc       <  2:61759/141‑1 (MQ=255)

                                                                                                                               cTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGc       <  1:138505/141‑1 (MQ=255)

                                                                                                                                tttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      >  1:211648/1‑141 (MQ=255)

                                                                                                                                   gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   <  1:29142/141‑1 (MQ=255)

TGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35864‑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: `‑`