BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,084,426	G→T	7.8%	intergenic (‑51/+24)	PP_t55 ← / ← PP_t54	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,084,426	0	G	T	7.8%	434.4 / 10.6	140	intergenic (‑51/+24)	PP_t55/PP_t54	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base G (84/42); new base T (3/10); total (88/52)
Fisher's exact test for biased strand distribution p-value = 4.52e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

AAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTT  >  NC_002947/5084277‑5084575

aaGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATAACCCCTGGTTGAGTGAGATGCCATTg                                                                                                                                                       >  2:2662382/1‑150 (MQ=11)

   cAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCaa                                                                                                                                                    <  2:7012663/150‑1 (MQ=255)

    aGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAg                                                                                                                                                   <  1:6597766/150‑1 (MQ=255)

         cACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGaa                                                                                                                                              <  1:2053011/150‑1 (MQ=255)

            ttCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTATACCAGATAACCGCTGGTTGACTGAGATGGCTTTGAAACCAGACCCa                                                                                                                                           >  2:4466065/1‑150 (MQ=14)

            ttCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGACATGGCTTTGCAAGCAGACCCa                                                                                                                                           >  2:6743059/1‑150 (MQ=255)

            ttCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATAACCGCTGGTTGACTGAAAAGGCTTTGAAACCAGAACCa                                                                                                                                           >  2:4466130/1‑150 (MQ=255)

             tCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCAc                                                                                                                                          >  1:1120212/1‑150 (MQ=255)

                tgattgTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTc                                                                                                                                       <  2:4293873/145‑1 (MQ=255)

                  aaaGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCaa                                                                                                                                     <  1:3814944/150‑1 (MQ=255)

                      tGGCGGACTCAGCAGGATTCGTACCTGCGACCGCTCGGTTCGTAGCCGATTACTCTATCCAGCTGAGCTATGTGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGaaa                                                                                                                                 <  1:3445247/150‑1 (MQ=255)

                          ggACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGTTTGAGTGAAAAGGCTTTGCAAGAAGAACCACTACAAGAAGGTgg                                                                                                                             >  2:402144/1‑150 (MQ=17)

                          ggACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTgg                                                                                                                             <  2:3432080/150‑1 (MQ=255)

                            aCTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCg                                                                                                                           <  1:4208155/150‑1 (MQ=255)

                                 gcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTCCCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTc                                                                                                                      >  2:5652046/1‑150 (MQ=255)

                                 gcagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGGTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGAGGCGGtctc                                                                                                                      >  2:2118194/1‑150 (MQ=255)

                                   agGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAGATGGCTTTGCAAGCAGAACCACTTCAAGAGTGTGGCGGACTcag                                                                                                                    >  2:1327081/1‑150 (MQ=255)

                                   agGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGACCTGGCTGTGCAAGCAGGACCACTTCAAGATGGTGGAGGACTcag                                                                                                                    >  2:4026980/1‑150 (MQ=17)

                                     gATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTcagca                                                                                                                  <  1:6102903/150‑1 (MQ=255)

                                       ttCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAgg                                                                                                                >  2:2073200/1‑150 (MQ=255)

                                           aaCCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTc                                                                                                            <  1:7042062/150‑1 (MQ=255)

                                           aaCCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTc                                                                                                            <  1:3303104/150‑1 (MQ=255)

                                             ccTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGa                                                                                                          >  2:6961190/1‑150 (MQ=255)

                                                 cGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTTTGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCt                                                                                                      <  2:7296970/150‑1 (MQ=255)

                                                     cGCTCGGTTCGTAGCCGAGTACTCTAGCCAGCGGAGCTATGAGGCCGCGATGGGCAGTTTTAAACCAGTTTACCGCTGGTGGACTGAAAAGGCGTGGCAAGCCGACCCCCTGCAAGATGGTGGCGGCCTCAGCAGGATTCGCTCCTGCgc                                                                                                  >  2:6793005/1‑149 (MQ=17)

                                                        tCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGAcca                                                                                               >  2:1343407/1‑149 (MQ=255)

                                                            ttCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTTACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCg                                                                                           <  1:6985738/150‑1 (MQ=255)

                                                             tCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCgg                                                                                          >  1:3483799/1‑150 (MQ=255)

                                                              cGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGt                                                                                         <  1:6995519/150‑1 (MQ=255)

                                                               gTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGtt                                                                                        >  2:874179/1‑150 (MQ=255)

                                                                tAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTc                                                                                       >  2:3769022/1‑150 (MQ=255)

                                                                     gAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGc                                                                                  >  1:4864342/1‑150 (MQ=255)

                                                                     gAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGc                                                                                  >  1:7275321/1‑150 (MQ=255)

                                                                     gAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGc                                                                                  >  1:822497/1‑150 (MQ=255)

                                                                      aGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGcc                                                                                 >  1:1750442/1‑150 (MQ=255)

                                                                        tACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                               <  1:2393548/150‑1 (MQ=255)

                                                                           tctATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTACACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTa                                                                            <  2:5296296/150‑1 (MQ=255)

                                                                           tctATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTa                                                                            >  1:2551519/1‑150 (MQ=255)

                                                                             tATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAct                                                                          >  1:511262/1‑150 (MQ=255)

                                                                             tATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAct                                                                          >  1:352543/1‑150 (MQ=255)

                                                                               tCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAAc                                                                                                                                             >  1:1254371/1‑81 (MQ=255)

                                                                               tCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAAc                                                                                                                                             <  2:1254371/81‑1 (MQ=255)

                                                                               tCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActct                                                                        >  1:6086399/1‑150 (MQ=255)

                                                                                       tgCTTTGAGGCCGCTAATGGCAGTTTTACACCAGATTACCGCTGGGTGACTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                                <  2:5173571/149‑1 (MQ=14)

                                                                                       aGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCt                                                                <  2:1120212/150‑1 (MQ=255)

                                                                                         cTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                              >  1:5890342/1‑150 (MQ=255)

                                                                                            tGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCt                                                           >  1:7344198/1‑150 (MQ=255)

                                                                                             gAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTa                                                          >  1:5123826/1‑150 (MQ=255)

                                                                                              agcccGCGTTTGGCAGTTTTATACCAGTTTACCTCTGGCTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                         <  2:4305220/147‑1 (MQ=11)

                                                                                              actCCGCGATTCGCAGTTTAAAACCAGATTACCTCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                         <  2:2016545/148‑1 (MQ=255)

                                                                                              aGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAt                                                         <  1:2835669/150‑1 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTa                                                                                    <  2:5269025/121‑1 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTa                                                                                    >  1:5269025/1‑121 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  1:4807173/1‑150 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  2:5503500/1‑150 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  1:6487381/1‑150 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  1:4404093/1‑150 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  1:3585942/1‑150 (MQ=255)

                                                                                                tCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  1:3697062/1‑150 (MQ=255)

                                                                                                 ccGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAg                                                                                                                                <  2:1560408/76‑1 (MQ=255)

                                                                                                  cgcgATTGGCAGTTTTAAACCAGATAACCGCTGGTTGAATGAAATGGCTTTCCAAGCAGAAACACTTCAAGAAAGAGGCCGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGACGAGTACTCAAACCAGCTGAGCTAagaga                                                     >  2:307292/1‑149 (MQ=17)

                                                                                                   gcgATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTc                                                    >  1:7423877/1‑150 (MQ=255)

                                                                                                    cgATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCa                                                                                                                                                 <  1:2942207/56‑1 (MQ=255)

                                                                                                    cgATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCa                                                                                                                                                 >  2:2942207/1‑56 (MQ=255)

                                                                                                          gCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAAtt                                             >  1:1001255/1‑150 (MQ=255)

                                                                                                          gCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAAtt                                             >  2:2121725/1‑150 (MQ=255)

                                                                                                           cAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTg                                            >  1:5042695/1‑150 (MQ=255)

                                                                                                            aGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCgagagtgg                                           >  2:1580892/1‑143 (MQ=255)

                                                                                                                   aaCCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGattt                                    >  1:7139669/1‑150 (MQ=255)

                                                                                                                    aCCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTa                                   >  2:3243046/1‑150 (MQ=255)

                                                                                                                    aCCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTa                                   >  2:5722373/1‑150 (MQ=255)

                                                                                                                     ccAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAg                                  >  1:6640380/1‑150 (MQ=255)

                                                                                                                                cTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCCTATTTAGACCAGATACCg                       >  1:1166825/1‑150 (MQ=255)

                                                                                                                                cTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCg                       >  1:164593/1‑150 (MQ=255)

                                                                                                                                      taCTGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGtt                 <  1:874094/149‑1 (MQ=25)

                                                                                                                                       aCTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       <  2:1410033/111‑1 (MQ=32)

                                                                                                                                       aCTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                       >  1:1410033/1‑111 (MQ=32)

                                                                                                                                       aCTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGttg                >  2:738935/1‑150 (MQ=40)

                                                                                                                                         tGAAATGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAg                                  <  1:245729/130‑1 (MQ=11)

                                                                                                                                         tGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGACACTGCGACCGCTAGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGAAAGTGGCAGGTTTAGACCAGATACCGCTGGtttat              >  2:659995/1‑150 (MQ=14)

                                                                                                                                         tGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGttgat              >  2:1428432/1‑150 (MQ=255)

                                                                                                                                         tGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGttgat              >  2:621914/1‑150 (MQ=255)

                                                                                                                                         tGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGttgat              >  1:1463046/1‑150 (MQ=255)

                                                                                                                                           aagatgCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGttgattg            <  1:7314979/145‑1 (MQ=255)

                                                                                                                                           aaaTGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGGGTTAGACCAGATACCGCTGGttgattg            >  2:7034357/1‑150 (MQ=255)

                                                                                                                                             aTGTCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGaa          <  1:5292582/150‑1 (MQ=255)

                                                                                                                                             aTGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGaa          >  2:1246042/1‑150 (MQ=38)

                                                                                                                                              tGGCTTTTCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAag         <  1:3118273/150‑1 (MQ=255)

                                                                                                                                                gCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCTGATACCGCTGGTTGATTGAagca       <  1:5652046/150‑1 (MQ=255)

                                                                                                                                                 cTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGa                                 >  2:7441414/1‑123 (MQ=37)

                                                                                                                                                 cTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGa                                 <  1:7441414/123‑1 (MQ=37)

                                                                                                                                                  ttttCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGaa          <  1:1067254/145‑1 (MQ=38)

                                                                                                                                                  ttttCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGaa          >  2:1067254/1‑145 (MQ=38)

                                                                                                                                                  tttGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                              <  1:2078491/93‑1 (MQ=33)

                                                                                                                                                  tttGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAagcagc     >  2:6852069/1‑150 (MQ=255)

                                                                                                                                                  tttGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAagcagc     >  1:4563800/1‑150 (MQ=255)

                                                                                                                                                  tttGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAAGAGTCCGCTATTGGCAGTGTTAGACCAGATACCGCTGGTTGATTGAagcagc     >  1:4563204/1‑150 (MQ=255)

                                                                                                                                                  ttagCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAagcagc     <  1:2662382/147‑1 (MQ=255)

                                                                                                                                                     gCATGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCTGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCtt  <  1:6280686/150‑1 (MQ=255)

                                                                                                                                                     gCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCtt  <  2:3347688/150‑1 (MQ=255)

AAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAAACCAGATTACCGCTGGTTGACTGAAATGGCTTTGCAAGCAGAACCACTTCAAGAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTT  >  NC_002947/5084277‑5084575

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