BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	2,291,888	G→A	8.7%	K359K (AAG→AAA)	USA300HOU_2155 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	2,291,888	0	G	A	8.7%	291.2 / 13.4	106	K359K (AAG→AAA)	USA300HOU_2155	transposase
Reads supporting (aligned to +/- strand): ref base G (58/38); new base A (2/7); total (61/45)
Fisher's exact test for biased strand distribution p-value = 3.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 106 total aligned reads displayed.

TTCAAAAGACGATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTG  >  CP000730/2291762‑2292024

ttCAAAAGACGATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGt                                                                                                                            >  2:513/1‑141 (MQ=11)

    aaaGACGATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGtt                                                                                                                        >  1:1034724/1‑141 (MQ=32)

     aaGACGATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTa                                                                                                                       >  2:834350/1‑141 (MQ=32)

         cGATATACTACGACGCTCCTACGAACTTGTCCAAGGATTAAGAAAAGATCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTTTGTGAGTAAAGGTGTATGGAAAGTGGTTAAata                                                                                                                   >  1:698250/1‑141 (MQ=11)

         cGATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAata                                                                                                                   <  1:847286/141‑1 (MQ=32)

          gATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAatat                                                                                                                  >  1:42462/1‑141 (MQ=32)

           atatACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAAtatt                                                                                                                 <  1:284930/141‑1 (MQ=32)

            tataCTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAACGGTGTATGGAAAGTGGTTAAAtatta                                                                                                                >  2:858881/1‑141 (MQ=17)

            tataCTACGACACTCCTACGAACTTGTCAAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAAtatta                                                                                                                >  1:1013823/1‑141 (MQ=17)

             ataCTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAAtattat                                                                                                               >  2:1023405/1‑141 (MQ=32)

                 tacGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGaa                                                                                                           <  2:477220/141‑1 (MQ=32)

                 tacGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGaa                                                                                                           <  1:404470/141‑1 (MQ=32)

                 tacGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGaa                                                                                                           <  1:404469/141‑1 (MQ=32)

                   cGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGaaaa                                                                                                         <  1:132553/141‑1 (MQ=32)

                     acacTCCTTCGAGCTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAACTTTATTAATCGTTTAAATTCATTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACa                                                                                                       <  1:102646/141‑1 (MQ=255)

                     acacTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACa                                                                                                       >  2:554746/1‑141 (MQ=32)

                     acacTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACa                                                                                                       >  2:1056287/1‑141 (MQ=32)

                        cTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCa                                                                                                    >  1:26340/1‑141 (MQ=32)

                        cTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCa                                                                                                    >  1:419580/1‑141 (MQ=32)

                           cTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAg                                                                                                 <  2:637519/141‑1 (MQ=32)

                           cTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAg                                                                                                 <  2:215176/141‑1 (MQ=32)

                             aCGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGa                                                                                               >  2:13767/1‑141 (MQ=32)

                                  cTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTa                                                                                          >  1:473054/1‑141 (MQ=32)

                                  cTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTa                                                                                          >  1:400721/1‑141 (MQ=32)

                                    tGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAg                                                                                        <  2:551352/141‑1 (MQ=32)

                                          aGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATAc                                                                                  <  1:110270/141‑1 (MQ=32)

                                            gATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACaa                                                                                >  1:726068/1‑141 (MQ=32)

                                              ttaagAAAAGATCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAAGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAAtt                                                                              <  2:652005/137‑1 (MQ=255)

                                              ttACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAAtt                                                                              >  1:263305/1‑141 (MQ=32)

                                              ttACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATC‑AAGGATGTTAAGAAATACAAttt                                                                             >  2:495051/1‑141 (MQ=17)

                                                 cGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATttat                                                                           >  2:602292/1‑141 (MQ=32)

                                                  gAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATttatt                                                                          >  2:759556/1‑141 (MQ=32)

                                                  gAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATttatt                                                                          >  2:793069/1‑141 (MQ=32)

                                                  gAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATttatt                                                                          >  1:21701/1‑141 (MQ=32)

                                                  gAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATttatt                                                                          >  1:201443/1‑141 (MQ=32)

                                                      aGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTAccc                                                                      <  2:733407/141‑1 (MQ=32)

                                                      aGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTAccc                                                                      <  2:120024/141‑1 (MQ=32)

                                                         ccTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGc                                                                   >  1:419514/1‑141 (MQ=32)

                                                          cTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAAGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCa                                                                  <  1:932288/141‑1 (MQ=17)

                                                           tAAGGTTATGTAATGGGCCTAAATTTATTAATCGTTTAAATTCAGGTAGTAAAAAGGCTGTGAGTGAGGTTGGATGTGAAGTGGGTGAATAGTATAGAGAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCAt                                                                 >  1:270318/1‑141 (MQ=11)

                                                                   tGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTaataat                                                         >  1:435789/1‑141 (MQ=32)

                                                                   tGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTaataat                                                         >  2:290984/1‑141 (MQ=32)

                                                                   tGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTaataat                                                         >  1:441217/1‑141 (MQ=32)

                                                                      aaTTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAAGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGt                                                      >  1:810426/1‑141 (MQ=17)

                                                                        ttGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGc                                                    <  2:45605/141‑1 (MQ=32)

                                                                        ttGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGc                                                    >  1:204957/1‑141 (MQ=32)

                                                                         tGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCt                                                   <  1:551070/141‑1 (MQ=32)

                                                                          ggCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCta                                                  <  2:261730/141‑1 (MQ=32)

                                                                          ggCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCta                                                  >  2:306359/1‑141 (MQ=32)

                                                                          ggCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCta                                                  <  2:154560/141‑1 (MQ=32)

                                                                           gCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCtat                                                 <  2:111282/141‑1 (MQ=32)

                                                                           gCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCtat                                                 >  2:301327/1‑141 (MQ=32)

                                                                           gCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCtat                                                 >  1:756584/1‑141 (MQ=32)

                                                                            ccTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCtata                                                >  2:333973/1‑141 (MQ=32)

                                                                             cTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAg                                               >  2:825648/1‑141 (MQ=32)

                                                                               aaaTTTCTTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGCGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATACTGGTGCTATAGaa                                             <  2:1038461/141‑1 (MQ=12)

                                                                                aaTTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAg                                            <  2:609458/141‑1 (MQ=32)

                                                                                aaTTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAAGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAg                                            <  1:799968/141‑1 (MQ=17)

                                                                                 aTTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAgg                                           >  1:636447/1‑141 (MQ=32)

                                                                                  tttATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAAGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGa                                          <  2:919939/141‑1 (MQ=17)

                                                                                   ttattaATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAAGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGaa                                         <  1:819324/141‑1 (MQ=17)

                                                                                    tattaATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAAt                                        <  1:50245/141‑1 (MQ=32)

                                                                                    tattaATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAAt                                        >  2:862363/1‑141 (MQ=32)

                                                                                    tattaATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAAt                                        <  1:447395/141‑1 (MQ=32)

                                                                                    tattaATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAAt                                        <  1:632049/141‑1 (MQ=32)

                                                                                         aTCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATtaata                                   >  1:841014/1‑141 (MQ=32)

                                                                                           cGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATtaataat                                 >  1:118406/1‑141 (MQ=32)

                                                                                            gTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATtaataata                                <  2:288812/141‑1 (MQ=32)

                                                                                              ttAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAg                              >  2:330944/1‑141 (MQ=32)

                                                                                              ttAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAg                              <  2:522539/141‑1 (MQ=32)

                                                                                              ttAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAg                              >  1:27592/1‑141 (MQ=32)

                                                                                                aaaTTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCAATAGAAGGAATTAATAATAAGAt                            >  2:355216/1‑141 (MQ=17)

                                                                                                aaaTTCAGTTAGTAAAAAGTCTGAGAGGAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGAt                            >  1:526335/1‑141 (MQ=17)

                                                                                                 aaTTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATa                           >  2:741407/1‑141 (MQ=32)

                                                                                                 aaTTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATa                           >  2:16162/1‑141 (MQ=32)

                                                                                                  aTTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATaa                          >  1:216970/1‑141 (MQ=32)

                                                                                                  aTTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATaa                          >  2:478981/1‑141 (MQ=32)

                                                                                                   ttCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATaaa                         <  1:440610/141‑1 (MQ=32)

                                                                                                   ttCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATaaa                         <  1:332196/141‑1 (MQ=32)

                                                                                                     cAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACAGCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAAt                       >  2:739543/1‑141 (MQ=17)

                                                                                                      aGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAAtt                      >  2:19332/1‑141 (MQ=32)

                                                                                                      aGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAAtt                      >  2:1062537/1‑141 (MQ=32)

                                                                                                        ttAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTaa                    <  1:958922/141‑1 (MQ=32)

                                                                                                           gTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCa                 >  2:568106/1‑141 (MQ=32)

                                                                                                           gTAAAAAGTCTGTGAGTAAGGGGGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATAGTGCTATAGAAGGAATCAATAATAAGATAAAATTAATCa                 >  1:662479/1‑141 (MQ=12)

                                                                                                               aaaGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTg             <  1:548066/141‑1 (MQ=34)

                                                                                                               aaaGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTg             <  1:1016766/141‑1 (MQ=34)

                                                                                                                aaGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCGATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGa            >  1:802396/1‑141 (MQ=21)

                                                                                                                 aGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGaa           <  1:363390/141‑1 (MQ=35)

                                                                                                                 aGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGaa           <  1:290984/141‑1 (MQ=34)

                                                                                                                   tCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAAtt         <  1:42539/141‑1 (MQ=35)

                                                                                                                    cTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAAttt        >  1:617325/1‑141 (MQ=34)

                                                                                                                     tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc       >  1:345215/1‑141 (MQ=35)

                                                                                                                     tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc       <  2:137369/141‑1 (MQ=35)

                                                                                                                     tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc       >  1:596515/1‑141 (MQ=35)

                                                                                                                     tgtgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTc       >  2:138654/1‑141 (MQ=35)

                                                                                                                       tgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCtt     >  2:635745/1‑141 (MQ=35)

                                                                                                                       tgAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCtt     >  1:803529/1‑141 (MQ=35)

                                                                                                                        gAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCttt    >  1:397960/1‑141 (MQ=35)

                                                                                                                          gTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTg  <  2:4998/141‑1 (MQ=35)

TTCAAAAGACGATATACTACGACACTCCTACGAACTTGTCCAAGGATTACGAAAAGACCTAAGGTTATGTAATTGGCCTAAATTTATTAATCGTTTAAATTCAGTTAGTAAAAAGTCTGTGAGTAAGGGTGTATGGAAAGTGGTTAAATATTATAGAAAACATCAAAGGATGTTAAGAAATACAATTTATTACCCAGCATTTAATAATGGTGCTATAGAAGGAATTAATAATAAGATAAAATTAATCAAGTGAATTTCTTTTG  >  CP000730/2291762‑2292024

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