BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,262,311	G→A	10.1%	T158T (ACC→ACT) ‡	PP_2857 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,262,311	0	G	A	10.1%	352.0 / 23.3	129	T158T (ACC→ACT) ‡	PP_2857	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (65/50); new base A (6/7); total (71/58)
Fisher's exact test for biased strand distribution p-value = 5.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 134 total aligned reads displayed.

TCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGC  >  NC_002947/3262163‑3262460

tctcCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTgg                                                                                                                                                      <  2:4187504/150‑1 (MQ=255)

  tcCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTa                                                                                                                                                    >  2:4421674/1‑150 (MQ=255)

      tACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGa                                                                                                                                                >  1:4737841/1‑150 (MQ=255)

      tACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCAACCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGa                                                                                                                                                >  2:604822/1‑150 (MQ=255)

          cTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTg                                                                                                                                            <  1:1572964/150‑1 (MQ=255)

           tATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTg                                                                                                                                                       <  1:3276869/138‑1 (MQ=255)

           tATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTg                                                                                                                                                       >  2:3276869/1‑138 (MQ=255)

           tATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCCGCTACCGCTGAAGGCCTGGACACTGGTAGTGACCTgg                                                                                                                                           >  2:3034387/1‑150 (MQ=255)

            aTTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAGGc                                                                                                                                          <  1:4151946/150‑1 (MQ=255)

             ttAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCAGGTCACTACCAGTGACCAGGcc                                                                                                                                         <  1:54830/150‑2 (MQ=255)

                  ttAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTAACGGTGAAGGCCAGGTCACTGGTAGTGACCTTGCGATAGc                                                                                                                                    >  2:222116/1‑150 (MQ=255)

                  ttAAGTGTATTTCTTAACTGCAGGGTCAATCCGCGAGACGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCACTCGGCTACCGGTGAAGGCCTGGTCACTTGTAGTTACCTGGCGATAGc                                                                                                                                    <  2:5388053/150‑1 (MQ=255)

                     aGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCctg                                                                                                                                 <  2:2368438/150‑1 (MQ=255)

                      gtgtCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCACGGCTACCGGTGAAGGCCAGGTCACTGGTAGTGACCTGGCGATAGCctgc                                                                                                                                >  2:3537806/1‑150 (MQ=255)

                       tgtCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCctgct                                                                                                                               <  2:3719312/150‑1 (MQ=255)

                           gTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGaa                                                                                                                           >  1:5098419/1‑150 (MQ=255)

                            ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGTCTACCGGTGAAGGCCTGGTCACTACCAGTGACCAGGcc                                                                                                                                         <  2:2817821/135‑2 (MQ=255)

                            ttCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTACCAGTGACCAGGcc                                                                                                                                         >  1:2817821/1‑134 (MQ=255)

                                  aCTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGt                                                                                                                    <  2:2593616/150‑1 (MQ=255)

                                    tGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTg                                                                                                                  >  1:4939309/1‑150 (MQ=255)

                                        ttgagcaTCCGGGTCGCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTcc                                                                                                              >  1:4140459/7‑150 (MQ=255)

                                             aaTCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGa                                                                                                           <  1:556717/148‑1 (MQ=255)

                                             aaTCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgt                                                                                                         >  2:3503235/1‑150 (MQ=255)

                                             aaTCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgt                                                                                                         >  1:1909227/1‑150 (MQ=255)

                                               tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGACTGTCAGTTGCTCCTGAgttg                                                                                                       >  2:3611225/1‑150 (MQ=255)

                                               tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttg                                                                                                       >  2:3611222/1‑150 (MQ=255)

                                               tCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTAGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttg                                                                                                       >  2:4725031/1‑150 (MQ=255)

                                                 cgcgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgtt                                                                                                     <  2:722628/150‑1 (MQ=255)

                                                  gcgACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGttgttg                                                                                                    >  2:2768352/1‑150 (MQ=255)

                                                       ccGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTa                                                                                                                                                    >  1:4993783/1‑97 (MQ=255)

                                                        cGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAgc                                                                                              <  1:2883219/150‑1 (MQ=255)

                                                              cTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACg                                                                                        >  1:2076915/1‑150 (MQ=255)

                                                               tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGt                                                                                       >  1:1253123/1‑150 (MQ=255)

                                                               tCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGt                                                                                       >  1:1253496/1‑150 (MQ=255)

                                                                 aaGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTc                                                                                     >  1:3713758/1‑150 (MQ=255)

                                                                 aaGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTc                                                                                     >  2:230947/1‑150 (MQ=255)

                                                                   ggTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTg                                                                                   <  2:4553066/150‑1 (MQ=255)

                                                                   ggTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTg                                                                                   >  1:706378/1‑150 (MQ=255)

                                                                    gtgtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGt                                                                                  >  2:4934761/1‑150 (MQ=255)

                                                                     tgtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGtt                                                                                 <  1:3878904/150‑1 (MQ=255)

                                                                      gtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTg                                                                                <  2:181268/150‑1 (MQ=255)

                                                                      gtTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTg                                                                                <  2:2910621/150‑1 (MQ=255)

                                                                        tGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCa                                                                              <  2:1516733/150‑1 (MQ=255)

                                                                         ggCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGt                                                                                                                                                  <  2:2817867/81‑1 (MQ=255)

                                                                         ggCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACAGGTGAAGGCCTGGTCACTGGTAGt                                                                                                                                                  >  1:2817867/1‑81 (MQ=255)

                                                                           ccATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                           <  2:2286356/150‑1 (MQ=255)

                                                                           ccATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                           <  2:2066005/150‑1 (MQ=255)

                                                                            cATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                          <  2:1056342/150‑1 (MQ=255)

                                                                             aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAgttgt                                                                                                      >  2:633069/1‑121 (MQ=255)

                                                                             aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAgg                                                                         >  2:4802980/1‑150 (MQ=255)

                                                                             aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAgg                                                                         >  2:1608386/1‑150 (MQ=255)

                                                                             aTTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAgg                                                                         <  2:1493838/150‑1 (MQ=255)

                                                                              ttCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAAt                                                                                                                          <  2:1902833/100‑1 (MQ=255)

                                                                              ttCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAAt                                                                                                                          >  1:1902833/1‑100 (MQ=255)

                                                                              ttCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGa                                                                        >  1:3074784/1‑150 (MQ=255)

                                                                                cGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACg                                                                      <  2:4950050/150‑1 (MQ=255)

                                                                                        cACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTgtcgtc                                                              >  2:127484/1‑150 (MQ=255)

                                                                                         aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGtcgtcg                                                             >  2:144453/1‑150 (MQ=255)

                                                                                         aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGtcgtcg                                                             <  1:1599702/150‑1 (MQ=255)

                                                                                         aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGtcgtcg                                                             >  2:3017451/1‑150 (MQ=255)

                                                                                         aCCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTGCTGTTGCATGAGGACGCTGtcgtcg                                                             >  2:1894030/1‑150 (MQ=255)

                                                                                          ccccGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGa                                                            <  1:4718848/150‑1 (MQ=255)

                                                                                                   tctGGTTCAACTGTATCACCCCTCGGCTACCGGTGAAGGCCTGGTCTCTGTTAGTGACCTGGCGATATCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCttttct                                                   <  1:3194164/150‑4 (MQ=255)

                                                                                                   tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCa                                                                                                                      >  2:4000077/1‑83 (MQ=255)

                                                                                                   tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCa                                                                                                                      <  1:4000077/83‑1 (MQ=255)

                                                                                                   tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                   <  2:4654288/150‑1 (MQ=255)

                                                                                                   tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                   >  1:4732256/1‑150 (MQ=255)

                                                                                                   tctGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCg                                                   >  1:3363343/1‑150 (MQ=255)

                                                                                                       gTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtg                                               <  2:2126308/150‑1 (MQ=255)

                                                                                                        ttCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTtgt                                              <  1:581931/150‑1 (MQ=255)

                                                                                                           aaCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                           <  1:4502673/118‑1 (MQ=255)

                                                                                                           aaCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGa                                                                           >  2:4502673/1‑118 (MQ=255)

                                                                                                           aaCTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCa                                           >  1:1055583/1‑150 (MQ=255)

                                                                                                                  tCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATc                                                                                                 <  1:514199/89‑1 (MQ=255)

                                                                                                                  tCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATc                                                                                                 >  2:514199/1‑89 (MQ=255)

                                                                                                                      cccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGAt                                >  1:2959325/1‑150 (MQ=255)

                                                                                                                      cccTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGAt                                <  2:5137101/150‑1 (MQ=255)

                                                                                                                           ggCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCa                           <  2:3594039/150‑1 (MQ=255)

                                                                                                                           ggCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCa                           <  2:4737841/150‑1 (MQ=255)

                                                                                                                            gCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAg                          >  2:2569473/1‑150 (MQ=255)

                                                                                                                                 cGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGt                     >  2:600748/1‑150 (MQ=255)

                                                                                                                                 cGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGt                     >  1:1765456/1‑150 (MQ=255)

                                                                                                                                   cataaGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTAGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTg                   >  2:3663189/4‑150 (MQ=255)

                                                                                                                                    tGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTgg                  <  2:2917450/150‑1 (MQ=255)

                                                                                                                                     gAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGc                 <  1:4421674/150‑1 (MQ=255)

                                                                                                                                         gCCAGGTCACTACCAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGc             >  2:3258483/1‑150 (MQ=255)

                                                                                                                                           caggTCACTACCAGTGACCAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                          >  1:3268983/3‑87 (MQ=38)

                                                                                                                                           caggTCACTACCAGTGACCAGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAg                                                                          <  2:3268983/85‑1 (MQ=38)

                                                                                                                                           cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCCGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg           >  2:1700209/1‑150 (MQ=255)

                                                                                                                                           cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg           >  2:2883391/1‑150 (MQ=255)

                                                                                                                                           cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg           >  1:4025432/1‑150 (MQ=255)

                                                                                                                                           cTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTg           >  1:4936565/1‑150 (MQ=255)

                                                                                                                                            tGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTgg          >  1:5178438/1‑150 (MQ=255)

                                                                                                                                            tGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTgg          <  1:4649028/150‑1 (MQ=255)

                                                                                                                                               tCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGgtg       <  2:2718085/150‑1 (MQ=255)

                                                                                                                                                  cTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTg    <  1:2245631/150‑1 (MQ=255)

                                                                                                                                                   tGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTgg   >  2:4040883/1‑150 (MQ=255)

                                                                                                                                                   tGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTgg   >  2:722427/1‑150 (MQ=255)

                                                                                                                                                    ggTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGc  <  1:222116/150‑1 (MQ=255)

                                                                                                                                                    ggTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGc  <  1:222196/150‑1 (MQ=255)

TCTCCGTACGCTATTAGGTTAAGTGTCGTTCTTAACTGCAGGGTCAATCCGCGACCCGGATGCTCAAGGTGTTGGCCATTCGGTTTCCCACCCCGGCACTCTGGTTCAACTGGATCACCCCTCGGCTACCGGTGAAGGCCTGGTCACTGGTAGTGACCTGGCGATAGCCTGCTGGAATGTCAGTTGCTCCTGAGTTGTTGATCAGCGCCACGTTCTGTTGCATGAGGACGCTGTCGTCGATGCTTTGCGGTTGTGCACCAATGCTGATGCTCAGTGCGTTGGCTTGCTGGGTGTTGGC  >  NC_002947/3262163‑3262460

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