BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,149,269	A→T	6.3%	A1072A (GCT→GCA)	rne ←	ribonuclease E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,149,269	0	A	T	6.3%	368.8 / 10.2	111	A1072A (GCT→GCA)	rne	ribonuclease E
Reads supporting (aligned to +/- strand): ref base A (58/46); new base T (4/3); total (62/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

CGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAG  >  NC_002947/2149120‑2149417

  ggTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagc                                                                                                                                                    <  1:921481/150‑1 (MQ=255)

   gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                                   >  1:285614/1‑150 (MQ=255)

   gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                                   <  1:5169473/150‑1 (MQ=255)

   gTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTgcagca                                                                                                                                                   >  2:351060/1‑150 (MQ=255)

     ttACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCagcagc                                                                                                                                                 <  2:4759030/150‑1 (MQ=255)

     ttACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCagcagc                                                                                                                                                 <  2:1443277/150‑1 (MQ=255)

     ttACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCagcagc                                                                                                                                                 <  1:460140/150‑1 (MQ=255)

        cccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTc                                                                                                                                              <  2:3896995/150‑1 (MQ=255)

         ccGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGGGGCTTCCAGGGCGGGTGCCTGGGTCGCTGCagcagcgtgg                                                                                                                                             >  2:4712388/1‑146 (MQ=255)

           gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTccca                                                                                                                                           <  1:4438921/150‑3 (MQ=255)

           gcgAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCtg                                                                                                                                           >  1:4282143/1‑150 (MQ=255)

               agagGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGcc                                                                                                                                       <  1:1089251/150‑1 (MQ=255)

                      cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                                <  1:5134006/150‑1 (MQ=255)

                      cGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTgg                                                                                                                                <  1:5134397/150‑1 (MQ=255)

                       gTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgc                                                                                                                               <  1:1833312/150‑1 (MQ=255)

                       gTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgc                                                                                                                               <  2:1653752/150‑1 (MQ=255)

                       gTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgc                                                                                                                               <  2:1953720/150‑1 (MQ=255)

                           ggCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgcagca                                                                                                                           >  1:3190175/1‑150 (MQ=255)

                            gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGcagcag                                                                                                                          >  2:1092256/1‑150 (MQ=255)

                            gCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGcagcag                                                                                                                          >  1:3678302/1‑150 (MQ=255)

                              aaaCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCAGTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTGTGGCAGCAGcc                                                                                                                        >  1:5404583/1‑150 (MQ=255)

                                aCCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTc                                                                                                                      <  1:2392978/150‑1 (MQ=255)

                                 ccAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCg                                                                                                                     >  1:448295/1‑150 (MQ=255)

                                 ccAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCg                                                                                                                     >  2:1101974/1‑150 (MQ=255)

                                   aaaaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGc                                                                                                                   >  1:3199237/1‑150 (MQ=255)

                                   aaaaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGc                                                                                                                   <  2:4316291/150‑1 (MQ=255)

                                    aaaaaaGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGcc                                                                                                                  >  2:310643/1‑150 (MQ=255)

                                          gCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGc                                                                                                            >  2:1236700/1‑150 (MQ=255)

                                          gCCCCGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGc                                                                                                            >  1:1193157/1‑150 (MQ=255)

                                            cccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                          >  1:1153974/1‑150 (MQ=255)

                                            cccGCCAGTGCAAACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCtt                                                                                                          >  2:786464/1‑150 (MQ=255)

                                             ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                                                                                                              >  2:2069061/1‑145 (MQ=255)

                                             ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTccc                                                                                                              <  1:2069061/145‑1 (MQ=255)

                                             ccGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTg                                                                                                         >  2:2434729/1‑150 (MQ=255)

                                              cGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgc                                                                                                        >  1:4805540/1‑150 (MQ=255)

                                                ccAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcgg                                                                                                      <  2:1305701/150‑1 (MQ=255)

                                                ccAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcgg                                                                                                      >  1:856917/1‑150 (MQ=255)

                                                ccAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTgcgg                                                                                                      >  2:3048299/1‑150 (MQ=255)

                                                      tGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGc                                                                                                <  2:2045119/150‑1 (MQ=255)

                                                        cACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCTGCAGCTTccc                                                                                                                                            >  2:5231904/1‑103 (MQ=255)

                                                         aCTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACt                                                                                             >  1:205827/1‑150 (MQ=255)

                                                               ggggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGt                                                                                       >  1:4202097/1‑150 (MQ=255)

                                                                 ggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                                     <  1:4547489/150‑1 (MQ=255)

                                                                 ggCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTgg                                                                                     <  2:3312843/150‑1 (MQ=255)

                                                                    tttttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTc                                                                                  >  1:4585792/1‑150 (MQ=255)

                                                                      tttCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGt                                                                                >  2:3983975/1‑150 (MQ=255)

                                                                         cattcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCg                                                                             <  2:1706103/150‑1 (MQ=255)

                                                                          attcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCgg                                                                            <  2:205827/150‑1 (MQ=255)

                                                                          attcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCgg                                                                            >  2:239113/1‑150 (MQ=255)

                                                                           ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                           <  2:842638/150‑1 (MQ=255)

                                                                           ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                           >  1:2429241/1‑150 (MQ=255)

                                                                           ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGt                                                                           >  1:2429359/1‑150 (MQ=255)

                                                                           ttcaGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTACACGTGGGTCGTTCGGt                                                                           <  2:3585591/150‑1 (MQ=255)

                                                                                  ttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgca                                                                                                                              <  2:3986818/92‑1 (MQ=255)

                                                                                  ttGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGgca                                                                                                                              >  1:3986818/1‑92 (MQ=255)

                                                                                   tGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGcc                                                                   <  1:3909809/150‑1 (MQ=255)

                                                                                   tGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGcc                                                                   <  1:3568224/150‑1 (MQ=255)

                                                                                   tGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGcc                                                                   <  1:3181630/150‑1 (MQ=255)

                                                                                   tGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGcc                                                                   <  1:3181643/150‑1 (MQ=255)

                                                                                       tATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTg                                                               >  1:2625069/1‑150 (MQ=255)

                                                                                           aGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCa                                                           <  2:3938731/150‑1 (MQ=255)

                                                                                             ccGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGc                                                         <  2:2625069/150‑1 (MQ=255)

                                                                                             ccGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGc                                                         >  2:5402655/1‑150 (MQ=255)

                                                                                                 gatgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                                     <  2:661749/150‑1 (MQ=255)

                                                                                                 gatgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGa                                                     <  2:661711/150‑1 (MQ=255)

                                                                                                  atgaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGAc                                                    >  1:53265/1‑150 (MQ=255)

                                                                                                    gaGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAg                                                  <  1:3434442/150‑1 (MQ=255)

                                                                                                      ggTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTg                                                <  1:3551429/150‑1 (MQ=255)

                                                                                                       gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTgg                                               <  2:1779205/150‑1 (MQ=255)

                                                                                                       gTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTgg                                               <  2:3199237/150‑1 (MQ=255)

                                                                                                          tgtgCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAg                                            >  1:1054287/1‑150 (MQ=255)

                                                                                                               tCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCg                                       >  1:3827981/1‑150 (MQ=255)

                                                                                                                cATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGa                                      <  1:2885889/150‑1 (MQ=255)

                                                                                                                cATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGa                                      <  1:4865657/150‑1 (MQ=255)

                                                                                                                  tCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGAcc                                    >  1:1746608/1‑150 (MQ=255)

                                                                                                                    ggcggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCAc                                  <  2:3190175/150‑1 (MQ=255)

                                                                                                                     gcggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACt                                 >  2:5068500/1‑150 (MQ=255)

                                                                                                                      cggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtg                                >  2:4176281/1‑150 (MQ=255)

                                                                                                                      cggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtg                                >  2:4176280/1‑150 (MQ=255)

                                                                                                                       ggcTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgc                               >  2:4693733/1‑150 (MQ=255)

                                                                                                                        gcTTCCAGGGCTGGTTCCTGGGAAGCGGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGTCAGGTTCGACCACtgca                              <  1:290944/150‑1 (MQ=255)

                                                                                                                              aGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACtgcaggtgca                        >  1:2067692/1‑150 (MQ=255)

                                                                                                                                  cTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGtt                    >  1:2834423/1‑150 (MQ=255)

                                                                                                                                  cTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGtt                    >  1:4276332/1‑150 (MQ=255)

                                                                                                                                   tGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTc                   <  1:3857316/150‑1 (MQ=255)

                                                                                                                                   tGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTc                   <  1:1907859/150‑1 (MQ=255)

                                                                                                                                   tGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTc                   >  2:280347/1‑150 (MQ=255)

                                                                                                                                     gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGa                 <  1:351060/150‑1 (MQ=255)

                                                                                                                                     gTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGa                 >  2:3349751/1‑150 (MQ=255)

                                                                                                                                         cTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccac             >  1:4739402/1‑150 (MQ=255)

                                                                                                                                          tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc            >  1:1643153/1‑150 (MQ=255)

                                                                                                                                          tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc            >  1:1239064/1‑150 (MQ=255)

                                                                                                                                          tGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTAGACCACTGCAGGTGCAGGTTCGaccacc            >  1:1643240/1‑150 (MQ=255)

                                                                                                                                          tGGGAAGCTGCAGCAGCTTCCGGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGaccacc            <  2:4102084/150‑1 (MQ=255)

                                                                                                                                           aggAAGCTGCTGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCg           >  2:586798/2‑150 (MQ=255)

                                                                                                                                              aaGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAg        <  2:3657660/150‑1 (MQ=255)

                                                                                                                                                 cTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCCGGGTCGACCACCGCAGGTg     >  1:5007290/1‑150 (MQ=255)

                                                                                                                                                  tGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGc    >  2:4781286/1‑150 (MQ=255)

                                                                                                                                                    cagcagCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAg  >  1:726103/1‑150 (MQ=255)

                                                                                                                                                    cagcagCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAg  >  1:5201932/1‑150 (MQ=255)

CGGGTTTACCCGCGAAGAGGCCCGTCAGGCAAACCAAAAAAAGCCCCGCCAGTTTGCACTGGCGGGGCTTTTTCATTCAGCTTTGGCTATCAGCCGTGATGAGGTTTGTGCTCATCGGCGGCTTCCAGGGCTGGTTCCTGGGAAGCTGCAGCAGCTTCCTGTGCCGCTTTGGCAGCAGCCTCGGCCTCCCGCTTGCGGCGACGCACTTCACGTGGGTCGTTCGGTGCACGGCCGTTGGCCAGCATGACAGTGGCAGGTTCGACCACTGCAGGTGCAGGTTCGACCACCGCAGGTGCAG  >  NC_002947/2149120‑2149417

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