BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	733,648	A→T	5.9%	intergenic (+19/‑28)	PP_t09 → / → PP_t10	tRNA‑Lys/tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	733,648	0	A	T	5.9%	530.8 / 14.3	177	intergenic (+19/‑28)	PP_t09/PP_t10	tRNA‑Lys/tRNA‑Pro
Reads supporting (aligned to +/- strand): ref base A (91/68); new base T (8/3); total (100/77)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 182 total aligned reads displayed.

CTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733499‑733797

cTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCa                                                                                                                                                       >  2:6071509/1‑150 (MQ=255)

  gCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCaaa                                                                                                                                                     <  1:4017762/150‑1 (MQ=255)

    tttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAaga                                                                                                                                                   <  2:3122198/150‑1 (MQ=255)

     ttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAagag                                                                                                                                                  <  1:4198094/150‑1 (MQ=255)

     ttGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAagag                                                                                                                                                  <  2:126455/150‑1 (MQ=255)

         gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGc                                                                                                                                              >  1:1610410/1‑150 (MQ=255)

         gAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATAACAAAGAGATGc                                                                                                                                              >  1:954928/1‑150 (MQ=255)

          aGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCTTAAGAACAGTTGGAAGAGCAGATGGCAATTAACCAATACGTCGAAGGATCGAAAACTACACTACCCACAATGTAGAAAAGGGCATATCAAAGAGATgac                                                                                                                                             >  1:1059974/1‑148 (MQ=17)

            cAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGACATGCCCt                                                                                                                                           >  1:3070735/1‑150 (MQ=255)

                 tGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCtttttt                                                                                                                                      <  1:1808913/150‑1 (MQ=255)

                  gTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttttttt                                                                                                                                     <  2:1347314/150‑1 (MQ=255)

                  gTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCttttttt                                                                                                                                     <  2:3015155/150‑1 (MQ=255)

                   tCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAATGGCAAGCACAAAGAGAAGACCtaggtttt                                                                                                                                    >  2:5566939/1‑143 (MQ=14)

                     aGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGa                                                                                                                                  >  2:840845/1‑150 (MQ=255)

                      ggAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGaa                                                                                                                                 <  2:6215013/150‑1 (MQ=255)

                       gAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTGTGaaa                                                                                                                                >  1:345500/1‑150 (MQ=255)

                          ttCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGTAacac                                                                                                                             >  1:4395332/1‑150 (MQ=255)

                             gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                          >  2:6120845/1‑150 (MQ=255)

                             gAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCg                                                                                                                          <  1:5896265/150‑1 (MQ=255)

                              aaTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCgg                                                                                                                         <  1:1245429/150‑1 (MQ=255)

                               aTCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        <  2:157126/150‑1 (MQ=255)

                                 ccccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgt                                                                                                                      <  2:5497946/150‑1 (MQ=255)

                                  cccTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtg                                                                                                                     >  1:3837432/1‑150 (MQ=255)

                                    cTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                   <  1:767579/150‑1 (MQ=255)

                                    cTCACCCCGACCATTTTCCGGGTCGTTAGATCAGTTGGTAGAGCAGTTTGCTTTTAACCCATTGGTCGTAGGTTAGAATCCTACTCTACCCACCATTTATAAAAGGGCATCTGAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTa                                                                                                                   <  2:2395491/150‑1 (MQ=255)

                                         cccGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCa                                                                                                              <  2:6618006/150‑1 (MQ=255)

                                              ccATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCg                                                                                                         >  1:3866038/1‑150 (MQ=255)

                                              ccATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCg                                                                                                         <  2:3253654/150‑1 (MQ=255)

                                               cATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCggg                                                                                                                        >  1:5684899/1‑134 (MQ=255)

                                                aTTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGt                                                                                                       >  2:2222109/1‑150 (MQ=255)

                                                   ttCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAACAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgc                                                                                                    <  2:2745590/150‑1 (MQ=255)

                                                       gggTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACCCGACCCACCCTGTAGAAAAGGGCATCTCCACGCGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgcgc                                                                                                <  2:1934725/150‑1 (MQ=21)

                                                       gggTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACAAGGAACAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAgcgcgc                                                                                                >  1:3604152/1‑150 (MQ=255)

                                                             ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCtt                                                                                          >  2:4501466/1‑150 (MQ=255)

                                                             ttAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTGTTGAAACACTAGGGGTGTAGAGCAGTCCGGTCGCGCACCTGCtt                                                                                          >  2:4693032/1‑150 (MQ=21)

                                                              tAGCTCAGTTGGTAGAGCAGTTGGCTTTTTACCTCTGGTTAGTAGGTTCGAATCCTTCACGACCCACCCTGTAGAAAAGGGCCTCTCAAAGAGATGCACTTTTTTTTGAAACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCttt                                                                                         <  1:6050168/150‑1 (MQ=25)

                                                               aGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                        <  1:6071509/150‑1 (MQ=255)

                                                               aGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCCCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTg                                                                                        <  2:2550143/150‑1 (MQ=255)

                                                                gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCACGTAGAACAGGGCATCTCCAAGCGCTGCAGTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  2:2723170/150‑1 (MQ=25)

                                                                gCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTAGTAGTTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGATATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTgg                                                                                       <  1:547506/150‑1 (MQ=255)

                                                                  tCAGTTGGTAGAGCAGTTGGCTTTTAACCACTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGACACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGa                                                                                     <  1:1011782/150‑1 (MQ=255)

                                                                  tCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAAAACTCGGGGTGTAGCGCAGTCCGGTCGCGCGCCTGCTTTGGGa                                                                                     >  2:3939042/1‑150 (MQ=255)

                                                                   cAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATATCAAAGAGATGCCCTTTTTGTTGAAACACACGGGGTGTAGCGCAGTCCGGGAGCGCGCCTGCGGTGGGAg                                                                                    >  2:4268266/1‑150 (MQ=25)

                                                                      tgggTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAg                                                                                 <  1:5312364/148‑1 (MQ=255)

                                                                             agCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCa                                                                          <  2:4113539/150‑1 (MQ=255)

                                                                             agCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCGTTGGGAGCAGGATGTCa                                                                          >  2:6489933/1‑150 (MQ=255)

                                                                               cAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCtt                                                                                          <  1:3686773/132‑1 (MQ=255)

                                                                               cAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCtt                                                                                          >  2:3686773/1‑132 (MQ=255)

                                                                                aGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCACTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGa                                                                       <  2:6371762/150‑1 (MQ=255)

                                                                                  ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTATCGCGCCTGCGTTGGGAGCAGGATGTCAGGAGt                                                                     >  2:457822/1‑150 (MQ=255)

                                                                                  ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGt                                                                     >  2:5263619/1‑150 (MQ=255)

                                                                                  ttGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTACAAAAGGCCAAAACAAAGAGATGCCCTTTTTATTGAAACACTCGCGTTGACGCCCACTCCGGAAGCGCGACAGCTTTGGGCGCAGGATGTCAGGAGt                                                                     >  2:1791503/1‑150 (MQ=21)

                                                                                  ttGGCTTTTAACCAATTGGTCGTAGGTACGAATCCTACACGACACACCATCTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTATTTGAAACAATCGGGGTGTAGCGCGGTACGGTATCGCGCATGCTTTGGGAGCAGGATGACAGGAGt                                                                     >  2:2894154/1‑150 (MQ=37)

                                                                                   tGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTACAAATGGCCATAACAAAGAGATGCCCTTTTTATAGAAACACACGGGTAGAACCGCACTCCGGAATCGCGCCTGATGTGGGAGCAGGATGTCAGGAgat                                                                    >  2:4892445/1‑148 (MQ=18)

                                                                                       ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGaa                                                                >  2:354724/1‑150 (MQ=255)

                                                                                       ttttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCC‑TTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAt                                                               >  2:1538877/1‑150 (MQ=255)

                                                                                        tttAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGTGGTGTAGAGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAAt                                                               >  2:4882989/1‑150 (MQ=255)

                                                                                          tAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATcc                                                             >  2:1019496/1‑150 (MQ=255)

                                                                                            aCCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATcccc                                                           >  1:6329127/1‑150 (MQ=255)

                                                                                             ccAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCt                                                          >  1:694816/1‑150 (MQ=255)

                                                                                              cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                         >  1:2236492/1‑150 (MQ=255)

                                                                                              cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                         >  1:2933108/1‑150 (MQ=255)

                                                                                              cAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTc                                                         >  2:3274936/1‑150 (MQ=255)

                                                                                               aaTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCa                                                        >  2:6936763/1‑150 (MQ=255)

                                                                                               aaTTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCa                                                        >  1:721741/1‑150 (MQ=255)

                                                                                                    gTCGTAGGTTCGAATCCTACACGACCCACCATGAACAAAAGGGAAAAACAAAGAGATGCCCTTTTTATTTAAACACTAGAGGTGAACCCAAGTCCGGAATCGCGCCTGCATAGAGAGAAGTATGTCAGGAGTACGAAACCCCTCAccaca                                                   >  2:6766357/1‑147 (MQ=16)

                                                                                                       gTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCTGGATGTCAGGAGTTCGAATCCCCTCACCCCGAcc                                                >  2:814555/1‑150 (MQ=255)

                                                                                                        tAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGgtgt                                                                                                                    <  2:2321260/81‑1 (MQ=255)

                                                                                                              cGAATCCTACACGAACCACCATGTAGAAATGGGCATCTCAAAGAGATTCCCTTTTTTTTGAAACTCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTg                                         <  2:2434367/150‑1 (MQ=255)

                                                                                                                  tCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTc                                     <  1:64540/150‑1 (MQ=255)

                                                                                                                    cTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGt                                   >  1:324474/1‑150 (MQ=255)

                                                                                                                     tACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGtt                                  >  1:4560108/1‑150 (MQ=255)

                                                                                                                        acGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATc                                                              >  2:1064882/1‑119 (MQ=255)

                                                                                                                          gACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTc                             >  2:5230785/1‑150 (MQ=255)

                                                                                                                           aCCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCa                            <  1:6120845/150‑1 (MQ=255)

                                                                                                                            cccACCATGTAGAAAAGGGCATCTCAAATAGATGACCTTTTTTTTGAAACACTAGGGGTGTATCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGATTTAGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAg                           >  1:2305493/1‑150 (MQ=21)

                                                                                                                             ccaccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGGGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGt                          >  2:7101148/1‑150 (MQ=255)

                                                                                                                              caccaTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGtt                         >  1:6486988/1‑150 (MQ=255)

                                                                                                                                 ttcaaaaaaaaGGGCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCGCAGTTGGt                      >  1:5008584/8‑150 (MQ=18)

                                                                                                                                  aTGTAGACAAGGGCATCTCAAAGAGTTGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGATGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCATGAGTTCGAATCCCCTCACCCCTAACATTTTTGGGTCTTTAGCTCAGTTGGTa                     <  1:6257119/150‑1 (MQ=21)

                                                                                                                                  aTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTa                     >  2:6805800/1‑150 (MQ=255)

                                                                                                                                  aTGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTa                     <  2:5846126/150‑1 (MQ=255)

                                                                                                                                   tGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTag                    >  1:2846580/1‑150 (MQ=255)

                                                                                                                                        aaaaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAg               >  1:2114334/1‑150 (MQ=255)

                                                                                                                                          aaGGGCATCTCAACGAGATGCCCTTTTTTTTGACACCCTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt             <  2:1610410/150‑1 (MQ=255)

                                                                                                                                          aaGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGtt             >  1:6655749/1‑150 (MQ=255)

                                                                                                                                           aGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTTGCTCAGTTGGTAGAGCAGTTg            >  1:1429167/1‑150 (MQ=255)

                                                                                                                                            gggCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  2:4370113/1‑150 (MQ=255)

                                                                                                                                            gggCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  1:3789789/1‑150 (MQ=255)

                                                                                                                                            gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTCGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  2:1367956/1‑150 (MQ=255)

                                                                                                                                            gggCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTgg           >  2:882803/1‑150 (MQ=255)

                                                                                                                                             ggCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGc          >  2:2187971/1‑150 (MQ=255)

                                                                                                                                             ggCATCTCTTTGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGc          >  1:5258774/1‑150 (MQ=255)

                                                                                                                                                 tctcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTCGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt      >  1:2199796/1‑150 (MQ=255)

                                                                                                                                                 tctcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCtttt      <  2:694816/150‑1 (MQ=255)

                                                                                                                                                   tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa    >  1:1779928/1‑150 (MQ=255)

                                                                                                                                                   tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa    >  2:4537285/1‑150 (MQ=255)

                                                                                                                                                   tcAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTaa    >  1:5431513/1‑150 (MQ=255)

                                                                                                                                                    cAAAGAGATGCCCTTTTTTTTGAAACACACGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAc   >  1:293985/1‑150 (MQ=255)

                                                                                                                                                     aaaGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAAcc  >  2:4412687/1‑150 (MQ=255)

CTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTCCGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACCAATTGGTCGTAGGTTCGAATCCTACACGACCCACCATGTAGAAAAGGGCATCTCAAAGAGATGCCCTTTTTTTTGAAACACTCGGGGTGTAGCGCAGTCCGGTAGCGCGCCTGCTTTGGGAGCAGGATGTCAGGAGTTCGAATCCCCTCACCCCGACCATTTTTGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGGCTTTTAACC  >  NC_002947/733499‑733797

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