BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,892,211	C→A	5.9%	G341G (GGC→GGA)	sdaA →	L‑serine dehydratase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,892,211	0	C	A	5.9%	1330.1 / 36.8	336	G341G (GGC→GGA)	sdaA	L‑serine dehydratase 1
Reads supporting (aligned to +/- strand): ref base C (180/136); new base A (10/10); total (190/146)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 345 total aligned reads displayed.

ATCGTTCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCC  >  CP009273/1892062‑1892360

aTCGTTCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGc                                                                                                                                                       <  1:3979812/150‑1 (MQ=255)

aTCGTTCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGc                                                                                                                                                       >  2:5074390/1‑150 (MQ=255)

aTCGTTCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGc                                                                                                                                                       >  1:314933/1‑150 (MQ=255)

     tCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGACGTTGGTTGCCAGGGCGAAGt                                                                                                                                                  >  2:1438885/1‑150 (MQ=255)

     tCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGACAGCCCGGACATCTATACCCGTTACTTTAGGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGt                                                                                                                                                  <  1:1831043/150‑1 (MQ=255)

      ccGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGAGTGGTGCCGTGTCTAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGACGTg                                                                                                                                                 >  2:1095931/1‑150 (MQ=255)

       cGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTCTGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTgg                                                                                                                                                <  2:4889195/150‑1 (MQ=255)

       cGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTgg                                                                                                                                                >  2:4562623/1‑150 (MQ=255)

          cAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGgtg                                                                                                                                             <  2:1596591/150‑1 (MQ=255)

            gTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTAGGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGtt                                                                                                                                           <  1:651409/150‑1 (MQ=255)

             tGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTg                                                                                                                                          >  2:35536/1‑150 (MQ=255)

                 ggCTTACTTTGACCAGTTTTTTTACTCGGTCAGCCCGGACATATTTACCCGTTTATTTATGGCTGCGGGCGCTATTGGTGAATTTTATATATTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTg                                                                                                                                      <  1:2737999/150‑1 (MQ=255)

                    ttACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTc                                                                                                                                   >  2:595005/1‑150 (MQ=255)

                    ttACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTc                                                                                                                                   <  1:1629901/150‑1 (MQ=255)

                     tACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCa                                                                                                                                  >  2:4114388/1‑150 (MQ=255)

                      aCTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCaa                                                                                                                                 >  1:3638202/1‑150 (MQ=255)

                       cTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAAt                                                                                                                                <  1:819553/150‑1 (MQ=255)

                       cTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAAt                                                                                                                                >  2:3261302/1‑150 (MQ=255)

                       cTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAAt                                                                                                                                >  1:3011428/1‑150 (MQ=255)

                          tGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGgc                                                                                                                             >  1:2993140/1‑150 (MQ=255)

                           gACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGgct                                                                                                                            <  1:3395576/150‑1 (MQ=255)

                            aCCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGgctg                                                                                                                           >  1:2929260/1‑150 (MQ=255)

                                  ttATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCgg                                                                                                                     >  1:4411666/1‑150 (MQ=255)

                                  ttATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCgg                                                                                                                     >  2:3868878/1‑150 (MQ=255)

                                  ttATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCgg                                                                                                                     >  1:4505621/1‑150 (MQ=255)

                                   tATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCggg                                                                                                                    >  1:985827/1‑150 (MQ=255)

                                      tGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCt                                                                                                                 >  1:750355/1‑150 (MQ=255)

                                       gAATCGGTCTGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCtt                                                                                                                <  2:506235/150‑1 (MQ=255)

                                       gAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCtt                                                                                                                >  1:3125482/1‑150 (MQ=255)

                                        aaTCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTg                                                                                                               >  2:1369951/1‑150 (MQ=255)

                                            ggTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGa                                                                                                           >  1:3481232/1‑150 (MQ=255)

                                            ggTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGa                                                                                                           >  1:950035/1‑150 (MQ=255)

                                              tCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAAc                                                                                                         >  2:4774833/1‑150 (MQ=255)

                                              tCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAAc                                                                                                         >  2:3745406/1‑150 (MQ=255)

                                               cAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAAct                                                                                                        >  2:2786945/1‑150 (MQ=255)

                                               cAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAAct                                                                                                        <  1:4451056/150‑1 (MQ=255)

                                                aGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAActg                                                                                                       >  2:2991524/1‑150 (MQ=255)

                                                 gCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGAATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAActgc                                                                                                      <  1:2729970/150‑1 (MQ=255)

                                                   ccGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAActgctg                                                                                                    <  2:1938611/150‑1 (MQ=255)

                                                      gACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGc                                                                                                 >  1:1891861/1‑150 (MQ=255)

                                                         aTCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGt                                                                                              >  2:826494/1‑150 (MQ=255)

                                                          tCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTa                                                                                             <  2:3684268/150‑1 (MQ=255)

                                                           cTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAg                                                                                            <  1:100776/150‑1 (MQ=255)

                                                             atccccGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGTGGTAGcc                                                                                          <  1:541133/147‑1 (MQ=255)

                                                             ataCCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGAGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGcc                                                                                          <  2:4193751/150‑1 (MQ=255)

                                                             ataCCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGAGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGcc                                                                                          <  2:1629604/150‑1 (MQ=255)

                                                              taCCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGccc                                                                                         >  1:179531/1‑150 (MQ=255)

                                                                    ttACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACa                                                                                   >  1:4930745/1‑150 (MQ=255)

                                                                    ttACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGAGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACa                                                                                   >  1:3868034/1‑150 (MQ=255)

                                                                         ttATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGttt                                                                              >  1:2141711/1‑150 (MQ=255)

                                                                          tATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTg                                                                             <  2:4716971/150‑1 (MQ=255)

                                                                          tATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTg                                                                             <  2:4716886/150‑1 (MQ=255)

                                                                           aTGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGc                                                                            >  2:1915989/1‑150 (MQ=255)

                                                                           aTGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGAGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGc                                                                            >  2:1999581/1‑150 (MQ=255)

                                                                            tGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCg                                                                           <  1:2071014/150‑1 (MQ=255)

                                                                            gggCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCg                                                                           <  1:4739809/149‑1 (MQ=255)

                                                                                  cGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTggcgg                                                                     >  1:1437069/1‑150 (MQ=255)

                                                                                       gcgATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGaa                                                                >  1:5074606/1‑150 (MQ=255)

                                                                                        cgATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGaaa                                                               <  2:1304588/150‑1 (MQ=255)

                                                                                          aTTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAAtt                                                             <  2:4876675/150‑1 (MQ=255)

                                                                                           ttGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACAGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTg                                                            <  1:3269061/150‑1 (MQ=255)

                                                                                            tGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTgg                                                           <  2:4360215/150‑1 (MQ=255)

                                                                                             ggTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGc                                                          >  1:4698796/1‑150 (MQ=255)

                                                                                              gTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCa                                                         <  1:4433917/150‑1 (MQ=255)

                                                                                                gCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATg                                                       >  1:481885/1‑150 (MQ=255)

                                                                                                gCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATg                                                       <  2:314931/150‑1 (MQ=255)

                                                                                                 cATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATgg                                                      >  2:4437775/1‑150 (MQ=255)

                                                                                                 cATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATgg                                                      <  2:672458/150‑1 (MQ=255)

                                                                                                  aTTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGa                                                     >  1:4138399/1‑150 (MQ=255)

                                                                                                    tGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAac                                                   >  2:2905666/1‑150 (MQ=255)

                                                                                                     gTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAaca                                                  <  2:2207360/150‑1 (MQ=255)

                                                                                                        taAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACaa                                               <  2:2539315/150‑1 (MQ=255)

                                                                                                           aaTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCt                                            >  1:1486934/1‑150 (MQ=255)

                                                                                                            aTGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCtt                                           >  2:3622171/1‑150 (MQ=255)

                                                                                                             tGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTg                                          <  1:1249324/150‑1 (MQ=255)

                                                                                                             tGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTg                                          >  2:2805178/1‑150 (MQ=255)

                                                                                                             tGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTg                                          <  1:2815488/150‑1 (MQ=255)

                                                                                                             tGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTg                                          >  2:3095940/1‑150 (MQ=255)

                                                                                                             tGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTg                                          <  1:1530922/150‑1 (MQ=255)

                                                                                                                aCGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGtt                                       <  1:25836/150‑1 (MQ=255)

                                                                                                                  gCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTa                                     <  1:2629762/150‑1 (MQ=255)

                                                                                                                       tATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCTATGGCTGCTGCTGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTg                                <  1:3324228/150‑1 (MQ=255)

                                                                                                                       tATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTg                                >  1:2106233/1‑150 (MQ=255)

                                                                                                                       tATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTg                                <  2:3996331/150‑1 (MQ=255)

                                                                                                                       tATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTGCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTg                                >  1:5127519/1‑150 (MQ=255)

                                                                                                                        aTTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGc                               >  1:1088238/1‑150 (MQ=255)

                                                                                                                           tCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGAc                            >  1:3920568/1‑150 (MQ=255)

                                                                                                                            ccGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGAcc                           >  1:3386393/1‑150 (MQ=255)

                                                                                                                              ggTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCg                         >  2:3183432/1‑150 (MQ=255)

                                                                                                                              ggTGCGGAAGTAGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCg                         >  2:4289503/1‑150 (MQ=255)

                                                                                                                                tGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGt                       >  1:229553/1‑150 (MQ=255)

                                                                                                                                tGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGt                       >  1:1783311/1‑150 (MQ=255)

                                                                                                                                 gCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGtt                      >  1:3771046/1‑150 (MQ=255)

                                                                                                                                      aGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTgcagg                 >  2:528717/1‑150 (MQ=255)

                                                                                                                                         tggttgCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTgcagggca              <  1:627402/150‑1 (MQ=255)

                                                                                                                                            ttgCCAGGGAGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGt           >  2:1187831/1‑150 (MQ=255)

                                                                                                                                             tgCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGtt          >  2:2616598/1‑150 (MQ=255)

                                                                                                                                              gCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTc         <  1:4562623/150‑1 (MQ=255)

                                                                                                                                                cAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAg       >  1:4251167/1‑150 (MQ=255)

                                                                                                                                                     cGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGcc  <  1:4080559/150‑1 (MQ=255)

ATCGTTCCGGCAGTGCTGGCTTACTATGACCACTTTATTGAATCGGTCAGCCCGGACATCTATACCCGTTACTTTATGGCAGCGGGCGCGATTGGTGCATTGTATAAAATGAACGCCTCTATTTCCGGTGCGGAAGTTGGTTGCCAGGGCGAAGTGGGTGTTGCCTGTTCAATGGCTGCTGCGGGTCTTGCAGAACTGCTGGGCGGTAGCCCGGAACAGGTTTGCGTGGCGGCGGAAATTGGCATGGAACACAACCTTGGTTTAACCTGCGACCCGGTTGCAGGGCAGGTTCAGGTGCC  >  CP009273/1892062‑1892360

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: `‑`