BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,256,887	A→T	7.7%	noncoding (434/451 nt)	PP_mr14 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,256,887	0	A	T	7.7%	751.0 / 36.0	230	noncoding (434/451 nt)	PP_mr14	cspA regulatory region
Reads supporting (aligned to +/- strand): ref base A (105/107); new base T (8/10); total (113/117)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 232 total aligned reads displayed.

CATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTT  >  NC_002947/1256738‑1257036

   tACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGaa                                                                                                                                                    <  1:1325301/150‑1 (MQ=255)

   tACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGaa                                                                                                                                                    <  1:657070/150‑1 (MQ=255)

      cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGc                                                                                                                                                 >  2:3720262/1‑150 (MQ=255)

       aTCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAAACCCGCCAGAAGGc                                                                                                                                                 >  1:2529312/1‑150 (MQ=255)

       aTCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCg                                                                                                                                                <  1:410124/150‑1 (MQ=255)

       aTCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCg                                                                                                                                                <  1:1977410/150‑1 (MQ=255)

          aGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGt                                                                                                                                             <  2:1172716/150‑1 (MQ=255)

           gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGtt                                                                                                                                            <  1:1157298/150‑1 (MQ=255)

           gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGtt                                                                                                                                            <  1:3914195/150‑1 (MQ=255)

           gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGtt                                                                                                                                            >  1:2131558/1‑150 (MQ=255)

            cGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCTTCTGGCGGGttt                                                                                                                                           <  2:476650/150‑1 (MQ=255)

            cGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGttt                                                                                                                                           >  2:719449/1‑150 (MQ=255)

             gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGtttt                                                                                                                                          <  2:77180/150‑1 (MQ=255)

             gAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGtttt                                                                                                                                          <  1:2573273/150‑1 (MQ=255)

             aaaGACTGACCAGCCAAACCGAGTTTGAACGCACACAGCGACTGAGCCGTACTTGAATCCTGTGCAGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATAACGAGCAGGCAT‑AAAAAAACACGCCAGAAGGCGGGattt                                                                                                                                          >  2:1244978/2‑150 (MQ=37)

              aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCTTCTGGCGGGttttt                                                                                                                                         <  1:3676547/150‑1 (MQ=255)

              aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGttttt                                                                                                                                         >  1:4456759/1‑150 (MQ=255)

              aaGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGGCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCTTCTGGCGGGttttt                                                                                                                                         <  1:2432375/150‑1 (MQ=255)

                gACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCTTCTGGCGGGttttttt                                                                                                                                       <  1:1496607/150‑1 (MQ=255)

                gACTGACCTGCCAGCCCGCGGTTGCCCGCCCACAGCGACTGCGCCGTCCTGGACGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGACGAGGCCAGCGCCTATTTCGAGCAGGCAT‑ACCCAAACCCGCCAGCAGGCGGGttttttt                                                                                                                                       <  2:184454/150‑1 (MQ=37)

                  cTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                     <  1:1294346/150‑3 (MQ=255)

                  cTGACCTGCCATACCGAGTTGGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                     <  1:121515/150‑3 (MQ=38)

                   tGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCGGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaa                                                                                                                                    >  1:2470142/1‑150 (MQ=255)

                   tGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaa                                                                                                                                    >  1:157737/1‑150 (MQ=255)

                   tGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaa                                                                                                                                    <  2:1483549/150‑1 (MQ=255)

                    gACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                   <  2:1740615/150‑1 (MQ=255)

                       cTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCt                                                                                                                                >  1:3136688/1‑150 (MQ=255)

                        tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                               >  1:1929294/1‑150 (MQ=255)

                          ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAgg                                                                                                                             >  2:458076/1‑150 (MQ=255)

                             tACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAg                                                                                                                          >  2:1080712/1‑150 (MQ=255)

                              aCCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                         >  2:1957123/1‑150 (MQ=255)

                              aCCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                         >  1:904124/1‑150 (MQ=255)

                              aCCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCCCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                         <  1:2677454/150‑1 (MQ=255)

                               ccGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCt                                                                                                                        >  2:2589187/1‑150 (MQ=255)

                                     ttccccGCACACAGCGACTGCGCCGTACATGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGACCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCg                                                                                                                  <  1:3663676/147‑1 (MQ=255)

                                       gCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAc                                                                                                                <  2:2176010/150‑1 (MQ=255)

                                           gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAg                                                                                                            >  2:2689272/1‑150 (MQ=255)

                                            cacacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAgg                                                                                                           <  1:163523/150‑1 (MQ=255)

                                             acacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGc                                                                                                          >  2:3755935/1‑150 (MQ=255)

                                              cacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCg                                                                                                         <  1:2308346/150‑1 (MQ=255)

                                               acaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGa                                                                                                        >  1:731371/1‑150 (MQ=255)

                                                 aGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAACAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATc                                                                                                      <  2:2662406/150‑1 (MQ=255)

                                                  gCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTGTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCt                                                                                                     >  2:2635079/1‑150 (MQ=255)

                                                    gACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGa                                                                                                   >  2:134248/1‑150 (MQ=255)

                                                      cTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                 <  2:4456759/150‑1 (MQ=255)

                                                          gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtct                                                                                             <  2:4372141/150‑1 (MQ=255)

                                                          gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtct                                                                                             <  1:1815080/150‑1 (MQ=255)

                                                          gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtct                                                                                             >  1:839422/1‑150 (MQ=255)

                                                           ccGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtctt                                                                                            <  2:2718411/150‑1 (MQ=255)

                                                                 ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCt                                                                                      >  1:1389730/1‑150 (MQ=255)

                                                                   gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AATAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAGTTAGGCGATCTGAACTTGTTCAGCGTGc                                                                                    >  1:2377374/1‑150 (MQ=255)

                                                                    aaGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCa                                                                                   <  2:1727583/150‑1 (MQ=255)

                                                                      gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATg                                                                                 >  1:2989549/1‑150 (MQ=255)

                                                                      gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATg                                                                                 >  1:3911605/1‑150 (MQ=255)

                                                                         tgtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCt                                                                              <  2:831153/150‑1 (MQ=255)

                                                                          gtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTAGTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCtt                                                                             >  1:3505445/1‑150 (MQ=255)

                                                                                tCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGc                                                                       >  2:2055109/1‑150 (MQ=255)

                                                                                 cccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGcc                                                                      <  1:3446352/150‑1 (MQ=255)

                                                                                  ccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAATACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcg                                                                     >  1:2831480/1‑150 (MQ=255)

                                                                                  ccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcg                                                                     >  1:2185830/1‑150 (MQ=255)

                                                                                  ccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcg                                                                     >  2:2863746/1‑150 (MQ=255)

                                                                                    tgtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgcg                                                                   <  2:4121981/150‑1 (MQ=255)

                                                                                    tgtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgcg                                                                   <  2:3230771/150‑1 (MQ=255)

                                                                                     gtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                  <  1:336755/150‑1 (MQ=255)

                                                                                     gtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                  <  2:535543/150‑1 (MQ=255)

                                                                                        ttcgCGGGCTCGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTa                                                               >  1:4401796/4‑150 (MQ=255)

                                                                                        ttcgCGGGCTCGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTa                                                               >  2:2507737/4‑150 (MQ=255)

                                                                                            cGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcga                                                           <  1:1771309/150‑1 (MQ=255)

                                                                                               gCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgacga                                                        >  1:3058173/1‑150 (MQ=255)

                                                                                                cGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGaa                                                       <  1:2211293/150‑1 (MQ=255)

                                                                                                    cccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTa                                                   >  1:4306109/1‑150 (MQ=255)

                                                                                                         gAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGc                                              <  1:165921/150‑1 (MQ=255)

                                                                                                           agagGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCt                                            <  2:4033690/150‑1 (MQ=255)

                                                                                                              ggCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGc                                         <  1:4352966/150‑1 (MQ=255)

                                                                                                                ccAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCct                                       >  1:4483344/1‑150 (MQ=255)

                                                                                                                ccAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCct                                       >  1:1695971/1‑150 (MQ=255)

                                                                                                                     gcCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAACGTAACAGCCTGGCCTTCttt                                  >  2:2944794/1‑150 (MQ=255)

                                                                                                                       cTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCa                                <  1:1080712/150‑1 (MQ=255)

                                                                                                                        tATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAg                               >  1:1964032/1‑150 (MQ=255)

                                                                                                                        tATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAg                               >  1:543596/1‑150 (MQ=255)

                                                                                                                         aTTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAgg                              >  2:2158708/1‑150 (MQ=255)

                                                                                                                              gAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGtttt                         >  2:3029935/1‑150 (MQ=255)

                                                                                                                              gAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGtttt                         >  2:4200402/1‑150 (MQ=255)

                                                                                                                                gCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGa                       <  1:2067052/150‑1 (MQ=255)

                                                                                                                                 cAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGaa                      <  2:1827802/150‑1 (MQ=255)

                                                                                                                                   ggCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGc                    <  2:1084238/150‑1 (MQ=255)

                                                                                                                                    gCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                   >  1:3806324/1‑150 (MQ=255)

                                                                                                                                    gCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                   >  1:2097901/1‑150 (MQ=255)

                                                                                                                                       t‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTc                <  2:3051791/150‑1 (MQ=255)

                                                                                                                                         aaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCa               >  2:3941024/1‑150 (MQ=255)

                                                                                                                                         aaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGACa               >  2:2928569/1‑150 (MQ=255)

                                                                                                                                         aaaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCa               <  1:4473656/150‑1 (MQ=255)

                                                                                                                                         aaaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCa               <  2:2529312/150‑1 (MQ=255)

                                                                                                                                          aaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAg              >  2:1059960/1‑150 (MQ=255)

                                                                                                                                           aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgac                                                          <  2:117023/105‑1 (MQ=255)

                                                                                                                                           aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGc             >  2:1006354/1‑150 (MQ=255)

                                                                                                                                            aaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCt            <  1:779262/150‑1 (MQ=255)

                                                                                                                                                cccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCGTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCGTGGa        >  2:1783158/1‑150 (MQ=255)

                                                                                                                                                 ccGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGAt       <  2:3523047/150‑1 (MQ=255)

                                                                                                                                                      aGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCtt  >  2:1131793/1‑150 (MQ=255)

CATTACCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCAT‑AAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTT  >  NC_002947/1256738‑1257036

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