BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,003,722	G→T	7.5%	E141* (GAA→TAA)	yedK →	DUF159 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,003,722	0	G	T	7.5%	1087.8 / 37.6	305	E141* (GAA→TAA)	yedK	DUF159 family protein
Reads supporting (aligned to +/- strand): ref base G (159/123); new base T (19/4); total (178/127)
Fisher's exact test for biased strand distribution p-value = 1.51e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

Only 100 of 308 total aligned reads displayed.

GGCAACATGGTCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGAAATCGCGGCC  >  CP009273/2003573‑2003871

   aaCATGGTCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAg                                                                                                                                                    >  1:3672855/1‑150 (MQ=255)

      aTGGTCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCg                                                                                                                                                 >  2:319379/1‑150 (MQ=255)

       tGGTCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGa                                                                                                                                                >  1:5018656/1‑150 (MQ=255)

         gTCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAg                                                                                                                                              <  1:3202775/150‑1 (MQ=255)

          tCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAgg                                                                                                                                             <  1:1062957/150‑1 (MQ=255)

          tCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAgg                                                                                                                                             <  1:5161423/150‑1 (MQ=255)

            gggCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGAt                                                                                                                                           <  1:2439211/150‑1 (MQ=255)

              gCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCCATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGAttt                                                                                                                                         <  1:1085152/150‑1 (MQ=255)

               cAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGAtttt                                                                                                                                        <  2:2457557/150‑1 (MQ=255)

                aaTCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGAttttt                                                                                                                                       <  1:1807919/150‑1 (MQ=255)

                    tGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATc                                                                                                                                   >  2:5153947/1‑150 (MQ=255)

                     gTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCg                                                                                                                                  >  2:815160/1‑150 (MQ=255)

                       tttGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTc                                                                                                                                >  1:1754723/1‑150 (MQ=255)

                       tttGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTc                                                                                                                                >  2:457457/1‑150 (MQ=255)

                         tGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCAc                                                                                                                              <  1:2369224/150‑1 (MQ=255)

                          gCCGATTGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACATCCTTTTTTTATCTATCGTTCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCAct                                                                                                                             <  1:4604678/150‑1 (MQ=255)

                          gCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCAct                                                                                                                             >  1:1204919/1‑150 (MQ=255)

                            cGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCActgc                                                                                                                           <  2:1492956/150‑1 (MQ=255)

                             gATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCActgct                                                                                                                          <  2:3739648/150‑1 (MQ=255)

                             gATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCActgct                                                                                                                          >  1:1428330/1‑150 (MQ=255)

                             gATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCActgct                                                                                                                          >  1:3167559/1‑150 (MQ=255)

                              aTGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTATTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCActgctg                                                                                                                         >  1:4276591/1‑150 (MQ=255)

                                ggCTGGTGTGAGTGGATATAAGATGGCGACATAAGACAGCCTTTTTTTATCTATCGTGCTGATGTTCAACCTTTTTTTCTGGCCGCGATAGGTAGCCCACCATTTGTTCGTGGTGTCGATGCCGAAGGATTTTTGATCGTCTCTGCTGCg                                                                                                                       <  2:1481498/150‑1 (MQ=255)

                                      tttGAGTGGAAAAAAGAAGGCGACAAAAAACAGACTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGAt                                                                                                                 <  1:640657/150‑1 (MQ=255)

                                      attGAGTGTAAAAAAGAAGGCGACAAAAATCAGCCTTTTTTTATCTATCGTGCTGATGGACAAACTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGAt                                                                                                                 <  1:2893107/149‑1 (MQ=255)

                                       ttGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATc                                                                                                                <  1:857981/150‑1 (MQ=255)

                                          aGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAg                                                                                                             >  1:307991/1‑150 (MQ=255)

                                            tGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGt                                                                                                           <  2:833566/150‑1 (MQ=255)

                                              gAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCt                                                                                                         <  1:3515107/150‑1 (MQ=255)

                                               aaaaaaGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTg                                                                                                        <  1:4027090/150‑1 (MQ=255)

                                                 aaaaGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGt                                                                                                      >  2:3029695/1‑150 (MQ=255)

                                                  aaaGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGGAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTa                                                                                                     >  1:1656450/1‑150 (MQ=255)

                                                   aagaagGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAg                                                                                                    >  2:2414983/1‑150 (MQ=255)

                                                     gaagGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGat                                                                                                  >  1:4660386/1‑150 (MQ=255)

                                                          cGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCa                                                                                             >  2:1232504/1‑150 (MQ=255)

                                                          cGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCa                                                                                             <  2:4736221/150‑1 (MQ=255)

                                                           gACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTAATTTTATGTCCGCGATAGGTAGCACACCATATGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCAGCGGCACATCAAGCTCTGGTAGATATTCAt                                                                                            >  2:4059726/1‑150 (MQ=255)

                                                           gACAAAAAACAGCCGTTTTTTATCTATCGTGCTGATTGACAAACTATTTTTATTTCCGCGATATGTAGCACACAATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCAt                                                                                            >  1:742254/1‑150 (MQ=255)

                                                            aCAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATg                                                                                           <  2:3550489/150‑1 (MQ=255)

                                                            aCAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATg                                                                                           >  2:2542111/1‑150 (MQ=255)

                                                            aCAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATg                                                                                           <  2:1085423/150‑1 (MQ=255)

                                                              aaaaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAc                                                                                         >  2:4326883/1‑150 (MQ=255)

                                                              aaaaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAc                                                                                         <  2:4977309/150‑1 (MQ=255)

                                                              aaaaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGAAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGACACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAc                                                                                         >  2:778757/1‑150 (MQ=255)

                                                               aaaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAcc                                                                                        <  1:2002010/150‑1 (MQ=255)

                                                               aaaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAcc                                                                                        <  2:4949616/150‑1 (MQ=255)

                                                               aaaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAcc                                                                                        >  1:1208040/1‑150 (MQ=255)

                                                                aaaaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAccg                                                                                       >  2:754796/1‑150 (MQ=255)

                                                                  aaCAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGAccgcc                                                                                     >  1:2654411/1‑150 (MQ=255)

                                                                    cAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACcgccgc                                                                                   >  2:3913773/1‑150 (MQ=255)

                                                                     aGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCgccgcc                                                                                  >  2:4423024/1‑150 (MQ=255)

                                                                         tttttttATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCTGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACt                                                                              <  2:345308/150‑1 (MQ=255)

                                                                            ttttATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGt                                                                           >  2:137954/1‑150 (MQ=255)

                                                                            ttttATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGt                                                                           <  1:4162697/150‑1 (MQ=255)

                                                                             tttATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTa                                                                          <  2:2929233/150‑1 (MQ=255)

                                                                              ttATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTAc                                                                         <  1:5206893/150‑1 (MQ=255)

                                                                               tatctatcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTTCt                                                                        >  1:194986/1‑150 (MQ=255)

                                                                               tatctatcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACt                                                                        <  2:2234780/150‑1 (MQ=255)

                                                                                 tctatcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGt                                                                      <  2:5032925/150‑1 (MQ=255)

                                                                                  ctatcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTc                                                                     <  2:4730244/150‑1 (MQ=255)

                                                                                  cgatcGTGCTGATGGACAACCTATTTTGATGGCCGCGAGAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGCTCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTTCTGTc                                                                     <  2:1848179/148‑1 (MQ=255)

                                                                                   tatcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCg                                                                    >  2:4002104/1‑150 (MQ=255)

                                                                                    atcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGc                                                                   >  2:3612120/1‑150 (MQ=255)

                                                                                    atcGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACAGCCGCCCACTGGTACTGTCGc                                                                   >  1:435736/1‑150 (MQ=255)

                                                                                      cGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCa                                                                 <  1:1725348/150‑1 (MQ=255)

                                                                                      cGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCa                                                                 >  2:2668610/1‑150 (MQ=255)

                                                                                        tGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGa                                                               <  1:4048921/150‑1 (MQ=255)

                                                                                           tGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGc                                                            >  1:479061/1‑150 (MQ=255)

                                                                                            gATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCt                                                           <  1:4538848/150‑1 (MQ=255)

                                                                                             aTGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTg                                                          <  2:915920/150‑1 (MQ=255)

                                                                                             aTGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTg                                                          <  1:2299290/150‑1 (MQ=255)

                                                                                              tGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTgc                                                         <  2:288200/150‑1 (MQ=255)

                                                                                              tGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTgc                                                         >  2:539732/1‑150 (MQ=255)

                                                                                               ggACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTgcg                                                        <  2:2546990/150‑1 (MQ=255)

                                                                                                    aCCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGa                                                   >  2:3082007/1‑150 (MQ=255)

                                                                                                    aCCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGa                                                   >  2:1597104/1‑150 (MQ=255)

                                                                                                      cTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAAt                                                 >  2:2016700/1‑150 (MQ=255)

                                                                                                      cTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAAt                                                 >  2:613803/1‑150 (MQ=255)

                                                                                                      cTATTTTAATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGAACGACACTTCTGCGGCATATAAAGGTCTGGTAGATATTCATGAACGCCGCACACTGGTACTGTCGCCAGAAGCTGCGCGGGAAt                                                 >  2:3332819/1‑150 (MQ=255)

                                                                                                        aTTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATgg                                               >  2:3062285/1‑150 (MQ=255)

                                                                                                        aTTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATgg                                               >  1:523856/1‑150 (MQ=255)

                                                                                                           tttATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATg                                            <  2:2290676/150‑1 (MQ=255)

                                                                                                                      gATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGAt                                 <  1:2714779/150‑1 (MQ=255)

                                                                                                                       aTAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGAtt                                <  1:1274214/150‑1 (MQ=255)

                                                                                                                        tAGGTAGCACACCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTa                               <  1:222240/150‑1 (MQ=255)

                                                                                                                          ggTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGc                             <  1:3976773/150‑1 (MQ=255)

                                                                                                                           gTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCg                            <  1:3939999/150‑1 (MQ=255)

                                                                                                                            tAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCgg                           >  2:1000843/1‑150 (MQ=255)

                                                                                                                             aGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGt                          >  1:1038609/1‑150 (MQ=255)

                                                                                                                             aGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGt                          >  2:4637541/1‑150 (MQ=255)

                                                                                                                               cacaCCATTTGAGCGTGGTGACTAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCTGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTaa                        <  1:1148465/150‑1 (MQ=255)

                                                                                                                                acacCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTaaa                       >  1:1608853/1‑150 (MQ=255)

                                                                                                                                     aTTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGc                  >  1:3638499/1‑150 (MQ=255)

                                                                                                                                     aTTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGc                  >  2:1805636/1‑150 (MQ=255)

                                                                                                                                            cGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGaa           >  2:1337803/1‑150 (MQ=255)

                                                                                                                                            cGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGaa           >  2:2382638/1‑150 (MQ=255)

                                                                                                                                               ggtgACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGAAATc        <  2:4851037/150‑1 (MQ=255)

                                                                                                                                                  gACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGAAATcgcg     >  1:3397281/1‑150 (MQ=255)

                                                                                                                                                   aCGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGAAATCGCgg    >  2:1120945/1‑150 (MQ=255)

                                                                                                                                                     gAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCAAGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGAAATCGCGGcc  >  2:1664042/1‑150 (MQ=255)

GGCAACATGGTCGGGCAATCTGTTTTGCCGATGGCTGGTTTGAGTGGAAAAAAGAAGGCGACAAAAAACAGCCTTTTTTTATCTATCGTGCTGATGGACAACCTATTTTTATGGCCGCGATAGGTAGCACACCATTTGAGCGTGGTGACGAAGCCGAAGGATTTTTGATCGTCACTGCTGCGGCAGATCAAGGTCTGGTAGATATTCATGACCGCCGCCCACTGGTACTGTCGCCAGAAGCTGCGCGGGAATGGATGCGGCAAGAGATTAGCGGTAAAGAAGCCTCAGAAATCGCGGCC  >  CP009273/2003573‑2003871

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