BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,453,165	T→A	5.2%	L47Q (CTG→CAG)	PP_1265 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,453,165	0	T	A	5.2%	1380.7 / 38.9	403	L47Q (CTG→CAG)	PP_1265	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (184/197); new base A (12/9); total (197/206)
Fisher's exact test for biased strand distribution p-value = 5.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 419 total aligned reads displayed.

GGATTGCACCATGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGCCAGTTGTGGCTG  >  NC_002947/1453016‑1453314

 gATTGCACCATGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCgctg                                                                                                                                                      >  2:620863/1‑150 (MQ=255)

    tGCACCATGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGctgctg                                                                                                                                                   >  1:6635317/1‑150 (MQ=255)

     gCACCATGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCtgctgc                                                                                                                                                  <  2:258114/150‑1 (MQ=255)

       cccATGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGATTGCCCTTGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTgcgc                                                                                                                                                <  2:7975119/149‑1 (MQ=255)

       aCCATGCCCCGTGAGATCGCCTTCCATGGCGTGTACAGGCCCACCATGACCTGGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTgcgc                                                                                                                                                <  2:8333637/150‑1 (MQ=255)

          aTGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTg                                                                                                                                             >  1:2541805/1‑150 (MQ=255)

           tGCCCCGTGAGATCGCCTTCCCTGGCGTGTACATGCCCACCATGACCTTGATGGTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGa                                                                                                                                            <  1:683283/150‑1 (MQ=255)

           tGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGa                                                                                                                                            <  1:6091234/150‑1 (MQ=255)

             ccccGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGTCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGc                                                                                                                                          <  2:1683163/150‑1 (MQ=255)

             ccccGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGc                                                                                                                                          <  2:8526909/150‑1 (MQ=255)

                   gagaTCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTc                                                                                                                                    <  2:2182439/150‑1 (MQ=255)

                   gagaTCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTc                                                                                                                                    <  1:7557800/150‑1 (MQ=255)

                     gaTCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGt                                                                                                                                  <  1:8087903/150‑1 (MQ=255)

                        cGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTg                                                                                                                               >  2:6959004/1‑150 (MQ=255)

                        cGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTg                                                                                                                               <  1:1038856/150‑1 (MQ=255)

                          ccTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCt                                                                                                                             <  2:5728364/150‑1 (MQ=255)

                           cTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCtt                                                                                                                            <  1:2065872/150‑1 (MQ=255)

                             tCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGt                                                                                                                          <  1:8531353/150‑1 (MQ=255)

                             tCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGATACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGt                                                                                                                          <  2:2224416/150‑1 (MQ=255)

                              ccATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGtt                                                                                                                         >  1:1105114/1‑150 (MQ=255)

                               cATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGCTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTc                                                                                                                        <  1:6404218/150‑1 (MQ=255)

                                    cGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGgcgc                                                                                                                   <  2:7471762/150‑1 (MQ=255)

                                        tACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATg                                                                                                               <  1:3114665/150‑1 (MQ=255)

                                          cATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGgc                                                                                                             >  1:7144082/1‑150 (MQ=255)

                                           aTGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCgctg                                                                                                                                                      <  1:5628105/108‑1 (MQ=255)

                                           aTGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGgcg                                                                                                            >  1:4901503/1‑150 (MQ=255)

                                               ccaccaTGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGctgt                                                                                                        <  1:7427871/150‑1 (MQ=255)

                                               ccaccaTGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGctgt                                                                                                        <  2:7268591/150‑1 (MQ=255)

                                                 accaTGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGctgtct                                                                                                      <  2:5036585/150‑1 (MQ=255)

                                                   caTGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGctgtctgt                                                                                                    >  1:5818167/1‑150 (MQ=255)

                                                      gACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTa                                                                                                 <  2:3854181/150‑1 (MQ=255)

                                                       aCCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTAc                                                                                                <  2:2683501/150‑1 (MQ=255)

                                                         cTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTg                                                                                              >  1:2628557/1‑150 (MQ=255)

                                                          ttGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTgg                                                                                             <  2:6735505/150‑1 (MQ=255)

                                                           tGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGt                                                                                            <  1:7257728/150‑1 (MQ=255)

                                                             aTGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTga                                                                                          <  2:4139139/150‑1 (MQ=255)

                                                             aTGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTga                                                                                          <  2:7664415/150‑1 (MQ=255)

                                                              tGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTgag                                                                                         >  2:489815/1‑150 (MQ=255)

                                                                  ccTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATc                                                                                     >  2:660113/1‑150 (MQ=255)

                                                                   cTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATcc                                                                                    >  1:3767277/1‑150 (MQ=255)

                                                                    tGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCt                                                                                   >  1:2181260/1‑150 (MQ=255)

                                                                      tttGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTc                                                                                 <  1:5195557/150‑1 (MQ=255)

                                                                           ccTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGa                                                                            >  1:4453630/1‑150 (MQ=255)

                                                                            cTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGaa                                                                           >  1:7654342/1‑150 (MQ=255)

                                                                             tGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGaaa                                                                          >  2:2790301/1‑150 (MQ=255)

                                                                               ggcctggcctGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGaaaaa                                                                        <  1:5067214/150‑1 (MQ=255)

                                                                                  ctggcctggGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGtt                                                                     >  2:4744758/1‑150 (MQ=255)

                                                                                   tggcctggGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGttc                                                                    >  2:4243807/1‑150 (MQ=255)

                                                                                   tggcctggGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGttc                                                                    >  1:1398316/1‑150 (MQ=255)

                                                                                   tggcctggGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGttc                                                                    <  2:5644616/150‑1 (MQ=255)

                                                                                       ctggggcctggACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCa                                                                >  2:2000943/1‑150 (MQ=255)

                                                                                       ctggggcctggACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCa                                                                <  2:7583575/150‑1 (MQ=255)

                                                                                       ctggggcctggACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCa                                                                <  2:3534650/150‑1 (MQ=255)

                                                                                          ggacctggACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGcc                                                             >  1:3504417/4‑150 (MQ=255)

                                                                                                 gACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGc                                                      >  1:2448003/1‑150 (MQ=255)

                                                                                                   ccGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCa                                                    <  2:5394560/150‑1 (MQ=255)

                                                                                                    cGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCAc                                                   <  1:7318005/150‑1 (MQ=255)

                                                                                                      gTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCAccc                                                 >  2:4311117/1‑150 (MQ=255)

                                                                                                      gTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCAccc                                                 >  2:5497742/1‑150 (MQ=255)

                                                                                                       ttCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCct                                                <  1:4723117/150‑1 (MQ=255)

                                                                                                         cATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCctgc                                              <  1:4184848/150‑1 (MQ=255)

                                                                                                          aTTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCctgct                                             <  1:8728503/150‑1 (MQ=255)

                                                                                                           ttGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCctgctg                                            >  1:2927834/1‑150 (MQ=255)

                                                                                                            tGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctgg                                           <  2:3968749/150‑1 (MQ=255)

                                                                                                              ccagccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtg                                         >  2:7002663/1‑150 (MQ=255)

                                                                                                                agccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtgct                                       >  1:5429527/1‑150 (MQ=255)

                                                                                                                 gccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtgctg                                      <  2:3767623/150‑1 (MQ=255)

                                                                                                                 gccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtgctg                                      <  1:2745932/150‑1 (MQ=255)

                                                                                                                 gccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtgctg                                      <  1:8202285/150‑1 (MQ=255)

                                                                                                                  ccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtgctgg                                     >  2:3454891/1‑150 (MQ=255)

                                                                                                                  ccaTGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCtgctggtgctgg                                     >  1:4292829/1‑150 (MQ=255)

                                                                                                                    atgatgGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTc                                   >  2:7198807/1‑150 (MQ=255)

                                                                                                                    atgatgGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTc                                   <  2:1214637/150‑1 (MQ=255)

                                                                                                                    atgatgGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCGACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTc                                   <  2:7517408/150‑1 (MQ=255)

                                                                                                                      gatgGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGc                                 <  2:887517/150‑1 (MQ=255)

                                                                                                                      gatgGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGc                                 >  2:4213383/1‑150 (MQ=255)

                                                                                                                      gatgGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGc                                 <  1:3363046/150‑1 (MQ=255)

                                                                                                                         ggCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGc                              <  1:2257794/150‑1 (MQ=255)

                                                                                                                            tACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGt                           <  2:8107092/150‑1 (MQ=255)

                                                                                                                              ccGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTgg                         <  2:7772396/150‑1 (MQ=255)

                                                                                                                              ccGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTgg                         <  1:481218/150‑1 (MQ=255)

                                                                                                                               cGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGc                        >  1:1316974/1‑150 (MQ=255)

                                                                                                                               cGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGc                        >  1:4551746/1‑150 (MQ=255)

                                                                                                                               cGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGc                        <  2:642072/150‑1 (MQ=255)

                                                                                                                                gCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGcc                       >  2:5901391/1‑150 (MQ=255)

                                                                                                                                 cttcttCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCa                      >  1:4997626/1‑150 (MQ=255)

                                                                                                                                  ttcttcTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCAt                     <  1:6712226/150‑1 (MQ=255)

                                                                                                                                    cttctGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCg                   <  2:4417267/150‑1 (MQ=255)

                                                                                                                                    cttctGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCg                   >  1:2912430/1‑150 (MQ=255)

                                                                                                                                    cttctGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCg                   <  2:1894153/150‑1 (MQ=255)

                                                                                                                                     ttctGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCgg                  >  1:3749547/1‑150 (MQ=255)

                                                                                                                                     ttctGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCGTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGcccgcgg                  >  2:156886/1‑145 (MQ=255)

                                                                                                                                       ctGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTc                <  1:7003573/150‑1 (MQ=255)

                                                                                                                                         ggCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGc              >  1:5948135/1‑150 (MQ=255)

                                                                                                                                         ggCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGc              <  2:2919473/150‑1 (MQ=255)

                                                                                                                                         ggCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGc              >  1:177171/1‑150 (MQ=255)

                                                                                                                                           cACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGCCa            <  1:888038/150‑1 (MQ=255)

                                                                                                                                             tcaggCGCAGCAGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGAGCTGGTCGCTGCCGTGGCCATCGGTCGCCAGt          >  2:2403620/4‑150 (MQ=255)

                                                                                                                                              ccGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGCCAGtt         >  1:7904004/1‑150 (MQ=255)

                                                                                                                                                 gcgcTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGCCAGTtgtg      <  2:3661952/150‑1 (MQ=255)

GGATTGCACCATGCCCCGTGAGATCGCCTTCCATGGCGTGTACATGCCCACCATGACCTTGATGTTCCTGTTTGCCCTGGGCCTGGCCTGGGGCCTGGACCGGTTCATTGCCAGCCATGATGGCTACCGCTTCTTCTGGCACCCGGCGCTGCTGCGCCTGAGCCTGTTCGTCTGCTTGTTCGGCGCCATGGCGCTGTCTGTCTACTGGTGAGAATCCTTCGATGAAAAAGTTCTTCAGCCTGATCGCCACCCTGCTGGTGCTGGTCGCTGCCGTGGCCATCGGTCGCCAGTTGTGGCTG  >  NC_002947/1453016‑1453314

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