BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,101,556	T→G	5.8%	P43P (CCA→CCC)	csgG ←	curli secretion channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,101,556	0	T	G	5.8%	390.3 / 12.2	135	P43P (CCA→CCC)	csgG	curli secretion channel
Reads supporting (aligned to +/- strand): ref base T (66/61); new base G (2/6); total (68/67)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 135 total aligned reads displayed.

CCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGC  >  NC_000913/1101409‑1101704

ccAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTgg                                                                                                                                                    >  2:476065/1‑150 (MQ=255)

 cAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGc                                                                                                                                                   >  1:350915/1‑150 (MQ=255)

 cAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGc                                                                                                                                                   <  1:31032/150‑1 (MQ=255)

 cAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGGGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGGTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGGCGGCGCTGGc                                                                                                                                                   <  2:1930604/150‑1 (MQ=255)

   gCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAg                                                                                                                                                 <  2:247650/150‑1 (MQ=255)

    cGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCACGGTTTCGTCCTGAATGTTGAATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAGa                                                                                                                                                >  2:1042953/1‑150 (MQ=255)

        aaTCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATggg                                                                                                                                            <  1:1699901/150‑1 (MQ=255)

            tttCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCaa                                                                                                                                        >  1:2135937/1‑150 (MQ=255)

                aGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCt                                                                                                                                    >  1:1032641/1‑150 (MQ=255)

                 gTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCtt                                                                                                                                   <  1:471444/150‑1 (MQ=255)

                  tGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCttt                                                                                                                                  >  2:74378/1‑150 (MQ=255)

                   gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                                 >  2:1956722/1‑150 (MQ=255)

                   gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTg                                                                                                                                 <  1:1820025/150‑1 (MQ=37)

                   gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGGTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCATCTGCCAGCGCCGACGGGTAAAATCTTTg                                                                                                                                 <  1:1595177/150‑1 (MQ=37)

                    ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCGGTTTTCGTGCTGAATGTTGTATACCGAAACAAATATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGt                                                                                                                                <  2:2282554/150‑1 (MQ=255)

                    ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGt                                                                                                                                <  2:1283753/150‑1 (MQ=255)

                    ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGt                                                                                                                                <  2:2352695/150‑1 (MQ=255)

                     cGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTa                                                                                                                               <  1:1461141/150‑1 (MQ=255)

                      gTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGttt                                                                                                                              <  2:1862460/150‑3 (MQ=37)

                        gACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGct                                                                                                                            >  1:363244/1‑150 (MQ=255)

                        gACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGct                                                                                                                            <  2:623024/150‑1 (MQ=255)

                            aGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGa                                                                                                                        <  1:1378268/150‑1 (MQ=255)

                              cATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGc                                                                                                                      <  1:404907/150‑1 (MQ=255)

                               aTTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCa                                                                                                                     >  2:551433/1‑150 (MQ=255)

                               aTTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCa                                                                                                                     >  1:2453543/1‑150 (MQ=255)

                                   ccGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                 <  2:705699/150‑1 (MQ=255)

                                    cGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAgg                                                                                                                <  1:64232/150‑1 (MQ=255)

                                        gcgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCAtt                                                                                                            <  2:1046903/150‑1 (MQ=255)

                                         cgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTa                                                                                                           >  1:1207038/1‑150 (MQ=255)

                                             ttGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGt                                                                                                       >  1:1700780/1‑150 (MQ=255)

                                              tGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTc                                                                                                      <  1:1020999/150‑1 (MQ=255)

                                              tGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTc                                                                                                      <  1:1160998/150‑1 (MQ=255)

                                                cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCgg                                                                                                    >  1:1501056/1‑150 (MQ=255)

                                                  gAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTc                                                                                                  >  2:1020630/1‑150 (MQ=255)

                                                       gcagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggcg                                                                                             <  1:1115540/150‑1 (MQ=255)

                                                       gcagTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggcg                                                                                             >  2:1126451/1‑150 (MQ=255)

                                                       gcagTGGAGAAGTTACTTGCCGGGTAGGGTTAAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggcg                                                                                             >  2:1910539/1‑150 (MQ=255)

                                                         agTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTggcggc                                                                                           <  1:540263/150‑1 (MQ=255)

                                                          gTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGAt                                                                                                                                               >  2:656980/1‑97 (MQ=255)

                                                          gTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGAt                                                                                                                                               <  1:656980/97‑1 (MQ=255)

                                                          gTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGct                                                                                          <  2:2326396/150‑1 (MQ=255)

                                                            ggAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttc                                                                                        <  2:2201557/150‑1 (MQ=255)

                                                             gagaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttct                                                                                       >  2:787969/1‑150 (MQ=255)

                                                               gaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                     >  2:2200137/1‑150 (MQ=255)

                                                               gaAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                     >  1:1247621/1‑150 (MQ=255)

                                                               gaAGTTACTTGCCGGGTACGGTTTAAATTGCCCGGTTTCGACCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGAACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                     >  2:2200144/1‑150 (MQ=255)

                                                                 aGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAg                                                                                   >  1:2186035/1‑150 (MQ=255)

                                                                  gTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAgg                                                                                  >  1:427632/1‑150 (MQ=255)

                                                                   ttACTTGCCGGGTAGGGGTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggc                                                                                 <  1:2093202/150‑1 (MQ=255)

                                                                      cTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggg                                                                              <  2:1727452/150‑1 (MQ=255)

                                                                          ccGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggggcgg                                                                          >  2:611940/1‑150 (MQ=255)

                                                                           cGGGGAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGGTGGATACCGACACAAAGAGTTTACGCGTCGGCGCGGGGAGAGGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTACTGTCGGGCTGGCGGCTGCGTTAGGCGGGGGGGt                                                                         <  2:1087607/150‑1 (MQ=255)

                                                                            gggTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGtt                                                                        >  1:1386711/1‑150 (MQ=255)

                                                                              gTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      >  1:1400709/1‑150 (MQ=255)

                                                                               tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAg                                                                     <  2:1885104/150‑1 (MQ=255)

                                                                                  ggTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTTGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                  >  1:2035425/1‑150 (MQ=255)

                                                                                    tttAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATATTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATcc                                                                <  1:620575/150‑1 (MQ=255)

                                                                                      tAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGGTTGTTTAGGCGGGGCGGTTAAGCATCCGc                                                              >  1:594587/1‑150 (MQ=255)

                                                                                      tAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGc                                                              <  2:458806/150‑1 (MQ=255)

                                                                                       aaaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCt                                                             <  1:469663/150‑1 (MQ=255)

                                                                                       aaaTTGCCCGGGTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCt                                                             <  1:1678881/150‑1 (MQ=255)

                                                                                        aaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTc                                                            <  2:1173398/150‑1 (MQ=255)

                                                                                        aaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTc                                                            <  1:116346/150‑1 (MQ=255)

                                                                                         aTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCa                                                           >  2:1114478/1‑150 (MQ=255)

                                                                                          ttGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAg                                                          <  2:817593/150‑1 (MQ=255)

                                                                                            gCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTa                                                        <  2:1518547/150‑1 (MQ=255)

                                                                                                 gTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCTCGAGGCATTAATGTCGGGCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATg                                                   <  2:1573300/150‑1 (MQ=255)

                                                                                                   ttCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGAc                                                 <  1:1415470/150‑1 (MQ=255)

                                                                                                     cGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACgg                                               <  2:2454987/150‑1 (MQ=255)

                                                                                                     cGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACgg                                               <  2:1032641/150‑1 (MQ=255)

                                                                                                         cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                           >  1:1801221/1‑150 (MQ=255)

                                                                                                         cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                           <  2:2037935/150‑1 (MQ=255)

                                                                                                         cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                           <  1:894493/150‑1 (MQ=255)

                                                                                                         cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCGGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTGAGTAACGTGACgcccac                                           >  2:425410/1‑145 (MQ=255)

                                                                                                                ttgtATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGa                                    <  1:74378/150‑1 (MQ=255)

                                                                                                                  gtATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGaaa                                  <  1:528409/150‑1 (MQ=255)

                                                                                                                    ataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTa                                                                       <  1:747648/111‑1 (MQ=255)

                                                                                                                     taCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATaa                               >  1:1435392/1‑150 (MQ=255)

                                                                                                                     taCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATaa                               >  1:423008/1‑150 (MQ=255)

                                                                                                                     taCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATaa                               >  2:424841/1‑150 (MQ=255)

                                                                                                                      aCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAg                              >  2:577386/1‑150 (MQ=255)

                                                                                                                      aCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAg                              >  2:401857/1‑150 (MQ=255)

                                                                                                                       ccGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                             >  1:104961/1‑150 (MQ=255)

                                                                                                                       ccGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                             >  2:564996/1‑150 (MQ=255)

                                                                                                                        cGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgcg                            >  2:1728869/1‑150 (MQ=255)

                                                                                                                             cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                       >  2:2397535/1‑150 (MQ=255)

                                                                                                                             cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                       <  1:114223/150‑1 (MQ=255)

                                                                                                                             cAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                       >  1:2271331/1‑150 (MQ=255)

                                                                                                                              aaaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCAt                      >  1:1031534/1‑150 (MQ=255)

                                                                                                                              aaaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCAt                      >  2:2459467/1‑150 (MQ=255)

                                                                                                                                aGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGa                    >  2:2251475/1‑150 (MQ=255)

                                                                                                                                aGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGa                    >  1:760508/1‑150 (MQ=255)

                                                                                                                                  aTTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGatt                  <  1:1070031/150‑1 (MQ=255)

                                                                                                                                     ttACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGattatt               >  2:1703723/1‑150 (MQ=255)

                                                                                                                                      tACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAttt              <  1:1892243/150‑1 (MQ=255)

                                                                                                                                       aCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAtttt             >  2:1160595/1‑150 (MQ=255)

                                                                                                                                        cccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATACGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc            >  2:1112740/1‑150 (MQ=255)

                                                                                                                                          cGTCGGCGCTGGCAGTTGGGTCAAATCTTTGTAGCTCTGAGCACGCGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATGCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCt          <  2:519584/150‑1 (MQ=255)

                                                                                                                                            tCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTa        <  1:1222763/150‑1 (MQ=255)

                                                                                                                                                  cTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGc  <  2:679189/150‑1 (MQ=255)

CCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGC  >  NC_000913/1101409‑1101704

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