BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,271,666	C→T	5.4%	intergenic (+50/‑50)	PP_t32 → / → PP_t33	tRNA‑Asp/tRNA‑Asp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,271,666	0	C	T	5.4%	348.4 / 5.0	110	intergenic (+50/‑50)	PP_t32/PP_t33	tRNA‑Asp/tRNA‑Asp
Reads supporting (aligned to +/- strand): ref base C (47/57); new base T (3/3); total (50/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 115 total aligned reads displayed.

GCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCC  >  NC_002947/2271520‑2271810

gCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACAt                                                                                                                                                  <  3:301313/149‑1 (MQ=255)

 ccAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATg                                                                                                                                                 >  4:222316/1‑149 (MQ=255)

  cAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCTGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGACA‑CAATGCACATgg                                                                                                                                                >  6:190681/1‑148 (MQ=255)

   agagagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGcc                                                                                                                                               <  1:277448/149‑1 (MQ=255)

   agagagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGcc                                                                                                                                               <  1:282201/149‑1 (MQ=255)

    gagagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCt                                                                                                                                              <  6:173085/149‑1 (MQ=255)

      gagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTATAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTACCGTCCGCTGCGCCATTTTTCAGTAACAGATGTGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGGCTTg                                                                                                                                            <  5:190681/149‑1 (MQ=255)

      gagGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTg                                                                                                                                            <  2:152306/149‑1 (MQ=255)

       agGAACGGACTGCGCTGCTGTAGTTCAGTCTGTTAGATTAGGGGCCTGTCACGGCGGGGGTCGCGGTTTCGAGTACCGTGCTCTTCGCCATTTTTCAGTAACAGATGGGTGCCTGGCAGCTGTCTGAAGCCA‑CAAGTCACATGCCTTg                                                                                                                                            <  4:176560/148‑1 (MQ=255)

       agGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCATGTGCCTTgt                                                                                                                                           <  6:180873/149‑2 (MQ=255)

        ggAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGGCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTg‑gg                                                                                                                                         <  8:172116/149‑1 (MQ=255)

        ggAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCATGTGCCTTgtg                                                                                                                                          <  7:226047/149‑3 (MQ=255)

        ggAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTGCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTg‑gg                                                                                                                                         <  3:1855/149‑1 (MQ=255)

        ggAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAAAACCGGCCTGTCACGCCGGGGGTCGCGCGATAGAGTACCGTCCGCTGCGACAATTTTCAGTAACAGATGGGTGCCTCGCACCGGTCTCAAGACA‑CAAGGAACATGCCTTg‑g                                                                                                                                          >  5:334199/1‑148 (MQ=21)

          aaCGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCt                                                                                                                                       >  6:217882/1‑149 (MQ=255)

          aaCGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCt                                                                                                                                       <  6:347471/149‑1 (MQ=255)

             gCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGAt                                                                                                                                    <  8:159303/149‑1 (MQ=255)

             gCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGAt                                                                                                                                    <  5:289794/149‑1 (MQ=255)

              cACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCATGTGCCTTgtg                                                                                                                                          <  4:248955/143‑3 (MQ=255)

              cACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCATGTGCCTTgtg                                                                                                                                          >  3:248955/1‑141 (MQ=255)

                 tGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATcc                                                                                                                                  >  1:26055/1‑147 (MQ=255)

                      aGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTc                                                                                                                           <  2:113586/149‑1 (MQ=255)

                         ggTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAAt                                                                                                                        <  2:82544/149‑1 (MQ=255)

                          gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                       <  6:310405/149‑1 (MQ=255)

                          gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                       <  8:162771/149‑1 (MQ=255)

                              tgcAGTCGGTTAGAATCCCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                                   <  8:255744/147‑1 (MQ=255)

                               tCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                                   <  5:217882/148‑1 (MQ=255)

                               tCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgc                                                                                                                   <  4:139250/148‑1 (MQ=255)

                                 aGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGgcgc                                                                                                                 <  5:58913/148‑1 (MQ=255)

                                      gTTATAATACCGGCCTGTCACGGCGGTGGTCTCGGGTTCGAGTCCCGTCCTCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCTGTCTTAAGCCA‑AAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCt                                                                                                           <  6:334199/149‑1 (MQ=255)

                                      gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAAGCACATCCCTTG‑CGCTCATCCCAGTTTCAATAGGTCGCAAGcc                                                                                                            >  3:15356/1‑148 (MQ=37)

                                          gAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                        <  7:273071/148‑1 (MQ=255)

                                           aaTACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAAc                                                                                                       >  5:52509/1‑148 (MQ=255)

                                            aTACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGa                                                                                                     <  6:88784/149‑1 (MQ=255)

                                             ccccGGCGTGACAGGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGa                                                                                                     >  5:373218/3‑148 (MQ=255)

                                               ccGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                   <  5:224178/148‑1 (MQ=255)

                                               ccGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                   >  6:224178/1‑148 (MQ=255)

                                                cGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                        >  5:88151/1‑142 (MQ=255)

                                                cGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGaa                                                                                                        <  6:88151/142‑1 (MQ=255)

                                                   ccTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAg                                                                                              <  6:52509/149‑1 (MQ=255)

                                                     tGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                            >  5:130739/1‑149 (MQ=255)

                                                       tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCg                                                                                            <  5:218833/147‑1 (MQ=255)

                                                         aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACg                                                                                                      >  4:291465/1‑135 (MQ=255)

                                                         aCGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACg                                                                                                      <  3:291465/135‑1 (MQ=255)

                                                                    cgcgGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGt                                                                              >  2:221849/1‑148 (MQ=255)

                                                                             gtgTCCCGTGCGCTGCGCGATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGTGACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATAcc                                                                     <  4:15356/146‑1 (MQ=255)

                                                                             gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                       >  1:5413/1‑98 (MQ=255)

                                                                             gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATc                                                                                                                       <  2:5413/98‑1 (MQ=255)

                                                                             gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                     >  2:265332/1‑100 (MQ=255)

                                                                             gAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCgg                                                                                                                     <  1:265332/100‑1 (MQ=255)

                                                                              aGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                   <  1:35159/117‑1 (MQ=255)

                                                                              aGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATa                                                                                                   >  2:35159/1‑117 (MQ=255)

                                                                                tCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGcc                                                                 >  4:13256/1‑149 (MQ=255)

                                                                                  ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCACCAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCaa                                                                                                               <  1:71547/102‑1 (MQ=255)

                                                                                  ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCaa                                                                                                               >  2:71547/1‑101 (MQ=255)

                                                                                  ccGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCt                                                                <  1:221849/148‑1 (MQ=255)

                                                                                  ccGTCCGCTGCGCCATTTATCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATACCAGTTTCAATCGGGCGCAAGACTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACAGGCCTg                                                               >  3:142148/1‑149 (MQ=255)

                                                                                   cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTa                                                                                         >  7:303670/1‑122 (MQ=255)

                                                                                   cGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTa                                                                                         <  8:303670/122‑1 (MQ=255)

                                                                                    gTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCGGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCg                                                                                >  1:81307/1‑130 (MQ=255)

                                                                                      ccGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTc                                                             <  6:130739/147‑1 (MQ=255)

                                                                                        gCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGc                                                         >  4:58065/1‑149 (MQ=255)

                                                                                         cTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTc                                                                                 >  5:250532/1‑124 (MQ=255)

                                                                                               cATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGAt                                                                                                    >  5:311859/1‑99 (MQ=255)

                                                                                               cATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGAt                                                                                                    <  6:311859/99‑1 (MQ=255)

                                                                                                       aGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGtt                                         >  5:307399/1‑149 (MQ=255)

                                                                                                         tAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCt                                                                <  8:127365/125‑1 (MQ=255)

                                                                                                         tAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTc                                        >  7:156544/1‑148 (MQ=255)

                                                                                                           aCAGATGGGTGCCTGGCACCGGTCTGAAGCCAGCAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTc                                        >  1:232102/1‑147 (MQ=255)

                                                                                                           aCAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAg                                     <  3:58065/149‑1 (MQ=255)

                                                                                                            cAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGaa                                                                         >  8:237847/1‑113 (MQ=255)

                                                                                                            cAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGt                                    >  2:267775/1‑149 (MQ=255)

                                                                                                             aGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTc                                        >  3:115289/1‑144 (MQ=255)

                                                                                                               aTGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTccc                                 <  3:13256/149‑1 (MQ=255)

                                                                                                                  ggTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTc                              >  6:63605/1‑149 (MQ=255)

                                                                                                                   gTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTcc                             >  8:27883/1‑149 (MQ=255)

                                                                                                                          gCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTgcgc                      >  5:252179/1‑149 (MQ=255)

                                                                                                                              cGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAATGATACGCAGCGGTAGTTCAGTCGGt                                                                              <  8:4204/90‑1 (MQ=255)

                                                                                                                                gTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTgcgc                      >  7:75510/1‑143 (MQ=255)

                                                                                                                                   tGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCgc             >  7:25439/1‑148 (MQ=255)

                                                                                                                                   tGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGt             >  8:264516/1‑149 (MQ=255)

                                                                                                                                    gAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg‑g                                                    <  8:29767/109‑1 (MQ=255)

                                                                                                                                    gAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCgg‑g                                                    >  7:29767/1‑109 (MQ=255)

                                                                                                                                    gAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGtt            >  6:320087/1‑149 (MQ=255)

                                                                                                                                        ccA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCaa        <  6:307399/149‑1 (MQ=255)

                                                                                                                                          cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTc          >  6:366571/3‑147 (MQ=255)

                                                                                                                                          cccaAGGCATGTGCCTTGTGGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTACCGTCCGCTGCGCCATCGTCGTTTc          >  3:69842/3‑147 (MQ=14)

                                                                                                                                          a‑caAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGTGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAg       <  1:267775/149‑1 (MQ=255)

                                                                                                                                          a‑caAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg      >  7:252675/1‑149 (MQ=255)

                                                                                                                                          a‑caAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAgg      >  2:342074/1‑149 (MQ=255)

                                                                                                                                          a‑caAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTCTCACGCCGG‑GCGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCTTCTTCGTTTCAAgg      >  7:252008/1‑149 (MQ=255)

                                                                                                                                               ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAg       >  3:288833/1‑144 (MQ=255)

                                                                                                                                               ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAg       <  4:288833/144‑1 (MQ=255)

                                                                                                                                               ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  <  8:156544/149‑1 (MQ=255)

                                                                                                                                               ggCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  >  5:132933/1‑149 (MQ=255)

                                                                                                                                                  acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGTCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  >  7:258617/1‑146 (MQ=255)

                                                                                                                                                  acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGtt                                         <  3:296023/107‑1 (MQ=255)

                                                                                                                                                  acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGtt                                         >  4:296023/1‑107 (MQ=255)

                                                                                                                                                  acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTccc  >  3:162245/1‑146 (MQ=255)

                                                                                                                                                  acaTGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTAGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTAAAGGTccc  >  4:275270/1‑146 (MQ=255)

GCCAGAGAGGAACGCACTGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATTTTTCAGTAACAGATGGGTGCCTGGCACCGGTCTGAAGCCA‑CAAGGCACATGCCTTG‑GGCTGATCCCAGTTTCAATCGGGCGCAAGCCTGAACGATACGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGG‑GGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATCTTCGTTTCAAGGTCCC  >  NC_002947/2271520‑2271810

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