BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	940,480	A→G	5.2%	A1715A (GCT→GCC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	940,480	0	A	G	5.2%	331.7 / 3.9	107	A1715A (GCT→GCC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand): ref base A (58/43); new base G (3/3); total (61/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 107 total aligned reads displayed.

CTGCTCGCCCGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGTGAGGCT  >  NC_002947/940333‑940626

cTGCTCGCCCGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAgc                                                                                                                                                  >  1:1370565/1‑150 (MQ=255)

cTGCTCGCCCGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAgc                                                                                                                                                  >  2:832999/1‑150 (MQ=255)

 tGCTCGCCCGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAgcg                                                                                                                                                 >  1:1372006/1‑150 (MQ=255)

      gCCCGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGAc                                                                                                                                            >  1:1972354/1‑150 (MQ=255)

        ccGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACgg                                                                                                                                          >  2:726616/1‑150 (MQ=255)

          ggCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCAGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCa                                                                                                                                        >  2:2083107/1‑150 (MQ=255)

           gCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCGGCGCTGACGGcct                                                                                                                                       <  1:1706390/150‑3 (MQ=255)

                 cgGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGc                                                                                                                                 >  1:483848/1‑150 (MQ=255)

                  ggCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCt                                                                                                                                <  1:1802700/150‑1 (MQ=255)

                  ggCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCt                                                                                                                                >  2:1912521/1‑150 (MQ=255)

                     tGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTcac                                                                                                                             <  2:1972354/150‑1 (MQ=255)

                                cGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggcg                                                                                                                  >  1:1130368/1‑150 (MQ=255)

                                cGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggcg                                                                                                                  <  1:1126748/150‑1 (MQ=255)

                                 gATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCGAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCAgcgggcgg                                                                                                                 >  2:1768562/1‑150 (MQ=255)

                                 gATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggcgg                                                                                                                 >  1:13745/1‑150 (MQ=255)

                                 gATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggcgg                                                                                                                 >  1:1757848/1‑150 (MQ=255)

                                  aTCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggcggg                                                                                                                <  1:1822771/150‑1 (MQ=255)

                                    cAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGt                                                                                                              >  2:1634202/1‑150 (MQ=255)

                                      ggAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTaa                                                                                                            >  1:357995/1‑150 (MQ=255)

                                       gAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAAt                                                                                                           <  2:1985329/150‑1 (MQ=255)

                                       gAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAAt                                                                                                           >  1:2045678/1‑150 (MQ=255)

                                       gAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACCGCAATGTTGCTCACCGCACCGGGCGGGGTAAt                                                                                                           >  2:2018702/1‑150 (MQ=255)

                                                gTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGgcgc                                                                                                  >  1:2233077/1‑150 (MQ=255)

                                                    gCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCg                                                                                              >  2:147690/1‑150 (MQ=255)

                                                     cTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACAGCACCGGGCGGGGTAATGTCCCGCGGATCgg                                                                                             >  2:1118684/1‑149 (MQ=255)

                                                      tCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGaa                                                                                            >  1:2397229/1‑150 (MQ=255)

                                                      tCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGaa                                                                                            <  1:784635/150‑1 (MQ=255)

                                                      tCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGaa                                                                                            <  1:858690/150‑1 (MQ=255)

                                                      tCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGaa                                                                                            >  1:1080796/1‑150 (MQ=255)

                                                         cacaACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAAtt                                                                                          <  2:2265670/149‑1 (MQ=255)

                                                         cacaACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAAtt                                                                                          >  1:2265670/1‑149 (MQ=255)

                                                         cacaACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTg                                                                                         >  2:1410494/1‑150 (MQ=255)

                                                             acacCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAgc                                                                                     <  2:2233077/150‑1 (MQ=255)

                                                                   gTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCa                                                                               <  2:641966/150‑1 (MQ=255)

                                                                        gATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGc                                                                          >  1:1684732/1‑150 (MQ=255)

                                                                             cGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGtttc                                                                     <  1:1410494/150‑3 (MQ=255)

                                                                              ggTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGcc                                                                    <  2:1757848/150‑1 (MQ=255)

                                                                              ggTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGcc                                                                    <  1:1588109/150‑1 (MQ=255)

                                                                              ggTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGcc                                                                    >  1:1851243/1‑150 (MQ=255)

                                                                                tGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGt                                                                  <  1:2239066/150‑1 (MQ=255)

                                                                                  ccACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCa                                                                <  1:715937/150‑1 (MQ=255)

                                                                                  ccACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCa                                                                <  2:13745/150‑1 (MQ=255)

                                                                                  ccACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCa                                                                <  1:1634202/150‑1 (MQ=255)

                                                                                    tcTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGc                                                              >  2:1652346/2‑150 (MQ=255)

                                                                                    aCTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCaga                                                                             <  1:2186786/135‑1 (MQ=255)

                                                                                    aCTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCaga                                                                             >  2:2186786/1‑135 (MQ=255)

                                                                                    aCTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCa                                                                <  1:1627827/148‑1 (MQ=255)

                                                                                    aCTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGc                                                              >  1:1657047/1‑150 (MQ=255)

                                                                                      tGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTcaccgcac                                                                                                                        >  2:1405228/1‑90 (MQ=255)

                                                                                      tGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTcaccgcac                                                                                                                        <  1:1405228/90‑1 (MQ=255)

                                                                                      tGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCAt                                                            >  1:2073137/1‑150 (MQ=255)

                                                                                        gCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATcc                                                          <  1:1118684/150‑1 (MQ=255)

                                                                                        gCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATcc                                                          >  2:1643463/1‑150 (MQ=255)

                                                                                               gcACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGc                                                   >  1:930280/1‑150 (MQ=255)

                                                                                               gcACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGc                                                   >  1:376428/1‑150 (MQ=255)

                                                                                                 aCCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGt                                                 >  1:462266/1‑150 (MQ=255)

                                                                                                 aCCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGt                                                 >  1:395708/1‑150 (MQ=255)

                                                                                                   cTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGGCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGcc                                                                    >  2:1391696/1‑129 (MQ=255)

                                                                                                   cTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGcc                                                                    <  1:1391696/129‑1 (MQ=255)

                                                                                                   cTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTAc                                               >  1:115060/1‑150 (MQ=255)

                                                                                                   cTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTAc                                               >  2:1203472/1‑150 (MQ=255)

                                                                                                   cTCAACAGTGGCGCCGGGTTCGCCGCGACCACTCAGGACCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTAc                                               >  1:1509537/1‑150 (MQ=255)

                                                                                                    tCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATccac                                                        >  1:2254723/1‑140 (MQ=255)

                                                                                                      aaCAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTc                                            <  2:357995/150‑1 (MQ=255)

                                                                                                        cAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGa                                          <  1:1203472/150‑1 (MQ=255)

                                                                                                        cAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGa                                          >  1:252212/1‑150 (MQ=255)

                                                                                                           tGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAggg                                       <  2:1370565/150‑1 (MQ=255)

                                                                                                               gcCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGc                                      >  1:955554/1‑147 (MQ=255)

                                                                                                                ccGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCCGCATCCACCAAGCGTACGTCGAGGGCTTCg                                  <  1:1643463/150‑1 (MQ=255)

                                                                                                                     ttCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGtt                             <  1:720153/150‑1 (MQ=255)

                                                                                                                      tCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTg                            <  2:930280/150‑1 (MQ=255)

                                                                                                                         ccgcGAACAATCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGcc                         >  2:102051/1‑150 (MQ=255)

                                                                                                                          cgcgACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCa                                                                                       >  1:2216652/1‑87 (MQ=255)

                                                                                                                          cgcgACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCa                                                                                       <  2:2216652/87‑1 (MQ=255)

                                                                                                                           gcgACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTg                       >  1:869229/1‑150 (MQ=255)

                                                                                                                             gACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTggg                     <  2:376428/150‑1 (MQ=255)

                                                                                                                              aCCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGc                    <  2:483848/150‑1 (MQ=255)

                                                                                                                                cGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAg                  <  2:115060/150‑1 (MQ=255)

                                                                                                                                   tCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCg               >  2:596984/1‑150 (MQ=255)

                                                                                                                                     aGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggc                                                                                                                   <  2:1616000/48‑1 (MQ=255)

                                                                                                                                     aGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACcgggc                                                                                                                   >  1:1616000/1‑48 (MQ=255)

                                                                                                                                      gggCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAg                                         >  2:1623570/1‑121 (MQ=255)

                                                                                                                                      gggCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAg                                         <  1:1623570/121‑1 (MQ=255)

                                                                                                                                      gggCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACa            <  1:102051/150‑1 (MQ=255)

                                                                                                                                        gCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAgg          <  1:2172901/150‑1 (MQ=255)

                                                                                                                                        gCCAGGCCGTCAGAGCCGACGGCAATGGTGCGCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAgg          <  2:631810/150‑1 (MQ=255)

                                                                                                                                        acattgCCGTCGGCGCTGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTgg                      >  1:1173080/6‑138 (MQ=255)

                                                                                                                                        acattgCCGTCGGCGCTGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTgg                      <  2:1173080/133‑1 (MQ=255)

                                                                                                                                        acattgCCGTCGGCGCTGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAgg          >  1:2081833/6‑150 (MQ=255)

                                                                                                                                        acattgCCGTCGGCGCTGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAgg          >  1:461262/6‑150 (MQ=255)

                                                                                                                                         ccAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGc                                                                          <  1:1201984/85‑1 (MQ=255)

                                                                                                                                         ccAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGc                                                                          >  2:1201984/1‑85 (MQ=255)

                                                                                                                                         ccAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAggg         >  1:1670131/1‑150 (MQ=255)

                                                                                                                                         ccAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGCCTTCGCCGTTGGCCTGGGCAGGCGACAggg         >  2:316538/1‑150 (MQ=255)

                                                                                                                                          cAGGCCGTCAGCGCCGACGGCTATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGt        >  2:2178994/1‑150 (MQ=255)

                                                                                                                                          cAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGTAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGt        >  2:1389720/1‑150 (MQ=255)

                                                                                                                                           aGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGTg       >  2:617383/1‑150 (MQ=255)

                                                                                                                                            ttccGTCAGCGCCGACGGCAATGTTGCTCATCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGTGa      <  2:395708/148‑1 (MQ=255)

                                                                                                                                            tgCCGTCGGCGCTGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGTGa      <  2:461262/149‑1 (MQ=255)

                                                                                                                                                gTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGTGAGGCt  <  2:1080796/150‑1 (MQ=255)

CTGCTCGCCCGGCTGCGCGGCTGGGGCGAGGTCGATCAGGAAGGTGCCGTTGGCTCCCACAACACCCGTGCCGATCACGGTGCCACTGGCATCGCGCACCTCAACAGTGGCGCCGGGTTCGCCGCGACCGCTCAGGGCCAGGCCGTCAGCGCCGACGGCAATGTTGCTCACCGCACCGGGCGGGGTAATGTCCGGCGCATCGAATTGCAGCGGCGCAGAGGCGTTGCCGTCAGCATCCACCAAGCGTACGTCGAGGGCTTCGCCGTTGGCCTGGGCAGGCGACAGGGTGAGGCT  >  NC_002947/940333‑940626

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