BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,372,779	A→T	5.6%	K102* (AAG→TAG) ‡	PP_4720 →	RNA‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,372,779	0	A	T	5.6%	1255.3 / 40.0	374	K102* (AAG→TAG) ‡	PP_4720	RNA‑binding protein
Reads supporting (aligned to +/- strand): ref base A (160/193); new base T (14/7); total (174/200)
Fisher's exact test for biased strand distribution p-value = 7.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 382 total aligned reads displayed.

CTCGGAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATA  >  NC_002947/5372631‑5372928

cTCGGAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACaa                                                                                                                                                      <  1:2556828/150‑1 (MQ=255)

    gAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTATCCAACATCCACCGTTACAAGTaa                                                                                                                                                  <  1:2407466/150‑1 (MQ=255)

     aaGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAAc                                                                                                                                                 <  2:8375640/150‑1 (MQ=255)

       gACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCTGACCTTCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACgg                                                                                                                                               <  2:4848487/150‑1 (MQ=255)

         ccGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACTTGTAACGGtg                                                                                                                                             <  2:1583033/150‑2 (MQ=255)

         ccGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTa                                                                                                                                             <  1:3279734/150‑1 (MQ=255)

          cgcgAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAg                                                                                                                                            <  2:7987579/150‑1 (MQ=255)

          cgcgAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAg                                                                                                                                            >  1:1459995/1‑150 (MQ=255)

              aagaaCGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTcc                                                                                                                                        <  1:5109603/150‑1 (MQ=255)

                gaaCGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGa                                                                                                                                      <  2:6888829/150‑1 (MQ=255)

                    gcgcCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAg                                                                                                                                  >  1:6776946/1‑150 (MQ=255)

                      gcCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTg                                                                                                                                <  1:1586879/150‑1 (MQ=255)

                        cGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGgc                                                                                                                              >  2:1960268/1‑150 (MQ=255)

                              cATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTg                                                                                                                        >  1:2747376/1‑150 (MQ=255)

                              cATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTg                                                                                                                        >  2:4594040/1‑150 (MQ=255)

                                tCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGaa                                                                                                                      >  2:743782/1‑150 (MQ=255)

                                 cGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAAc                                                                                                                     <  1:6142836/150‑1 (MQ=255)

                                    cGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAAcgcg                                                                                                                  >  2:56030/1‑150 (MQ=255)

                                          gtgCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGa                                                                                                            <  1:5479127/150‑1 (MQ=255)

                                             cAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGTTCACCGCGccactgatc                                                                                                         <  1:3036510/150‑8 (MQ=255)

                                             cAAGGCCGGCCGGGCCGAGATGGTGCAGACCATCGGCAAGAAAGCCCTGAGCTACCGCAAGAACCCACAGCCGACCACGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCTATCAGTGGCGCGGTGAACGCGCCCTGAcct                                                                                                         <  1:6979587/150‑2 (MQ=255)

                                                ggccggccGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTa                                                                                                      >  1:2880394/1‑150 (MQ=255)

                                                 gccggccgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTAc                                                                                                     <  2:5964174/150‑1 (MQ=255)

                                                  ccggccgTGCCGGGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTAcc                                                                                                    >  1:1248299/1‑150 (MQ=255)

                                                  ccggccgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTAcc                                                                                                    >  1:5981500/1‑150 (MQ=255)

                                                  caggacttGCCGAGCTGTTGCAGACCATCGTCAATAAAGCCCTGATCTACCTCAATAACCCACAGACGAAAAAGCAGCTTTCCAACATCCACCTTTACAAGTAACGGTAGCTCCTATCAGTGGCGCTGTGAACGCGCCATGAACGGTAcc                                                                                                    <  2:310531/143‑1 (MQ=255)

                                                       cgTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCTGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTg                                                                                               <  1:5162180/150‑1 (MQ=255)

                                                        gTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTgg                                                                                              <  2:7859856/150‑1 (MQ=255)

                                                        gTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTgg                                                                                              <  1:5773817/150‑1 (MQ=255)

                                                          gCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGc                                                                                            <  1:8484768/150‑1 (MQ=255)

                                                          gCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGc                                                                                            >  1:5539208/1‑150 (MQ=255)

                                                             gAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACCGCCGAACAAGCAGCTGCCCAACATCCACCGTTACAAGTACCGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGccac                                                                                         <  1:18646/150‑1 (MQ=255)

                                                               gCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGccacca                                                                                       <  1:7805975/150‑1 (MQ=255)

                                                                      cAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGAt                                                                                <  2:6473194/150‑1 (MQ=255)

                                                                      cAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGAt                                                                                <  2:7289731/150‑1 (MQ=255)

                                                                         aCCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGCCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGcc                                                                             <  1:3410727/150‑1 (MQ=255)

                                                                           cATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCgg                                                                           >  2:8796196/1‑150 (MQ=255)

                                                                            aTCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGa                                                                          >  2:1259869/1‑150 (MQ=255)

                                                                             tCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGaa                                                                         >  1:4283358/1‑150 (MQ=255)

                                                                             tCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGaa                                                                         >  2:5003401/1‑150 (MQ=255)

                                                                               ggCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTACCGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGaaaa                                                                       <  1:785109/150‑1 (MQ=255)

                                                                                gCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAg                                                                      >  2:6979260/1‑150 (MQ=255)

                                                                                gCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAg                                                                      <  1:6677298/150‑1 (MQ=255)

                                                                                 cAATAAAGCCCTGATCTACCGCAAGAACCCCCAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCggtcaggg                                                                     <  1:2182245/150‑7 (MQ=255)

                                                                                  aaGAAAGCCCTGATTTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGcc                                                                    >  1:5448893/1‑150 (MQ=255)

                                                                                  aaGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGcc                                                                    >  1:8084199/1‑150 (MQ=255)

                                                                                    gAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCa                                                                  <  2:7738540/150‑1 (MQ=255)

                                                                                         cccTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCCGCAcc                                                             >  1:3618269/1‑150 (MQ=255)

                                                                                          ccTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCa                                                            <  2:6624571/150‑1 (MQ=255)

                                                                                          ccTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCa                                                            <  2:2964223/150‑1 (MQ=255)

                                                                                           cTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAg                                                           >  1:1120946/1‑150 (MQ=255)

                                                                                           cTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAg                                                           >  2:5495472/1‑150 (MQ=255)

                                                                                              aTCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGaa                                                        <  2:2783910/150‑1 (MQ=255)

                                                                                               tCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCTTTACAAGTAACGGTAGCTCCTATCTGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAg                                                       <  1:1565388/150‑1 (MQ=255)

                                                                                                cTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGc                                                      >  1:2003997/1‑150 (MQ=255)

                                                                                                cTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGc                                                      >  1:4409227/1‑150 (MQ=255)

                                                                                                cTACCGCAAGAACCCACAGCCGAACAAGCAGCTGACCAACAACCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGc                                                      >  1:4409645/1‑150 (MQ=255)

                                                                                                 tACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCa                                                     <  1:6294931/150‑1 (MQ=255)

                                                                                                  aCCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGACAAGCCCAGCACCAGGAAGCAg                                                    >  1:4634849/1‑150 (MQ=255)

                                                                                                      cAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACa                                                >  2:4334892/1‑150 (MQ=255)

                                                                                                      cAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACa                                                >  2:135426/1‑150 (MQ=255)

                                                                                                          aaCCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGACCATCTg                                            >  1:4828883/1‑150 (MQ=255)

                                                                                                           aCCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGc                                           <  2:4674826/150‑1 (MQ=255)

                                                                                                           aCCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGc                                           <  2:2034298/150‑1 (MQ=255)

                                                                                                             ccACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCa                                         >  2:3498000/1‑150 (MQ=255)

                                                                                                               acaGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAgc                                       >  1:93749/1‑150 (MQ=255)

                                                                                                                 aGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAgcgc                                     >  2:7116748/1‑150 (MQ=255)

                                                                                                                   ccGAACAAGCAGCTGTCCAACATCCCCCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTc                                   <  2:5292165/150‑1 (MQ=255)

                                                                                                                    cGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCg                                  <  2:5213280/150‑1 (MQ=255)

                                                                                                                     gAACAAGCAGCTGTCCAACATCCCCCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGc                                 <  2:2212100/150‑1 (MQ=255)

                                                                                                                     gAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGc                                 <  2:8198841/150‑1 (MQ=255)

                                                                                                                      aaCAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGcc                                >  1:6431030/1‑150 (MQ=255)

                                                                                                                          agcagcTGTCAAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACg                            <  1:5394151/150‑1 (MQ=255)

                                                                                                                           gcagcTGCCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACgg                           >  1:7122378/1‑150 (MQ=255)

                                                                                                                           gcagcGGTCCAACCTCCACCGTTACAAGTACCGGTAGCTCCGATCAGGGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACgg                           <  1:8274657/150‑1 (MQ=255)

                                                                                                                             agcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCTCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGACGCAGAACATCTGCCAGCGCTCGCCCACGGaa                         >  1:8420410/1‑150 (MQ=255)

                                                                                                                             agcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGaa                         <  2:1817575/150‑1 (MQ=255)

                                                                                                                             agcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGaa                         >  2:1612477/1‑150 (MQ=255)

                                                                                                                             agcTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGaa                         >  1:1754810/1‑150 (MQ=255)

                                                                                                                                tGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATa                      <  2:4487983/150‑1 (MQ=255)

                                                                                                                                tGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATa                      <  2:6306120/150‑1 (MQ=255)

                                                                                                                                    cAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGt                  >  1:1020242/1‑150 (MQ=255)

                                                                                                                                     aaCATCCCCCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTa                 <  1:7466752/150‑1 (MQ=255)

                                                                                                                                     aaCATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGTCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTa                 <  2:3886392/150‑1 (MQ=255)

                                                                                                                                     aaCATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGACATACCGTa                 >  2:3962543/1‑150 (MQ=255)

                                                                                                                                     aaCATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTa                 <  1:6561132/150‑1 (MQ=255)

                                                                                                                                         tCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAgcgc             >  1:313518/1‑150 (MQ=255)

                                                                                                                                         tCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAgcgc             >  1:104530/1‑150 (MQ=255)

                                                                                                                                         tCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAgcgc             >  2:7583606/1‑150 (MQ=255)

                                                                                                                                           taCCGTTACTTGTAACGGTGGCTCCGATCAGTGGCGCGGTGAACACGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCaa           >  1:4085946/2‑150 (MQ=255)

                                                                                                                                           cACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCaa           <  1:1975762/150‑1 (MQ=255)

                                                                                                                                            aCCGTTACTTGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGaa                         >  1:5712832/1‑135 (MQ=255)

                                                                                                                                            aCCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAg          >  2:5915435/1‑150 (MQ=255)

                                                                                                                                              cGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGgt        <  2:1144370/150‑1 (MQ=255)

                                                                                                                                              cGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGgt        <  2:4925051/150‑1 (MQ=255)

                                                                                                                                              cGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGgt        >  2:8635516/1‑150 (MQ=255)

                                                                                                                                               gTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTTGCGCAAGgtg       <  2:1420100/150‑1 (MQ=255)

                                                                                                                                                ttACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGgtgt      <  2:3708806/150‑1 (MQ=255)

                                                                                                                                                    aaGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATa  >  2:3317642/1‑150 (MQ=255)

                                                                                                                                                    aaGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATa  >  1:8750722/1‑150 (MQ=255)

CTCGGAAGACCGCGAAGAACGCGCCGAGACCATCGCCGAGCTGTGCAAGGCCGGCCGTGCCGAGCTGGTGCAGACCATCGGCAAGAAAGCCCTGATCTACCGCAAGAACCCACAGCCGAACAAGCAGCTGTCCAACATCCACCGTTACAAGTAACGGTAGCTCCGATCAGTGGCGCGGTGAACGCGCCCTGACCGGTACCGGCTGGGCCACCAGCACGATGCCGGAAAAGCCCAGCACCAGGAAGCAGAACATCTGCCAGCGCTCGCCCACGGAAATACCGTAGCGCAAGGTGTAATA  >  NC_002947/5372631‑5372928

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