BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,559,170	C→A	6.3%	D258Y (GAT→TAT)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,559,170	0	C	A	6.3%	712.9 / 2.1	237	D258Y (GAT→TAT)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand): ref base C (74/134); new base A (26/2); total (101/136)
Fisher's exact test for biased strand distribution p-value = 3.91e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

GTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGG  >  CP009273/3559022‑3559319

gTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATca                                                                                                                                                      <  2:3114240/150‑2 (MQ=255)

gTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCt                                                                                                                                                      >  2:1594386/1‑150 (MQ=255)

gTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCt                                                                                                                                                      >  2:1594512/1‑150 (MQ=255)

 ttttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATAtt                                                                                                                                                     >  1:2769798/1‑150 (MQ=255)

  tttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTc                                                                                                                                                    <  2:1709333/150‑1 (MQ=255)

   ttATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTCGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTcc                                                                                                                                                   >  2:319775/1‑150 (MQ=255)

    tATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCa                                                                                                                                                  <  1:2657182/150‑1 (MQ=255)

    tATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCa                                                                                                                                                  <  1:1354339/150‑1 (MQ=255)

       gCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCAcat                                                                                                                                               >  2:2070288/1‑148 (MQ=255)

         aaCTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGcc                                                                                                                                             <  1:1916657/150‑1 (MQ=255)

             gATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACTGCCGCGa                                                                                                                                         >  2:1786862/1‑150 (MQ=255)

             gATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTATTGCACAGCCGCGa                                                                                                                                         >  1:4409509/1‑150 (MQ=255)

               tAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAag                                                                                                                                       <  1:1551101/150‑1 (MQ=255)

                aaTGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAagt                                                                                                                                      <  1:2617465/150‑1 (MQ=255)

                  tGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTATTAAAATGCCGCGGagtag                                                                                                                                    >  1:2807313/1‑150 (MQ=255)

                   ggcggcTTACAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAagtagt                                                                                                                                   <  1:3873227/150‑1 (MQ=255)

                    gcggcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTc                                                                                                                                  <  2:3093793/150‑1 (MQ=255)

                     cggcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCa                                                                                                                                 <  1:4284679/150‑1 (MQ=255)

                     cggcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTTGTCa                                                                                                                                 >  2:1303776/1‑150 (MQ=255)

                       gcTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCAcc                                                                                                                               <  2:2917309/150‑1 (MQ=255)

                         ttAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCt                                                                                                                             <  1:2104737/150‑1 (MQ=255)

                         ttAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCt                                                                                                                             <  2:3649433/150‑1 (MQ=255)

                         ttAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCt                                                                                                                             <  2:703991/150‑1 (MQ=255)

                         ttAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCAcccg                                                                                                                             >  1:4484293/1‑149 (MQ=255)

                            aaaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGt                                                                                                                          <  2:4034284/150‑1 (MQ=255)

                             aaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGtt                                                                                                                         >  2:2535723/1‑150 (MQ=255)

                             aaTGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGtt                                                                                                                         <  2:1917833/150‑1 (MQ=255)

                                 tCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACACTTCCGCGAAGTAGTCACCCGGGTTGAAt                                                                                                                     >  2:2667335/1‑150 (MQ=255)

                                  ccTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAtt                                                                                                                    <  1:1582577/150‑1 (MQ=255)

                                  ccTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAAtt                                                                                                                    <  2:186008/150‑1 (MQ=255)

                                       aTGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCg                                                                                                               >  2:1605769/1‑150 (MQ=255)

                                          gTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAgg                                                                                                            <  1:3819464/150‑1 (MQ=255)

                                           tCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGt                                                                                                           <  2:2312200/150‑1 (MQ=255)

                                           tCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGt                                                                                                           <  1:3256229/150‑1 (MQ=255)

                                              aGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTATGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGAGGGTaa                                                                                                        >  2:1263032/1‑150 (MQ=255)

                                                gAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAtt                                                                                                      <  1:3768745/150‑1 (MQ=255)

                                                gAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAtt                                                                                                      <  2:3871824/150‑1 (MQ=255)

                                                gAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATATTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCGAGGTTaaat                                                                                                      >  2:802548/1‑148 (MQ=255)

                                                 aaaCCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAAttt                                                                                                     >  2:3965239/1‑150 (MQ=255)

                                                    ccAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTGAATTTCGc                                                                                                  <  1:3965278/150‑1 (MQ=255)

                                                     cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCt                                                                                                 <  1:3988083/150‑1 (MQ=255)

                                                     cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGACGTAGTCACCCGGGTTGTATTTACCGAGGTTAATTTAGCt                                                                                                 >  1:3238678/1‑150 (MQ=255)

                                                     cAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATATTACCCAACCGCGAAGTATTCAACCGGGTTTAATTTAACGAGGTTAATTTTGCt                                                                                                 >  1:455433/1‑150 (MQ=255)

                                                       ggAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTAACGAGGTTAATTTCGCTAc                                                                                               >  2:1959919/1‑150 (MQ=255)

                                                        gAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACt                                                                                              <  2:561779/150‑1 (MQ=255)

                                                        gAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACt                                                                                              <  2:966573/150‑1 (MQ=255)

                                                         aaGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTg                                                                                             <  2:878975/150‑1 (MQ=255)

                                                           gTATTCCTGACGCAGGCGCATCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGATTGGTTTTTATATTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGATGTTAATTTCGCTACTGGc                                                                                           <  2:2130620/150‑1 (MQ=255)

                                                           gTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGc                                                                                           <  2:951468/150‑1 (MQ=255)

                                                           gTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGc                                                                                           <  2:4521127/150‑1 (MQ=255)

                                                             agtCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCtt                                                                                         <  2:3916862/148‑1 (MQ=255)

                                                             aTTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTTTATTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGAGGTTAATTTAGCTACTGGcgt                                                                                         >  1:122450/1‑148 (MQ=255)

                                                                ccTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTTGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGTTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTgcg                                                                                      <  1:4382865/150‑1 (MQ=255)

                                                                    aCGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACt                                                                                  <  1:401218/150‑1 (MQ=255)

                                                                        aGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTACACAGCCGCGAAGTAGTCACCCTGGTTGAATTTAACGAGGTTAATTTAGCTACTGGCTTGCGCACTCcac                                                                              >  1:297603/1‑150 (MQ=255)

                                                                         ggCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCcaca                                                                             <  1:353607/150‑1 (MQ=255)

                                                                            gcAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAAc                                                                          >  2:1724062/1‑150 (MQ=255)

                                                                             cAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACg                                                                         <  1:2701234/150‑1 (MQ=255)

                                                                                 tCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTAACGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGc                                                                     >  1:2263684/1‑150 (MQ=255)

                                                                                      cccGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGtt                                                                <  1:2981218/150‑1 (MQ=255)

                                                                                        cGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTgg                                                              <  1:1357841/150‑1 (MQ=255)

                                                                                        cGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTgg                                                              <  1:3851144/150‑1 (MQ=255)

                                                                                          gAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTc                                                            <  1:2151141/150‑1 (MQ=255)

                                                                                             tAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTcgcg                                                         <  2:46164/150‑1 (MQ=255)

                                                                                             tAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTcgcg                                                         <  2:121778/150‑1 (MQ=255)

                                                                                               ggTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTc                                                       <  2:2769798/150‑1 (MQ=255)

                                                                                                   gAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtg                                                   <  2:2573292/150‑1 (MQ=255)

                                                                                                   gAGTCATCCGGATACAGGACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGgtg                                                   <  2:1418921/150‑1 (MQ=255)

                                                                                                     gTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTc                                                 <  1:594583/150‑1 (MQ=255)

                                                                                                      tCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCg                                                <  1:2543228/150‑1 (MQ=255)

                                                                                                        aTCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTTTCTTCCACAGCCCCGAACTAGTCACCCAGGTTGAAATTAACGAGGTTAATTTCGCTAATTGCTTTCGCACTACACACACACCCCGTGTTGGTAGCGTCGGTGTCGTa                                              >  2:109858/1‑150 (MQ=255)

                                                                                                             gATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAg                                         <  1:687495/150‑1 (MQ=255)

                                                                                                             gATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTACACAGCCGCGGAGTATTAACCCTGGTTGGATTTACCGAGGGTAATTTCGCTACTGGCTTGCGCCCTCCACCAACCCAGCGTGTTGGTCGCGTCGGTGTCGTTAcacg                                         >  2:3068437/1‑147 (MQ=255)

                                                                                                              aTACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAgg                                        <  2:1572338/150‑1 (MQ=255)

                                                                                                               tACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAggg                                       <  1:601164/150‑1 (MQ=255)

                                                                                                                 cAGTACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGCAGTAGTAACCCGGGTTGAATTTAACGAGGTTAATTTCGCTACTGGCTTGCGCACTACACAAACGCAGCGGGTTTGTCGCGTCGGTGTCGTTACCAGGGAt                                     >  2:1201891/1‑150 (MQ=255)

                                                                                                                  aGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGAtt                                    <  1:1724062/150‑1 (MQ=255)

                                                                                                                   gTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTa                                   <  2:4409509/150‑1 (MQ=255)

                                                                                                                    tACGCGAGATACGTTCTCGGAGTGGTGTTTATCTTCCACAGCCGCGAAGTAGTCACCCGGGTTGAATTTAACGGGGTTAATTTCGCTACGGGCTTGCGCACTCCACAAACGCCGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTAt                                  >  1:730570/1‑150 (MQ=255)

                                                                                                                     aCGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATc                                 <  2:2490152/150‑1 (MQ=255)

                                                                                                                         gagaTACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCGCGAAGTAGTCACCCAGGTTGAATTTACCCAGGTTAATTTCGCTACTGGCGTGCGCACTCCACAAACGCAGCGTGTTGGGCGCGGCGGTGGCGTAACCAGGGATTAACTGAt                             >  1:1848468/1‑150 (MQ=255)

                                                                                                                         gagaTACGTTCTCGGAGTGGTGTTTATCTTACACAGCCGCGCAGTAGTCACCCGGGTTTAATTTACCGAGGTTAATTTCGCTAATGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCACGGATTATCTGAt                             >  1:4124295/1‑150 (MQ=255)

                                                                                                                          agaTACGTTCTCGGAGTGGTTTTTTTATTCCACTTCCGCGGAGTAGGAACCCAGGTTTAATTTAACGAGGTTAATTTCGCTACTGGCTGGCGCACACCACAAACCCAGCGTGTTGGTCGCGGCGGTGGCGTAACCACGGCTTATCTGAta                            >  2:3435259/1‑149 (MQ=255)

                                                                                                                          agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGAta                            >  2:1078163/1‑149 (MQ=255)

                                                                                                                          agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATc                            <  1:1776594/150‑1 (MQ=255)

                                                                                                                          agaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATc                            <  1:1776383/150‑1 (MQ=255)

                                                                                                                           gaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCa                                          >  2:3994896/1‑135 (MQ=255)

                                                                                                                           gaTACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCg                           <  2:893651/150‑1 (MQ=255)

                                                                                                                             tACGTTCTCGGAGTGGTTTTTATCTTCCACAGCCCCGAAGTTTTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTAAAGGCGTTCGCACACCACAAAAGCAGCGGGTTGGTCGCGTCGGTGTCGTTACCAGGGATTATCTGATCGTa                         >  2:430259/1‑150 (MQ=255)

                                                                                                                               cGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAg                       <  2:851001/150‑1 (MQ=255)

                                                                                                                               cGTTCTCGGAGTGGTTTTTATATTACACATCCGCGAAGTAGTCACCCTGGTTTAAATTAACGAGGTTAATTTAGCTACTTGCTTTCGCACACCACAAACCCACCCTGTTTGGGGCGGCGGTGTCGTAACCAGGGATTATATTATCGTAAg                       >  2:4180922/1‑150 (MQ=255)

                                                                                                                                gTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGc                      >  2:3007256/1‑150 (MQ=255)

                                                                                                                                          tGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTAGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGa            >  2:617000/1‑150 (MQ=255)

                                                                                                                                            gTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAt          >  1:2065747/1‑150 (MQ=255)

                                                                                                                                               tttATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAAtctc       >  1:784465/1‑150 (MQ=255)

                                                                                                                                                ttATCTTCCACTGCCGCGAAGTAGTCACCCTGTTTGAATTTACCGTGTTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGGGTGTTGGTCGCGTCGGTGTCTTAACCTGGTATTATCTGATCGTAAGCGACTCCCAGAATctct      <  2:1283492/150‑1 (MQ=255)

                                                                                                                                                ttATCTACCACTGCCGCGTTGTAGTCACACTGGTTGAATTTACCGAGTTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGAAGGTTTTTGGTCGCGTCGGTGTCTTAACCAGGTATTATCTGTTCGTAAGCGACTCCCAGAATctct      <  2:2355687/150‑1 (MQ=255)

                                                                                                                                                   tCTTCCACTGCCGCTATGTAGACACCCTGGTTGAATTTACCGAGGTTATTTTCGATACTGGCTTGCGCACTCCACAAACGCCGCGTGTTGGTCGCGTGGGTGTCTTACCCTGGGATTATCTTATCTTAAGCGACTCCCAGAATCTCTTCg   <  1:1473264/150‑1 (MQ=255)

                                                                                                                                                    cTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGCCTCCCAGAATCTCTTCgg  >  1:1005354/1‑150 (MQ=255)

                                                                                                                                                    cTTCCACTGCCGCCACGTAGTCCCCCTGTTTTAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACACCAAAAACGCAGCGTGTGGGTCGCGGGGGTGTCGTAACCAGGGATTATCTGATCGTACGCGACTCCCAGAATCTCTTCgg  <  1:319775/150‑1 (MQ=255)

GTTTTATGCAACTGATAATGGCGGCTTAAAATGTCCTGAATGGTCGAGGAAACCAGGAAGTATTCCTGACGCAGGCGCAGCTCACGCCCGGAGTAGGTGGAGTCATCCGGATACAGTACGCGAGATACGTTCTCGGAGTGGTTTTTATCTTCCACTGCCGCGAAGTAGTCACCCTGGTTGAATTTACCGAGGTTAATTTCGCTACTGGCTTGCGCACTCCACAAACGCAGCGTGTTGGTCGCGTCGGTGTCGTAACCAGGGATTATCTGATCGTAAGCGACTCCCAGAATCTCTTCGG  >  CP009273/3559022‑3559319

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