BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,084,593	A→T	6.1%	intergenic (‑67/+5)	PP_t54 ← / ← PP_t53	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,084,593	0	A	T	6.1%	752.9 / 17.0	241	intergenic (‑67/+5)	PP_t54/PP_t53	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand): ref base A (128/94); new base T (3/14); total (131/110)
Fisher's exact test for biased strand distribution p-value = 1.77e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 247 total aligned reads displayed.

GAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAAACAACTA  >  NC_002947/5084445‑5084741

  aaGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCaaaa                                                                                                                                                   <  1:1165222/150‑1 (MQ=255)

   aGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAg                                                                                                                                                  >  2:4306458/1‑150 (MQ=255)

       gCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                              <  1:3359155/150‑1 (MQ=255)

       gCGGACGCTGCCGGATCCGACCCCGCCACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGc                                                                                                                                              <  1:88823/150‑1 (MQ=18)

        cGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCg                                                                                                                                             <  2:159256/150‑1 (MQ=255)

          gACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTTTGAAGTGGcttt                                                                                                                                           <  1:3394329/150‑4 (MQ=255)

          gACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGcttt                                                                                                                                           <  1:1438068/150‑4 (MQ=255)

            cTCAGCAGGATTCGACCCTGCGACCGCTCGGTTCTTAGCCTAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACt                                                                                                                                         <  2:7186362/150‑1 (MQ=255)

             tCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCGATTGGCAGTTTTAGACCAGATACCCCGGGGTGAATGGAGCAGCCTTGCCAACAAAGCCACGTCAAAAGTGGCGGAcag                                                                                                                                        >  2:862871/1‑148 (MQ=17)

             tCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTTTGAAGTGGctttgtt                                                                                                                                        <  1:3773510/150‑7 (MQ=18)

                 cagGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTcagca                                                                                                                                    <  1:1674643/150‑1 (MQ=255)

                  agGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCCAGTGGCAGATTTAGACCAGATACCGCTGGATGATTGAAGCAGCCTTGCAAACAAAGCCACATCAAAACTGGAGGACTcagcag                                                                                                                                   >  2:7508427/1‑150 (MQ=18)

                      ttCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGAtt                                                                                                                               <  1:1938070/150‑1 (MQ=255)

                         gAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTTTGAAGTGGctt                                                                                                                                            >  2:2080821/1‑132 (MQ=255)

                          aaCCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGaa                                                                                                                           <  1:2342479/150‑1 (MQ=255)

                           aCCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAAc                                                                                                                          >  2:7752908/1‑150 (MQ=255)

                           aCCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAAc                                                                                                                          <  1:2198133/150‑1 (MQ=255)

                             cTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCt                                                                                                                        >  1:3878289/1‑150 (MQ=255)

                              tGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTg                                                                                                                       <  1:8239836/150‑1 (MQ=255)

                               gCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCCATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGc                                                                                                                      >  1:2717233/1‑150 (MQ=255)

                               gCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGc                                                                                                                      <  1:7536365/150‑1 (MQ=255)

                                cGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCg                                                                                                                     <  1:8594946/150‑1 (MQ=255)

                                cGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAAGTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCACAAGTGGCGGACTCAGCAGGATTCGAACCTGCg                                                                                                                     >  2:6433206/1‑150 (MQ=255)

                                 gACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGa                                                                                                                    <  2:1743410/150‑1 (MQ=255)

                                 gACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGa                                                                                                                    <  1:8290054/150‑1 (MQ=255)

                                  cccGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCAATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGAGGGTAGATTGAAGCAGCCTTGCAAACACAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGAc                                                                                                                   <  2:6906408/149‑1 (MQ=18)

                                      cTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTAAGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCt                                                                                                               <  2:5078460/150‑1 (MQ=255)

                                       tCGGTTCGTAGCCGAGTACTCTATCCAGCTGATCTATGAGTCCGCTATTGGCAGTTTTAGACCAGATCCCGCTGGTTGATTGAAGCAGCCTTGCAAACAACGCCACGTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTc                                                                                                              <  2:3165722/150‑1 (MQ=18)

                                       tCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTc                                                                                                              >  2:2238327/1‑150 (MQ=255)

                                           ttCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGtt                                                                                                          >  2:4711212/1‑150 (MQ=255)

                                              gTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGt                                                                                                       >  2:2471625/1‑150 (MQ=255)

                                                    gAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGa                                                                                                 <  1:8562802/150‑1 (MQ=255)

                                                        aCTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAc                                                                                             >  2:5959730/1‑150 (MQ=255)

                                                         ctctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTAct                                                                                            >  2:6194425/1‑150 (MQ=255)

                                                          tctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActc                                                                                           >  1:444896/1‑150 (MQ=255)

                                                          tctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActc                                                                                           >  1:7153064/1‑150 (MQ=255)

                                                           ctATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTActct                                                                                          <  1:8157347/150‑1 (MQ=255)

                                                            tATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                         >  1:7503247/1‑150 (MQ=255)

                                                            tATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTa                                                                                         >  1:1693530/1‑150 (MQ=255)

                                                             aTCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                                                        >  1:6881548/1‑150 (MQ=255)

                                                             aTCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTAt                                                                                        >  1:4864118/1‑150 (MQ=255)

                                                              tCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                       <  1:4089670/150‑1 (MQ=255)

                                                               ccAGCTGAGCTATGAGTCCGCATTTGGAAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATcc                                                                                      <  2:2922704/150‑1 (MQ=255)

                                                                 aGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                                                    >  2:391083/1‑150 (MQ=255)

                                                                 aGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAg                                                                                    >  1:3499025/1‑150 (MQ=255)

                                                                    tGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTg                                                                                 >  1:7407144/1‑150 (MQ=255)

                                                                     gAGCTAGGAGTCCGCTCGTGGCATTTTTAGAGCAGACCCCGCTGGTTGATAGAAGCAGCCTCGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGa                                                                                <  2:4487134/150‑1 (MQ=18)

                                                                         tATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTa                                                                            >  2:5835333/1‑150 (MQ=255)

                                                                         tATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTa                                                                            >  2:465220/1‑150 (MQ=255)

                                                                            gAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGa                                                                         >  2:1355208/1‑150 (MQ=255)

                                                                               tCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAATCAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTc                                                                      >  2:4919436/1‑150 (MQ=255)

                                                                               tCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTc                                                                      >  2:4355312/1‑150 (MQ=255)

                                                                                 cGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCg                                                                    >  1:5998710/1‑150 (MQ=255)

                                                                                 cGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCg                                                                    >  2:5351764/1‑150 (MQ=255)

                                                                                  gCACTTGTCAGTTATAGACAAGATACCTCTGGTTGATTGAAGCATCCTTTAAAACAAAGACACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCgt                                                                   <  2:5410329/150‑1 (MQ=18)

                                                                                         gCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGc                                                            >  2:4561998/1‑150 (MQ=255)

                                                                                         gCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGc                                                            >  2:3964510/1‑150 (MQ=255)

                                                                                         gCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGc                                                            >  1:6855784/1‑150 (MQ=255)

                                                                                         gCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGc                                                            >  2:1949086/1‑150 (MQ=255)

                                                                                               ttAGACCAGATCCCGCTGGTTGATTGAAGCATACTCTATAAGAAAGACACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCa                                                      <  2:2193818/150‑1 (MQ=18)

                                                                                                  gACCAGATACCGCGGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATta                                                   >  2:7827336/1‑150 (MQ=255)

                                                                                                   aCCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtat                                                  >  1:6526140/1‑150 (MQ=255)

                                                                                                    cgagctACCGCTGTTTGTTTTATGGAGCCTATCTAAGAACGACACTTCTAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtata                                                 <  1:5314161/145‑1 (MQ=17)

                                                                                                    ccAGATACCTCTGGTTGATTGAAGCAGCATTTCTAACAAAGCCTCGTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtata                                                 <  1:1788658/150‑1 (MQ=21)

                                                                                                    ccAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtata                                                 <  1:7034307/150‑1 (MQ=255)

                                                                                                    ccAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtata                                                 <  1:8062283/150‑1 (MQ=255)

                                                                                                    ccAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtata                                                 <  1:923732/150‑1 (MQ=255)

                                                                                                    ccAGAGCCCGCGGGGTGAGTGAAGCCGCCTGGCCCACAACGCCACGTCCAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATtata                                                 <  2:1471598/150‑1 (MQ=17)

                                                                                                     cAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAg                                                >  2:2702108/1‑150 (MQ=255)

                                                                                                     cAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAg                                                >  1:5222155/1‑150 (MQ=255)

                                                                                                     cAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAg                                                >  1:5058580/1‑150 (MQ=255)

                                                                                                       gATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGt                                              <  1:1501976/150‑1 (MQ=255)

                                                                                                             gCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGa                                        >  2:8287611/1‑150 (MQ=255)

                                                                                                               tGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGaaa                                      <  1:4668191/150‑1 (MQ=255)

                                                                                                                     attgaAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAtttt                                <  2:1240612/150‑1 (MQ=255)

                                                                                                                     attgaAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAtttt                                >  1:2983094/1‑150 (MQ=255)

                                                                                                                      ttgaAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACAttttt                               <  1:4689692/150‑1 (MQ=255)

                                                                                                                       tgaAGGAGCATTGCAAACAAAGCCTCTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTg                              <  1:5476045/150‑1 (MQ=255)

                                                                                                                       tgaAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATc                                                                                       <  1:2582475/93‑1 (MQ=255)

                                                                                                                         aaGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCa                            >  2:8189702/1‑150 (MQ=255)

                                                                                                                          agcagcCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCaa                           <  1:8343902/150‑1 (MQ=255)

                                                                                                                          agcagcCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCaa                           >  1:6168186/1‑150 (MQ=255)

                                                                                                                           gcagcCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCaag                          >  1:5581599/1‑150 (MQ=255)

                                                                                                                           gcagcCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCaag                          >  1:3261262/1‑150 (MQ=255)

                                                                                                                                 ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                    >  2:7279855/1‑150 (MQ=255)

                                                                                                                                 ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                    >  1:2644893/1‑150 (MQ=255)

                                                                                                                                 ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                    >  2:5120923/1‑150 (MQ=255)

                                                                                                                                 ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                    >  1:8661259/1‑150 (MQ=255)

                                                                                                                                 ttGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCtt                    >  1:3380321/1‑150 (MQ=255)

                                                                                                                                  tGCAAAGAAAGCCTTTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCttt                   <  2:8081105/150‑1 (MQ=255)

                                                                                                                                    caaacaaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTc                 >  2:6686484/1‑150 (MQ=255)

                                                                                                                                     aaacaaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCg                >  2:1141500/1‑150 (MQ=255)

                                                                                                                                      aacaaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTTGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGt               <  2:4643714/150‑1 (MQ=255)

                                                                                                                                      aacaaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGt               <  1:4306458/150‑1 (MQ=255)

                                                                                                                                         aaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtat            <  1:7210464/150‑1 (MQ=255)

                                                                                                                                         aaaGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtat            <  1:4355312/150‑1 (MQ=255)

                                                                                                                                          ccgCCACTTTTGAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTtata           >  2:2992060/3‑150 (MQ=255)

                                                                                                                                            gCCACTTTTGAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATaaa         >  2:8367779/1‑150 (MQ=255)

                                                                                                                                            gCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATaaa         >  2:658872/1‑150 (MQ=255)

                                                                                                                                                 ttCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAaacaac    <  1:778723/150‑1 (MQ=255)

GAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGCAATTGGCAGTTTTAGACCAGATACCGCTGGTTGATTGAAGCAGCCTTGCAAACAAAGCCACTTCAAAAGTGGCGGACTCAGCAGGATTCGAACCTGCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTATGAGTCCGTGTCTTGCCGCGCATTATAGGTCGCTGAAACATTTTTGCAAGAAGCTTTTCGTTATAAAACAACTA  >  NC_002947/5084445‑5084741

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