BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,571,551	C→A	12.3%	D439Y (GAC→TAC)	hsdM ←	DNA methyltransferase M

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,571,551	0	C	A	12.3%	1135.1 / 90.8	334	D439Y (GAC→TAC)	hsdM	DNA methyltransferase M
Reads supporting (aligned to +/- strand): ref base C (141/149); new base A (21/22); total (162/172)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 335 total aligned reads displayed.

TCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCT  >  CP009273/4571403‑4571700

tCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCg                                                                                                                                                      <  1:100516/150‑1 (MQ=255)

 cTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCgg                                                                                                                                                     >  1:601927/1‑150 (MQ=255)

  tGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGgt                                                                                                                                                    <  2:694667/150‑1 (MQ=255)

       cTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGtttt                                                                                                                                               <  1:101828/150‑1 (MQ=255)

        tAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGttttt                                                                                                                                              <  1:1126021/150‑1 (MQ=255)

        tAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGttttt                                                                                                                                              >  1:97322/1‑150 (MQ=255)

        tAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGttttt                                                                                                                                              >  1:3891605/1‑150 (MQ=255)

          aTACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGt                                                                                                                                            >  1:1585200/1‑150 (MQ=255)

           tACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGtt                                                                                                                                           >  1:4079772/1‑150 (MQ=255)

           tACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGtt                                                                                                                                           >  2:950467/1‑150 (MQ=255)

               tccggctccggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTctct                                                                                                                                       >  1:3326491/1‑150 (MQ=255)

                 cggctccggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTc                                                                                                                                     >  1:4331352/1‑150 (MQ=255)

                 cggctccggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTc                                                                                                                                     <  2:2812598/150‑1 (MQ=255)

                  ggctccggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCg                                                                                                                                    >  2:2962761/1‑150 (MQ=255)

                      ccggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGt                                                                                                                                <  1:3818589/150‑1 (MQ=255)

                       cggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTc                                                                                                                               <  2:1178900/150‑1 (MQ=255)

                        ggcAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCg                                                                                                                              <  1:4022603/150‑1 (MQ=255)

                           aGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCa                                                                                                                           <  1:2198849/150‑1 (MQ=255)

                            ggCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCaa                                                                                                                          <  2:4350387/150‑1 (MQ=255)

                                  cGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCg                                                                                                                    <  2:3108130/150‑1 (MQ=255)

                                   ggCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGt                                                                                                                   >  2:2761598/1‑150 (MQ=255)

                                   ggCATAAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGt                                                                                                                   <  1:3245742/150‑1 (MQ=255)

                                    gCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGtc                                                                                                                  <  1:865729/150‑1 (MQ=255)

                                       tCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTctct                                                                                                               <  2:2742943/150‑1 (MQ=255)

                                       tCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTctct                                                                                                               <  1:1960613/150‑1 (MQ=255)

                                        cAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCtt                                                                                                              >  2:3443313/1‑150 (MQ=255)

                                        cAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCtt                                                                                                              >  1:3503562/1‑150 (MQ=255)

                                         aaTACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTc                                                                                                             >  1:1423971/1‑150 (MQ=255)

                                          aTACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCg                                                                                                            <  2:4033908/150‑1 (MQ=255)

                                          aTACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCg                                                                                                            <  1:1444579/150‑1 (MQ=255)

                                          aTACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTAGCTGTCGGCAACTTCCGTCTCTTCg                                                                                                            >  1:779473/1‑150 (MQ=255)

                                           tACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCgg                                                                                                           <  1:2898381/150‑1 (MQ=255)

                                            aCTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGc                                                                                                          >  2:4317338/1‑150 (MQ=255)

                                             cTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCg                                                                                                         >  2:569801/1‑150 (MQ=255)

                                                tCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTa                                                                                                      <  2:1875880/150‑1 (MQ=255)

                                                  tttATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTaaa                                                                                                    >  2:4164758/1‑150 (MQ=255)

                                                   ttATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTaaaa                                                                                                   <  2:328274/150‑1 (MQ=255)

                                                   ttATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTaaaa                                                                                                   >  1:2322418/1‑150 (MQ=255)

                                                    tATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAAc                                                                                                  <  2:113856/150‑1 (MQ=255)

                                                    tATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCATCGCTGTCGGCAACATCCGTCTCTTAGGCGTTaaaaa                                                                                                  >  1:2800109/1‑149 (MQ=255)

                                                     aTCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACt                                                                                                 >  1:1307583/1‑150 (MQ=255)

                                                       cTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTcc                                                                                               >  2:4437094/1‑150 (MQ=255)

                                                        tttCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCa                                                                                              >  2:4389836/1‑150 (MQ=255)

                                                          tCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAtt                                                                                            >  1:1665126/1‑150 (MQ=255)

                                                           cAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTc                                                                                           >  2:2493509/1‑150 (MQ=255)

                                                            aGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCa                                                                                          >  1:4628619/1‑150 (MQ=255)

                                                             gCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCAc                                                                                         >  1:2950638/1‑150 (MQ=255)

                                                             gCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCAc                                                                                         <  2:1644070/150‑1 (MQ=255)

                                                               cAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCt                                                                                       <  1:2078171/150‑1 (MQ=255)

                                                                aGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCtt                                                                                      >  2:4422976/1‑150 (MQ=255)

                                                                aGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCtt                                                                                      >  2:366609/1‑150 (MQ=255)

                                                                aGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCtt                                                                                      >  2:2275231/1‑150 (MQ=255)

                                                                  gagaTATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCa                                                                                    >  2:1971296/1‑150 (MQ=255)

                                                                  gagaTATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCa                                                                                    <  1:48958/150‑1 (MQ=255)

                                                                    gaTATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGt                                                                                  <  2:1771102/150‑1 (MQ=255)

                                                                    gaTATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGt                                                                                  >  2:1384934/1‑150 (MQ=255)

                                                                      tatCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTgc                                                                                <  2:1833188/150‑1 (MQ=255)

                                                                          cAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCgg                                                                            >  1:3564320/1‑150 (MQ=255)

                                                                          cAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCgg                                                                            >  1:4482258/1‑150 (MQ=255)

                                                                           aGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCggg                                                                           <  1:2596129/150‑1 (MQ=255)

                                                                           aGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCGGTCGGCAACTTCCGTCTCTTAGGCGTTAAAACACCAATCACCATCAGTGCGCggg                                                                           >  1:1196076/1‑150 (MQ=255)

                                                                             cGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCt                                                                         >  1:1409290/1‑150 (MQ=255)

                                                                              gAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCtt                                                                        <  2:1657436/150‑1 (MQ=255)

                                                                              gAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCtt                                                                        <  2:3538801/150‑1 (MQ=255)

                                                                              gAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCtt                                                                        >  2:3774460/1‑150 (MQ=255)

                                                                              gAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCAGTCACCTTCAGTGCGCGGGCtt                                                                        <  2:2042582/150‑1 (MQ=255)

                                                                               aaTCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTa                                                                       >  2:4139208/1‑150 (MQ=255)

                                                                                aTCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTaa                                                                      >  2:766956/1‑150 (MQ=255)

                                                                                aTCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTaa                                                                      <  2:4711598/150‑1 (MQ=255)

                                                                                       tttGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGc                                                               <  2:3881069/150‑1 (MQ=255)

                                                                                           cggtgcggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGt                                                           >  2:2210768/1‑150 (MQ=255)

                                                                                            ggtgcggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTc                                                          <  1:3941391/150‑1 (MQ=255)

                                                                                               gcggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTc                                                       <  2:4733985/150‑1 (MQ=255)

                                                                                                cggATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCg                                                      >  1:547291/1‑150 (MQ=255)

                                                                                                  gATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTCGCAAGATGCTGGTCGGGGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTGCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGcc                                                    >  2:4568826/1‑150 (MQ=255)

                                                                                                   aTCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCa                                                   <  1:4253852/150‑1 (MQ=255)

                                                                                                      cACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATac                                                >  2:3881720/1‑150 (MQ=255)

                                                                                                      cACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATac                                                >  2:3679636/1‑150 (MQ=255)

                                                                                                        cTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATacac                                              <  2:820722/150‑1 (MQ=255)

                                                                                                              gCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTc                                        <  1:3633449/150‑1 (MQ=255)

                                                                                                                tGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCaa                                      >  1:1051601/1‑150 (MQ=255)

                                                                                                                  aaCTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCacac                                    >  2:4659222/1‑150 (MQ=255)

                                                                                                                  aaCTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAAc                                    <  2:2322418/150‑1 (MQ=255)

                                                                                                                   aCTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACg                                   >  2:3040200/1‑150 (MQ=255)

                                                                                                                     ttGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGc                                 <  2:2683215/150‑1 (MQ=255)

                                                                                                                      tGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCt                                <  1:2626173/150‑1 (MQ=255)

                                                                                                                       gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTg                               <  1:3129492/150‑1 (MQ=255)

                                                                                                                       gcgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTg                               <  1:911526/150‑1 (MQ=255)

                                                                                                                        cgcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGc                              <  2:3360538/150‑1 (MQ=255)

                                                                                                                         gcCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCa                             <  2:630366/150‑1 (MQ=255)

                                                                                                                            aGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAAt                          <  2:235531/150‑1 (MQ=255)

                                                                                                                             gCGGCTGGTAGCAAGATGCTGGTGGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATg                         <  1:65522/150‑1 (MQ=255)

                                                                                                                             gCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATg                         <  1:2509223/150‑1 (MQ=255)

                                                                                                                                     tAGCAAGTTGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCgg                 <  1:3585703/150‑1 (MQ=255)

                                                                                                                                      aGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGt                <  1:4170080/150‑1 (MQ=255)

                                                                                                                                       gCAAGATGCTGGTCTGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTATACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTa               <  2:916206/150‑1 (MQ=255)

                                                                                                                                         aaGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTaaa             >  2:1238414/1‑150 (MQ=255)

                                                                                                                                           gATGCTGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACg           <  2:2886498/150‑1 (MQ=255)

                                                                                                                                             tGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGgt         <  1:3040774/150‑1 (MQ=255)

                                                                                                                                                tGGTAGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGtgtg      >  2:778805/1‑150 (MQ=255)

TCTGCCGCTAATACATCCGGCTCCGGCAGGCTGTCGGCATCAATACTGTCTTTATCTTTCAGCCAGGAGATATCCAGCGAATCGGATTTTGCGGTGCGGATCCACTCACGGCTGAACTTGCGCCAGCGGCTGGTAGCAAGATGCTGGTCGGTGTTTTTGTTCTCTTCGCTGTCGGCAACTTCCGTCTCTTCGGCGTTAAAACTCCATTCACCTTCAGTGCGCGGGCTTAAACCGTGCGGGTCTTCGCCATACACGCGCTCAAACGGCTGCAAATGCTCGTCGGTAAACGGTGTGCGCT  >  CP009273/4571403‑4571700

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